Archaeoseismology: Methodological issues and procedure

Archaeoseismic research contributes important data on past earthquakes. A limitation of the usefulness of archaeoseismology is due to the lack of continuous discussion about the methodology. The methodological issues are particularly important because archaeoseismological investigations of past earthquakes make use of a large variety of methods. Typical in situ investigations include: (1) reconstruction of the local archaeological stratigraphy aimed at defining the correct position and chronology of a destruction layer, presumably related to an earthquake; (2) analysis of the deformations potentially due to seismic shaking or secondary earthquake effects, detectable on walls; (3) analysis of the depositional characteristics of the collapsed material; (4) investigations of the local geology and geomorphology to define possible natural cause(s) of the destruction; (5) investigations of the local factors affecting the ground motion amplifications; and (6) estimation of the dynamic excitation, which affected the site under investigation. Subsequently, a ‘territorial’ approach testing evidence of synchronous destruction in a certain region may delineate the extent of the area struck by the earthquake. The most reliable results of an archaeoseismological investigation are obtained by application of modern geoarchaeological practice (archaeological stratigraphy plus geological–geomorphological data), with the addition of a geophysical-engineering quantitative approach and (if available) historical information. This gives a basic dataset necessary to perform quantitative analyses which, in turn, corroborate the archaeoseismic hypothesis. Since archaeoseismological investigations can reveal the possible natural causes of destruction at a site, they contribute to the wider field of environmental archaeology, that seeks to define the history of the relationship between humans and the environment. Finally, through the improvement of the knowledge on the past seismicity, these studies can contribute to the regional estimation of seismic hazard.     

Lunar noise correlation,imaging and monitoring

Passive seismic techniques have revolutionarised seismology,leading for example to increased resolution in surface wave tomography,to the possibility to monitor changes in the propagation medium,and to many new processing strategies in seismic exploration.Here we review applications of the new techniques to a very particular dataset,namely data from the Apollo 17 lunar network.The special conditions of the lunar noise environment are investigated,illustrating the interplay between the properties of the noise and the ability to reconstruct Green's functions.With a dispersion analysis of reconstructed Rayleigh waves new information about the shallow shear velocity structure of the Moon are obtained.Passive image interferometry is used to study the effect of temperature changes in the subsurface on the seismic velocities providing direct observation of a dynamic process in the lunar environment.These applications highlight the potential of passive techniques for terrestrial and planetary seismology.     

Predicting the effect of biofouling on ship resistance using CFD

This paper proposes a Computational Fluid Dynamics (CFD) based unsteady RANS model which enables the prediction of the effect of marine coatings and biofouling on ship resistance and presents CFD simulations of the roughness effects on the resistance and effective power of the full-scale 3D KRISO Container Ship (KCS) hull.Initially, a roughness function model representing a typical coating and different fouling conditions was developed by using the roughness functions given in the literature. This model then was employed in the wall-function of the CFD software and the effects of a typical as applied coating and different fouling conditions on the frictional resistance of flat plates representing the KCS were predicted for a design speed of 24 knots and a slow steaming speed of 19 knots using the proposed CFD model. The roughness effects of such conditions on the resistance components and effective power of the full-scale 3D KCS model were then predicted at the same speeds. The resulting frictional resistance values of the present study were then compared with each other and with results obtained using the similarity law analysis. The increase in the effective power of the full-scale KCS hull was predicted to be 18.1% for a deteriorated coating or light slime whereas that due to heavy slime was predicted to be 38% at a ship speed of 24 knots. In addition, it was observed that the wave resistance and wave systems are significantly affected by the hull roughness and hence viscosity.     

Micro-fracture instabilities in granular solids

We study the aspect of unstable behavior (like strain localization bands) in elastic solids as a consequence of micro-fracturing. A two-scale approach of computational homogenization is considered. The macroscopic behavior is investigated by finite element computations on a unit cell. At the micro-level, we consider a granular structure with elastic grains. The inter-granular boundaries are modeled with cohesive laws, friction and unilateral contact. We show that decohesion between grains gives rise to macro-instabilities, indicated by the loss of ellipticity, typical for deformation localization bands. The relation between the microscopic softening on inter-granular boundaries and the onset of macro-instabilities is studied through numerical examples. The influence of the cohesive law and friction parameters is analyzed. For periodic distributions of granular structures, we prove the loss of periodicity by failure and the corresponding size dependence effect in the homogenized response. We present numerical examples of bifurcation of solutions for granular cell structures and of particular solutions specific to elementary volumes with periodic cell distribution. Size dependence appears in the unstable regime and is strongly influenced by cohesion and friction parameters.     

关于控制迭代过程的结束条件

迭代法属于逐次逼近法,应用很广。它是求解函数方程、方程组和微分方程等的一种基本而重要的数值方法。在具体求解时,需要给出控制迭代过程结束的条件。现有的某些计算方法(或数值分析)的教材和书,在论述求解函数方程的计算方案、流程图(或框图)及算法语言程序时,常常将此条件取为相邻两次迭代值xk、xk-1的差小于给定的准确度本文指出这个条件是不完全的,讨论了该条件的正确表达式,由误差估计式给出了控制迭代过程结束条件的五种形式。     

地形追随坐标系下埃特尔位涡的计算

随着数值天气预报模式的日益普及以及位涡在不同天气、气候时、空尺度上的广泛应用,对采用模式输出资料计算的埃特尔（Ertel）位涡的精度提出了更高要求。将气象要素从中尺度模式普遍采用的地形追随坐标系向定义位涡的z或p坐标系的插值过程（插值法）,是导致埃特尔位涡计算误差的主要原因。文中利用地形追随坐标系与z（p）坐标系的变换关系,导出了z坐标系定义的埃特尔位涡在地形追随坐标系中的表达形式,实现了直接在模式格点上计算埃特尔位涡（直接法）,在提高计算精度的同时,保留了其在z坐标系中的物理意义,方便了分析与应用。为了进一步分析直接法在减少计算误差方面的作用,使用WRF模式对2016年7月发生在华北地区复杂地形背景下的一次黄淮气旋爆发过程进行了模拟,利用模式输出资料对此次天气过程中强降水时段插值法计算误差的分布进行了分析。结果表明,利用直接法计算位涡可有效减少插值法引入的计算误差,中低层均方根误差可达0.5 PVU,中高层可达0.3 PVU。将其应用到中小尺度对流性天气精细化结构的分析及气候统计研究中,可以有效减少插值法误差对结果产生的不良影响,提高计算和分析准确度。      

SNREI地球模型负荷勒夫数数值计算的新进展

基于Ｌｏｎｇｍａｎ和Ｆａｒｒｅｌｌ的表面负荷形变理论，深入研究计算ＳＮＲＥＩ弹性地球模型表面负荷勒夫数的数值计算方法，取得新的进展：采用了一种新的代换方法─ｒ￣ｎ因子法，比传统的无量纲化法简单易行，且能达到很高的计算精度；确立了检验表面边界条件、微分方程积分线性性和计算结果的收敛性，用于评价数值计算结果精度的有效准则；给出了ＰＲＥＭ地球模型的一组高精度计算结果．     

差分法计算地震动旋转分量

旋转地震学是一门研究由地震、爆破以及工程振动等引起的地球介质运动的新兴学科.旋转运动的研究由来已久,但是由于缺少高精度的旋转分量地震仪,所以旋转运动的研究大多仅限于理论方面.差分法作为利用平动分量获取旋转分量的一种计算方法,在理论研究方面较为成熟,但是缺乏实际数据的验证.通过在对旋转运动研究现状充分调研的基础上,利用模...     

A study of microearthquake seismicity and focal mechanisms within the Sea of Marmara (NW Turkey) using ocean bottom seismometers (OBSs)

We have carried out seismological observations within the Sea of Marmara (NW Turkey) in order to investigate the seismicity induced after Gölcük–İzmit (Kocaeli) earthquake (M_w 7.4) of August 17, 1999, using ocean bottom seismometers (OBSs). High-resolution hypocenters and focal mechanisms of microearthquakes have been investigated during this Marmara Sea OBS project involving deployment of 10 OBSs within the Çınarcık (eastern Marmara Sea) and Central-Tekirdağ (western Marmara Sea) basins during April–July 2000. Little was known about microearthquake activity and their source mechanisms in the Marmara Sea. We have detected numerous microearthquakes within the main basins of the Sea of Marmara along the imaged strands of the North Anatolian Fault (NAF). We obtained more than 350 well-constrained hypocenters and nine composite focal mechanisms during 70 days of observation. Microseismicity mainly occurred along the Main Marmara Fault (MMF) in the Marmara Sea. There are a few events along the Southern Shelf. Seismic activity along the Main Marmara Fault is quite high, and focal depth distribution was shallower than 20 km along the western part of this fault, and shallower than 15 km along its eastern part. From high-resolution relative relocation studies of some of the microearthquake clusters, we suggest that the western Main Marmara Fault is subvertical and the eastern Main Marmara Fault dips to south at 45°. Composite focal mechanisms show a strike-slip regime on the western Main Marmara Fault and complex faulting (strike-slip and normal faulting) on the eastern Main Marmara Fault.