Near-field and far-field simulation of accelerograms of Sikkim earthquake of September 18, 2011 using modified semi-empirical approach

The semi-empirical approach for modeling of strong ground motion given by Midorikawa (Tectonophysics 218:287?C295, 1993) has been modified in the present paper for component wise simulation of strong ground motion. The modified approach uses seismic moment in place of attenuation relation for scaling of acceleration envelope. Various strong motion properties like directivity effect and dependence of peak ground acceleration with respect to surface projection of source model have been studied in detail in the present work. Recently, Sikkim earthquake of magnitude 6.9 (M _w ) that occurred on September 18, 2011 has been recorded at various near-field and far-field strong motion stations. The modified semi-empirical technique has been used to confirm the location and parameters of rupture responsible for this earthquake. Strong motion record obtained from the iterative modeling of the rupture plane has been compared with available strong motion records from near as well as far-field stations in terms of root mean square error between observed and simulated records. Several possibilities of nucleation point, rupture velocity, and dip of rupture plane have been considered in the present work and records have been simulated at near-field stations. Final selection of model parameters is based on root mean square error of waveform comparison. Final model confirms southward propagating rupture. Simulations at three near-field and twelve far-field stations have been made using final model. Comparison of simulated and observed record has been made in terms of peak ground acceleration and response spectra at 5?% damping. Comparison of simulated and observed record suggests that the method is capable of simulating record which bears realistic appearance in terms of shape and strong motion parameters. Present work shows that this technique gives records which matches in a wide frequency range for Sikkim earthquake and that too from simple and easily accessible parameters of the rupture plane.     

Modeling of strong motion generation areas of the 2011 Tohoku,Japan earthquake using modified semi-empirical technique

Modification in the semi-empirical technique for the simulation of strong ground motion has been introduced to incorporate the strong motion generation areas (SMGA) in the modeled rupture plane. Strong motion generation areas identified within the rupture plane of the Tohoku earthquake of March 11, 2011 (M _w = 9.0), have been modeled using this modified technique. Two different source models having four and five SMGAs, respectively, are considered for modeling purpose. Strong motion records using modified semi-empirical technique have been simulated at two near-field stations located at epicentral distance of 137 and 140 km, respectively, using two different source models. Comparison of the observed and simulated acceleration waveforms is made in terms of root mean square error (RMSE) at both stations. Minimum root mean square error of the waveform comparison has been obtained at both the stations for source model having five SMGAs. Simulations from same rupture model have been made at other four stations lying at epicentral distance between 154 and 249 km. Comparison of observed and simulated records has been made in terms of RMSE in acceleration records, velocity records and response spectra at each six station. Simulations have been made at six other stations to obtain distribution of peak ground acceleration and peak ground velocity with hypocentral distance. Peak ground acceleration and velocity from simulated and observed records are compared at twelve stations surrounding the source of Tohoku earthquake. Comparison of waveforms and parameters extracted from observed and simulated strong motion records confirms the efficacy of the developed modified technique to model earthquake characterized by SMGAs.     

Flood risk assessment in a coastal lagoon under present and future scenarios: Ria de Aveiro case study

Iranian strong motion records as well as detailed conditions of their instrument sites and the characteristics of their causative seismic sources are compiled and processed. The dataset consists of 2286 three-component records from 461 Iranian earthquakes with at least two high-quality records having moment magnitude from 3.9 to 7.3. These records are about 20% of the Iranian database and are suitable for seismic hazard analysis and engineering applications. Perhaps for the first time in the literature, the distance to the surface projection of the fault is reported for a great number of records corresponding to earthquakes with M > 6.0. The raw accelerations are processed using the wavelet de-noising method. Having corrected and filtered these raw data, the pseudospectral accelerations are calculated for each of the three components of time series, separately. In addition to the ground motion parameters, a large and comprehensive list of metadata characterizing the recording conditions of each record is also developed. Moreover, careful revision of the characteristics of the earthquakes such as location, magnitude, style of faulting and fault rupture plane geometry, if available, is carried out using the best available information in a scientifically sound manner. Finally, we also focus on special ground motion records including records with peak ground acceleration (PGA) >300 cm/s² and distances less than 30 km. These are “exceptional” records in the Iranian dataset and include less than 2% of the selected dataset.     

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.     

Deterministic seismic scenario for Gujarat region, India

In this study, the modified stochastic method based on dynamic corner frequency has been used for the simulation of strong ground motions in Gujarat region. The earthquake-generating faults have been identified in the Gujarat region on the basis of past seismicity of the region. In all, 19 probable faults have been identified with 12 in Kachchh region, 5 in Saurashtra and 2 in Mainland Gujarat region. The maximum magnitude has been assigned to each fault based on the regional tectonic environment and past seismicity. The strong ground motions from these identified sources have been estimated at numerous points distributed all over Gujarat region on a grid. The peak ground acceleration (PGA) values have been extracted from the accelerograms and contoured. The spatial distribution of maximum of 19 PGA values at every grid point have been described and discussed. The ground motions at the surface of 32 important cities of the Gujarat have been estimated by incorporating the site amplification functions. The site amplification functions are obtained using the local earthquake data. These cities are located on various types of geological formations. We note that the site amplification functions have modified the character of the records and amplified the acceleration values at almost all the sites. The Kachchh region can expect surface accelerations between 400 and 800 cm/s², Saurashtra between 100 and 200 cm/s² and Mainland less than 50 cm/s² from a future large earthquake. The obtained results are useful for disaster mitigation measures, strengthening the existing built environment and design of structures in the region.     

A neural network approach for attenuation relationships: An application using strong ground motion data from Turkey

This paper presents an application of neural network approach for the prediction of peak ground acceleration (PGA) using the strong motion data from Turkey, as a soft computing technique to remove uncertainties in attenuation equations. A training algorithm based on the Fletcher–Reeves conjugate gradient back-propagation was developed and employed for three sample sets of strong ground motion. The input variables in the constructed artificial neural network (ANN) model were the magnitude, the source-to-site distance and the site conditions, and the output was the PGA. The generalization capability of ANN algorithms was tested with the same training data. To demonstrate the authenticity of this approach, the network predictions were compared with the ones from regressions for the corresponding attenuation equations. The results indicated that the fitting between the predicted PGA values by the networks and the observed ones yielded high correlation coefficients (R²). In addition, comparisons of the correlations by the ANN and the regression method showed that the ANN approach performed better than the regression. Even though the developed ANN models suffered from optimal configuration about the generalization capability, they can be conservatively used to well understand the influence of input parameters for the PGA predictions.     

Use of site amplification and anelastic attenuation for the determination of source parameters of the Sikkim earthquake of September 18, 2011, using far-field strong-motion data

An earthquake of magnitude 6.9 (M _w) occurred in the Sikkim region of India on September 18, 2011. This earthquake is recorded on strong-motion network in Uttarakhand Himalaya located about 900 km away from the epicenter of this earthquake. In this paper acceleration record from six far-field stations has been used to compute the source parameters of this earthquake. The acceleration spectra of ground motion at these far-field stations are strongly affected by both local site effects and near-site anelastic attenuation. In the present work the spectrum of S-phase recorded at these far-field stations has been corrected for anelastic attenuation at both source and site and the site amplification terms. Site amplifications at different stations and near-site shear wave attenuation factor have been computed by the technique of inversion of acceleration spectra given by Joshi et al. (Pure Appl Geophys 169:1821–1845, 2012a). For estimation of site amplification and shear wave quality factor [Q _β (f)] at the recording sites, ten local events recorded at various stations between July 2011 and December 2011 have been used. The obtained source spectrum from acceleration records is compared with the theoretical source spectrum defined by Brune (J Geophys Res 76:5002, 1970) at each station for both horizontal components of the records. Iterative forward modeling of theoretical source spectrum gives the average estimate of seismic moment (M _o), source radius (r _o) and stress drop (Δσ) as (3.2 ± 0.8) × 10²⁶ dyne cm, 13.3 ± 0.8 km and 59.2 ± 8.8 bars, respectively, for the Sikkim earthquake of September 18, 2011.     

Energy flux of strong earthquakes

In this paper, the energy flux of strong earthquakes at a station is determined considering the progressive rupture of a fault as the source of earthquakes. It is found that the motion of the source and the relative position of the station with respect to the fault are important in determining the energy density, the energy flux and the duration of the earthquake at this station. There is a “sphere of influence” beyond which the source may be assumed to be stationary. The analytical results are in good agreement with those of the 5 strong motion records obtained very near the fault from the Parkfield event of 27th June, 1966. 21 strong motion records are studied for energy densities at the stations from which a magnitude-energy relationship is obtained which agrees closely with other existing relationships.     

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.     

Site response zones and short-period earthquake ground motion projections for the Las Vegas Basin

A deterministic seismic hazard analysis was conducted to address the effect of local soil conditions on earthquake-induced strong ground motion in the Las Vegas Basin, Nevada (US). Using a large geological and geotechnical database, two response units were defined: a fine-grained unit, predominantly clay; and a coarse-grained unit, predominantly gravel. A moderate number of high-quality shallow shear wave velocity measurements were collected from which characteristic shear wave velocity profiles were developed for each response unit. An equivalent-linear one-dimensional site response model was used. The model was calibrated using a basin-wide, small-strain ground motion database. Calibration tests showed that ground motion projections become increasingly conservative with increasing ground-motion amplitude. Projections were overconservative for the coarse-grained response unit, likely due to the sparseness of the velocity database. For the earthquake response analyses, historical ground motions were used to model characteristic ‘bedrock’ motion for earthquakes on 10 faults judged to be critical. Response spectral envelopes were generated for each unit through Monte-Carlo simulations. For the fine-grained response unit, 95th percentile peak ground acceleration, peak spectral acceleration and predominant period were 310 cm/s², 1100 cm/s², and 0.29 s, respectively. With respect to codified design spectra, projections are lower at short periods and higher at long periods. Projections of peak spectral accelerations for the coarse-grained response unit, were more than double that of codified spectra; however, they are believed to be overconservative. Near-fault effects and basin-edge effects, though potentially important, were not considered in these analyses.     

City of Bingöl in May 2003: Assessment of strong ground motion records

The main shock of Bingöl earthquake (M _W = 6.4) recorded by six accelerometers in the area occurred at 03:27 local time on May 1, 2003. The largest acceleration value of north–south component was recorded as 545.5 cm/s² at the nearest station which it is 12 km away from the epicenter of earthquake. Especially, 0.15 s short period was observed when high spectral acceleration value occurred. An acceleration value greater than 50 gal was recorded at the BNG (Bingöl) station and structural damage occurred within 6.5 s was very important for the near source and strong ground motion seismology. The recorded peak acceleration values were greater than the estimated empirical acceleration values. However, the structural damage was not as high and widespread as expected. This occurrence was explained by considering the factors of earthquake source, frequency content, effective duration, effective acceleration value, local soil conditions, rupture direction and attenuation.     

Free field ground motions at various depths derived from actual seismograph records by using elastodynamic infinite elements

This paper presents a method for the analysis of free field ground motions at various depths from actual seismograph records. The proposed method replaces the propagating rupture on the fault surface by a fictitious focal point and a seismograph station in the vicinity of the given soil site. A fictitious focal point and seismograph station scheme is used to calibrate the free field ground motion of the soil site. The harmonic vibration analysis of the half space, which elastodynamic infinite elements are adopted in the far field and finite elements in the near field, yields the Fourier transform of any constituent wave in terms of the amplitude of the forces applied at the focal point, thus in terms of the Fourier transforms of the accelerations in three orthogonal directions recorded at a seismograph station. If a seismograph station is in the far field (as it is usually), the recorded values there can be transformed into the values at the infinite element node by means of appropriate shape functions. The results can be more refined in case of the availability of simultaneous records at more seismograph stations. In case of N stations, we assume N fictitious focal points, and solve a set of 3N simultaneous equations of 3N forces, three orthogonal forces at each focal point. Naturally, only the results by the present method were ‘exact’ at the calibrated points, and should be accurate in-between and near those points. The accelerations at the points near seismograph stations at various depths can also be obtained. An appropriate inverse Fourier transform algorithm will properly yields all results as time functions.     

Effect of frequency-dependent radiation pattern in the strong motion simulation of the 2011 Tohoku earthquake,Japan, using modified semi-empirical method

We perform a strong ground motion simulation using a modified semi-empirical technique (Midorikawa in Tectonophysics 218:287–295, 1993), with frequency-dependent radiation pattern model. Joshi et al. (Nat Hazards 71:587–609, 2014) have modified the semi-empirical technique to incorporate the modeling of strong motion generation areas (SMGAs). A frequency-dependent radiation pattern model is applied to simulate high-frequency ground motion more precisely. Identified SMGAs (Kurahashi and Irikura in Earth Planets Space 63:571–576, 2011) of the 2011 off the Pacific coast of Tohoku earthquake (M _w  = 9.0) were modeled using this modified technique. We analyzed the effect of changing seismic moment values of SMGAs on the simulated acceleration time series. Final selection of the moment values of SMGAs is based on the root-mean-square error (RMSE) of waveform comparison. Records are simulated for both frequency-dependent and constant radiation pattern function. Simulated records for both cases are compared with observed records in terms of peak ground acceleration, peak ground velocity and pseudo-acceleration response spectra at different stations. Comparison of simulated and observed records in terms of RMSE suggests that the method is capable of simulating record, which matches in a wide frequency range for this earthquake and bears realistic appearance in terms of shape and strong motion parameters. The results confirm the efficacy and suitability of rupture model defined by five SMGAs for the developed modified technique.     

Model parameter optimization for site-dependent simulation of ground motion by simulated annealing: Re-evaluation of the Ashigara Valley prediction experiment

More than 40 groups from 10 different countries participated in a weak- and strong-motion prediction experiment at Ashigara Valley which required the blind prediction of time series, spectra and spectral ratios for selected and instrumented sedimentary sites with well-known geotechnical properties. The wide scatter of the results of this experiment have raised a number of questions as to how to model high-frequency ground motion in the presence of available geotechnical and geophysical data. Using a simulated annealing waveform inversion method, we have tried to optimize and automate the model construction for ID site-dependent ground-motion simulation. We found a whole set of successfull models which provide good waveform fit (r > 0.8) for the observed displacement records at site KS2 but also yield sufficiently accurate response spectra and peak value predictions for both surface and downhole site. This shows that ID models are fully adequate to model the site conditions at least for the weak motion data. The resulting successful layer models consistently show a slower, less dense, and slightly thicker low velocity coverage with stronger damping than the official geotechnical model. Furthermore, their statistical properties directly measure the sensitivity of the individual parameters for the simulations. The critical re-evaluation of our own prediction which was based on stochastic simulation shows that although this approach has its greatest merits in situations where little information is available, it can also be successfully applied to model individual records if sufficient care is taken to determine the source parameters. Simulated annealing waveform inversion has shown to be a powerful tool to optimize that process.     

Estimation of model parameter of Sumatra earthquake using empirical Green’s function technique and generation of hypothetical earthquake scenario for Andaman Island, India

Empirical Green??s function (EGF) technique is considered to be most effective technique for simulation of ground motions due to a finite earthquake source. In the present paper, this technique has been used to simulate ground motion due to a great earthquake. The coastal region of Sumatra Island has been visited by a great earthquake on December 26, 2004. This earthquake has been recorded at several broadband stations including a nearest broadband station PSI in Indonesia. The shear wave contributions in both horizontal components have been simulated at PSI station using EGF technique. The comparison of simulated and observed waveform has been made for various possibilities of rupture parameters in terms of root mean square error. The final rupture model supports rupture velocity of 3.0?km/s with nucleation point supporting northward propagating rupture that coincide with high-slip asperity defined by Sorensen et al. (Bull Seism Soc Am 97:S139?CS151, 2007). The final modeling parameters have been used to simulate record at MDRS station in coastal state of Tamilnadu, India. In an attempt to model a scenario of great earthquake in the Andaman Island, a hypothetical rupture plane is modeled in this region. The event occurred on August 10, 2008 of magnitude 6.2 (M _w ) recorded on strong motion array at Port Blair has been used as EGF to simulate records due to the hypothetical great earthquake. Possibilities of earthquake due to the oblique strike-slip and thrust mechanism have been modeled in the present paper. Several possibilities of nucleation point for both cases has been considered, and it is seen that variation of peak ground acceleration at Port Blair station for strike-slip and thrust mechanism is 126?C738 gals and 647?C2,571 gals, respectively, which indicate high seismic hazard potential of Andaman Island.     

汶川地震紫坪铺面板堆石坝地震波输入研究

对国家强震动台网中心紫坪铺面板堆石坝区域台站实测主震记录以及大坝台网实测余震记录进行分析,研究主震与强余震地震动的基本特征。分别选取茂县地办、郫县走石山、成都中和这3组基岩台站实测主震地震动,紫坪铺台站2008年11月6日实测余震地震动以及按水工抗震规范人工生成地震动作为数值计算的地震动输入,对紫坪铺大坝进行三维动力有限元分析,并与实测结果进行对比。研究表明,郫县走石山与成都中和2个远场台站位于断层下盘,其实测地震动的加速度反应谱长周期（0.65 s以后）分量过于丰富,不宜作为断裂带附近紫坪铺大坝的地震动输入;紫坪铺大坝台站实测的余震地震动1 Hz附近（大坝基频）的频率成分相对较少,且持续时间较短,以至于难以激发大坝响应;对比坝顶实测地震动加速度反应谱和数值计算反应谱,建议汶川地震中紫坪铺大坝动力计算时可采用茂县地办台站实测地震动或按抗震设计规范反应谱人工生成地震动。     

Investigation of the variability of strong ground motions from Vrancea earthquakes

A systematic investigation of the applicability of several ground motion prediction models for Vrancea intermediate-depth seismic source is conducted in this research. Two ground motion prediction models recommended by previous evaluations (Vacareanu et al. in Bull Earthq Eng 11(6):1867–1884, 2013a; Pavel et al. in Earthq Struct 6(1):1–18, 2014), as well as two new state-of-the-art ground motion prediction equations (Vacareanu et al. in J Earthq Eng, 2013b; Earthq Struct 6(2):141–161, 2014) are tested using an increased strong ground motion database consisting of 150 recordings from Vrancea subcrustal earthquakes. The evaluation is performed by using several goodness-of-fit parameters from the literature. Moreover, the applicability of the single-station sigma method is also investigated by using the same strong ground motion database recorded in 30 seismic stations from southern and eastern Romania. The influence of the soil conditions on the numerical results obtained in this study is investigated and discussed using the results provided by the analysis of variance method. The impact of the single-station standard deviation on the levels of seismic hazard is also assessed in this study, and the results show, in the analyzed cases, significant reductions of the hazard levels.     

Seismic Site Classification of Recording Stations in Tarai Region of Uttarakhand,from Multiple Approaches

Repeated earthquakes (EQs) are clear indication of alarming seismicity which can be witnessed across Indian subcontinent. Increase in population density with inappropriate construction practice repeatedly rise alarm that in comparison to damage scenarios experienced during previous major to great EQs in India, future catastrophes would be manifold. Performing regional seismic hazard as well as site response studies can possibly help in accurate estimation of probable future seismic scenario. Site class (SC) of EQ recording stations is an important part of both seismic hazard as well as site response analyses. In seismic hazard analysis, suitable attenuation relations are often selected based on comparison of recorded ground motion with proposed ground motion as per selected attenuation relation for the same SC. Thus, unless SC of recorded ground motions is known, suitability of selected attenuation relation cannot be validated. In addition, recent studies suggest that for same soil column, ground motion may amplify at the surface from minimal to very high depending upon input motion characteristics. Thus again, unless SC of recording station is not known, recorded ground motion cannot be considered with confidence as outcrop or base motion for region specific site response studies. In the present work, SC of eight recording stations located in Tarai region of Uttarakhand, India located adjacent to the Himalayan belt and which are part of PESMOS database, are established by three different methods namely; equivalent linear ground response analysis, generalized inversion technique and horizontal to vertical spectral ratio method. Collectively all these three methods suggest same SC for each of the eight recording stations including Roorkee, Rishikesh, Dehradun etc. Further, obtained SC based on the present study is considerably different from available SC as per PESMOS database. However, present findings are matching with recent published work. Obtained results can be very helpful in developing surface seismic hazard using regional ground motion records towards minimizing future EQ induced damages.     

Seismic hazard evaluation for dams in northern Colorado, U.S.A.

A seismic hazard evaluation for three dams in the Rocky Mountains of northern Colorado is based on a study of the historical seismicity. To model earthquake occurrence as a random process utilizing a maximum likelihood method, the catalog must exhibit random space-time characteristics. This was achieved using a declustering procedure and correction for completeness of recording. On the basis of the resulting a- and b-values, probabilistic epicentral distances for a 2 × 10^–5 annual probability were calculated. For a random earthquake of magnitude M _L 6.0–6.5, this distance is 15 km. Suggested ground motion parameters were estimated using a probabilistic seismic hazard analysis. Critical peak horizontal accelerations at the dams are 0.22g if median values are assumed and 0.39g if variable attenuation and seismicity rates are taken into account. For structural analysis of the dams, synthetic acceleration time series were calculated to match the empirical response spectra. In addition, existing horizontal strong motion records from two Mammoth Lakes, California earthquakes were selected and scaled to fit the target horizontal acceleration response spectra.     

芦山地震冶勒大坝强震监测资料分析

冶勒沥青混凝土心墙堆石坝最大坝高为124.5 m,坝址区地震烈度高,地质条件复杂,两岸坝基条件严重不对称。大坝上布设了9台强震仪组成的强震监测台阵,曾获得2008年汶川地震和攀枝花地震的大坝强震监测记录。2013年4月20日四川省雅安市芦山县发生里氏7.0级地震,冶勒大坝距震中约212.5 km,坝址区震感较为强烈,强震监测台阵获得了此次地震较为完整的有效记录。对芦山地震主震记录进行时域分析和频谱分析,总结冶勒大坝在芦山地震中的动力反应规律,并与汶川地震时坝体动力反应进行对比分析。研究表明,芦山地震主震时冶勒大坝最大加速度记录为47.043 cm/s2,最长持续时间为76.98 s,坝顶动力放大效应明显;芦山和汶川地震时大坝动力反应规律的差异与地震波频谱特性及大坝自振特性等密切相关。总体而言,冶勒大坝在震后运行安全稳定,芦山地震未对冶勒大坝造成明显不利影响。