Hybrid attenuation model for estimation of peak ground accelerations in the Kutch region, India

The Kutch region of Gujarat in India is the locale of one of the most devastating earthquake of magnitude (M _w) 7.7, which occurred on January 26, 2001. Though, the region is considered as seismically active region, very few strong motion records are available in this region. First part of this paper uses available data of strong motion earthquakes recorded in this region between 2006 and 2008 years to prepare attenuation relation. The developed attenuation relation is further used to prepare synthetic strong motion records of large magnitude earthquakes using semiempirical simulation technique. Semiempirical simulation technique uses attenuation relation to simulate strong ground motion records of any target earthquake. The database of peak ground acceleration obtained from simulated records is used together with database of peak ground acceleration obtained from observed record to develop following hybrid attenuation model of wide applicability in the Kutch region: $$ \begin{aligned} \ln \left( {\text{PGA}} \right) & = - 2.56 + 1.17 \, M_{\text{w}} - \, 0.015R - 0.0001\ln \left( {E + 15} \right) \\ &\quad 3.0 \le M_{\text{w}} \le 8.2;\quad 12 \le R \le 120;\quad {\text{std}} . {\text{ dev}}.(\sigma ): \pm 0.5 \\ \end{aligned} $$ ln ( PGA ) = ? 2.56 + 1.17 M w ? 0.015 R ? 0.0001 ln ( E + 15 ) 3.0 ≤ M w ≤ 8.2 ; 12 ≤ R ≤ 120 ; std . dev . ( σ ) : ± 0.5 In the above equation, PGA is maximum horizontal ground acceleration in gal, M _w is moment magnitude of earthquake, R is hypocentral distance, and E is epicentral distance in km. The standard deviation of residual of error in this relation is 0.5. This relation is compared with other available relations in this region, and it is seen that developed relation gives minimum root mean square error in comparison with observed and calculated peak ground acceleration from same data set. The applicability of developed relation is further checked by testing it with the observed peak ground acceleration from earthquakes of magnitude (M _w), 3.6, 4.0, 4.4, and 7.7, respectively, which are not included in the database used for regression analysis. The comparison demonstrates the efficacy of developed hybrid attenuation model for calculating peak ground acceleration values in the Kutch region.     

Attenuation relations of peak ground acceleration and velocity in the Southern Dead Sea Transform region

Using the recorded earthquake strong ground motion, the attenuation of peak ground acceleration (PGA) and peak ground velocity (PGV) are derived in the southern Dead Sea Transform region. The expected values of strong motion parameters from future earthquakes are estimated from attenuation equations, which are determined by regression analysis on real accelerograms. In this study, the method of Joyner and Boor [Bull Seismol Soc Am 71(6):2011–2038, 1981] was selected to produce the attenuation model for the southern Dead Sea Transform region. The dataset for PGA consists of 57 recordings from 30 earthquakes and for PGV 26 recordings from 19 earthquakes. The attenuation relations developed in this study are proposed as replacement for former probabilistic relations that have been used for a variety of earthquake engineering applications. The comparison between the derived PGA relations from this study with the former relations clearly shows significant lower values than the other relations.     

Strong ground motion attenuation in Aswan area, Egypt

Recent and paleo seismicity indicate that moderate seismic activity is relatively large for Aswan area. This is a warning on the possibility of occurrence of earthquakes in the future too. No strong motion records are available in Aswan area for engineers to rely upon. Consequently, the seismological modeling is an alternative approach till sufficient instrumental records around Aswan become available. In the present study, we have developed new ground motion attenuation relationship for events spanning 4.0?≤? M _w?≤?7.0 and distance to the surface projection of the fault up to 100 km for Aswan based on a statistically simulated seismological model. We generated suites of ground motion time histories using stochastic technique. The ground motion attenuation relation describes the dependence of the strength of the ground motions on the earthquake magnitude and distance from the earthquake. The proposed equation for peak ground acceleration (PGA) for the bed rock is in the form of: $ {\mathbf{log}}{\text{ }}\left( {{\mathbf{PGA}}/{\mathbf{gal}}} \right){\text{ }} = {\mathbf{1}}.{\mathbf{24}} + {\mathbf{0}}.{\mathbf{358}}{M_{\mathbf{w}}} - {\text{ }}{\mathbf{log}}\left( {\mathbf{R}} \right){\text{ }}-{\text{ }}{\mathbf{0}}.{\mathbf{008}}{\text{ }}{\mathbf{R}}{\text{ }} + {\text{ }}{\mathbf{0}}.{\mathbf{22}}{\text{ }}{\mathbf{P}} $ . Where PGA is the peak ground acceleration in gal (cm/s²); M_w, its moment magnitude; R is the closest distance between the rupture projection and the site of interest; and the factor P is a dummy variable. It is observed that attenuation of strong motion in Aswan is correlated with those used before in Egypt.     

Robust attenuation relations for peak time-domain parameters of strong ground motions

This study presents new attenuation models for the estimation of peak ground acceleration (PGA), peak ground velocity (PGV), and peak ground displacement (PGD) using a hybrid method coupling genetic programming and simulated annealing, called GP/SA. The PGA, PGV, and PGD were formulated in terms of earthquake magnitude, earthquake source to site distance, average shear-wave velocity, and faulting mechanisms. A worldwide database of strong ground motions released by Pacific Earthquake Engineering Research Center (PEER) was employed to establish the models. A traditional genetic programming analysis was performed to benchmark the proposed models. For more validity verification, the GP/SA models were employed to predict the ground-motion parameters of the Iranian plateau earthquakes. Sensitivity and parametric analyses were carried out and discussed. The results show that the GP/SA attenuation models can offer precise and efficient solutions for the prediction of estimates of the peak time-domain characteristics of strong ground motions. The performance of the proposed models is better than or comparable with the attenuation relationships found in the literature.     

Modeling of strong ground motions for 1991 Uttarkashi, 1999 Chamoli earthquakes, and a hypothetical great earthquake in Garhwal–Kumaun Himalaya

In this study, stochastic finite fault modeling is used to simulate Uttarkashi (1991) and Chamoli (1999) earthquakes using all available source, path, and site parameters available for the region. These two moderate earthquakes are recorded at number of stations of a strong motion network. The predicted peak ground accelerations at these stations are compared with the observed data and the ground motion parameters are constrained. The stress drop of Uttarkashi and Chamoli earthquakes is constrained at 77 and 65?bars, respectively, whereas the quality factor Q _C is 112 $ f^{0.97} $ and 149 $ f^{0.95} $ for these two regions. The high-frequency attenuation parameter Kappa is in the range 0.04?C0.05. The constrained ground motion parameters are then used to simulate Mw 8.5 earthquake in central seismic gap region of Himalaya. Two scenarios are considered with epicenter of future great earthquake at locations of Uttarkashi and Chamoli earthquakes using above constrained parameters. The most vulnerable towns are the towns of Dehradun and Almora where expected PGA is in excess of 600?cm/s² at V_S30 520?m/s when the epicenter of the great earthquake is at the location of Uttarkashi (1991) earthquake. The towns of Shimla and Chandigarh can expect PGA close to 200?cm/s². Whereas when the epicenter of the great earthquake is at the location of Chamoli (1999) earthquake, the towns of Dehradun and Almora can expect PGA of around 500 and 400?cm/s², respectively, at V_S30 620?m/s. The National Capital Region, Delhi can expect accelerations of around 80?cm/s² in both the cases. The PGA contour maps obtained in this study can be used to assess the seismic hazard of the region and identify vulnerable areas in and around central Himalaya from a future great earthquake.     

Generation of isoseismal maps for Mexico using velocity and acceleration data

Acceleration and velocity data of large earthquakes recorded by the Mexican National Broadband Seismological Network were used to generate isoseismal maps for Mexico. The seismic data consist of 99 events recorded at 27 seismic stations located in the southern and northern regions of Mexico from 2004 to 2009. The magnitude (Mw) of these events ranged from 4.1 to 7.1. Peak ground velocity values (PGV) and peak ground accelerations (PGA) were estimated, and velocity-derived peak ground accelerations (PGAv) were calculated. No important differences between PGAv and PGA values were found; thus, both parameters were used in the isoseismal determination. The generated synthetic isoseismal maps were compared with those from an existing catalog for large earthquakes in Mexico. Using empiric relations between PGA and MMI (modified Mercalli intensity) and PGV and MMI, the obtained simulated intensity maps showed similar behavior to those reported in the catalog. The results indicate that the PGAv values can be used to determine intensities when acceleration records (PGA) are unavailable.     

Calcite dissolution kinetics in Na-Ca-Mg-Cl brines

Sedimentary basins can contain close to 20% by volume of pore fluids commonly classified as brines. These fluids can become undersaturated with respect to calcite as a result of migration, dispersive mixing, or anthropogenic injection of CO₂. This study measured calcite dissolution rates in geologically relevant Na-Ca-Mg-Cl synthetic brines (50-200 g L⁻¹ TDS). The dissolution rate dependency on brine composition, pCO₂ (0.1-1 bar), and temperature (25.0-82.5 °C) was modeled using the empirical rate equation

R=k(1-Ω)ⁿ     

Site response of the Ganges basin inferred from re-evaluated macroseismic observations from the 1897 Shillong, 1905 Kangra,and 1934 Nepal earthquakes

We analyze previously published geodetic data and intensity values for the M _s = 8.1 Shillong (1897), M _s = 7.8 Kangra (1905), and M _s = 8.2 Nepal/Bihar (1934) earthquakes to investigate the rupture zones of these earthquakes as well as the amplification of ground motions throughout the Punjab, Ganges and Brahmaputra valleys. For each earthquake we subtract the observed MSK intensities from a synthetic intensity derived from an inferred planar rupture model of the earthquake, combined with an attenuation function derived from instrumentally recorded earthquakes. The resulting residuals are contoured to identify regions of anomalous intensity caused primarily by local site effects. Observations indicative of liquefaction are treated separately from other indications of shaking severity lest they inflate inferred residual shaking estimates. Despite this precaution we find that intensites are 1–3 units higher near the major rivers, as well as at the edges of the Ganges basin. We find evidence for a post-critical Moho reflection from the 1897 and 1905 earthquakes that raises intensities 1–2 units at distances of the order of 150 km from the rupture zone, and we find that the 1905 earthquake triggered a substantial subsequent earthquake at Dehra Dun, at a distance of approximately 150 km. Four or more M = 8 earthquakes are apparently overdue in the region based on seismic moment summation in the past 500 years. Results from the current study permit anticipated intensities in these future earthquakes to be refined to incorporate site effects derived from dense macroseismic data.     

Thermodynamic studies of pyrrhotite-pyrite equilibria in the Ag-Fe-S system by solid-state galvanic cell technique at 518-723 K and total pressure of 1 atm

The reaction FeS₂(cr) + 2Ag(cr) = ‘FeS’(cr) + Ag₂S(cr) was studied by measuring the temperature dependence of the electromotive force (EMF) of the all-solid-state galvanic cell with common gas space:

(-)Pt|Ag|AgI|Ag₂S,‘FeS’,FeS₂|Pt(+)     

Numerical evaluation of seismic response of shallow foundation on loose silt and silty sand

This study includes the results of a set of numerical simulations carried out for sands containing plastic/non-plastic fines, and silts with relative densities of approximately 30?40% under different surcharges on the shallow foundation using FLAC 2D. Each model was subjected to three ground motion events, obtained by scaling the amplitude of the El Centro (1940), Kobe (1995) and Kocaeli (1999) Q12earthquakes. Dynamic behaviour of loose deposits underlying shallow foundations is evaluated through fully coupled nonlinear effective stress dynamic analyses. Effects of nonlinear soil structure interaction (SSI) were also considered by using interface elements. This parametric study evaluates the effects of soil type, structure weight, liquefiable soil layer thickness, event parameters (e.g., moment magnitude of earthquake (M _w ), peak ground acceleration PGA, PGV/PGA ratio and the duration of strong motion (D _5?95) and their interactions on the seismic responses. Investigation on the effects of these parameters and their complex interactions can be a valuable tool to gain new insights for improved seismic design and construction.     

Near-source strong motion database catalog for Iran

This paper discusses a newly developed high-quality integrated dataset of shallow earthquake ground motions that occurred in Iran, from 1976 to 2013. A total of 860 three-component strong motion records are processed from 183 earthquake events, moment magnitudes 5.0?≤?M _w ?≤?7.4, and rupture distances of R _RUP  ≤ 120 km. Strong motion data from Iran having special tectonic features and shallow earthquakes with depths less than 35 km are included. This paper presents a thorough procedure used to collect and to generate a database following the Next-Generation Attenuation-West research projects. This database can be used in the development and ranking of ground motion models and for seismological and engineering hazard and risk analyses. Unprocessed strong motion records are obtained from the Iranian Strong Motion Network (ISMN). The time series collected were thoroughly examined through several rounds of quality reviews. The newly generated database includes the peak ground acceleration, peak ground velocity, and pseudo-spectral acceleration for the 5% damped with periods ranging from 0.01 to 10 s. The database also includes ground motion information and source characterization and parameters. This study is the near-source compiled ground motion database that can be used for Iran, and it is consistent with standard worldwide databases.     

Zinc complexation in aqueous sulfide solutions: Determination of the stoichiometry and stability of complexes via ZnS(cr) solubility measurements at 100 °C and 150 bars

The solubility of ZnS_(cr) was measured at 100 °C, 150 bars in sulfide solutions as a function of sulfur concentration (m(S_total) = 0.02-0.15) and acidity (pH_t = 2-11). The experiments were conducted using a Ti flow-through hydrothermal reactor enabling the sampling of large volumes of solutions at experimental conditions, with the subsequent concentration and determination of trace quantities of Zn. Prior to the experiments, a long-term in situ conditioning of the solid phase was performed in order to attain the reproducible Zn concentrations (i.e. solubilities). The ZnS_(cr) solubility product was monitored in the course of the experiment. The following species were found to account for Zn speciation in solution: Zn²⁺ (pH_t < 3), (pH_t 3-4.5), (pH_t 5-8), and ZnS(HS)⁻ (pH_t > 8) (pH_t predominance regions are given for m(S_total) = 0.1). Solubility data collected in this study at pH_t > 3 were combined with the ZnS_(cr) solubility product determined at lower pH to yield the following equilibrium constants (t = 100 °C, P = 150 bars):

     

Sr/Ca and Ca/Ca fractionation during inorganic calcite formation: I. Sr incorporation

Partitioning of strontium during spontaneous calcite formation was experimentally studied using an advanced CO₂-diffusion technique. Results at different precipitation rates and T = 5, 25, and 40 °C show that at constant temperature Sr incorporation into calcite is controlled by the precipitation rate (R in μmol/m²/h) according to the individual expressions

     

Probabilistic seismic hazard zonation of Syria

Earthquake hazard maps for Syria are presented in this paper. The Peak Ground Acceleration (PGA) and the Modified Mercalli Intensity (MMI) on bedrock, both with 90% probability of not being exceeded during a life time of 50, 100 and 200 years, respectively are developed. The probabilistic PGA and MMI values are evaluated assuming linear sources (faults) as potential sources of future earthquakes. A new attenuation relationship for this region is developed. Ten distinctive faults of potential earthquakes are identified in and around Syria. The pertinent parameters of each fault, such as theb-parameter in the Gutenberg-Richter formula, the annual rate ₄ and the upper bound magnitudem ₁ are determined from two sets of seismic data: the historical earthquakes and the instrumentally recorded earthquake data (AD 1900–1992). The seismic hazard maps developed are intended for preliminary analysis of new designs and seismic check of existing civil engineering structures.     

Effect of pH and organic ligands on the kinetics of smectite dissolution at 25 °C

Forward dissolution rates of Na-Montmorillonite (Wyoming) SWy-2 smectite (Ca_0.06Na_0.56)[Al_3.08Fe(III)_0.38Mg_0.54] [Si_7.93 Al_0.07]O₂₀(OH)₄ were measured at 25 °C in a mixed-flow reactor equipped with interior dialysis compartment (6-8 kDa membrane) as a function of pH (1-12), dissolved carbonate (0.5-10 mM), phosphate (10⁻⁵ to 0.03 M), and nine organic ligands (acetate, oxalate, citrate, EDTA, alginate, glucuronic acid, 3,4-dihydroxybenzoic acid, gluconate, and glucosamine) in the concentration range from 10⁻⁵ to 0.03 M. In organic-free solutions, the Si-based rates decrease with increasing pH at 1 ? pH ? 8 with a slope close to −0.2. At 9 ? pH ? 12, the Si-based rates increase with a slope of ∼0.3. In contrast, non-stoichiometric Mg release weakly depends on pH at 1 ? pH ? 12 and decreases with increasing pH. The empirical expression describing Si-release rates [R, mol/cm²/s] obtained in the present study at 25 °C, I = 0.01 M is given by

     

Precipitation and mixing properties of the “disordered” (Mn,Ca)CO3 solid solution

The enthalpy of mixing of the calcite-rhodochrosite (Ca,Mn)CO₃ solid solution was determined at 25 °C from calorimetric measurements of the enthalpy of precipitation of solids with different compositions. A detailed study of the broadening of powder X-ray diffraction peaks shows that most of the precipitates are compositionally homogeneous. All the experimental enthalpy of mixing (ΔH^m) values are positive and fit reasonably well (R² = 0.86) to a Guggenheim function of three terms:

     

基于复合方法的天水盆地宽频带地震动模拟

本文采用有限差分和随机振动合成结合的复合方法,模拟了当礼县—罗家堡断裂发生矩震级Mw7.7级大地震时,在天水盆地产生的宽频带地震动场,分析了在设定地震条件下盆地内的地震动分布特征,为该区黄土地震滑坡分析提供了地震动参数结果。结果显示:(1)有限差分法和随机振动合成法可以很好地互补,得到盆地内地表宽频带地震动;(2)地震在盆地区域产生了强烈地震动,PGA(峰值加速度)介于150~900 gal,离断层较近的区域东南角的PGA最大,随着断层距的增加,PGA逐渐减小。河谷南侧的PGA值相比北侧较大,具备诱发滑坡的强大动力条件;(3)盆地区域PGV(峰值速度)最大为120 cm/s。受第四系覆盖层放大效应和地形放大效应共同影响,水平向地震动在盆地区域东侧和中部具有较大PGV,而西侧PGV相对较小。竖向地震动在盆地区域东侧较弱,而在中部和西侧较强,特别是最西侧陡峭的山坡上,PGV达到了最大值。此外,竖向地震动明显受到覆盖层厚度的影响,譬如在盆地区域南侧的中间部位,也具有较大的PGV。     

Seawater nutrient and carbonate ion concentrations recorded as P/Ca, Ba/Ca, and U/Ca in the deep-sea coral Desmophyllum dianthus

As paleoceanographic archives, deep sea coral skeletons offer the potential for high temporal resolution and precise absolute dating, but have not been fully investigated for geochemical reconstructions of past ocean conditions. Here we assess the utility of skeletal P/Ca, Ba/Ca and U/Ca in the deep sea coral D. dianthus as proxies of dissolved phosphate (remineralized at shallow depths), dissolved barium (trace element with silicate-type distribution) and carbonate ion concentrations, respectively. Measurements of these proxies in globally distributed D. dianthus specimens show clear dependence on corresponding seawater properties. Linear regression fits of mean coral Element/Ca ratios against seawater properties yield the equations: P/Ca_coral (μmol/mol) = (0.6 ± 0.1) P/Ca_sw(μmol/mol) - (23 ± 18), R² = 0.6, n = 16 and Ba/Ca_coral(μmol/mol) = (1.4 ± 0.3) Ba/Ca_sw(μmol/mol) + (0 ± 2), R² = 0.6, n = 17; no significant relationship is observed between the residuals of each regression and seawater temperature, salinity, pressure, pH or carbonate ion concentrations, suggesting that these variables were not significant secondary dependencies of these proxies. Four D. dianthus specimens growing at locations with Ω_arag ? 0.6 displayed markedly depleted P/Ca compared to the regression based on the remaining samples, a behavior attributed to an undersaturation effect. These corals were excluded from the calibration. Coral U/Ca correlates with seawater carbonate ion: U/Ca_coral(μmol/mol) = (−0.016 ± 0.003) (μmol/kg) + (3.2 ± 0.3), R² = 0.6, n = 17. The residuals of the U/Ca calibration are not significantly related to temperature, salinity, or pressure. Scatter about the linear calibration lines is attributed to imperfect spatial-temporal matches between the selected globally distributed specimens and available water column chemical data, and potentially to unresolved additional effects. The uncertainties of these initial proxy calibration regressions predict that dissolved phosphate could be reconstructed to ±0.4 μmol/kg (for 1.3-1.9 μmol/kg phosphate), and dissolved Ba to ±19 nmol/kg (for 41-82 nmol/kg Ba_sw). Carbonate ion concentration derived from U/Ca has an uncertainty of ±31μmol/kg (for ). The effect of microskeletal variability on P/Ca, Ba/Ca, and U/Ca was also assessed, with emphasis on centers of calcification, Fe-Mn phases, and external contaminants. Overall, the results show strong potential for reconstructing aspects of water mass mixing and biogeochemical processes in intermediate and deep waters using fossil deep-sea corals.     

Mechanisms of equilibrium and kinetic oxygen isotope effects in synthetic aragonite at 25 °C

Aragonite was precipitated in the laboratory at 25 °C in isotopic equilibrium with Na-Ca-Mg-Cl-CO₃ solutions at two different pH values (i.e., pH = ∼8.2 and ∼10.8) by the constant addition method. On the basis of the oxygen isotope composition of the aragonite precipitates, it was demonstrated that the equilibrium aragonite-water fractionation factor is independent of the pH of the parent solution and equal to:

1000lnα_{(aragonite-H2O)}=29.12±0.09     

Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar

Here we report on an experimental investigation of the relation between the dissolution rate of albite feldspar and the Gibbs free energy of reaction, ΔG_r. The experiments were carried out in a continuously stirred flow-through reactor at 150 °C and pH_(150 °C) 9.2. The dissolution rates R are based on steady-state Si and Al concentrations and sample mass loss. The overall relation between ΔG_r and R was determined over a free energy range of −150 < ΔG_r < −15.6 kJ mol⁻¹. The data define a continuous and highly non-linear, sigmoidal relation between R and ΔG_r that is characterized by three distinct free energy regions. The region furthest from equilibrium, delimited by −150 < ΔG_r < −70 kJ mol⁻¹, represents an extensive dissolution rate plateau with an average rate . In this free energy range the rates of dissolution are constant and independent of ΔG_r, as well as [Si] and [Al]. The free energy range delimited by −70 ? ΔG_r ? −25 kJ mol⁻¹, referred to as the ‘transition equilibrium’ region, is characterized by a sharp decrease in dissolution rates with increasing ΔG_r, indicating a very strong inverse dependence of the rates on free energy. Dissolution nearest equilibrium, defined by ΔG_r > −25 kJ mol⁻¹, represents the ‘near equilibrium’ region where the rates decrease as chemical equilibrium is approached, but with a much weaker dependence on ΔG_r. The lowest rate measured in this study, R = 6.2 × 10⁻¹¹ mol m⁻² s⁻¹ at ΔG_r = −16.3 kJ mol⁻¹, is more than two orders of magnitude slower than the plateau rate. The data have been fitted to a rate equation (adapted from Burch et al. [Burch, T. E., Nagy, K. L., Lasaga, A. C., 1993. Free energy dependence of albite dissolution kinetics at 80 °C and pH 8.8. Chem. Geol.105, 137-162]) that represents the sum of two parallel reactions

R=k₁[1-exp(-ng^m₁)]+k₂[1-exp(-g)]^m₂,