Structures exhibit coupled torsional and translational responses to earthquake ground motion input if their centres of floor mass and their centres of resistance do not coincide. However, torsional motions may occur even in nominally symmetric structures due to accidental eccentricity and torsional ground motions. The sources giving rise to accidental eccentricity include the difference between the assumed and actual distributions of mass and stiffness, asymmetric yielding strength, non-linear patterns of force-deformation relationships, and differences in coupling of the structural foundation with the supporting soil. Symmetric and regular buildings that are properly designed have a much higher ability to survive a strong earthquake event than asymmetric buildings and their response to earthquake loading is far more straightforward to predict and design for. On the other hand, even though the response of asymmetric buildings is more unpredictable, designers still have to compromise structural regularity to accommodate functional and aesthetic needs. As a result, serious and widespread damage associated with structural asymmetry has been observed repeatedly in past major earthquakes. In the first studies examining torsional effects in buildings, attention was focused on the elastic structural behaviour of single-storey buildings and the main purpose was to achieve a complete understanding of the effects of mass and stiffness eccentricities and to evaluate them by simple static models. However, as the response of real structures is mainly inelastic, these studies gave poor information on torsional behaviour and interest has moved towards non-linear response studies. In an effort to clarify some of the issues influencing the inelastic torsional response of multistorey asymmetric structures, this thesis presents a series of coherent parametric investigations. These investigations include comparing the response of various reference models to the performance of code-designed torsionally unbalanced structures. An extensive parametric investigation of torsionally responding structures designed as stipulated by a selection of major earthquake building codes is presented and the adequacy of the static torsional provisions is assessed for a wide range of structural configurations and parameters. Detailed investigations of torsionally asymmetric structures incorporating frame elements oriented along both orthogonal axes of the structure are also conducted and the effect of including the second earthquake component to simultaneously excite the structural models is quantified. The relative merits and deficiencies of each code provision are discussed and a new proposed optimised method is tested. All fundamental conclusions from the investigations conducted are presented and various topics for further research are proposed, which are considered to be both necessary and pertinent for increasing and refining the knowledge and understanding the complex behaviour of multistorey torsionally asymmetric buildings.

The main focus of this thesis is to effectively estimate levels of peak ground acceleration (PGA) during a seismic event for a given site. This will be achieved by applying regression analysis via a mixed model methodology to data collected from previously recorded seismic events collected from the PEER Strong Motion Database using a program written in MATLAB. The basic mixed model combines both fixed and random effect terms. Two models are analyzed and compared based on varying combinations of predictor variables, such as magnitude, distance, shear wave velocity, and site class. While the primary objective of this thesis solely examines the modeling of PGA, the same methodology can be applied in predicting other ground motion intensity parameters such as Peak Ground Velocity (PGV) or the spectral ordinate at a given vibration period.

Properly accounting for the uncertainties in predicting ground motion parameters is critical for Probabilistic Seismic Hazard Analyses (PSHA). This is particularly important for critical facilities that are designed for long return period motions. Non-ergodic PSHA is a framework that allows for this proper accounting of uncertainties. This, in turn, allows for more informed decisions by designers, owners and regulating agencies. The ergodic assumption implies that the standard deviation applicable to a specific source-path-site combination is equal to the standard deviation estimated using a database with multiple source-path-site combinations. The removal of the ergodic assumption requires dense instrumental networks operating in seismically active zones so that a sufficient number of recordings are made. Only recently, with the advent of networks such as the Japanese KiK-net network has this become possible. This study contributes to the state of the art in earthquake engineering and engineering seismology in general and in non-ergodic seismic hazard analysis in particular. The study is divided in for parts. First, an automated protocol was developed and implemented to process a large database of strong ground motions for GMPE development. A comparison was conducted between the common records in the database processed within this study and other studies. The comparison showed the viability of using the automated algorithm to process strong ground motions. On the other hand, the automated algorithm resulted in narrower usable frequency bandwidths because of the strict criteria adopted for processing the data. Second, an approach to include path-specific attenuation rates in GMPEs was proposed. This approach was applied to a subset of the KiK-net database. The attenuation rates across regions that contains volcanoes was found to be higher than other regions which is in line with the observations of other researchers. Moreover, accounting for the path-specific attenuation rates reduced the aleatoric variability associated with predicting pseudo-spectral accelerations. Third, two GMPEs were developed for active crustal earthquakes in Japan. The two GMPEs followed the ergodic and site-specific formulations, respectively. Finally, a comprehensive residual analysis was conducted to find potential biases in the residuals and propose models to predict some components of variability as a function of some input parameters. / Ph. D.

The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.

The characteristics of earthquake shaking are affected by the local site conditions. The effects of the local soil conditions are often quantified via an amplification factor (AF), which is defined as the ratio of the ground motion at the soil surface to the ground motion at a rock site at the same location. Amplification factors can be defined for any ground motion parameter, but most commonly are assessed for acceleration response spectral values at different oscillator periods. Site amplification can be evaluated for a site by conducting seismic site response analysis, which models the wave propagation from the base rock through the site-specific soil layers to the ground surface. An alternative to site-specific seismic response analysis is site amplification models. Site amplification models are empirical equations that predict the site amplification based on general characteristics of the site. Most of the site amplification models that already used in ground motion prediction equations characterize a site with two parameters: the average shear wave velocity in the top 30 m (VS30) and the depth to bedrock. However, additional site parameters influence site amplification and should be included in site amplification models. To identify the site parameters that help explain the variation in site amplification, ninety nine manually generated velocity profiles are analyzed using seismic site response analysis. The generated profiles have the same VS30 and depth to bedrock but a different velocity structure in the top 30 m. Different site parameters are investigated to explain the variability in the computed amplification. The parameter Vratio, which is the ratio of the average shear wave velocity between 20 m and 30 m to the average shear wave velocity in the top 10 m, is identified as the site parameter that most affects the computed amplification for sites with the same VS30 and depth to bedrock. To generalize the findings from the analyses in which only the top 30 m of the velocity profile are varied, a suite of fully randomized velocity profiles are generated and site response analysis is used to compute the amplification for each site for a range of input motion intensities. The results of the site response analyses conducted on these four hundred fully randomized velocity profiles confirm the influence of Vratio on site amplification. The computed amplification factors are used to develop an empirical site amplification model that incorporates the effect of Vratio, as well as VS30 and the depth to bedrock. The empirical site amplification model includes the effects of soil nonlinearity, such that the predicted amplification is a function of the intensity of shaking. The developed model can be incorporated into the development of future ground motion prediction equations. / text

This dissertation details work aimed at better understanding topographic effects in earthquake ground motions. The experiment, conducted in Central-Eastern Utah, used frequent and predictable seismicity produced by underground longwall coal mining as a source of low-intensity ground motions. Locally-dense arrays of seismometers deployed over various topographic features were used to passively monitor seismic energy produced by mining-induced implosions and/or stress redistribution in the subsurface. The research consisted of two separate studies: an initial feasibility experiment (Phase I) followed by a larger-scale main study (Phase II). Over 50 distinct, small-magnitude (M[subscript 'L'] < 1.6) seismic events were identified in each phase. These events were analyzed for topographic effects in the time domain using the Peak Ground Velocity (PGV), and in the frequency domain using the Standard Spectral Ratio (SSR) method, the Median Reference Method (MRM), and the Horizontal-to-Vertical Spectral Ratio (HVSR) method. The polarities of the horizontal ground motions were also visualized using directional analyses. The various analysis methods were compared to assess their ability to estimate amplification factors and determine the topographic frequencies of interest for each feature instrumented. The MRM was found to provide the most consistent, and presumably accurate, estimates of the amplification factor and frequency range for topographic effects. Results from this study clearly indicated that topographic amplification of ground motions does in fact occur. These amplifications were very frequency dependent, and the frequency range was correctly estimated in many, but not all, cases using simplified, analytical methods based on the geotechnical and geometrical properties of the topography. Amplifications in this study were found to generally range from 2 to 3 times a reference/baseline site condition, with some complex 3D features experiencing amplifications as high as 10. Maximum amplifications occurred near the crest of topographic features with slope angles greater than approximately 15 degrees, and the amplifications were generally oriented in the direction of steepest topographic relief, with some dependency on wave propagation direction. / text

Elevated water tanks are widely used to store water for drinking as well as for fire extinguishing purposes. After a severe earthquake, the need of water for drinking as well as fire control will increase dramatically. To ensure that water tanks remain functional after an earthquake, proper analysis method should be followed in order to calculate the response of a structure for earthquake. In this study, the lateral forces developed during earthquake are investigated from commercially available SAP2000 software and the results are compared with the 2006 edition of the ACI standard "Seismic Design of Liquid-Containing Concrete Structures and Commentary" (ACI 350.3-06). The elevated concrete tank is modeled for full, half-full and empty conditions. Linear modal time history analysis is performed using scaled ground motions. Three-directional ground motion records from five different earthquakes have been scaled to the design level and applied to the structure. Sloshing behavior of water inside the tank and the effect of vertical ground motion on the columns have been investigated. It is found that, vertical ground motions can increase the axial forces in columns by up to 20 %, and the ACI 350.3-06 design method is not always conservative. As seismic response depends on both the dynamic properties of the structure and the spectral characteristics of ground motions, more research is needed to understand and model the seismic response of elevated water tanks.

Kyoto University (京都大学) / 0048 / 新制・課程博士 / 博士(工学) / 甲第16788号 / 工博第3509号 / 新制||工||1531(附属図書館) / 29463 / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 清野 純史, 教授 小池 武, 准教授 古川 愛子 / 学位規則第4条第1項該当

This dissertation work investigates the effects of earthquake rupture complexity and heterogeneities in Earth structure on near-field ground-motions. More specifically, we address two key issues in seismology: (1) near-field ground-shaking variability as function of distance and azimuth for unilateral directive ruptures, and (2) impact of rupture complexity and seismic scattering on Mach wave coherence associated with supershear rupture propagation. We examine earthquake ground-motion variability associated with unilateral ruptures based on ground-motion simulations of the MW 7.3 1992 Landers earthquake, eight simplified source models, and a MW 7.8 rupture simulation (ShakeOut) for the San Andreas fault. Our numerical modeling reveals that the ground-shaking variability in near-fault distances (< 20 km) is larger than that given by empirical ground motion prediction equations. In addition, the variability decreases with increasing distance from the source, exhibiting a power-law decay. The high near-field variability can be explained by strong directivity effects whose influence weaken as we move away from the fault. At the same time, the slope of the power-law decay is found to be dominantly controlled by slip heterogeneity. Furthermore, the ground-shaking variability is high in the rupture propagation direction whereas low in the directions perpendicular to it. However, the variability expressed as a function of azimuth is not only sensitive to slip heterogeneity, but also to rupture velocity. To study Mach wave coherence for supershear ruptures, we consider heterogeneities in rupture parameters (variations in slip, rise time and rupture speed) and 3D scattering media having small-scale random heterogeneities. The Mach wave coherence is reduced at near-fault distances (< 10 km) by the source heterogeneities. At the larger distances from the source, medium scattering plays the dominant role in reducing the Mach wave coherence. Combined effect of the source and medium heterogeneities on the supershear ruptures produce peak ground accelerations consistent with the estimates from empirical ground motion prediction equations. Therefore, we suggest that supershear ruptures may be more common in nature than detected.

Genetic Expression Programming (GEP) is used to create Ground Motion Predicting Equations (GMPEs) for the average peak ground acceleration using 12,854 ground motion records obtained from the NGA-WEST2 project. The predictor set considered in this research consists of the moment magnitude, dip angle, rake angle, depth to the top of fault rupture, Joyner Boore distance, closest distance to the ruptured fault area, and the shear wave velocity in the top 30 m of the site. Four out of 23 candidate models were able to fairly predict the PGA for magnitudes larger than 4.5 and compared well with existing GMPEs in literature. GEP was capable of reasonably predicting the physical importance of the magnitude and distance parameters. However, other parameters often were either not fitted, or fitted as regression coefficients. The results illustrate GEP’s potential as a viable alternative to regression methods currently used in developing GMPEs.