Ensuring the accuracy of free-field inversion is crucial in determining seismic excitation for soil-structure interaction (SSI) systems. Due to the spherical and cylindrical diffusion properties of body waves and surface waves, the near-fault zone presents distinct free-field responses compared to the far-fault zone. Consequently, existing far-fault free-field inversion techniques are insufficient for providing accurate seismic excitation for SSI systems within the near-fault zone. To address this limitation, a tailored near-fault free-field inversion method based on a multi-objective optimization algorithm is proposed in this study. The proposed method establishes an inversion framework for both spherical body waves and cylindrical surface waves and then transforms the overdetermined problem in inversion process into an optimization problem. Within the multi-objective optimization model, objective functions are formulated by minimizing the three-component waveform differences between the observation point and the delayed reference point. Additionally, constraint conditions are determined based on the attenuation property of propagating seismic waves. The accuracy of the proposed method is then verified through near-fault wave motion characteristics and validated against real downhole recordings. Finally, the application of the proposed method is investigated, with emphasis on examining the impulsive property of underground motions and analyzing the seismic responses of SSI systems. The results show that the proposed method refines the theoretical framework of near-fault inversion and accurately restores the free-field characteristics, particularly the impulsive features of near-fault motions, thereby providing reliable excitation for seismic response assessments of SSI systems.
A utility tunnel is an infrastructure that consolidates multiple municipal pipeline systems into a shared underground passage. As long linear structures inevitably cross different soils, this paper aims to accurately assess the seismic damage to a shallow-buried utility tunnel in a non-homogeneous zone by employing a viscous-spring artificial boundary and deriving the corresponding nodal force equations. The three-dimensional model of the utility tunnel-soil system is established using finite element software, and a plug-in is developed to simulate the three-dimensional oblique incidence of SV waves with a horizontal non-homogeneous field. In this study, the maximum interstory displacement angle of the utility tunnel is used as the damage indicator. Analysis of structural vulnerability based on IDA method using PGA as an indicator of seismic wave intensity, which considers the angle of oblique incidence of SV waves, the type of seismic waves, and the influence of the nonhomogeneous field on the seismic performance of the utility tunnel. The results indicate that the failure probability of the utility tunnel in different soil types increases with the incident angle and PGA. Additionally, the failure probability under the pulse wave is higher than that under the non-pulse wave; Particular attention is given to the states of severe damage (LS) and collapse (CP), particularly when the angle of incidence is 30 degrees and the PGA exceeds 0.6g, conditions under which the probability of failure is higher. Additionally, the failure probability of the non-homogeneous zone is greater than that of sand and clay; the maximum interlayer displacement angle increases with the incident angle, accompanied by greater PGA dispersion, indicating the seismic wave intensity. The maximum inter-layer displacement angle increases with the incident angle, and the dispersion of the seismic wave intensity indicator (PGA) becomes greater. This paper proposes vulnerability curves for different working conditions, which can serve as a reference for the seismic design of underground structures.
This study investigates the influence of wood pellet fly ash blended binder (WABB) on the mechanical properties of typical weathered granite soils (WS) under a field and laboratory tests. WABB, composed of 50 % wood pellet fly ash (WA), 30 % ground granulated blast furnace slag (GGBS), and 20% cement by dry mass, was applied at dosages of 200-400 kg/m3 to four soil columns were constructed at a field site deposited with WS. After 28 days, field tests, including coring, standard penetration tests (SPT), and permeability tests, revealed enhanced soil cementation and reduced permeability, indicating a denser soil matrix. Unconfined compressive tests (UCT) and free-free resonant column (FFRC) tests on field cores at 28 and 56 days, compared with laboratory specimens and previously published data, demonstrated strength gains 1.2-2.1 times higher due to field-induced stress. The presence of clay minerals influenced the WABB's interaction and microstructure development. Correlations between seismic waves, small-strain moduli, and strength were developed to monitor in-situ static and dynamic stiffness gain of WABB-stabilized weathered granite soils.
The influence of surface Rayleigh waves (SRWs) on the seismic behavior of three archetype nonconforming reinforced concrete (RC) buildings including weak first story with four, six, and eight stories when subjected to earthquake ground motions (EQGMs) recorded during the strong September 19, 2017 Mw7.1 earthquake in Mexico City, is discussed in this paper. For this purpose, ground acceleration time histories corresponding to the retrograde and prograde components of SRWs were extracted from EQGMs collected at the accelerographic stations placed at the transition and soft soil sites. It was found that the SWRs contribute to about 50% of the median maximum IDR demand (IDRmax) triggered by the as-recorded earthquake ground motions at the ground level of the four- and six-story building models, while their contribution is about 30% of IDRmax for the eight-story building model. It should be noted that that SRWs induce median IDRmax demands to the four-story building model larger than about 11% and 49% than those to the six- and eight-story models, respectively, for soft soil sites. Moreover, the prograde component can trigger IDRmax demands in the four-story building model larger than 73% and 45% than those for the six- and eight-story models, respectively, for the transition sites. Particularly, it was shown that SRWs induce median IDRmax demands in excess of 0.35% at the first level of the archetype building models, which is associated to the light cracking damage state of nonductile RC columns, and even in excess of IDRmax of 0.71% associated to the severe cracking damage state when the record-to-record variability is considered in the IDRmax demand (i.e. the 84th percentile of IDRmax). Although the earthquake ground motion component of the surface Rayleigh waves was negligible in the median IDRmax, this study showed that the effect of the directionality of IDRmax is important for the CH84 station, where significant polarity of spectral ordinates was identified in previous studies.
With the rapid development of infrastructure in western China, numerous arch bridges have been constructed as vital transportation hubs spanning river canyons. Understanding the impact of canyon topography on the seismic response of long-span half-through arch bridges crossing canyons is essential. This study first establishes a seismic input method for oblique P-wave and SV-wave incidence, based on the viscous-spring artificial boundary theory, which transforms ground motions into equivalent nodal loads on artificial boundaries. The feasibility of this proposed method is systematically validated. Subsequently, parametric investigations are carried out to explore the effects of seismic wave incidence angle, canyon depth-to-breadth ratio and soil elastic modulus on the ground motion amplification characteristics in V-shaped canyons under oblique P-wave and SV-wave excitations. Finally, dynamic response patterns of the arch ribs and the stress-strain relationships at critical structural components are thoroughly analyzed. Key findings reveal that SV-waves induce significantly different ground motion amplification effects compared to P-waves, with the wave incidence angle and canyon width-to-depth ratio being crucial influencing factors. The connection between the arch footings and the concrete cross braces constitutes the most vulnerable region, frequently exhibiting maximum stresses that exceed the yield strength of C40 concrete under multiple scenarios. Notably, when the depth-to-breadth ratio (D/B) is 0.75, the peak stress at the arch footings reaches 5.18 x 10(7)kPa, surpassing the yield stress threshold of C40 concrete. These findings highlight the need for special seismic fortification measures at these critical connections during bridge design. This research offers valuable insights into the seismic design of long-span arch bridges in complex topographic conditions.
This study investigates the liquefaction characteristics of deep soil layers and their subsequent effects on the seismic response of subway station structures, utilizing shaking table tests and inputting seismic waves of varying principal frequencies. Macroscopically, the liquefaction of deep soil strata does not result in surface manifestations such as water spraying and sand bubbling. However, it still induces cracking and damage to the soil surrounding the structure. Analyzing from the perspective of the pore pressure ratio reveals that the ratio under free-field conditions is significantly lower than under structural conditions. Additionally, the pore pressure ratio caused by the Beijing Hotel wave is greater, followed by the Beijing artificial wave, while the Ming Shan wave results in the smallest ratio. In terms of the station structure, the structural acceleration and tensile strain increment induced by the Beijing Hotel wave are the most significant, followed by the Beijing artificial wave, with the least effect from the Ming Shan wave. This indicates that the liquefaction behavior of deep soil layers is primarily influenced by the overlying load and the frequency characteristics of seismic waves. The construction of subway stations reduces the overlying loads on soil layers, increasing the likelihood of soil layer liquefaction. Meanwhile, a lower main frequency of the seismic wave results in a higher degree of liquefaction in the deep soil layers. The seismic response of the station structure is contingent on the frequency characteristics of the seismic wave. The lower the primary frequency of the seismic wave, the higher the seismic response of the station structure. Furthermore, the liquefaction behavior of the deep soil layers also impacts the seismic response of the station structure, particularly the tensile strain response of the top and bottom slabs of the station structure.
Resonance can significantly amplify a structure's response to seismic loads, leading to extended damage, especially in critical infrastructure like nuclear power plants. Thus, this study focuses on the resonance effects of the dynamic interaction between layered soil, pile foundations, and nuclear island structures, which is particularly important given the limited availability of bedrock sites for such facilities. Specifically, this study explores the resonance behavior of nuclear islands under various seismic conditions through large-scale shaking table tests by developing a dynamic interaction model for layered soil-pile-nuclear island systems. The proposed model comprises a 3 x 3 pile group supporting the upper structure of a nuclear island embedded within a three-layer soil profile. Sinusoidal waves of varying frequencies identify the factors influencing the system's resonance response. Besides, the resonance effects are validated by inputting seismic motions based on compressed acceleration time histories. Furthermore, the impact of non-primary frequency components on structural resonance is assessed by comparing sinusoidal wave components. The findings reveal that resonance effects increase as the amplitude of the input seismic motion increases to a certain threshold, after which the effect stabilizes. This trend is particularly pronounced in the bending moment response at the pile head. Additionally, an independent resonance phenomenon is observed in the superstructure, suggesting that its resonance effects should be considered separately in nuclear island design. Similar resonance effects are observed when the predominant frequency of sinusoidal waves closely matches the compressed seismic motions, suggesting that sinusoidal inputs effectively simulate structural resonance during seismic design testing.
Liquefaction and shear sliding (i.e., slides) are common failure modes of cohesionless seafloor under ocean waves. However, existing research has rarely focused on shear-sliding failure, especially when considering wave- induced residual pore water pressure. Additionally, the relationship between shear-sliding failure and liquefaction is not well understood. In this study, a slice method is developed to assess the shear-sliding failure in cohesionless seafloor under non-linear waves, incorporating the effect of both oscillatory and residual seabed responses. The applicability of various liquefaction criteria is discussed, based on the interrelation between the shear-sliding and liquefaction zones. The results indicate that the seabed soil is more prone to shear-sliding failure than liquefaction under wave-induced pore water pressure. When only oscillatory pore water pressure is considered, the liquefaction criteria, assuming the initial vertical effective stress vanishes due to the excess pore water pressure, better identify the liquefaction zone, which is enveloped by and overlaps with the shear- sliding zone at a factor of safety of zero. In cases where both oscillatory and residual pore pressure coexist, the unified liquefaction criterion, which also assumes onset of liquefaction at zero vertical effective stress, provides more reliable predictions of the liquefaction zone. As residual pore pressure accumulates, the difference between shear sliding and liquefaction depths becomes more pronounced. A sensitivity analysis of shear-sliding depth with varying soil parameters indicates that relative density exerts the most significant influence, followed by the effective internal friction angle, while the shear modulus has the least effect. The effect of variations in soil parameters on shear-sliding depth diminishes to some extent with prolonged wave action.
The Tibetan Plateau, a critical region influencing both local and global atmospheric circulation, climate dynamics, hydrology and terrestrial ecosystems, is undergoing climate-driven changes, including glacial retreat, permafrost thaw and groundwater changes. Despite its importance, implementing continuous and systematic observations has been challenging due to the area's high altitude and extreme climate conditions. In this context, seismic interferometry emerges as a cost-effective method for the continuous monitoring of subsurface structural changes driven by environmental factors and internal geophysical processes. We investigate subsurface evolution using four years of seismic data from nine stations on the northeastern Tibetan Plateau, by applying coda wave interferometry across multiple frequency bands. Our findings highlight seismic velocity changes within the frequency bands 5-10, 0.77-1.54, and 0.25-0.51 Hz, revealing depth-dependent seasonal and long-term changes. Near-surface and deeper strata exhibit similar seasonal patterns, with velocities increasing in winter and decreasing in summer driven by changes in hydrological processes, while intermediate ice-water phase strata show contrasting behaviour due to thermal elastic strain. Long-term trends suggest that the upper subsurface layer is affected by melting water and precipitation originating from Kunlun Mountains, whereas deeper layer reflect groundwater level variations influenced by climate change and human activities. This study provides insights into the environmental evolution of the Tibetan Plateau and its impact on managing local groundwater resources.
Estimating the landscape and soil freeze-thaw (FT) dynamics in the Northern Hemisphere (NH) is crucial for understanding permafrost response to global warming and changes in regional and global carbon budgets. A new framework for surface FT-cycle retrievals using L-band microwave radiometry based on a deep convolutional autoencoder neural network is presented. This framework defines the landscape FT-cycle retrieval as a time-series anomaly detection problem, considering the frozen states as normal and the thawed states as anomalies. The autoencoder retrieves the FT-cycle probabilistically through supervised reconstruction of the brightness temperature (TB) time series using a contrastive loss function that minimizes (maximizes) the reconstruction error for the peak winter (summer). Using the data provided by the Soil Moisture Active Passive (SMAP) satellite, it is demonstrated that the framework learns to isolate the landscape FT states over different land surface types with varying complexities related to the radiometric characteristics of snow cover, lake-ice phenology, and vegetation canopy. The consistency of the retrievals is assessed over Alaska using in situ observations, demonstrating an 11% improvement in accuracy and reduced uncertainties compared to traditional methods that rely on thresholding the normalized polarization ratio (NPR).