Uncertainties of landslide susceptibility prediction: Influences of random errors in landslide conditioning factors and errors reduction by low pass filter method

In the existing landslide susceptibility prediction(LSP)models,the influences of random errors in landslide conditioning factors on LSP are not considered,instead the original conditioning factors are directly taken as the model inputs,which brings uncertainties to LSP results.This study aims to reveal the influence rules of the different proportional random errors in conditioning factors on the LSP un-certainties,and further explore a method which can effectively reduce the random errors in conditioning factors.The original conditioning factors are firstly used to construct original factors-based LSP models,and then different random errors of 5%,10%,15% and 20%are added to these original factors for con-structing relevant errors-based LSP models.Secondly,low-pass filter-based LSP models are constructed by eliminating the random errors using low-pass filter method.Thirdly,the Ruijin County of China with 370 landslides and 16 conditioning factors are used as study case.Three typical machine learning models,i.e.multilayer perceptron(MLP),support vector machine(SVM)and random forest(RF),are selected as LSP models.Finally,the LSP uncertainties are discussed and results show that:(1)The low-pass filter can effectively reduce the random errors in conditioning factors to decrease the LSP uncertainties.(2)With the proportions of random errors increasing from 5%to 20%,the LSP uncertainty increases continuously.(3)The original factors-based models are feasible for LSP in the absence of more accurate conditioning factors.(4)The influence degrees of two uncertainty issues,machine learning models and different proportions of random errors,on the LSP modeling are large and basically the same.(5)The Shapley values effectively explain the internal mechanism of machine learning model predicting landslide sus-ceptibility.In conclusion,greater proportion of random errors in conditioning factors results in higher LSP uncertainty,and low-pass filter can effectively reduce these random errors.

The Factors of the Failure in the Treatment of Congenital Equinus Clubfoot Varus by the Ponseti Method at the Reference Health Care Center (CSREF) of Commune III of Bamako

Clubfoot malformation is the most common serious congenital anomaly affecting the foot in children. Its treatment by the Ponseti method is simple, profitable and widely used in the world. Objective: The objective of this study was to identify the factors of the failure of its treatment by the Ponseti method. Material and Method: We carried out a retrospective and descriptive study of cases of congenital equinus clubfoot varus at the Reference Health Care Center of Commune III of Bamako over 26 months from September 2020 to November 2022. Data were treated with the utmost anonymity. Result: This study was performed on 44 children seen for clubfoot: male (68%) and female (32%), with a sex ratio of 2.1. We obtained 13 cases of recidivism including 7 boys and 6 girls. We found 21 cases of unilateral and 23 bilateral;among which 9 recurrences were found against 4 in the unilateral forms. 85% of recurrences did not have good adherence to the splint and 62% did not come regularly for follow-up consultation. We obtained 33 children with idiopathic clubfeet (75%) with a recurrence of 24%, and 7 children with secondary clubfeet with 71 % recurrence. There was no recurrence in the postural type. Among the recurrences, 38.5% started treatment between 1 and 6 months, 23.1% from 0 to 1 month and 15.4% at 2 years and more. 85% of recurrences had a Pirani score between 4.5 to 6 at the start of treatment and 15% at a score of 2.5 to 4. Conclusion: The factors of the failure of the Ponseti method are mainly non-compliance with treatment, secondary clubfeet, and a high Pirani score at the start of treatment.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

Quantifying the impact of earthquakes and geological factors on spatial heterogeneity of debris-flow prone areas:A case study in the Hengduan Mountains

Understanding the spatial heterogeneity of debris-flow-prone areas holds significant implications for regional risk management, particularly in seismically active regions with geological faults. Despite the significance of this knowledge, a comprehensive quantification of the influence of regional topographical and geological factors on the spatial heterogeneity of debris-flow-prone areas has been lacking. This study selected the Hengduan Mountains, an earthquake-prone region characterized by diverse surface conditions and complex landforms, as a representative study area. An improved units zoning and objective factors identification methodology was employed in earthquake and fault analysis to assess the impact of seismic activity and geological factors on spatial heterogeneity of debrisflow prone areas. Results showed that the application of GIS technology with hydrodynamic intensity and geographical units analysis can effectively analyze debris-flow prone areas. Meanwhile, earthquake and fault zones obviously increase the density of debrisflow prone catchments and make them unevenly distributed. The number of debris-flow prone areas shows a nonlinear variation with the gradual increase of geomorphic factor value. Specifically, the area with 1000 m-2500 m elevation difference, 25°-30° average slope, and 0.13-0.15 land use index is the most favorable conditions for debris-flow occurrence;The average annual rainfall from 600 to 1150 mm and landslides gradient from 16° to 35° are the main causal factors to trigger debris flow. Our study sheds light on the quantification of spatial heterogeneity in debris flow-prone areas in earthquake-prone regions, which can offer crucial support for post-debris flow risk management strategies.

Evaluation of soil erosion vulnerability in Hubei Province of China using RUSLE model and combination weighting method

Soil erosion has been recognized as a critical environmental issue worldwide.While previous studies have primarily focused on watershed-scale soil erosion vulnerability from a natural factor perspective,there is a notable gap in understanding the intricate interplay between natural and socio-economic factors,especially in the context of spatial heterogeneity and nonlinear impacts of human-land interactions.To address this,our study evaluates the soil erosion vulnerability at a provincial scale,taking Hubei Province as a case study to explore the combined effects of natural and socio-economic factors.We developed an evaluation index system based on 15 indicators of soil erosion vulnerability:exposure,sensitivity,and adaptability.In addition,the combination weighting method was applied to determine index weights,and the spatial interaction was analyzed using spatial autocorrelation,geographical temporally weighted regression and geographical detector.The results showed an overall decreasing soil erosion intensity in Hubei Province during 2000 and 2020.The soil erosion vulnerability increased before 2000 and then.The areas with high soil erosion vulnerability were mainly confined in the central and southern regions of Hubei Province(Xiantao,Tianmen,Qianjiang and Ezhou)with obvious spatial aggregation that intensified over time.Natural factors(habitat quality index)had negative impacts on soil erosion vulnerability,whereas socio-economic factors(population density)showed substantial spatial variability in their influences.There was a positive correlation between soil erosion vulnerability and erosion intensity,with the correlation coefficients ranging from-0.41 and 0.93.The increase of slope was found to enhance the positive correlation between soil erosion vulnerability and intensity.

Failure evaluation and control factor analysis of slope block instability along traffic corridor in Southeastern Tibet

The instability of slope blocks occurred frequently along traffic corridor in Southeastern Tibet(TCST),which was primarily controlled by the rock mass structures.A rapid method evaluating the control effects of rock mass structures was proposed through field statistics of the slopes and rock mass structures along TCST,which combined the stereographic projection method,modified M-JCS model,and limit equilibrium theory.The instabilities of slope blocks along TCST were then evaluated rapidly,and the different control factors of instability were analyzed.Results showed that the probabilities of toppling(5.31%),planar(16.15%),and wedge(35.37%)failure of slope blocks along TCST increased sequentially.These instability modes were respectively controlled by the anti-dip joint,the joint parallel to slope surface with a dip angle smaller than the slope angle(singlejoint),and two groups of joints inclined out of the slope(double-joints).Regarding the control effects on slope block instability,the stabilization ability of doublejoints(72.7%),anti-dip joint(67.4%),and single-joint(57.6%)decreased sequentially,resulting in different probabilities of slope block instability.Additionally,nearby regional faults significantly influenced the joints,leading to spatial heterogeneity and segmental clustering in the stabilization ability provided by joints to the slope blocks.Consequently,the stability of slope blocks gradually weakened as they approached the fault zones.This paper can provide guidance and assistance for investigating the development characteristics of rock mass structures and the stability of slope blocks.

Accelerating Factor Xa inhibitor discovery with a de novo drug design pipeline

Small-molecule drugs are essential for maintaining human health. The objective of this study is to identify a molecule that can inhibit the Factor Xa protein and be easily procured. An optimization-based de novo drug design framework, Drug CAMD, that integrates a deep learning model with a mixed-integer nonlinear programming model is used for designing drug candidates. Within this framework, a virtual chemical library is specifically tailored to inhibit Factor Xa. To further filter and narrow down the lead compounds from the designed compounds, comprehensive approaches involving molecular docking,binding pose metadynamics(BPMD), binding free energy calculations, and enzyme activity inhibition analysis are utilized. To maximize efficiency in terms of time and resources, molecules for in vitro activity testing are initially selected from commercially available portions of customized virtual chemical libraries. In vitro studies assessing inhibitor activities have confirmed that the compound EN300-331859shows potential Factor Xa inhibition, with an IC_(50)value of 34.57 μmol·L^(-1). Through in silico molecular docking and BPMD, the most plausible binding pose for the EN300-331859-Factor Xa complex are identified. The estimated binding free energy values correlate well with the results obtained from biological assays. Consequently, EN300-331859 is identified as a novel and effective sub-micromolar inhibitor of Factor Xa.

A Calculation Method of Double Strength Reduction for Layered Slope Based on the Reduction of Water Content Intensity

The calculation of the factor of safety(FOS)is an important means of slope evaluation.This paper proposed an improved double strength reductionmethod(DRM)to analyze the safety of layered slopes.The physical properties of different soil layers of the slopes are different,so the single coefficient strength reduction method(SRM)is not enough to reflect the actual critical state of the slopes.Considering that the water content of the soil in the natural state is the main factor for the strength of the soil,the attenuation law of shear strength of clayey soil changing with water content is fitted.This paper also establishes the functional relationship between different reduction coefficients.Then,a USDFLD subroutine is programmed using the secondary development function of finite element software.Controlling the relationship between field variables and calculation time realizes double strength reduction applicable to the layered slope.Finally,by comparing the calculation results of different examples,it is proved that the stress and displacement distribution of the critical slope state obtained by the improved method is more realistic,and the calculated safety factor is more reliable.The newly proposedmethod considers the difference of intensity attenuation between different soil layers under natural conditions and avoids the disadvantage of the strength reduction method with uniform parameters,which provides a new idea and method for stability analysis of layered and complex slopes.

Factors affecting land dissection density in geomorphological regions of China

Land dissection density(LDD)provides morphological evidence regarding prior intense soil erosion and quantifies the distribution of land dissections.A comprehensive understanding of the potential factors influencing the spatial pattern and value of the LDD is vital in geological disasters,soil erosion,and other related domains.Land dissection phenomena in China affects large areas with different morphological,pedological,and climatic characteristics.Prior studies have focused on the potential factors influencing the LDD at a watershed scale.However,these results are insufficient to reflect the status quo of dissection development and its primary influencing factors on a national scale.LDD’s spatial patterns and the dominant factors at a regional scale in millions of square kilometers remain to be ascertained.This study used the geomorphon-based method and the geographical detector model to quantify the spatial pattern of LDD over China and identify the dominant factors affecting this pattern in China’s six first-order geomorphological regions(GR1~GR6).The results yield the following findings:(1)LDD in China ranges from 0~4.55 km/km^(2),which is larger in central and eastern regions than in other regions of China;(2)dominant factors and their dominant risk subcategories vary with each geomorphological region’s primary internal and external forces;(3)the influence of natural factors is more significant on the large regional scale in millions of square kilometers compared to anthropogenic factors;relief degree of land surface(RDLS)is dominant in GR1,GR2,and GR5;the slope is dominant in GR6,soil type is dominant in GR3 and GR4,and lithology plays a critical role in the dominant interactions of GR3,GR4,and GR6;(4)the interactions between factors on LDD’s spatial pattern have a more significant effect than individual factors.

Optimal Shape Factor and Fictitious Radius in the MQ-RBF:Solving Ill-Posed Laplacian Problems

To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).

Effects of Different Grinding Methods on the Quality of Soybean Bean Milk

[Objectives] This study was conducted to improve the nutritional value of soybean milk, enrich the variety and taste of soybean milk, and find healthy food that is more conducive to people s nutritional needs. [Methods] Whole soybean milk was prepared by grinding with a grinding wheel at a low concentration (low-concentration grinding) and a stainless steel mill at a high concentration (high-concentration grinding). The sensory, physical and chemical characteristics and anti-nutritional factors of whole soybean milk produced by different grinding methods were studied. [Results] Compared with low-concentration grinding, the protein content in soybean milk prepared by high-concentration grinding increased by 24%, and the dietary fiber content increased by 74.7%. Before and after high-pressure homogenization, the particle size D(4, 3) of soybean milk prepared by low-concentration grinding was 212.1 and 93.59 μm, respectively, and the particle size D(4, 3) of soybean milk prepared by high-concentration grinding was 134.0 and 64.64 μm, respectively. The trypsin inhibitor activity and phytic acid content of soybean milk prepared by high-concentration grinding were significantly lower than those of soybean milk prepared by low-concentration grinding. [Conclusions] This study improves the diet structure of the broad masses of people, strengthens people s physique, and provides a new idea for the implementation and development of China s "Soybean Action Programme".

Enhancing Stress Intensity Factor Reduction in Cracks Originating from a Circular Hole in a Rectangular Plate under Uniaxial Stress through Piezoelectric Actuation

Circular holes are commonly employed in engineering designs;however, they often serve as locations where cracks initiate and propagate. This paper explores a novel approach to structural repair by utilizing piezoelectric actuators. The primary focus of this study is to investigate the influence of an adhesively bonded piezoelectric actuator patch placed above a circular hole on the stress intensity factor (SIF) in an aluminium plate. The plate is subjected to uniaxial tensile stress, while the piezoelectric actuator is excited with varying voltage levels. The analysis is conducted using the finite element method (FEM), a powerful numerical technique for simulating complex structures. The study assesses the stress distribution and employs the SIF as an adequate criterion for evaluating the impact of different patch configurations. The results indicate a strong correlation between the applied voltage and the SIF. Whether the SIF increases or decreases depends on the polarization of the piezoelectric actuator. Particularly noteworthy is the finding that rectangular patches in a horizontal orientation significantly reduce the SIF compared to other patch geometries. Moreover, double-sided patches exhibit a pronounced decrease in the SIF compared to single-sided patches. In summary, this research underscores the potential of piezoelectric actuators in mitigating stress intensity in structures with circular hole with crack initiation. It offers valuable insights into the influence of applied voltage, patch geometry, and patch placement on the SIF, thereby contributing to developing effective strategies for enhancing structural integrity.

Integrated interpretation of dual frequency induced polarization measurement based on wavelet analysis and metal factor methods

In mineral exploration, the apparent resistivity and apparent frequency （or apparent polarizability） parameters of induced polarization method are commonly utilized to describe the induced polarization anomaly. When the target geology structure is significantly complicated, these parameters would fail to reflect the nature of the anomaly source, and wrong conclusions may be obtained. A wavelet approach and a metal factor method were used to comprehensively interpret the induced polarization anomaly of complex geologic bodies in the Adi Bladia mine. Db5 wavelet basis was used to conduct two-scale decomposition and reconstruction, which effectively suppress the noise interference of greenschist facies regional metamorphism and magma intrusion, making energy concentrated and boundary problem unobservable. On the basis of that, the ore-induced anomaly was effectively extracted by the metal factor method.

Estimation of quality factors by energy ratio method

The quality factor Q, which reflects the energy attenuation of seismic waves in subsurface media, is a diagnostic tool for hydrocarbon detection and reservoir characterization. In this paper, we propose a new Q extraction method based on the energy ratio before and after the wavelet attenuation, named the energy-ratio method（ERM）. The proposed method uses multipoint signal data in the time domain to estimate the wavelet energy without invoking the source wavelet spectrum, which is necessary in conventional Q extraction methods, and is applicable to any source wavelet spectrum; however, it requires high-precision seismic data. Forward zero-offset VSP modeling suggests that the ERM can be used for reliable Q inversion after nonintrinsic attenuation（geometric dispersion, refl ection, and transmission loss） compensation. The application to real zero-offset VSP data shows that the Q values extracted by the ERM and spectral ratio methods are identical, which proves the reliability of the new method.

METHOD TO CALCULATE BENDING CENTER AND STRESS INTENSITY FACTORS OF CRACKED CYLINDER UNDER SAINT_VENANT BENDING

Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross_section is not thin_walled, but of small torsion rigidity is proposed. Some numerical examples are given.

Q-factor estimation in CMP gather and the continuous spectral ratio slope method

The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.

CFD-Based Method of Determining Form Factor k for Different Ship Types and Different Drafts

The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low-speed model test. However, this method is problematic for ships with bulbous bows and transom. In this article, a Computational Fluid Dynamics(CFD)-based method is introduced to obtain k for different type of ships at different drafts, and a comparison is made between the CFD method and the model test. The results show that the CFD method produces reasonable k values. A grid generating method and turbulence model are briefly discussed in the context of obtaining a consistent k using CFD.

METHOD TO CALCULATE STRESS INTENSITY FACTOR OF V-NOTCH IN BI-MATERIALS

Based on Zak＇s stress function, the eigen-equation of stress singularity ofbi-materials with a V-notch was obtained. A new definition of stress intensity factor for a perpendicular interfacial V-notch of bi-material was put forward. The effects of shear modulus and Poisson＇s ratio of the matrix material and attaching material on eigen-values were analyzed. A generalized expression for calculating/（i of the perpendicular V-notch of bi-materials was obtained by means of stress extrapolation. Effects of notch depth, notch angle and Poisson＇s ratio of materials on the singular stress field near the tip of the V-notch were analyzed systematically with numerical simulations. As an example, a finite plate with double edge notches under uniaxial uniform tension was calculated by the method presented and the influence of the notch angle and Poisson＇s ratio on the stress singularity near the tip of notch was obtained.

DYNAMIC CHARACTERISTIC ANALYSIS OF FUZZY-STOCHASTIC TRUSS STRUCTURES BASED ON FUZZY FACTOR METHOD AND RANDOM FACTOR METHOD

A new fuzzy stochastic finite element method based on the fuzzy factor method and random factor method is given and the analysis of structural dynamic characteristic for fuzzy stochastic truss structures is presented. Considering the fuzzy randomness of the structural physical parameters and geometric dimensions simultaneously, the structural stiffness and mass matrices axe constructed based on the fuzzy factor method and random factor method; from the Rayleigh＇s quotient of structural vibration, the structural fuzzy random dynamic characteristic is obtained by means of the interval arithmetic; the fuzzy numeric characteristics of dynamic characteristic axe then derived by using the random variable＇s moment function method and algebra synthesis method. Two examples axe used to illustrate the validity and rationality of the method given. The advantage of this method is that the effect of the fuzzy randomness of one of the structural parameters on the fuzzy randomness of the dynamic characteristic can be reflected expediently and objectively.

Rutting influencing factors and prediction model for asphalt pavements based on the factor analysis method

To clarify the importance of various influencing factors on asphalt pavement rutting deformation and determine a screening method of model indicators,the data of the RIOHTrack full-scale track were examined using the factor analysis method(FAM).Taking the standard test pavement structure of RIOHTrack as an example,four rutting influencing factors from different aspects were determined through statistical analysis.Furthermore,the common influencing factors among the rutting influencing factors were studied based on FAM.Results show that the common factor can well characterize accumulative ESALs,center-point deflection,and temperature,besides humidity,which indicates that these three influencing factors can have an important impact on rutting.Moreover,an empirical rutting prediction model was established based on the selected influencing factors,which proved to exhibit high prediction accuracy.These analysis results demonstrate that the FAM is an effective screening method for rutting prediction model indicators,which provides a reference for the selection of independent model indicators in other rutting prediction model research when used in other areas and is of great significance for the prediction and control of rutting distress.