Geodynamics of the deep sedimentary basin of the Caspian Sea region: paragenetic correlation of seismicity and mud volcanism

In this paper we analyze seismic regime and earthquake depth distribution and correlation of seismicity and mud volcanism in the Azerbaijan and the Caspian Sea region. For the present region we have calculated accurate source locations, seismic activity, earthquake repetition and released earthquake energy parameters. It is shown that the active tectonic processes in the region are concentrated within the thick sedimentary cover that we consider as a general source of contemporary stress and a main structural element responsible for the origin of regional earthquakes. The correlation of seismicity and mud volcanism is of paragenetic character.     

A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluation of seismic ground motion induced by topographic irregularity

Results of an extensive numerical study on the 2D scattering of seismic waves by local topography are presented. The investigation has been conducted using the direct boundary element method. Several types of topography (slopes, canyons and ridges) are considered. The influences of some key parameters, such as exciting frequency and geometry of the irregular feature, on surface ground motion are studied in detail. It is found that local topographic conditions play an important role in the modification of seismic ground motion at the irregular feature itself and its neighbourhood. The present results can be considered to be useful from the viewpoint of earthquake engineering and seismology.     

Two formulations for dynamic response of a cylindrical cavity in cross‐anisotropic porous media

Two formulations for calculating dynamic response of a cylindrical cavity in cross‐anisotropic porous media based on complex functions theory are presented. The basis of the method is the solution of Biot's consolidation equations in the complex plane. Employing two groups of potential functions for solid skeleton and pore fluid (each group includes three functions), the u–w formulation of Biot's equations are solved. Difference of these two solutions refers to use of two various potential functions. Equations for calculating stress, displacement and pore pressure fields of the medium are mentioned based on each two formulations. Copyright © 2009 John Wiley & Sons, Ltd.     

Wave-induced stresses and pore pressures in sloping seabeds

By using the Mellin transform, a general solution is obtained for the total stresses and pore pressure distribution in an infinite wedge under fairly general conditions of loading. The results for the particular case in which each surface is subjected to static sinusoidal pressures are given in an explicit form. The infinite integrals can be evaluated by numerical integration methods. This analytical method is verified by different cases. The results of this study could be used to determine the total stresses and pore pressures induced by a wave on a sloping seabed.     

Numerical Simulation of Solute Transport in Saturated Porous Media with Bounded Domains

This work presents the first attempt to develop unconditionally stable, implicit finite difference solutions of one-sided spatial fractional advection-dispersion equation (s-FADE) by imposing the nonzero Dirichlet boundary condition (ND BC) or the nonzero fractional Robin boundary condition (NFR BC) at inlet boundary and the zero fractional Neumann boundary condition (ZFN BC) at outlet boundary. The results of the numerical studies performed using artificial solute transport parameters demonstrated that the numerical solution with the NFR BC as the inlet boundary produced much more realistic concentration values. The numerical solution with the NFR BC at the inlet boundary was capable of correctly describing the Fickian and non-Fickian behaviors of the solute transport at different α values, and it had the relatively same accuracy at different numbers of the spatial nodes. Also, the practical application of the numerical solution with the NFR BC as the inlet boundary was investigated by conducting tracer experiments in homogeneous and heterogeneous soil columns. According to the obtained results, this numerical solution described well solute transport in the homogenous and heterogeneous soils. The α values of the homogeneous and heterogeneous soils were obtained in the ranges of 1.849–1.999 and 1.248–1.570, respectively, which were in excellent agreement with the physical properties of the soils. In a nutshell, the numerical solution of the s-FADE with the NFR BC as the inlet boundary can be successfully applied to describe the solute transport in the homogeneous and heterogeneous soils with bounded spatial domains.     

Experimental vulnerability curves for the residential buildings of Iran

Iran is one of the most seismically active countries of the world located on the Alpine-Himalayan earthquake belt. More than 180,000 people were killed due to earthquakes in Iran during the last five decades. Considering the fact that most Iranians live in masonry and non-engineered houses, having a comprehensive program for decreasing the vulnerability of society holds considerable importance. For this reason, loss estimation should be done before an earthquake strikes to prepare proper information for designing and selection of emergency plans and the retrofitting strategies prior to occurrence of earthquake. The loss estimation process consists of two principal steps of hazard analysis and vulnerability assessment. After identifying the earthquake hazard, the first step is to evaluate the vulnerability of residential buildings and lifelines and also the social and economic impacts of the earthquake scenarios. Among these, residential buildings have specific importance, because their destruction will disturb the daily life and result in casualties. Consequently, the vulnerability assessment of the buildings in Iran is important to identify the weak points in the built environment structure. The aim of this research is to prepare vulnerability curves for the residential buildings of Iran to provide a proper base for estimating probable damage features by future earthquakes. The estimation may contribute fundamentally for better seismic performance of Iranian societies. After a brief review of the vulnerability assessment methods in Iran and other countries, through the use of the European Macroseismic method, a model for evaluating the vulnerability of the Iranian buildings is proposed. This method allows the vulnerability assessment for numerous sets of buildings by defining the vulnerability curves for each building type based on the damage observations of previous earthquakes. For defining the vulnerability curves, a building typology classification is presented in this article, which is representative of Iranian building characteristics. The hazard is described in terms of the macroseismic intensity and the EMS-98 damage grades have been considered for classifying the physical damage to the buildings. The calculated vulnerability indexes and vulnerability curves show that for engineered houses there is not any notable difference between the vulnerability of Iranian and Risk-UE building types. For the non-engineered houses, the vulnerability index of brick and steel structures is less than the corresponding values of the other unreinforced masonry buildings of Iran. The vulnerability index of unreinforced and masonry buildings of Iran are larger than the values of the similar types in Risk-UE and so the Iranian buildings are more vulnerable in this regard.     

Boundary integral formulation and two-dimensional fundamental solutions for dynamic behavior analysis of unsaturated soils

In this paper the coupled equations governing the dynamic behavior of unsaturated soils are derived based on the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) [Gatmiri B. Analysis of fully coupled behavior of unsaturated porous medium under stress, suction and temperature gradient. Final report of CERMES-EDF, 1997] and Gatmiri et al. (1998) [Gatmiri B, Delage P, Cerrolaza M, UDAM: a powerful finite element software for the analysis of unsaturated porous media. Adv Eng Software 1998; 29(1): 29–43]. In this formulation, the solid skeleton displacements, water pressure and air pressure are presumed to be independent variables. The Boundary Integral formulations as well as fundamental solutions for such a dynamic u–p_w–p_a theory are presented in this paper for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated two dimensional (2D) fundamental solutions for such deformable porous medium with linear elastic behavior are derived in Laplace transform domain using the method of Hörmander. Finally, some numerical results are presented to show the accuracy of the proposed solutions. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case.     

Predicting river water temperatures using stochastic models: case study of the Moisie River (Québec,Canada)

Successful applications of stochastic models for simulating and predicting daily stream temperature have been reported in the literature. These stochastic models have been generally tested on small rivers and have used only air temperature as an exogenous variable. This study investigates the stochastic modelling of daily mean stream water temperatures on the Moisie River, a relatively large unregulated river located in Québec, Canada. The objective of the study is to compare different stochastic approaches previously used on small streams to relate mean daily water temperatures to air temperatures and streamflow indices. Various stochastic approaches are used to model the water temperature residuals, representing short‐term variations, which were obtained by subtracting the seasonal components from water temperature time‐series. The first three models, a multiple regression, a second‐order autoregressive model, and a Box and Jenkins model, used only lagged air temperature residuals as exogenous variables. The root‐mean‐square error (RMSE) for these models varied between 0·53 and 1·70 °C and the second‐order autoregressive model provided the best results. A statistical methodology using best subsets regression is proposed to model the combined effect of discharge and air temperature on stream temperatures. Various streamflow indices were considered as additional independent variables, and models with different number of variables were tested. The results indicated that the best model included relative change in flow as the most important streamflow index. The RMSE for this model was of the order of 0·51 °C, which shows a small improvement over the first three models that did not include streamflow indices. The ridge regression was applied to this model to alleviate the potential statistical inadequacies associated with multicollinearity. The amplitude and sign of the ridge regression coefficients seem to be more in agreement with prior expectations (e.g. positive correlation between water temperature residuals of different lags) and make more physical sense. Copyright © 2006 John Wiley & Sons, Ltd.     

Seismic site effects by an optimized 2D BE/FE method I. Theory, numerical optimization and application to topographical irregularities

This paper deals with the evaluation of seismic site effects due to the local topographical and geotechnical characteristics. The amplification of surface motions is calculated by a numerical method combining finite elements in the near field and boundary elements in the far field (FEM/BEM). The numerical technique is improved by time truncation. In the first part of this article, the accuracy and the relevance of this optimized method are presented. Moreover, parametric studies are done on slopes, ridges and canyons to characterize topographical site effects. The second part deals with sedimentary valleys. The complexity of the combination of geometrical and sedimentary effects is underlined. Extensive parametrical studies are done to discriminate the topographical and geotechnical effects on seismic ground movement amplifications in two-dimensional irregular configurations. Characteristic coefficients are defined to predict the amplifications of horizontal displacements. The accuracy of this quantitative evaluation technique is tested and discussed.