Estimation of seismic hazard in Gujarat region, India

The seismic hazard in the Gujarat region has been evaluated. The scenario hazard maps showing the spatial distribution of various parameters like peak ground acceleration, characteristics site frequency and spectral acceleration for different periods have been presented. These parameters have been extracted from the simulated earthquake strong ground motions. The expected damage to buildings from future large earthquakes in Gujarat region has been estimated. It has been observed that the seismic hazard of Kachchh region is more in comparison with Saurashtra and mainland. All the cities of Kachchh can expect peak acceleration in excess of 500?cm/s² at surface in case of future large earthquakes from major faults in Kachchh region. The cities of Saurashtra can expect accelerations of less than 200?cm/s² at surface. The mainland Gujarat is having the lowest seismic hazard as compared with other two regions of Gujarat. The expected accelerations are less than 50?cm/s² at most of the places. The single- and double-story buildings in Kachchh region are at highest risk as they can expect large accelerations corresponding to natural periods of such small structures. Such structures are relatively safe in mainland region. The buildings of 3?C4 stories and tall structures that exist mostly in cities of Saurashtra and mainland can expect accelerations in excess of 100?cm/s² during a large earthquake in Kachchh region. It has been found that a total of 0.11 million buildings in Rajkot taluka of Saurashtra are vulnerable to total damage. In Kachchh region, 0.37 million buildings are vulnerable. Most vulnerable talukas are Bhuj, Anjar, Rapar, Bhachau, and Mandvi in Kachchh district and Rajkot, Junagadh, Jamnagar, Surendernagar and Porbandar in Saurashtra. In mainland region, buildings in Bharuch taluka are more vulnerable due to proximity to active Narmada-Son geo-fracture. The scenario hazard maps presented in this study for moderate as well as large earthquakes in the region may be used to augment the information available in the probabilistic seismic hazard maps of the region.     

Estimating the vulnerability and loss functions of residential buildings

This paper is concerned with an investigation of the damage to residential buildings in two areas within Gilan and Zanjan provinces, Iran, caused by the Manjil-Rudbar earthquake of 20 June 1990. A statistical correlation between the observed ground motion and the damage to the residential buildings is derived for overall damaged buildings and expressed as the vulnerability function. The loss function is calculated by combining the seismic hazard with the vulnerability function.The study of vulnerability and annual seismic hazard shows that the specific annual risk for the range of motion of 0.18 to 0.5 g is equal to 0.02. This indicates that the specific risk for semi-engineered residential buildings with a lifetime of 20 years is about 33%. This study also shows that in large cities, such as Tehran, located in seismic areas, the extent of damage according to the vulnerability function will be 45 and 70% for expected maximum accelerations of 0.3 and 0.4 g, respectively.     

Mapping seismic risk: the current crisis

The seismic risks to which populations are exposed should be estimated reliably for mitigation and preparation of response to disastrous earthquakes. Three parameters need to be known: Population numbers, properties of the built environment, and the seismic hazard. If we focus on large cities, we can say that at least one of these is known satisfactorily, namely the population, but not the other two. In the developing world, the numbers of buildings in a city are known only approximately, their distribution into building types (resistance to shaking) has to be assumed, and the distribution of types throughout the city is unknown. Recent verification of the world seismic hazard map has shown that it is grossly misleading: Instrumental measurements of accelerations due to six earthquakes were about three times larger, on average, than the maximum likely accelerations shown on the map; the macroseismic intensities reported for the last 60 earthquakes with M ≥ 7.5 were all significantly larger than expected, based on the hazard map (by 2.3 intensity units for the 12 deadliest earthquakes); and calculations of losses of life based on the hazard map underestimate the losses sustained in the 12 recent earthquakes with more than 1,000 fatalities by two to three orders of magnitude. This means that the seismic risk in most of the approximately 1,000 large cities at risk in the developing world is unknown. To remedy this intolerable situation, models for the built environment in cities need to be constructed, using cost-effective analyses of satellite images, and worst case scenario estimates of the losses in case of the nearest maximum credible earthquake.     

Assessment of seismic hazards for Hurghada,Red Sea,Egypt

The entrance of the southern Suez Gulf of the Red Sea is known to be an area of high seismic activity in Egypt. The high rate of seismic activity in this area is mainly related to the adjustment in motion at the triple junction between the African plate, the Arabian plate, and the Sinai microplate. The present study attempts to estimate the Probabilistic Seismic Hazard Analysis (PSHA) for Hurghada site. This was done in two steps; the first one is by estimating specific parameters for the site, such as the mean seismic activity, λ, the Gutenberg-Richter parameter, b, and the maximum regional magnitude, m _max. The second step is by selecting a ground motion parameter that is applicable to Hurghada site. The procedure permits the combination of both historical and recent instrumental data. The results of the hazard assessment, expressed as the worst case scenario, detect that Hurghada is exposed to the maximum credible earthquake event of magnitude m _max = 7.1 ± 0.31, at hypocentral distance of 31.6 ± 10 km. The possibility of the maximum Peak Ground Acceleration (PGA), which occurred in relation to this event at Hurghada site, is equal to 0.29 g. The mean return periods with the selected accelerations for Hurghada, a horizontal acceleration of 0.1 g, is expected to occur once every 74–106 years, with an average of one every 90 years. This result which obtained from the hazard assessment can be used as an input data for a seismic risk assessment.     

Use of backpropagation neural network for landslide monitoring: a case study in the higher Himalaya

Static and dynamic rock slope stability analyses were performed using a numerical discontinuum modelling technique for a 700-m high rock slope in western Norway. The rock slope has been investigated by the Geological Survey of Norway (NGU), which has been carrying out rock slide studies for the county Møre and Romsdal in western Norway. The purpose of numerical modelling was to estimate the volume of the rock mass that could potentially slide under static and dynamic forces. This estimation was required to assess the run-up heights (tsunami) in a fjord that could potentially be caused by the rockslide. Three cases have been simulated for predicting the behaviour of the rock slope. First, an initial static loading is applied in the numerical model to simulate the prevailing rock mass conditions at the site. Second, saturated and weathered joint conditions are modelled by reducing the residual friction angle along the discontinuities of the rock mass. In doing so, the model simulates the effect of degradation of discontinuities in the rock slope. Third, a dynamic loading, based on peak ground accelerations expected in the area, is applied to simulate dynamic earthquake conditions.

These numerical studies have provided some useful insights into the deformation mechanisms in the rock slope. Both sliding and rotation of blocks start to occur once the residual friction angle along the discontinuities is reduced and when the region is shaken by a strong earthquake. The results indicate that, due to variations in the inclination of discontinuities, the entire slope does not become unstable and that down-slope sliding and rotation of blocks occur mainly on the top layers of the slope. Within the range of parameter values considered for this study, it is unlikely that the whole rock slope can be destabilised. The study provides an illustration of how the geo-mechanical properties of a rock mass can be integrated in a discontinuum rock slope model, which is used for predicting the behaviour of the slope under existing environmental and earthquake conditions. This model has helped not only to better understand the dynamics of the rockslide but also to estimate the potential rock volume that can become unstable when subjected to static and dynamic loads.     

Probabilistic seismic hazard analysis for Kakrapar atomic power station, Gujarat, India

Seismic ground motion caused by earthquakes mainly affects the constructions and structures around its area of influence. In this context, the probabilistic seismic hazard analysis (PSHA) is a scientific step towards the safety analysis of any major construction such as nuclear power plant. Thus, the present study focused to estimate seismic hazard level at different probabilities for Kakrapar nuclear power plant located in the Western India. The hazard curves for the study area are developed following the procedure of PSHA suggested by Cornell–McGuire. Three source zones, Narmada-Tapti zone (NTZ), Rann of Kuchchh (ROK), and west passive margin (WPM), are classified on the basis of seismicity and tectonic setting of the study area. The estimated maximum magnitude (m _max) for NTZ, ROK, and WPM are 6.9 ± 0.57, 6.5 ± 0.64, and 6.1 ± 0.64, respectively. Logic tree approach has been used for the development of hazard curves to account the epistemic uncertainties associated with the analysis. For maximum credible earthquake [MCE, i.e., the probability of exceedance of 2 % in 50 years (return period of ~2,500 years)], the peak spectral acceleration (i.e., PSA at 0.2 s) expected around 5 km of the Kakrapar nuclear power plant (site) is 0.23 g from all source zones; however, at exact site location, it is 0.18 g. The PSA values due to NTZ, ROK, and WPM based on MCE are 0.22, 0.065, and 0.052 g, respectively. In case of design-based earthquake (DBE, i.e., 50 % probability in 50 years (return period of ~110 years)), the calculated maximum spectral acceleration (SA) from all source zones is about 0.045 g. The PSA distribution for the DBE from the NTZ has reached a maximum value of 0.042 g; however, PSA for ROK and WPM is considerably low with a maximum value of 0.022 and 0.021 g, respectively. Considering the MCE and DBE, the estimated PSA at 0.2 s has a highest value of ~0.23 g from all source zones. Spectral accelerations (SAs) correspond to different periods are presented, and SA plots for NTZ zone can be considered as response spectra for the KAPS site. Deaggregation of PSHA in the present study is also discussed. PGA values reported in seismic zonation map and global seismic hazard analysis program around the present study area range from 0.05 to 0.2 g which is slightly lower than the peak acceleration obtained in this study. The results of this study would facilitate in the performance of the site-specific seismic probabilistic safety analysis.     

Regional seismic hazard study for the eastern Mediterranean (Trans-Jordan,Levant and Antakia) and Sinai region

The aim of this study is to assess the seismic hazard in the eastern Mediterranean and Sinai region using a probabilistic approach. An updated earthquake catalogue for the period 1 to 1993 AD that covers the area between latitude 27°–37°N and longitude 32°–39°E, has been used. Using the new seismic-tectonic map for the area, 10 line-sources are delineated. These lines or fault zones are thought to represent the main sources for the seismic potential in the area. The results are demonstrated as iso-contour lines of the peak-ground acceleration. The iso-acceleration contours represent 90% probability that these peak values will not be exceeded over periods of 50, 100 and 200 years, respectively. This study concludes that the seismic hazard severity is highest along the Jordan Dead Sea transform fault system, namely from south of the Gulf of Aqaba, Dead Sea-Jordan River, Tiberia Lake, Rachaya, Ed Damur, Yammuneh Fault, and Ghab Fault in the north. For the 50 year iso-contour map, the major cities of Amman, Damascus, and Beirut lay around the 2 m s⁻² contour line, while Jerusalem lies along the 3 m s⁻² line. Antakia in Turkey has the highest seismic potential severity (around 5 m s⁻²) while in Cyprus the maximum hazard is expected to reach 4 m s⁻² for the coming 50 years.     

Structure of the earth's crust in Jordan from potential field data

Gravity and magnetic data were collected and used to study the crustal structure of Jordan. Three new geophysical maps of Jordan were created: a Moho discontinuity map, a crystalline basement surface map, and a map showing the lowest limit of magnetic blocks. Depths of the Curie Isotherm were also calculated. Results indicate that the depth to the Moho discontinuity in Jordan varies from 32 to 33 km in the northwest to 38 km in the southeast. The basement complex rocks outcrop on the surface in the southwest but lie at about 8 km in the northeast. The Curie Isotherm (585 °C) lies at a depth of about 10 km in the area east of the Dead Sea and dips southeastward towards the Al-Sirhan (Wadi Sirhan), southeast Jordan, where it is located at 35 km depth. Local isostasy of rock masses (blocks) in Jordan does not occur. Nevertheless, this does not rule out the possible existence of isostasy in a regional scale at greater depths within the mantle.     

Blast Induced Damage and Dynamic Behaviour of Hangingwalls in Bench Stoping

This paper presents the results of a comprehensive monitoring program designed to investigate the extent of blast induced damage experienced by rock masses extracted by bench stoping methods. An array of triaxial geophones and extensometers were used to monitor blast vibration attenuation and measure hangingwall deformations during stope extraction. In addition, pre and post surveys of the hangingwall rock mass were conducted using a TV borehole camera and cavity survey instrumentation. These surveys were later used to calibrate damage profiles into the stope hangingwalls.

Peak particle velocity, hangingwall deformation measurements and stope surveys were used to develop a site specific damage model that allowed engineers to asses drilling and blasting configurations to minimise the extent of pre-conditioning and damage. In addition the study included the analysis of the frequency response, displacements and accelerations experienced by the excavation as extraction and mine filling progressed. This work aimed at improving our understanding of the influence of blasting on the dynamic behaviour of stope hangingwalls.

The study demonstrated that estimates of the maximum extent of rock mass pre-conditioning and/or damage made through the application of the Holmberg-Persson approach compared well with measured results. In addition, the study found that dynamic loading imparted on an exposed hangingwall from subsequent stope blasting was also expected to contribute to rock mass weakening and that mine filling was crucial to arrest further deterioration. Hangingwall accelerations were used to demonstrate that larger openings may be more susceptible to dynamic loading.     

Flow and sediment dynamics in a low sinuosity, braided river: Calamus River, Nebraska Sandhills

The interaction between channel geometry, flow, sediment transport and deposition associated with a midstream island was studied in a braided to meandering reach of the Calamus River, Nebraska Sandhills. Hydraulic and sediment transport measurements were made over a large discharge range using equipment operated from catwalk bridges. The relatively low sinuosity channel on the right-hand side of the island carries over 70% of the water discharge at high flow stages and 50–60% at low flow stages. As a result, mean velocity, depth, bed shear stress and sediment transport rate tend to be greater here than in the more strongly curved left-hand channel. The loci of maximum flow velocity, depth and bed shear stress are near the centre of the channel upstream of the island, but then split and move towards the outer banks of both channels downstream. Variations in these loci depend on the flow stage. Topographically induced across-stream flows are generally stronger than the weak, curvature-induced secondary circulations. Water surface topography is controlled mainly by centrifugal accelerations and local changes in downstream flow velocity. The averaged water surface slope of the study reach varies very little with discharge, having values between 0·00075 and 0·00090. As bed shear stress generally varies in a similar way to mean velocity, friction coefficients vary little, normally being in the range 0·07–0·13. These values are similar to those in straight channels with sandy dune-covered beds. Bedload is moved mainly as dunes at all flow stages. Grain size is mainly medium sand with coarse sand moved in thalwegs adjacent to the cut banks, and with fine sand at the downstream end of the island. These patterns of flow velocity, depth, water surface topography, bed shear stress, bedload transport rate and mean grain size can be accurately predicted using theoretical models of flow, bed topography and sediment transport rate in single river bends, applied separately to the left and right channels. During high flow stages deposition occurs persistently near the downstream end of the island, and cut banks are eroded. Otherwise, erosion and deposition occurs only locally within the channel as discharge varies. Abandonment and filling of a strongly curved channel segment may occur by migration of an upstream bar into the channel entrance at a high flow stage.     

Geochemical investigation of phosphorus and nitrogen in the hypersaline Dead Sea

Dissolved and particulate phosphorus, ammonia, and nitrate concentrations were measured in the Dead Sea water column (330 m depth, salinity ca. 340 g/L), in the Lower Jordan River and in springs along its western coast. Dissolved phosphate in the water column is fairly homogeneous, at about 35 μg P/L. Particulate P shows larger variability, 30–50 μg/L. The main inputs of dissolved phosphate into the lake are diffusion from bottom sediments (58%) and the Jordan River inflow (41%). Input from springs is negligible (1%). Biological uptake is a minor removal pathway because in the present Dead Sea, primary production takes place only when major floods occur and dilute the upper layer, about once in 10 years. When this happens, only about 200 ton P, out of a whole-lake reservoir of ca. 5000 ton P, are removed from the biologically active upper layer of about 7 m. Chemical removal pathways, not yet investigated, such as coprecipitation with authigenic aragonite, could be of importance. The average ammonia concentration in the water column has gradually increased from 5.9 mg/L in 1960 to 8.9 mg/L, in 1991. Diffusion from the bottom sediments is a major source of ammonia to the Dead Sea. The annual input from freshwater inflows has been of lesser importance in the 1960s and 1970s. A pollution-derived increase in the ammonia concentration of the Jordan River in the 1980s may partially explain the concomitant rise in the ammonia load of the lake. It is also suggested that following the 1979 overturn, and the yearly turnovers of the 1980s, ammonia might have been produced within the oxygenated water column by mineralization. Nitrate concentration was very low, 20 μg N/L, in the 1960s and increased to 0.2 to 0.5 mg/L in 1981. This increase is shown to be due to human pollution of the Jordan River. We conclude that nutrient concentration in the Dead Sea water column is controlled by physical and chemical factors, whereas biological effects are minimal.     

Modelling groundwater over-extraction in the southern Jordan Valley with scarce data

To deal with the challenge of groundwater over-extraction in arid and semi-arid environments, it is necessary to establish management strategies based on the knowledge of hydrogeological conditions, which can be difficult in places where hydrogeological data are dispersed, scarce or present potential misinformation. Groundwater levels in the southern Jordan Valley (Jordan) have decreased drastically in the last three decades, caused by over-extraction of groundwater for irrigation purposes. This study presents a local, two-dimensional and transient numerical groundwater model, using MODFLOW, to characterise the groundwater system and the water balance in the southern Jordan Valley. Furthermore, scenarios are simulated regarding hydrological conditions and management options, like extension of arable land and closure of illegal wells, influencing the projection of groundwater extraction. A limited dataset, literature values, field surveys, and the ‘crop water-requirement method’ are combined to determine boundary conditions, aquifer parameters, and sources and sinks. The model results show good agreement between predicted and observed values; groundwater-level contours agree with the conceptual model and expected flow direction, and, in terms of water balance, flow volumes are in accordance with literature values. Average annual water consumption for irrigation is estimated to be 29 million m³ and simulation results show that a reduction of groundwater pumping by 40% could recover groundwater heads, reducing the water taken from storage. This study presents an example of how to develop a local numerical groundwater model to support management strategies under the condition of data scarcity.     

Site-specific probabilistic seismic hazard analysis for northern part of the Qeshm Island,Iran

This paper presents results of a site-specific probabilistic seismic hazard analysis for northern part of the Qeshm Island, one the most seismic prone areas of Iran. Seismotectonic and seismicity properties of seismic sources in the study area were characterized and used for evaluation of various strong ground motion parameters implementing the classical Cornell’s PSHA approach. The results show that peak rock accelerations for 475-year return period are 0.4 and 0.27 g, respectively, for 84th and 50th percentiles while being about 0.37 and 0.61 g for 2475-year return period. These values are slightly smaller than those read from national seismic zonation maps which can be attributed to the considered conservatism for development of such design maps. In order to incorporate local site conditions, a series of dynamic site response analyses based on the equivalent linear approach were also employed. The results indicate that the presence of soft subsurface deposits at the site significantly alters the fundamental characteristics of the response spectra. The obtained median (50th percentile) peak ground accelerations for 975-year return period range between 0.49 and 0.54 g at different locations in the study site showing minor amplifications relative to their corresponding bedrock acceleration of 0.48 g. Finally, the obtained site-specific spectrum was compared with the standard spectrum mandated by the design codes. In this regard, the agreement was found to be reasonable at period ranges shorter than about 0.5 s, while the differences were more obvious at longer periods. This reveals the need for implementation of site-specific design spectrum to avoid underestimation or overestimation of seismic forces for designing critically important structures especially when softer deposits are encountered.     

Seismic response of pile foundations in liquefiable soil: parametric study

The performance of pile foundations in liquefiable soil subjected to earthquake loading is a very complex process. The strength and stiffness of the soil decrease due to the increase in pore pressure. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. This paper presents the response of vertical piles in liquefiable soil under seismic loads. A finite difference model, known as fast Lagrangian analysis of continua, is used to study the pile behavior considering a nonlinear constitutive model for soil liquefaction and pile?Csoil interaction. The maximum lateral displacement and maximum pile bending moment are obtained for different pile diameters, earthquake predominant frequencies, Arias intensities, and peak accelerations. It is found that the maximum lateral displacement and the maximum pile bending moment increase when the predominant earthquake frequency value decreases for a given peak acceleration value.     

Recharge of groundwater through multi-stage reservoirs in a desert basin

Artificial recharge of groundwater through geologic formations is of interest in Jordan as it is a proven water resource management technique that increases the efficiency of water use. This paper documents the research conducted into designing a series of small recharge dams on three connected valleys. With 50 potential sites previously screened, the focus was on the final selected location in the Jordanian desert, in which detailed topographical and hydrogeological surveys were used to determine the optimal locations for construction of surface water retention structures. The research produced designs for an integrated system of seven 2.5-m dikes and one 5-m dam with a total potential recharge volume of 500,000 m³. Once constructed, these multi-stage reservoirs will conserve part of the water being lost and will provide additional groundwater for future use. The research fills a gap in the lack of standard designs that can be expanded and replicated in other areas in Jordan and elsewhere.     

Dynamic Response of Sandy Slope Under Coupling of Earthquake and Groundwater

Based on dynamic elastic–plastic time history analysis method, dynamic responses of sandy slope in different groundwater levels are studied under near field earthquakes and far field earthquakes. The results show that the accelerations of the slope decrease with the increasing of groundwater levels, which demonstrates that the existence of groundwater damps vibration. The accelerations have different laws along the height of the slope under far field and near field earthquakes, however, the maximum values all appear at the top of the slope. The maximum horizontal displacement of the slope increases with the increasing of groundwater levels. The influence of far field earthquake on the deformation of slope toe is greater, while the influence of near field earthquake on the deformation of the slope top is greater. Thus, the top and the toe of the slope should be protected as the key position under earthquakes. The influence of high groundwater level on stability of the slope is more serious, and the safety factor calculated by pseudo-static method in slope seismic code will be lower for the slope in the high groundwater level.     

The Major Gulf of the Aqaba Earthquake, 22 November 1995 – Maximum Intensity Distribution

The major earthquake of 22 November 1995, with a moment magnitude M_W = 7.1 and a local magnitude M_L 6.2, was the beginning of a seismic swarm that occurred in the central part of the Gulf of Aqaba. During this swarm, thousands of earthquakes occurred with local magnitude ranging between 2 and 6.2 from 22 November 1996 to 31 December 1997, when 2089 earthquakes were detected and/or analyzed by the Jordan Seismological Observatory (JSO). The major earthquake strongly affected the near shoreline cities (Figure 1). The maximum observed intensity on these cities was VIII on the modified Mercalli intensity (MMI) scale. A questionnaire was distributed in the main cities of Jordan one week after the major earthquake. The results of this investigation, which demonstrated the observed intensity distribution for Aqaba city, shows a relationship between local conditions, such as geological foundations and topographical characters, and the extent of the destruction. This conclusion was supported by the maximum peak ground acceleration (PGA) measurements inside Aqaba and Eilat cities. From the results of this questionnaire compiled inside Jordan, and other reports and readings compiled from other nearby countries outside Jordan, a preliminary regional iso-intensity map for this major earthquake of 22 November 1995 is presented in this study.     

A simplified approach for studying the influence of vertical accelerations on seismic response of slopes

ABSTRACT

A simplified approach is presented for estimating permanent displacements in slopes as a result of both vertical and horizontal seismic accelerations. A study of 52 earthquake records showed that the time difference between maximum horizontal and vertical accelerations varied between 0 and 10.3 s. The approach is illustrated for an earth dam embankment by analysing the effects of five of the above earthquake records. The approach combines a pseudo-static slope stability analysis for estimation of the critical (or yield) horizontal-vertical acceleration combinations, and a Newmark type displacement analysis. Guidelines are presented for conservative choice of soil strength parameters of saturated clays for use in the stability analysis. While permanent displacements of up to 40 cm were predicted without considering the vertical acceleration component, no additional displacement above 3.5 cm resulted when this component was included. The predicted additional displacement was consistently less than 10%, and in 50% of the analyses, vertical acceleration led to smaller predicted displacements. The simple approach may be applied in analysis for any slope using real earthquake records. Using existing, empirical expressions for permanent displacement, based only on horizontal accelerations, the effect of the vertical accelerations may be conservatively estimated by increasing the displacement by 10%.     

Geochemical studies on ozokerite from the Dead Sea area

Ozokerite, a natural mixture of long-chain paraffins, is found at Ein-Humar (Jordan) on the eastern escarpment of the Dead Sea. The Ein-Humar ozokerite is mainly composed of saturated n-paraffins from C₃₀ to presumably well into the C_50′s with a maximum at C₃₈. IR, NMR and chromatographic data indicate the absence of cycloparaffinic and aromatic molecules. The ozokerite is assumed to be related to the heavy asphalts and oils in the Dead Sea area. It probably forms by separation of high-molecular paraffins from non-paraffinaceous heavy oils or asphalts by geothermal heat and precipitation in the conduits due to lowering of temperature.     

Effect of coupling excess pore pressure and deformation on nonlinear seismic soil response

The excess pore pressure (\(\Delta p_w\)) generation and consequent reduction in effective stress lead to the softening of a liquefiable soil deposit that can alter ground motions in terms of amplitude, frequency content and duration. However, total stress models, which are the most currently used, do not take into account coupling of excess pore pressures and soil deformations. To assess this effect, two analyses were made: (1) a Biot hydraulic and mechanical computation of a saturated soil deposit with coupling pore pressures and soil deformations and (2) a mechanical computation of a decoupled model with same initial behaviour. Both analyses were performed with a fully nonlinear elastoplastic multi-mechanism model. As \(\Delta p_w\) depends on the soil properties, two soils were analysed: loose-to-medium and medium-to-dense sand. The results regarding the profile of maximum accelerations and shear strains, the surface accelerations and their corresponding response spectra are analysed. The mean values of the normalized response spectra ratio of surface accelerations between the coupled and decoupled model show a deamplification of low and high frequencies (i.e. at frequencies lower than 1.0 Hz and higher than 10 Hz) that tend to increase with the liquefaction zone size. Coupling of \(\Delta p_w\) and soil deformation is therefore of great importance to accurately model the ground motion response. On the contrary, while peak acceleration predictions could be conservative, the amplification on the low frequencies could be largely underestimated which could be highly prejudicial for flexible buildings.