Support vector machine for evaluating seismic-liquefaction potential using shear wave velocity

The use of the shear wave velocity data as a field index for evaluating the liquefaction potential of sands is receiving increased attention because both shear wave velocity and liquefaction resistance are similarly influenced by many of the same factors such as void ratio, state of stress, stress history and geologic age. In this paper, the potential of support vector machine (SVM) based classification approach has been used to assess the liquefaction potential from actual shear wave velocity data. In this approach, an approximate implementation of a structural risk minimization (SRM) induction principle is done, which aims at minimizing a bound on the generalization error of a model rather than minimizing only the mean square error over the data set. Here SVM has been used as a classification tool to predict liquefaction potential of a soil based on shear wave velocity. The dataset consists the information of soil characteristics such as effective vertical stress (sigma'(v0)), soil type, shear wave velocity (V-s) and earthquake parameters such as peak horizontal acceleration (a(max)) and earthquake magnitude (M). Out of the available 186 datasets, 130 are considered for training and remaining 56 are used for testing the model. The study indicated that SVM can successfully model the complex relationship between seismic parameters, soil parameters and the liquefaction potential. In the model based on soil characteristics, the input parameters used are sigma'(v0), soil type. V-s, a(max) and M. In the other model based on shear wave velocity alone uses V-s, a(max) and M as input parameters. In this paper, it has been demonstrated that Vs alone can be used to predict the liquefaction potential of a soil using a support vector machine model. (C) 2010 Elsevier B.V. All rights reserved.

Exact free-surface flows for shallow-water equations .1. - the incompressible case

A comprehensive exact treatment of free surface flows governed by shallow water equations (in sigma variables) is given. Several new families of exact solutions of the governing PDEs are found and are shown to embed the well-known self-similar or traveling wave solutions which themselves are governed by reduced ODEs. The classes of solutions found here are explicit in contrast to those found earlier in an implicit form. The height of the free surface for each family of solutions is found explicitly. For the traveling or simple wave, the free surface is governed by a nonlinear wave equation, but is arbitrary otherwise. For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed; in another case, the free surface is a horizontal plane while the flow underneath is a sine wave. The existence of simple waves on shear flows is analytically proved. The interaction of large amplitude progressive waves with shear flow is also studied.

Field investigations of the 30 September 1993 earthquake in Maharashtra

An engineering field study of the 30 September Marathwada earthquake is reported. The study covered Osmanabad, Latur, Sholapur, Bijapur, Gulbarga and Bidar districts. The level of ground acceleration was estimated based on tilting of free standing objects. The study shows that the epicentral intensity has been VIII on the UNESCO scale. The horizontal ground acceleration near the epicentre has been estimated to be about 0.2 g.

Seepage effects on sand-bed channels

ExperimentS were conducted to find the effects of seepage on flow over a sand bed in a straight rectangular flume under two conditions: (1) When the channel bed is plane, horizontal, and nontransporting; and (2) when the bed is transporting at a constant sediment concentration. Effects of both injection and suction, caused by seepage flow into and out of the channel bed, are studied for condition 1; and only suction effects are studied for condition 2, Three sands, sizes 0.34 mm, 0.53 mm, and 0.80 mm, are used in the study. It is found that seepage can cause an increase or decrease in the bed shear stress relative to no seepage for the two conditions. The change in bed shear stress depends on the relative magnitudes of the bed shear stress and the critical shear stress of particles under the no-seepage condition, sediment concentration, and the seepage rate. Quantitative relationships giving the ratio of bed shear stresses with and without seepage are presented for both conditions of the bed. A procedure to estimate the changes in bed shear stress, friction factor, Manning's n, and stream power due to seepage for known initial conditions of the channel and the amount of applied seepage is presented

Natural Convection Heat Transfer From an Isothermal Vertical Surface to a Stable Thermally Stratified Fluid

Natural convection from an isothermal vertical surface to a thermally stratified fluid is studied numerically. A wide range of stratification levels is considered. It is shown that at high levels of ambient thermal stratification, a portion at the top of the plate absorbs heat, while a horizontal plume forms around a location where the plate temperature equals the ambient temperature. The plume is shown to be inherently unsteady, and its transient nature is investigated in detail. The effect of the temperature defect in striating the plume is discussed. Average Nusselt number data are presented for Pr = 6.0 and 0.7.

Exact free-surface flows for shallow-water equations .2. the compressible case

Exact free surface flows with shear in a compressible barotropic medium are found, extending the authors' earlier work for the incompressible medium. The barotropic medium is of finite extent in the vertical direction, while it is infinite in the horizontal direction. The ''shallow water'' equations for a compressible barotropic medium, subject to boundary conditions at the free surface and at the bottom, are solved in terms of double psi-series, Simple wave and time-dependent solutions are found; for the former the free surface is of arbitrary shape while for the latter it is a damping traveling wave in the horizontal direction, For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed. In the case of an isothermal medium, when gamma = 1, we again find simple wave and time-dependent solutions.

Simulation of ENSO-Related Surface Winds in the Tropical Pacific by an Atmospheric General Circulation Model Forced by Observed Sea Surface Temperatures

The authors present the simulation of the tropical Pacific surface wind variability by a low-resolution (R15 horizontal resolution and 18 vertical levels) version of the Center for Ocean-Land-Atmosphere Interactions, Maryland, general circulation model (GCM) when forced by observed global sea surface temperature. The authors have examined the monthly mean surface winds acid precipitation simulated by the model that was integrated from January 1979 to March 1992. Analyses of the climatological annual cycle and interannual variability over the Pacific are presented. The annual means of the simulated zonal and meridional winds agree well with observations. The only appreciable difference is in the region of strong trade winds where the simulated zonal winds are about 15%-20% weaker than observed, The amplitude of the annual harmonics are weaker than observed over the intertropical convergence zone and the South Pacific convergence zone regions. The amplitudes of the interannual variation of the simulated zonal and meridional winds are close to those of the observed variation. The first few dominant empirical orthogonal functions (EOF) of the simulated, as well as the observed, monthly mean winds are found to contain a targe amount of high-frequency intraseasonal variations, While the statistical properties of the high-frequency modes, such as their amplitude and geographical locations, agree with observations, their detailed time evolution does not. When the data are subjected to a 5-month running-mean filter, the first two dominant EOFs of the simulated winds representing the low-frequency EI Nino-Southern Oscillation fluctuations compare quite well with observations. However, the location of the center of the westerly anomalies associated with the warm episodes is simulated about 15 degrees west of the observed locations. The model simulates well the progress of the westerly anomalies toward the eastern Pacific during the evolution of a warm event. The simulated equatorial wind anomalies are comparable in magnitude to the observed anomalies. An intercomparison of the simulation of the interannual variability by a few other GCMs with comparable resolution is also presented. The success in simulation of the large-scale low-frequency part of the tropical surface winds by the atmospheric GCM seems to be related to the model's ability to simulate the large-scale low-frequency part of the precipitation. Good correspondence between the simulated precipitation and the highly reflective cloud anomalies is seen in the first two EOFs of the 5-month running means. Moreover, the strong correlation found between the simulated precipitation and the simulated winds in the first two principal components indicates the primary role of model precipitation in driving the surface winds. The surface winds simulated by a linear model forced by the GCM-simulated precipitation show good resemblance to the GCM-simulated winds in the equatorial region. This result supports the recent findings that the large-scale part of the tropical surface winds is primarily linear.

Geometrically Simple Logarithmic Weir

This paper discusses the design and experimental verification of a geometrically simple logarithmic weir. The weir consists of an inward trapezoidal weir of slope 1 horizontal to n vertical, or 1 in n, over two sectors of a circle of radius R and depth d, separated by a distance 2t. The weir parameters are optimized using a numerical optimization algorithm. The discharge through this weir is proportional to the logarithm of head measured above a fixed reference plane for all heads in the range 0.23R less than or equal to h less than or equal to 3.65R within a maximum deviation of +/-2% from the theoretical discharge. Experiments with two weirs show excellent agreement with the theory by giving a constant average coefficient of discharge of 0.62. The application of this weir to the field of irrigation, environmental, and chemical engineering is highlighted.

Profiling of late Tertiary-early Quaternary surface in the lower reaches of Narmada valley using microtremors

In this paper, we propose the first approximation for thickness of Quaternary sediment and late Quaternary early Tertiary topography for the part of lower reaches of Narmada valley in a systematic way using the shallow seismic method, that records both horizontal and vertical components of the microtremor (ambient noise) caused by natural processes. The measurements of microtremors were carried out at 31 sites spaced at a grid interval of 5 km s using Lennartz seismometer (5 s period) and City shark-II data acquisition system. The signals recorded were analysed for horizontal to the vertical (H/V) spectral ratio using GEOPSY software. For the present study, we concentrate on frequency range between 0.2 Hz and 10 Hz. The thickness of unconsolidated sediments at various sites is calculated based on non-linear regression equations proposed by Ibs-von Seht and Wohlenberg (1999) and Parolai et al. (2002). The estimated thickness is used to plot digital elevation model and cross profiles correlating with geomorphology and geology of the study area. (C) 2011 Elsevier Ltd. All rights reserved.

Airborne measurements of aerosol scattering properties above the MABL over Bay of Bengal during W_ICARB - characteristics and spatial gradients

Realizing the importance of aerosol characterization and addressing its spatio-temporal heterogeneities over Bay of Bengal (BoB), campaign mode observations of aerosol parameters were carried out using simultaneous cruise, aircraft and land-based measurements during the Winter Integrated Campaign for Aerosols gases and Radiation Budget (W_ICARB). Under this, airborne measurements of total and hemispheric backscatter coefficients were made over several regions of coastal India and eastern BoB using a three wavelength integrating nephelometer. The measurements include high resolution multi-level (ML) sorties for altitude profiles and bi-level (BL) sorties for spatial gradients within and above the Marine Atmospheric Boundary Layer (MABL) over BoB. The vertical profiles of the scattering coefficients are investigated in light of the information on the vertical structure of the atmospheric stability, which was derived from the collocated GPS (Global Positioning System) aided radiosonde ascents. In general, the altitude profiles revealed that the scattering coefficient remained steady in the convectively well-mixed regions and dropped off above the MABL. This decrease was quite rapid off the Indian mainland, while it was more gradual in the eastern BoB. Investigation on horizontal gradients revealed that the scattering coefficients over northern BoB are 3 to 4 times higher compared to that of central BoB within and above the MABL. A north-south gradient in scattering coefficients is observed over Port Blair in the eastern BoB, with values decreasing from south to north, which is attributed to the similar gradient in the surface wind speed, which can be replicated in the sea salt abundance. The gradients are parameterized using best-fit analytical functions.

Levitation of a drop over a film flow

A vertical jet of water impinging on a horizontal surface produces a radial film flow followed by a circular hydraulic jump. We report a phenomenon where fairly large (1 mi) drops of liquid levitate just upstream of the jump on a thin air layer between the drop and the film flow. We explain the phenomenon using lubrication theory. Bearing action both in the air film and the water film seems to be necessary to support large drops. Horizontal support is given to the drop by the hydraulic jump. A variety of drop shapes is observed depending on the volume of the drop and liquid properties. We show that interaction of the forces due to gravity, surface tension, viscosity and inertia produces these various shapes.

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

Seismic bearing capacity factors for spread foundations

The effect of horizontal earthquake body forces on the bearing capacity of foundations has been examined computationally in a rigorous manner by employing the method of stress characteristics. The bearing capacity factors N-c, N-q and N-y, due to the components of soil cohesion, ground surcharge pressure and soil unit weight respectively, have been plotted as a function of earthquake acceleration coefficient (a(h)) for different values of soil friction angle (phi). The inclusion of earthquake body forces causes a considerable reduction in the bearing capacity factors. The bearing capacity factors N-c and N-q are seen to be approximately of the same magnitude as those reported in the literature on the basis of different solution methods. However, the obtained values of N-y are found to be significantly smaller than the available results. The nature of the pressure distribution along the footing base and the geometry of the observed failure patterns vary with the consideration of earthquake body forces.

Seismic passive earth pressure coefficients using the method of characteristics

The method of characteristics was used to generate passive earth pressure coefficients for an inclined wall retaining cohesionless backfill material in the presence of pseudostatic horizontal earthquake body forces. The variation of the passive earth pressure coefficients K-pq and K-pgamma with changes in horizontal earthquake acceleration coefficient due to the components of soil unit weight and surcharge pressure, respectively, has been obtained; a closed-form solution for K-pq is also provided. The passive earth resistance has been found to decrease sharply with an increase in the magnitude of horizontal earthquake acceleration. The computed passive earth pressure coefficients were found to be the lowest when compared to all of the previous solutions available in the literature.

Seismic passive resistance in soils for negative wall friction

: In the presence of pseudo-static seismic forces, passive earth pressure coefficients behind retaining walls were generated using the limit equilibrium method of analysis for the negative wall friction angle case (i.e., the wall moves upwards relative to the backfill) with logarithmic spirals as rupture surfaces. Individual density, surcharge, and cohesion components were computed to obtain the total minimum seismic passive resistance in soils by adding together the individual minimum components. The effect of variation in wall batter angle, ground slope, wall friction angle, soil friction angle, and horizontal and vertical seismic accelerations on seismic passive earth pressures are considered in the analysis. The seismic passive earth pressure coefficients are found to be highly sensitive to the seismic acceleration coefficients both in the horizontal and the vertical directions. The results are presented in graphical and tabular formats.