Analysis of faulting destruction and water supply pipeline damage from the first mainshock of the February 6,2023 Türkiye earthquake doublet

In 2023,two consecutive earthquakes exceeding a magnitude of 7 occurred in Türkiye,causing severe casualties and economic losses.The damage to critical urban infrastructure and building structures,including highways,railroads,and water supply pipelines,was particularly severe in areas where these structures intersected the seismogenic fault.Critical infrastructure projects that traverse active faults are susceptible to the influence of fault movement,pulse velocity,and ground motions.In this study,we used a unique approach to analyze the acceleration records obtained from the seismic station array(9 strong ground motion stations)located along the East Anatolian Fault(the seismogenic fault of the MW7.8 mainshock of the 2023 Türkiye earthquake doublet).The acceleration records were filtered and integrated to obtain the velocity and displacement time histories.We used the results of an on-site investigation,jointly conducted by China Earthquake Administration and Türkiye’s AFAD,to analyze the distribution of PGA,PGV,and PGD recorded by the strong motion array of the East Anatolian Fault.We found that the maximum horizontal PGA in this earthquake was 3.0 g,and the maximum co-seismic surface displacement caused by the East Anatolian Fault rupture was 6.50 m.As the fault rupture propagated southwest,the velocity pulse caused by the directional effect of the rupture increased gradually,with the maximum PGA reaching 162.3 cm/s.We also discussed the seismic safety of critical infrastructure projects traversing active faults,using two case studies of water supply pipelines in Türkiye that were damaged by earthquakes.We used a three-dimensional finite element model of the PE(polyethylene)water pipeline at the Islahiye State Hospital and fault displacement observations obtained through on-site investigation to analyze pipeline failure mechanisms.We further investigated the effect of the fault-crossing angle on seismic safety of a pipeline,based on our analysis and the failure performance of the large-diameter Thames Water pipeline during the 1999 Kocaeli earthquake.The seismic method of buried pipelines crossing the fault was summarized.

Gravity Fault Subsidence and Beach Ridges Progradation in Quinta-Cassino (RS) Coastal Plain, Brazil

Ground penetrating radar (GPR) surveys have being applied to investigate very near-surface stratification of sedimentary units in coastal plains and to define their depositional conditions. This paper presents, however, low-frequency GPR survey to investigate fault-related depositional systems at greater depths. The Quinta-Cassino area in the Rio Grande do Sul Coastal Plain (RGSCP, Brazil) shows a wide strandplain that is made off by very long, continuous, and linear geomorphic features (beach ridges). This strandplain extends for ~70 km southward. The beach ridges show low-angle truncations against the Quinta escarpment, and also truncations in the strandplain. The traditional approach points that RGSCP was developed by juxtaposition of four lagoons/barrier systems as consequence of sea level changes;previous model assumes that no deformational episode occurred in RGSCP. The geophysical and geological surveys carried out in this area showed the existence of listric fault controlling the beach ridges in the escarpments and hanging-wall blocks. The radargrams could distinguish Pleistocene basement unit anticlockwise rotation, thickening of beach ridges radarfacies close to listric normal faults, and horst structures. These deformational features indicate that the extensional zone of a large-scale gravity-driven structure controlled the mechanical subsidence, the Holocene sedimentation and its stratigraphic and geomorphic features in the Quinta-Cassino area to build up an asymmetric delta. The results point to a new approach in dealing with RGSCP Holocene evolution.

Gravity-Driven Listric Growth Fault and Sedimentation in the Lagoa do Peixe, Rio Grande do Sul Coastal Plain, Brazil

High frequency, high resolution GPR surveys are successfully applied to investigate near-surface stratification architecture of sedimentary units in coastal plains and to define their depositional conditions. However, low frequency GPR surveys to investigate fault-related depositional systems at greater depths are scarce. This survey was designed investigate a > 100 km long linear escarpment that controls the northwest margin of the Lagoa do Peixe, an important lagoon in Rio Grande do Sul Coastal Plain (RGSCP, Brazil). The traditional approach points that RGSCP was developed by juxtaposition of four lagoons/barrier systems as consequence of sea level changes;no deformational structure is admitted to exist before. The low frequency GPR (50 MHz, RTA antenna) and geological surveys carried out in the RGSCP showed the existence of a large, gravity-driven listric growth fault controlling the Lagoa do Peixe escarpment and hangingwall sedimentation. The radargrams in four subareas along the Lagoa do Peixe Growth Fault could be interpreted following the seismic expression of rift-related depositional systems. The radargrams enabled to distinguish three main lagoonal deposition radarfacies. The lower lagoonal radarfacies is a convex upward unit, thicker close to growth fault;the radarfacies geometry indicates that fault displacement rate surpasses the sedimentation rate, and its upper stratum is aged ~3500 l4C years BP. The second lagoonal radarfacies is a triangular wedge restricted to the lagoon depocenter, whose geometry indicates that fault displacement and the sedimentation rates kept pace. The upper lagoonal radarfacies is being deposited since 1060 ± 70 l4C years BP, under sedimentation rate higher than fault displacement rate. The results indicate that low frequency GPR surveys can help in investigating fault-related depositional systems in coastal zones. They also point to a new approach in dealing with RGSCP stratigraphy.

Spatial distribution of near-fault ground motion

Near-fault strong ground motion of strike-slip and dip-slip of vertical and inclined rectangular fault in half-space and layered half-space is analyzed by dislocation source model. The Fourier spectra ratio of ground motion is adopted to study the characteristics of near-fault ground motion. For both slip models, near-fault strong ground motion with high amplitude is located in a narrow belt area along the projection of the fault on the ground and mainly controlled by the sub-faults nearby. Directivity of strike-slip fault is more dominant in long period for components perpendicular to the fault, and more dominant in long period for components parallel to the fault for dip-slip fault. The deeper the location of the source is, the more slowly the amplitude of ground motion attenuates. There is obvious hanging wall effect in ground motion of inclined fault, and the spatial distribution of ground motion is asymmetric which coincides with observational data. Finally, a fitting function of spatial distribution for near-fault ground motion is proposed and compared with near source factors of the 1997 Uniform Building Code of USA.

Dominant pulse simulation of near fault ground motions

In this study, a new mathematical model is developed composed of two parts, including harmonic and polynomial expressions for simulating the dominant velocity pulse of near fault ground motions. Based on a proposed velocity function, the corresponding expressions for the ground acceleration and displacement time histories are also derived. The proposed model is then fitted using some selected pulse-like near fault ground motions in the Next Generation Attenuation （NGA） project library. The new model is not only simple in form but also simulates the long-period portion of actual velocity near fault records with a high level of precision. It is shown that the proposed model-based elastic response spectra are compatible with the near fault records in the neighborhood of the prevailing frequency of the pulse. The results indicate that the proposed model adequately simulates the components of the time histories. Finally, the energy of the proposed pulse was compared with the energy of the actual record to confirm the compatibility.

Traveling wave effect on the seismic response of a steel arch bridge subjected to near fault ground motions

In the 1990s, several major earthquakes occurred throughout the world, with a common observation that near fault ground motion （NFGM） characteristics had a distinct impact on causing damage to civil engineering structures that could not be predicted by using far field ground motions. Since then, seismic responses of structures under NFGMs have been extensively examined, with most of the studies focusing on structures with relatively short fundamental periods, where the traveling wave effect does not need to be considered. However, for long span bridges, especially arch bridges, the traveling wave （only time delay considered） effect may be very distinct and is therefore important. In this paper, the results from a case study on the seismic response of a steel arch bridge under selected NFGMs is presented by considering the traveling wave effect with variable apparent velocities. The effects of fling step and long period pulses of NFGMs on the seismic responses of the arch bridge are also discussed.

Direct use of peak ground motion parameters for the estimation of inelastic displacement ratio of SDOF systems subjected to repeated far fault ground motions

This study is aimed at developing statistical equations to estimate the inelastic displacement ratio of singledegree-of-freedom systems subjected to far fault repeated earthquakes. In the study, peak ground motion parameters are used to define the scatter of the original data. The ratio of peak ground acceleration to peak ground velocity, and peak ground velocity of the ground motion records and structural parameters such as period of vibration and lateral strength ratio are used in the proposed equations. For the development of the equations, nonlinear time history analyses of single-degree-offreedom systems are conducted. Then, the results are used in a multivariate regression procedure. The equations are verified by comparing the estimated results with the calculated results. The average error and coefficient of variation of the proposed equations are presented. The analyses results revealed that the direct use of peak ground motion parameters for the estimation of inelastic displacement ratio significantly reduced the scatter in the original data and yielded accurate results. From the comparative results it is also observed that results obtained using equations specific to peak ground velocity or peak ground acceleration to peak ground velocity ratio are similar.

Influences of pre-existing fracture on ground deformation induced by normal faulting in mixed ground conditions

Physical model tests have been conducted by various researchers to investigate fault rupture propagation and ground deformation induced by bedrock faulting. However, the effects of pre-existing fracture on ground deformation are not fully understood. In this work, six centrifuge tests are reported to investigate the influence of pre-existing fracture on ground deformation induced by normal faulting in sand, clay and nine-layered soil with interbedded sand and clay layers. Shear box tests were conducted to develop a filter paper technique, which was adopted in soil model preparation to simulate the effects of pre-existing fracture in centrifuge tests. Centrifuge test results show that ground deformation mechanism in clay, sand and nine-layered soil strata is classified as a stationary zone, a shearing zone and a rigid body zone. Inclination of the strain localization is governed by the dilatancy of soil material. Moreover, the pre-existing fracture provides a preferential path for ground deformation and results in a scarp at the ground surface in sand. On the contrary, fault ruptures are observed at the ground surface in clay and nine-layered soil strata.

Single-phase-to-ground fault protection based on zero-sequence current ratio coefficient for low-resistance grounding distribution network

Definite-time zero-sequence over-current protection is presently used in systems whose neutral point is grounded by a low resistance(low-resistance grounding systems).These systems frequently malfunction owing to their high settings of the action value when a high-impedance grounding fault occurs.In this study,the relationship between the zero-sequence currents of each feeder and the neutral branch was analyzed.Then,a grounding protection method was proposed on the basis of the zero-sequence current ratio coefficient.It is defined as the ratio of the zero-sequence current of the feeder to that of the neutral branch.Nonetheless,both zero-sequence voltage and zero-sequence current are affected by the transition resistance,The influence of transition resistance can be eliminated by calculating this coefficient.Therefore,a method based on the zero-sequence current ratio coefficient was proposed considering the significant difference between the faulty feeder and healthy feeder.Furthermore,unbalanced current can be prevented by setting the starting current.PSCAD simulation results reveal that the proposed method shows high reliability and sensitivity when a high-resistance grounding fault occurs.

3D simulation of near-fault strong ground motion: comparison between surface rupture fault and buried fault

In this paper, near-fault strong ground motions caused by a surface rupture fault （SRF） and a buried fault （BF） are numerically simulated and compared by using a time-space-decoupled, explicit finite element method combined with a multi-transmitting formula （MTF） of an artificial boundary. Prior to the comparison, verification of the explicit element method and the MTF is conducted. The comparison results show that the final dislocation of the SRF is larger than the BF for the same stress drop on the fault plane. The maximum final dislocation occurs on the fault upper line for the SRF; however, for the BE the maximum final dislocation is located on the fault central part. Meanwhile, the PGA, PGV and PGD of long period ground motions （≤ 1 Hz） generated by the SRF are much higher than those of the BF in the near-fault region. The peak value of the velocity pulse generated by the SRF is also higher than the BE Furthermore, it is found that in a very narrow region along the fault trace, ground motions caused by the SRF are much higher than by the BF. These results may explain why SRFs almost always cause heavy damage in near-fault regions compared to buried faults.

An Ellipsoidal Set-Membership Approach to Distributed Joint State and Sensor Fault Estimation of Autonomous Ground Vehicles

This paper is concerned with the problem of distributed joint state and sensor fault estimation for autonomous ground vehicles subject to unknown-but-bounded(UBB)external disturbance and measurement noise.In order to improve the estimation reliability and performance in cases of poor data collection and potential communication interruption,a multisensor network configuration is presented to cooperatively measure the vehicular yaw rate,and further compute local state and fault estimates.Toward this aim,an augmented descriptor vehicle model is first established,where the unknown sensor fault is modeled as an auxiliary state of the system model.Then,a new distributed ellipsoidal set-membership estimation approach is developed so as to construct an optimized bounding ellipsoidal set which guarantees to contain the vehicle’s true state and the sensor fault at each time step despite the existence of UBB disturbance and measurement noises.Furthermore,a convex optimization algorithm is put forward such that the gain matrix of each distributed estimator can be recursively obtained.Finally,simulation results are provided to validate the effectiveness of the proposed approach.

Fault Line Selection Method Considering Grounding Fault Angle for Distribution Network

In the distribution network system with its neutral point grounding via arc suppression coil, when single-phase grounding fault occurred near zero-crossing point of the phase voltage, the inaccuracy of the line selection always existed in existing methods. According to the characteristics that transient current was different between the fault feeder and other faultless feeders, wavelet transformation was performed on data of the transient current within a power frequency cycle after the fault occurred. Based on different fault angles, wavelet energy in corresponding frequency band was chosen to compare. The result was that wavelet energy in fault feeder was the largest of all, and it was larger than sum of those in other faultless feeders, when the bus broke down, the disparity between each wavelet energy was not significant. Fault line could be selected out by the criterion above. The results of MATLAB/simulink simulation experiment indicated that this method had anti-interference capacity and was feasible.

Application of Atomic Sparse Decomposition to Feature Extraction of the Fault Signal in Small Current Grounding System

Applying the atomic sparse decomposition in the distribution network with harmonics and small current grounding to decompose the transient zero sequence current that appears after the single phase to ground fault occurred. Based on dictionary of Gabor atoms and matching pursuit algorithm, the method extracts the atomic components iteratively from the feature signals and translated them to damped sinusoidal components. Then we can obtain the parametrical and analytical representation of atomic components. The termination condition of decomposing iteration is determined by the threshold of the initial residual energy with the purpose of extract the features more effectively. Accordingly, the proposed method can extract the starting and ending moment of disturbances precisely as well as their magnitudes, frequencies and other features. The numerical examples demonstrate its effectiveness.

Small-current grounding fault location method based on transient main resonance frequency analysis

The small-current grounding fault in distribution network is hard to be located because of its weak fault features.To accurately locate the faults,the transient process is analyzed in this paper.Through the study we take that the main resonant frequency and its corresponding component is related to the fault distance.Based on this,a fault location method based on double-end wavelet energy ratio at the scale corresponding to the main resonant frequency is proposed.And back propagation neural network(BPNN)is selected to fit the non-linear relationship between the wavelet energy ratio and fault distance.The performance of this proposed method has been verified in different scenarios of a simulation model in PSCAD/EMTDC.

Representation of near-fault pulse-type ground motions

Near-fault ground motions with long-period pulses have been identified as critical in the design of structures. To aid in the representation of this special type of motion, eight simple pulses that characterize the effects of either the flingstep or forward-directivity are considered. Relationships between pulse amplitudes and velocity pulse period for different pulses are discussed. Representative ratios and peak acceleration amplification can exhibit distinctive features depending on variations in pulse duration, amplitude and the selected acceleration pulse shape. Additionally, response spectral characteristics for the equivalent pulses are identified and compared in terms of fixed PGA and PGV, respectively. Response spectra are strongly affected by the duration of pulses and the shape of the basic pulses. Finally, dynamic time history response features of a damped SDOF system subjected to pulse excitations are examined. These special aspects of pulse waveforms and their response spectra should be taken into account in the estimation of ground motions for a project site close to a fault.

An investigation of the influence of near-fault ground motion parameters on the pile's response in liquefiable soil

The performance of a soil-pile system can be significantly influenced by ground motion parameters. However, few research efforts have been performed to provide a complete description of the influence of key ground motion parameters on the pile’s behavior in liquefiable soil. In this study, a three-dimensional finite element（FE） model, incorporating a multisurface plasticity solid-fluid fully coupled formulation soil constitutive model, is developed and calibrated based on centrifuge test data. Seventy-two near-fault non-pulse-like（NF-NP） and seventy-two near-fault pulse-like（NF-P） ground motion records are studied with the calibrated FE model to distinguish the effects of several common ground motion parameters soon afterwards. Base on the parametric study results, a simple index, RPGV/PGA（i.e., the ratio of peak ground velocity（PGV） to peak ground acceleration（PGA））, shows its capability on characterizing the pile behavior under both NF-NP and NF-P ground motions. Furthermore, two equations are developed to characterize the relationships between the RPGV/PGA as well as the maximum pile’s moments and displacements. In general, this study can be helpful to gain new insights on the influence of typical index parameters for near-field ground motions on the response of the pile foundation in liquefiable soil.

Study on equivalent velocity pulse of nearfault ground motions

Near-fault strong ground motions that resulted in serious structural damage are characterized by directivity effect and pulse-type motion. Large-amplitude and long-period pulses are contained in the velocity time-history traces of near-fault pulse-type records. A reasonable model of equivalent velocity pulse is proposed on the basis of the ex- isted models in this paper to simplify the calculation and analysis. Based on the large amount of collected near-fault strong earthquakes records, the parameters describing equivalent velocity pulse model such as pulse period, pulse intensity and number of predominant pulses are studied, and comparison is made with the results obtained by others models. The proposed model is contributive to the seismic design for structures in near-fault areas.

Near-fault ground motions with prominent acceleration pulses:pulse characteristics and ductility demand

Major earthquakes of last 15 years （e.g., Northridge 1994, Kobe 1995 and Chi-Chi 1999） have shown that many near-fault ground motions possess prominent acceleration pulses. Some of the prominent ground acceleration pulses are related to large ground velocity pulses, others are caused by mechanisms that are totally different from those causing the velocity pulses or fling steps. Various efforts to model acceleration pulses have been reported in the literature. In this paper, research results from a recent study of acceleration pulse prominent ground motions and an analysis of structural damage induced by acceleration pulses are summarized. The main results of the study include： （1） temporal characteristics of acceleration pulses; （2） ductility demand spectrum of simple acceleration pulses with respect to equivalent classes of dynamic systems and pulse characteristic parameters; and （3） estimation of fundamental period change under the excitation of strong acceleration pulses. By using the acceleration pulse induced linear acceleration spectrum and the ductility demand spectrum, a simple procedure has been developed to estimate the ductility demand and the fundamental period change of a reinforced concrete （RC） structure under the impact of a strong acceleration pulse.

A coherency function model of ground motion at base rock corresponding to strike-slip fault

At present, the method to study spatial variation of ground motions is statistic analysis based on dense array re-cords such as SMART-1 array, etc. For lacking of information of ground motions, there is no coherency function model of base rock and different style site. Spatial variation of ground motions in elastic media is analyzed by deterministic method in this paper. Taking elastic half-space model with dislocation source of fault, near-field ground motions are simulated. This model takes strike-slip fault and earth media into account. A coherency func-tion is proposed for base rock site.

Study on inelastic displacement ratio spectra for near-fault pulse-type ground motions

In displacement-based seismic design, inelastic displacement ratio spectra （IDRS） are particularly useful for estimating the maximum lateral inelastic displacement demand of a nonlinear SDOF system from the maximum elastic displacement demand of its counterpart linear elastic SDOF system. In this study, the characteristics of IDRS for near-fault pulse-type ground motions are investigated based on a great number of earthquake ground motions. The influence of site conditions, ratio of peak ground velocity （PGV） to peak ground acceleration （PGA）, the PGV, and the maximum incremental velocity （MIV） on IDRS are also evaluated. The results indicate that the effect of near-fault ground motions on IDRS are significant only at periods between 0.2 s - 1.5 s, where the amplification can approach 20%. The PGV/PGA ratio has the most significant influence on IDRS among the parameters considered. It is also found that site conditions only slightly affect the IDRS.