A computational workflow for rupture‐to‐structural‐response simulation and its application to Istanbul

Scenario‐based earthquake simulations at regional scales hold the promise in advancing the state‐of‐the‐art in seismic risk assessment studies. In this study, a computational workflow is presented that combines (i) a broadband Green's function‐based fault‐rupture and ground motion simulation—herein carried out using the “UCSB (University of California at Santa Barbara) method”, (ii) a three‐dimensional physics‐based regional‐scale wave propagation simulation that is resolved at  Hz, and (iii) a local soil‐foundation‐structure finite element analysis model. These models are interfaced with each other using the domain reduction method. The innermost local model—implemented in ABAQUS—is additionally enveloped with perfectly matched layer boundaries that absorb outbound waves scattered by the structures contained within it. The intermediate wave propagation simulation is carried out using Hercules , which is an explicit time‐stepping finite element code that is developed and licensed by the CMU‐QUAKE group. The devised workflow is applied to a  km region on the European side of Istanbul, which was modeled using detailed soil stratigraphy data and realistic fault rupture properties, which are available from prior microzonation surveys and earthquake scenario studies. The innermost local model comprises a chevron‐braced steel frame building supported by a shallow foundation slab, which, in turn, rests atop a three‐dimensional soil domain. To demonstrate the utility of the workflow, results obtained using various simplified soil‐structure interaction analysis techniques are compared with those from the detailed direct model. While the aforementioned demonstration has a limited scope, the devised workflow can be used in a multitude of ways, for example, to examine the effects of shallow‐layer soil nonlinearities and surface topography, to devise site‐ and structure‐specific seismic fragilities, and for calibrating regional loss models, to name a few.     

用近场记录研究唐山地震的震源过程

通过对唐山地震近场记录的波形拟合,研究了震源短周期运动过程的基本特征,推断了主要破裂事件的震源参数.结果表明,高频地震波中的主要部分来自几个局部区域的破裂过程.它们的错距分布相当集中,上升时间为1.1-1.5s.由5次破裂事件组成的震源模型中引人随机扰动量,较好地描绘了地面运动的高频状态.     

Estimation of Winds at Different Isobaric Levels Based on the Observed Winds at 850 hPa Level Using Double Fourier Series

ＥｓｔｉｍａｔｉｏｎｏｆＷｉｎｄｓａｔＤｉｆｆｅｒｅｎｔlｓｏｂａｒｉｃＬｅｖｅｌｓＢａｓｅｄｏｎｔｈｅＯｂｓｅｒｖｅｄ　Ｗｉｎｄｓａｔ８５０ｈＰａＬｅｖｅｌＵｓｉｎｇＤｏｕｂｌｅＦｏｕｒｉｅｒＳｅｒｉｅｓＳ．Ｎ．ＢａｖａｄｅｋａｒａｎｄＲ．Ｍ．Ｋｈ?..     

Application of numerical wave-propagation techniques to study local soil effects: The case of Benevento (Italy)

Many of the numerical techniques used for seismic zonation studies treat one-dimensional structural models and/or the incidence of plane polarized body waves. These techniques are often not adequate for laterally heterogeneous structures and for sources that are not located beneath the site of interest. In such cases a more rigorous treatment of the combined effects of the source, the path and the site response is needed. This can be accomplished with a hybrid approach combining modal summation and the finite-difference technique. To demonstrate the differences between these techniques, the ground motion in the city of Benevento (Italy) is modelled. We first compare the results obtained with one-and two-dimensional structural models for vertical incidence of plane polarized body waves. These results are then compared with those obtained with the hybrid approach for two-dimensional structural models.The comparisons have allowed us to find important differences in the response obtained with the different modelling techniques. For the same site, these differences consist of strong variations in amplitude and in the shape of the spectral amplifications. For a seismic source which is not located beneath the site, vertical incidence of waves significantly overestimates the local hazard in a laterally homogeneous structure. For a laterally heterogeneous area, we can conclude that one-dimensional modelling fails to estimate the seismic hazard, whereas for a seismic source which is not located beneath the site of interest, two-dimensional modelling with vertical incidence of plane polarized body waves may not allow reliable estimates to be made of the frequency bands at which amplifications occur. The results obtained for two-dimensional structural models are used for a zonation of the city of Benevento.     

Seismic ground motion in a layered half-space due to a Haskell-type source. I: Theory

In this paper, closed-form analytic expressions for the frequency-wave number domain Fourier amplitudes of the displacement field at the free surface of a layered, anelastic half-space are established. The displacement field is caused by a seismic source described by a shear dislocation propagating with constant velocity over a rectangular fault (Haskell's model). Three-dimensional plane wave propagation is considered in the layered half-space using a propagator-based formalism. The wave radiation from the source is decoupled into P-SV and SH motions and the two problems are treated separately. First, analytic expressions are calculated for the displacement field at the free surface due to unidirectional unit impulses. Then, these expressions are used to compute solutions for the displacement field due to effective point sources associated with a pure strike slip and a pure dip slip. Finally, these solutions are combined and integrated over the rectangular fault area to establish closed-form analytic expressions of the total displacement field at the free surface.     

板块构造理论研究新进展

本文评述了80年代以来板块构造理论的重要研究成果,就板块划分与板块运动学研究,大洋中脊与板块的扩张,活动边缘与板块的俯冲,板块再造和造山研究等方面的进展,展示板块构造的许多概念已更趋复杂化,板块构造理论变得更加成熟了。     

Complex Site Effects in Thessaloniki (Greece): II. 2D SH Modelling and Engineering Insights

This paper presents results of numerical modelling of site response for Thessaloniki, obtained with two different 2D methods; a finite difference and a finite element method. Ground motion across a 2D model of the subsoil of the city has been simulated for vertically incident SH waves. The predominance of locally generated surface waves is very clear in the synthetic seismograms of a weak event and of stronger ones. These results are then compared with the observations in time domain and frequency domain. The role of the soil formations with high attenuation in the lateral propagation and the effect of the differential motion close to the lateral variations are also pinpointed. The stronger events were finally used to compute strong ground motion in order to reveal and to discuss practical engineering aspects such as peak ground acceleration value, the most familiar indicator in seismic norms, the soil to rock spectral coefficients for the period bandwidth of interest, and the aggravation factor in terms of 2D to 1D response spectra as a useful ruler to account for complex site effects.     

One to One Resonance at High Inclination

We report results from long term numerical integrations and analytical studies of particular orbits in the circular restricted three-body problem. These are mostly high-inclination trajectories in 1 : 1 resonance starting at or near the triangular Lagrangian L₅ point. In some intervals of inclination these orbits have short Lyapunov times, from 100 to a few hundred periods, yet the osculating semi-major axis shows only relatively small fluctuations and there are no escapes from the 1 : 1 resonance. The eccentricity of these chaotic orbits varies in an erratic manner, some of the orbits becoming temporarily almost rectilinear. Similarly the inclination experiences large variations due to the conservation of the Jacobi constant. We studied such orbits for up to 10⁹ periods in two cases and for 10⁶ periods in all others, for inclinations varying from 0° to 180°. Thus our integrations extend from thousands to 10 million Lyapunov times without escapes of the massless particle. Since there are no zero-velocity curves restricting the motion this opens questions concerning the reason for the persistence of the 1 : 1 resonant motion. In the theory sections we consider the mechanism responsible for the observed behavior. We derive the averaged 1 : 1 resonance disturbing function, to second order in eccentricity, to calculate a critical inclination found in the numerical experiment, and examine motion close to this inclination. The cause of the chaos observed in the numerical experiments appears to be the emergence of saddle points in the averaged disturbing potential. We determine the location of several such saddle points in the (, ) plane, with being the mean longitude difference and the argument of pericentre. Some of the saddle points are illustrated with the aid of contour plots of the disturbing function. Motion close to these saddles is sensitive to initial conditions, thus causing chaos.This revised version was published online in October 2005 with corrections to the Cover Date.