Modeling of co- and post-seismic surface deformation and gravity changes of M W 6.9 Yushu, Qinghai,earthquake

Based on the elastic dislocation theory, multilayered crustal model, and rupture model obtained by seismic waveform inversion, we calculated the coand post-seismic surface deformation and gravity changes caused by the Yushu M W 6.9 earthquake occurred on April 14, 2010. The observed GPS velocity field and gravity field in Yushu areas are disturbed by the coand post-seismic effects induced by Yushu earthquake, thus the theoretical coand post-seismic deformation and gravity changes will provide important modification for the background tectonic movement of Yushu and surrounding regions. The time relaxation results show that the influences of Yushu earthquake on Yushu and surrounding areas will last as long as 30 to 50 years. The maximum horizontal displacement, vertical uplift and settlement are about 1.96, 0.27 and 0.16 m, respectively, the maximal positive and negative value of gravity changes are 8.892×10-7 m·s-2 and -4.861×10-7 m·s-2 , respectively. Significant spatial variations can be found on the coand post-seismic effects： The co-seismic effect mainly concentrates in the region near the rupture fault, while viscoelastic relaxation mostly acts on the far field. Therefore, when using the geodetic data to research tectonic motion, we should not only consider the effect of co-seismic caused by earthquake, but also pay attention to the effect of viscoelastic relaxation.

Near-field surface deformation during the April 20,2013,Ms7.0 Lushan earthquake measured by 1-Hz GNSS

The April 20, 2013, Ms7.0 Lushan earthquake was successfully recorded by closely spaced Continuous Global Positioning System （CGPS） stations owned by the Crustal Movement Observation Network of Chi- na （CMONC）. The 1-Hz GNSS data from eight CGPS stations, which are located between 30 km and 200 km from the hypocenter, were processed within quasi-real-time. The near-field surface deformation indicated the following characteristics ： the near-field movements were limited to several centimeters ; the peak of the deformation wave was significantly larger than the static permanent offset; at the beginning of the event, the north wall of the fault moved to the southeast as the south wall moved to the southwest ; station SCTQ, which was the closest station to the hypocenter at 30 km, had the largest static permanent displacement of 2 cm; the peaks of the deformation waves were 1.5 cm, 5 cm and 3 cm, to the east, the south and vertically upward, respectively ; and the peaks of velocity and acceleration, derived from the deformation, were 3.4 cm/s and 5.3 cm/s^2,respectively.

Specular reflection based sensing surface deformation of gas tungsten arc weld pool

A sensing system is developed to measure the weld pool boundary and pool suoface deformation in gas tungsten arc welding. LaserStrobe technique is used to eliminate the strong arc light inteoference, and specular reflection from the pool suoface is sensed to describe the relation between the deformed stripes and pool surface depression. Clear images of both the pool boundary and the deformed stripes edges are obtained during gas tungsten arc welding process, which lays foundation for realtime monitoring the pool suoface depression and weld penetration.

Co-and post-seismic surface deformation and gravity changes of MS 7.0 Lushan,earthquake

On April 20, 2013, an earthquake with mag- nitude 7.0 occurred in the southwest of the Longmenshan fault system in and around Lushan County, Sichuan Province, China. This devastating earthquake killed hun- dreds of people, injured 10 thousand others, and collapsed countless buildings. In order to analyze the potential risk of this big earthquake, we calculate the co- and post-seismic surface deformation and gravity changes of this event. In this work, a multilayered crustal model is designed, and the elastic dislocation theory is utilized to calculate the co- and post-seismic deformations and gravity changes. During the process, a rupture model obtained by seismic waveform inversion （Liu et al. Sci China Earth Sci 56（7）： 1187-1192, 2013） is applied. The time-dependent relaxation results show that the influences on Lushan and its surrounding areas caused by the Ms7.0 Lushan earthquake will last as long as 10 years. The maximum horizontal displacement, vertical uplift, and settlement are about 5 cm, 21.24 cm, and 0.16 m, respectively; the maximal positive and nega- tive values of gravity changes are 45 and -0.47 μGal, respectively. These results may be applied to evaluate the long-term potential risk caused by this earthquake and to provide necessary information for post-earthquake reconstruction.

Reconstruction of surface deformation characteristics in alpine canyons under shadow conditions

Monitoring deformation in high undulating mountainous environments is critical for surface process research and disaster prevention studies. Although observations based on interferometric Synthetic Aperture Radar(InSAR) are an excellent tool for monitoring deformation, the shadow phenomena can limit its application. Based on a series of geomorphic parameters and limited InSAR observation data, surface deformations were reconstructed in areas with missing observations by constructing a random forest model to compensate for the shadow phenomenon at the grid-scale. The findings suggest that this method can be used to rebuild landscape variation characteristics in places where observation data is lacking. The dominant slope direction in the observation area corresponded to a more significant correlation between the reconstructed topography deformation characteristics and the observation. In addition, when building this model, consideration was given to the geomorphic parameter selection, elevation variation, hypsometric integral value, slope form, lithology, slope variation,and aspect variation;these parameters can significantly affect the surface deformation, which is closely related to their spatial autocorrelation. These findings are significant for eliminating the shadow phenomenon, which often occurs in In SAR observations taken over alpine canyon regions. The terrain and lithology of the underlying surface should be considered when reconstructing the surface deformation characteristics of the shadow region by using satellite observation data.

Approximate Tangent Plane Method for Calculating Surface Deformation in Elastic Contact Problems

A flew numerical method for constructing a pressure distribution to calculate surface elastic deformationcaused by normal contact pressure is developed in this paper. The pressure distribution over one of nonequidistantrectangles is fitted by an approximate tangent plane(ATP), which is formed by five pressure samples. Because thepressure distribution could be expressed as an one order linear polynomial, the iterative expression of elasticdeformation deduced by this method is simple, and the numerical accuracy is higher.

Rainfall-triggered waste dump instability analysis based on surface 3D deformation in physical model test

Landslide is the second largest natural disaster after earthquake. It is of significance to study the evolution laws and failure mechanism of landslides based on its surface 3D deformation information. Based on the rainfall-triggered waste dump instability model test, we studied the failure mechanisms of the waste dump by integrating surface deformation and internal slope stress and proposed novel parameters for identifying landslide stability. We developed a noncontact measurement device, which can obtain millimeter-level 3D deformation data for surface scene in physical model test;Then we developed the similar materials and established a test model for a waste dump. Based on the failure characteristics of slope surface, internal stress of slope body and displacement contours during the whole process, we divided the slope instability process in model test into four stages: rainfall infiltration and surface erosion, shallow sliding, deep sliding, and overall instability. Based on the obtained surface deformation data, we calculated the volume change during slope instability process and compared it with the point displacement on slope surface. The results showed that the volume change can not only reflect the slow-ultra acceleration process of slope failure, but also fully reflect the above four stages and reduce the fluctuations caused by random factors. Finally, this paper proposed two stability identification parameters: the volume change rate above the slip surface and the relative velocity of volume change rate. According to the calculation of these two parameters in model test, they can be used for study the deformation and failure mechanism of slope stability.

Numerical assessment of the influence of former mining activities and plasticity of rock mass on deformations of terrain surface

This paper presents the results of numerical simulations carried out to confirm the influence of former mining activities on deformation of the mining terrain.The assessment of deformation changes was carried out with the use of FLAC3 D program based on the finite difference method.Numerical calculations were carried out for the example of actual mining operations in seams 703/1-2 and 707/2 of‘‘Marcel"Coal Mine.Taking into account the influence of the model’s plastic features and the so-called activation of a higher occurring seam in conducted simulations enabled obtaining a very good description of the measured subsidence.Based on the results one may state that numerical model can be used to assess the influence of former mining activities and the direction of conducted exploitation on deformations of the mining terrain.These factors are not recognized by geometric-integral theories commonly used for predicting the influence of mining operations on the surface.The results presented in this paper confirm that the applied method of simulating the phenomenon of reactivation of post-mining goafs is correct.

Lower Crustal Flow and Its Relation to the Surface Deformation and Stress Distribution in Western Sichuan Region, China

The channel flow model was gradually being accepted with the more important multidisci- plinary evidences from geology and geophysics, but how the lower crustal flow influenced the surface deformation quantitatively was unknown. Here, we develop a three-dimensional viscoelastic model to explore the mechanical relations between the lower crustal flow and the surface deformation in western Sichuan. Based on numerous tests, our results show that the modeled results fit well with the observed GPS data when the lower crust flows faster than the upper crust about 11 mm/a in the rhombic block, which can be useful to understand the possible mechanism of the surface deformation in western Siehuan. Moreover, taking the Xianshuihe fault as an example, we preliminarily analyze the relation between the active fault and stress field, according to the boundary constraints that deduced from the best model. The results show that the maximum shear stress on the Xianshuihe fault zone is mainly located in the fault terminal, intersections and the bend of the fault geometry, the stress level on the northwestern segment that has the high slip rate is relatively high. Additionally, with the reduction of the Young＇s modulus in the fault zone, it＇s conducive to generate the greater strain distribution, hence forming the high stress level.

Determination of convergence of underground gas storage caverns using non-invasive methodology based on land surface subsidence measurement

Undergroundgas storage caverns aremonitoredfor environmental safety in termsof equipmentandpotential emissions,particularly methane emissions from the underground and above-ground parts of the storage facility.Periodical measurements of land surface deformations and costly echometric measurements of convergence of individual storage facilities are carried out.The aims of environmental monitoring are:(1)to eliminate potential hazards in the shortest time,(2)assess the overall impact of intensive operation of storage facilities on the environment,(3)developmonitoringmethods relevant to environmental protection,and(4)take actions in case of failure.The paper presents a solution to the problem of determination of the convergence of underground caverns in a salt rock mass based on the results of land surface subsidence measurements carried out using the Gauss-Markov equalization algorithm.Themethod makes it possible for ongoing control of cavern volume convergence after each subsidence measurement on the ground surface and determining the actual impact of the use frequency(injection-mediumconsumption)on the convergence in time.The presentedmethodology is universal and verified on caverns located in a salt rockmass.The Gauss-Markov inversion model is the first used in this area,hence its application is significant.

Noise reduction and periodic signal extraction for GNSS height data in the study of vertical deformation

Global navigation satellite system(GNSS)technique has irreplaceable advantages in the continuous monitoring of surface deformation.Reducing noise to improve the signal-to-noise ratio(SNR)and extract the concerned signals is of great significance.As an improved algorithm of empirical mode decomposition(EMD),complete ensemble empirical mode decomposition with adaptive noise(CEEMDAN)algorithm has better signal processing ability.Using the CEEMDAN algorithm,the height time series of 29GNSS stations in Chinese mainland were analyzed,and good denoising effects and extraction from periodic signals were achieved.The numerical results showed that the annual signal obtained with the CEEMDAN algorithm was significantly based on Lomb_Scargle spectrum analysis,and large differences in the long-term signals were found between the stations at different locations in Chinese mainland.With respect to data denoising,compared with the EMD and wavelet denoising algorithms,the CEEMDAN algorithm respectively improved the SNR by 29.35% and 36.54%,increased the correlation coefficient by 8.67% and 11.96%,and reduced root mean square error(RMSE)by 44.68% and 43.48%,indicating that the CEEMDAN algorithm had better denoising behavior than the other two algorithms.In addition,the results demonstrated that different denoising methods had little influence on estimating the annual vertical deformation velocity.The extraction of periodic signals showed that more components were retained by using the CEEMDAN algorithm than the EMD algorithm,which indicated that the CEEMDAN algorithm had advantages over frequency aliasing.In conclusion,the CEEMDAN algorithm was recommended for processing the GNSS height time series to analyze the vertical deformation due to its excellent features of denoising and the extraction of periodic signals.

Atomistic simulations of the surface severe plastic deformation-induced grain refinement in polycrystalline magnesium:The effect of processing parameters

Magnesium(Mg)based alloys are promising candidates for many applications,but their untreated surfaces usually have low strength and hardness.In this study,a single point diamond turning(SPDT)technique was applied to refine the grain size and improve the mechanical properties of the surface layers of Mg-Li alloys.By refining grains in the topmost layer to the nanometer scale(~60 nm),the surface hardness was found to be enhanced by approximately 60%.The atomic plastic deformation process during the SPDT was then studied by the real-time atomistic molecular dynamics(MD)simulations.A series of MD simulations with different combinations of parameters,including rake angle,cutting speed and cutting depth,were conducted to understand their influences on the microstructural evolution and associated plastic deformation mechanisms on the surface layer of the workpieces.The MD simulation results suggest that using increased rake angle,cutting speed and cutting depth can help to achieve better grain refinement.These simulation results,which provide atomic-level details of the deformation mechanism,can assist the parameter design for the SPDT techniques to achieve the high-performance heterogeneous nanostructured materials.

A comprehensive review on the processing-property relationships of laser strengthened magnesium

Among the existing series of softer metals,magnesium(Mg)has attracted much attention due to its impressive strength-to-weight ratio.However,due to its ease of deformability,Mg tends to suffer from rapid degradation in a wide variety of abrasive and electrochemical environments.One method of improving its surface properties is through surface modification techniques.Among the existing techniques,laser shock peening(LSP)has been one of the most widely utilized processes due to its surface-hardening-like effects.Despite this understanding,a comprehensive review has yet to exist that encapsulates the strengthening mechanism of LSP for Mg and its influence in degradation environments.This review aims to encapsulate the existing research around the LSP field for Mg.Specifically,an understanding of the surface-strengthening effects in relation to its mechanical,tribological,corrosion,and tribo-corrosion characteristics is elucidated.Additionally,the feasibility of LSP for Mg materials in critical industries is also discussed.Through this work,a novel understanding of LSP for Mg can be understood,which can provide a future direction for research in this field.

Experimental study of evaporation of sessile water droplet on PDMS surfaces

Evaporation of sessile water droplet on polydimethylsiloxane （PDMS） surfaces with three different curing ratios （5：1, 10：1, and 20：1） was experimentally investigated in this paper. We show that the constant contact radius （CCR） evaporation on surface with high curing ratio lasts longer than that with low curing ratio. We also measured Young＇s moduli of PDMS films by using atomic force microscopy （AFM） and simulated surface deformation of PDMS films induced by sessile water droplet. With increasing curing ratio of PDMS film, Young＇s modulus of PDMS film is getting lower, and then there will be larger surface deformation and more elastic stored energy. Since such energy acts as a barrier to keep the three-phase contact line pinned, thus it will result in longer CCR evaporation on PDMS surface with higher curing ratio.

Substrate elastic deformation due to vertical component of liquid-vapor interfacial tension

Young＇s equation is a fundamental equation in capillarity and wetting, which reflects the balance of the horizontal components of the three interracial tensions with the contact angle （CA）. However, it does not consider the vertical component of the liquid-vapor interracial tension （VCLVIT）. It is now well understood that the VCLVIT causes the elastic deformation of the solid substrate, which plays a significant role in the fabrication of the microfluidic devices because of the wide use of the soft materials. In this paper, the theoretical, experimental, and numerical aspects of the problem are reviewed. The effects of the VCLVIT-induced surface deformation on the wetting and spreading, the deflection of the microcantilever, and the elasto.capillarity and electro- elasto.capillarity are discussed. Besides a brief review on the historical development and the recent advances, some suggestions on the future research are also provided. Key words

Numerical Simulation of Sloshing with Large Deforming Free Surface by MPS-LES Method

Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

Quick Surface Treatment of AZ31B by AC Micro-arc Oxidation

In order to explore an effective way to shorten treatment time and enhance the quality of treatment coating, AC micro-arc oxidation was conducted to treat the surface of AZ31 B deformation magnesium alloy in KF＋KOH treatment solution. The infl uences of micro-arc oxidation parameters such as concentration of KF, concentration of KOH, output voltage of booster, temperature of treatment solution, and treatment time on treatment coating thickness were raveled out under different conditions. The structure and composition of treatment coating were determined, the growth mechanism of treatment coating was discussed, and the quick surface treatment technology for compact treatment coating with maximum thickness was developed. The experimental results show that： A maximum 33 μm-thick compact treatment coating, consisting of MgF2 and MgO mainly, can be formed on AZ31 B in 112 s under the conditions of 1 132 g/L KF, 382 g/L KOH, 66 V for output voltage of booster and 34 ℃ of treatment solution which were optimized by a genetic algorithm from the model established by artifi cial neural networks. There are no ＂crater-shaped＂ pores in this treatment coating as the heat shock resulting from the smooth variation of AC sinusoidal voltage is far smaller than that of the rigidly varied DC or pulse current. The treatment time is only one sixth of that adopted in the other surface treatment technology at best, principally for the reason that the coating can always grow irrespective of the electric potential of AZ31 B. This investigation lays a fi rm foundation for the extensive application of magnesium alloy.

Development and Deformation Characteristics of Large Ancient Landslides in the Intensely Hazardous Xiongba-Sela Section of the Jinsha River, Eastern Tibetan Plateau, China

The upstream Jinsha River,located in the eastern Tibetan Plateau,has been experiencing intense geological hazards characterized by a high density of ancient landslides,significant deformation and reactivation challenges.In this study,remote sensing interpretation,field investigations,and Small Baseline Subset Interferometric Synthetic Aperture Radar(SBAS-InSAR)technologies have been employed.Along a 17 km stretch of the Jinsha River,specifically in the Xiongba-Sela segment,16 large-scale ancient landslides were identified,9 of which are currently undergoing creeping deformation.Notably,the Sela and Xiongba ancient landslides exhibit significant deformation,with a maximum deformation rate of-192 mm/yr,indicating a high level of sliding activity.The volume of the Sela ancient landslide is estimated to be 1.8×108 to 4.5×108 m3,and characterized by extensive fissures and long-term creeping deformation.The SBAS-InSAR results revealed significant spatial variations in the deformation of the Sela ancient landslide,generally displaying two secondary zones of intense deformation,and landslide deformation exhibits nonlinear behavior with time.Between January 2016 and February 2022,Zone III1 on the southwest side of the Sela ancient landslide,experienced a maximum cumulative deformation of-857 mm,with a maximum deformation rate of-108 mm/yr.Zone III2,on the northeast side of the Sela ancient landslide,the maximum cumulative deformation was-456 mm,with a maximum deformation rate of-74 mm/yr;among these,the H2 and H4 secondary bodies on the south side of III1 are in the accelerative deformation stage and at the Warn warning level.We propose that the large-scale flood and debris flow disasters triggered by the Baige landslide-dammed lake-dam broken disaster chain in Tibetan Plateau during October and November 2018 caused severe erosion at the foot of downstream slopes.This far-field triggering effect accelerated the creep of the downstream ancient landslides.Consequently,the deformation rate of Zone III2 of the Sela ancient landslide increased by 6 to 8 times,exhibiting traction-type style reactivation.This heightened activity raises concerns about the potential for large-scale or overall reactivation of the landslide,posing a risk of damming the Jinsha River and initiating a dam-break disaster chain.Our research on the reactivation characteristics and mechanisms of large ancient landslides in high deep-cut valleys provides valuable guidance for geological hazard investigation and risk prevention.

Numerical Simulation of Current Density Distribution in Keyhole Double-Sided Arc Welding

In the double-sided arc welding system (DSAW) composing of PAW+TIG arcs, the PAW arc is guided by the TIG arc so that the current mostly flows through the direction of the workpiece thickness and the penetration is greatly improved. To analyze the current density distribution in DSAW is beneficial to understanding of this process. Considering all kinds of dynamic factors acting on the weldpool, this paper discusses firstly the surface deformation of the weldpool and the keyhole formation in PAW+TIG DSAW process on the basis of the magnetohydrodynamic theory and variation principles. Hence, a model of the current density distribution is developed. Through numerical simulation, the current density distribution in PAW+TIG DSAW process is quantitatively analyzed. It shows that the minimal radius of keyhole formed in PAW+TIG DSAW process is 0.5 mm and 89.5 percent of current flows through the keyhole.