A method to interpret fracture aperture of rock slope using adaptive shape and unmanned aerial vehicle multi-angle nap-of-the-object photogrammetry

The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods are inadequate for collecting data on high-steep rock slopes in complex mountainous regions.This study establishes a high-resolution three-dimensional model of a rock slope using unmanned aerial vehicle(UAV)multi-angle nap-of-the-object photogrammetry to obtain edge feature points of fractures.Fracture opening morphology is characterized using coordinate projection and transformation.Fracture central axis is determined using vertical measuring lines,allowing for the interpretation of aperture of adaptive fracture shape.The feasibility and reliability of the new method are verified at a construction site of a railway in southeast Tibet,China.The study shows that the fracture aperture has a significant interval effect and size effect.The optimal sampling length for fractures is approximately 0.5e1 m,and the optimal aperture interpretation results can be achieved when the measuring line spacing is 1%of the sampling length.Tensile fractures in the study area generally have larger apertures than shear fractures,and their tendency to increase with slope height is also greater than that of shear fractures.The aperture of tensile fractures is generally positively correlated with their trace length,while the correlation between the aperture of shear fractures and their trace length appears to be weak.Fractures of different orientations exhibit certain differences in their distribution of aperture,but generally follow the forms of normal,log-normal,and gamma distributions.This study provides essential data support for rock and slope stability evaluation,which is of significant practical importance.

Investigation of the block toppling evolution of a layered model slope by centrifuge test and discrete element modeling

Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.

Dynamic response and failure process of horizontal-layered fractured structure rock slope under strong earthquake

Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the dynamic response,failure mode,and spectral characteristics of rock slope with HLFS under strong earthquake conditions were investigated based on the large-scale shaking table model test.On this basis,multiple sets of numerical calculation models were further established by UDEC discrete element program.Five influencing factors were considered in the parametric study of numerical simulations,including slope height,slope angle,bedding-plane spacing and secondary joint spacing as well as bedrock dip angle.The results showed that the failure process of rock slope with HLFS under earthquake action is mainly divided into four phases,i.e.,the tensile crack of the slope shoulder joints and shear dislocation at the top bedding plane,the extension of vertical joint cracks and increase of shear displacement,the formation of step-through sliding surfaces and the instability,and finally collapse of fractured rock mass.The acceleration response of slopes exhibits elevation amplification effect and surface effect.Numerical simulations indicate that the seismic stability of slopes with HLFS exhibits a negative correlation with slope height and angle,but a positive correlation with bedding-plane spacing,joint spacing,and bedrock dip angle.The results of this study can provide a reference for seismic stability evaluation of weathered rock slopes.

Effect of slope angle on fractured rock masses under combined influence of variable rainfall infiltration and excavation unloading

Intense precipitation infiltration and intricate excavation processes are crucial factors that impact the stability and security of towering and steep rock slopes within mining sites.The primary aim of this research was to investigate the progression of cumulative failure within a cracked rock formation,considering the combined effects of precipitation and excavation activities.The study was conducted in the Huangniuqian eastern mining area of the Dexing Copper Mine in Jiangxi Province,China.An engineering geological investigation was conducted,a physical model experiment was performed,numerical calculations and theoretical analysis were conducted using the matrix discrete element method(Mat-DEM),and the deformation characteristics and the effect of the slope angle of a fractured rock mass under different scenarios were examined.The failure and instability mechanisms of the fractured rock mass under three slope angle models were analyzed.The experimental results indicate that as the slope angle increases,the combined effect of rainfall infiltration and excavation unloading is reduced.A novel approach to simulating unsaturated seepage in a rock mass,based on the van Genuchten model(VGM),has been developed.Compared to the vertical displacement observed in a similar physical experiment,the average relative errors associated with the slope angles of 45,50,and 55were 2.094%,1.916%,and 2.328%,respectively.Accordingly,the combined effect of rainfall and excavation was determined using the proposed method.Moreover,the accuracy of the numerical simulation was validated.The findings contribute to the seepage field in a meaningful way,offering insight that can inform and enhance existing methods and theories for research on the underlying mechanism of ultra-high and steep rock slope instability,which can inform the development of more effective risk management strategies.

Evaluation of slope stability through rock mass classification and kinematic analysis of some major slopes along NH-1A from Ramban to Banihal, North Western Himalayas

The network of Himalayan roadways and highways connects some remote regions of valleys or hill slopes,which is vital for India’s socio-economic growth.Due to natural and artificial factors,frequency of slope instabilities along the networks has been increasing over last few decades.Assessment of stability of natural and artificial slopes due to construction of these connecting road networks is significant in safely executing these roads throughout the year.Several rock mass classification methods are generally used to assess the strength and deformability of rock mass.This study assesses slope stability along the NH-1A of Ramban district of North Western Himalayas.Various structurally and non-structurally controlled rock mass classification systems have been applied to assess the stability conditions of 14 slopes.For evaluating the stability of these slopes,kinematic analysis was performed along with geological strength index(GSI),rock mass rating(RMR),continuous slope mass rating(CoSMR),slope mass rating(SMR),and Q-slope in the present study.The SMR gives three slopes as completely unstable while CoSMR suggests four slopes as completely unstable.The stability of all slopes was also analyzed using a design chart under dynamic and static conditions by slope stability rating(SSR)for the factor of safety(FoS)of 1.2 and 1 respectively.Q-slope with probability of failure(PoF)1%gives two slopes as stable slopes.Stable slope angle has been determined based on the Q-slope safe angle equation and SSR design chart based on the FoS.The value ranges given by different empirical classifications were RMR(37-74),GSI(27.3-58.5),SMR(11-59),and CoSMR(3.39-74.56).Good relationship was found among RMR&SSR and RMR&GSI with correlation coefficient(R 2)value of 0.815 and 0.6866,respectively.Lastly,a comparative stability of all these slopes based on the above classification has been performed to identify the most critical slope along this road.

Efficient slope reliability and sensitivity analysis using quantile-based first-order second-moment method

This paper introduces a novel approach for parameter sensitivity evaluation and efficient slope reliability analysis based on quantile-based first-order second-moment method(QFOSM).The core principles of the QFOSM are elucidated geometrically from the perspective of expanding ellipsoids.Based on this geometric interpretation,the QFOSM is further extended to estimate sensitivity indices and assess the significance of various uncertain parameters involved in the slope system.The proposed method has the advantage of computational simplicity,akin to the conventional first-order second-moment method(FOSM),while providing estimation accuracy close to that of the first-order reliability method(FORM).Its performance is demonstrated with a numerical example and three slope examples.The results show that the proposed method can efficiently estimate the slope reliability and simultaneously evaluate the sensitivity of the uncertain parameters.The proposed method does not involve complex optimization or iteration required by the FORM.It can provide a valuable complement to the existing approximate reliability analysis methods,offering rapid sensitivity evaluation and slope reliability analysis.

Reliability analysis of slope stability by neural network,principal component analysis,and transfer learning techniques

The prediction of slope stability is considered as one of the critical concerns in geotechnical engineering.Conventional stochastic analysis with spatially variable slopes is time-consuming and highly computation-demanding.To assess the slope stability problems with a more desirable computational effort,many machine learning(ML)algorithms have been proposed.However,most ML-based techniques require that the training data must be in the same feature space and have the same distribution,and the model may need to be rebuilt when the spatial distribution changes.This paper presents a new ML-based algorithm,which combines the principal component analysis(PCA)-based neural network(NN)and transfer learning(TL)techniques(i.e.PCAeNNeTL)to conduct the stability analysis of slopes with different spatial distributions.The Monte Carlo coupled with finite element simulation is first conducted for data acquisition considering the spatial variability of cohesive strength or friction angle of soils from eight slopes with the same geometry.The PCA method is incorporated into the neural network algorithm(i.e.PCA-NN)to increase the computational efficiency by reducing the input variables.It is found that the PCA-NN algorithm performs well in improving the prediction of slope stability for a given slope in terms of the computational accuracy and computational effort when compared with the other two algorithms(i.e.NN and decision trees,DT).Furthermore,the PCAeNNeTL algorithm shows great potential in assessing the stability of slope even with fewer training data.

Numerical analysis of downward progressive landslides in long natural slopes with sensitive clay

Landslides occurring in sensitive clay often result in widespread destruction,posing a significant risk to human lives and property due to the substantial decrease in undrained shear strength during deformation.Assessing the consequences of these landslides is challenging and necessitates robust numerical methods to comprehensively investigate their failure mechanisms.While studies have extensively explored upward progressive landslides in sensitive clays,understanding downward progressive cases remains limited.In this study,we utilised the nodal integration-based particle finite element method(NPFEM)with a nonlinear strain-softening model to analyse downward progressive landslides in sensitive clay on elongated slopes,induced by surcharge loads near the crest.We focused on elucidating the underlying failure mechanisms and evaluating the effects of different soil parameters and strainsoftening characteristics.The simulation results revealed the typical pattern for downward landslides,which typically start with a localised failure in proximity to the surcharge loads,followed by a combination of different types of failure mechanisms,including single flow slides,translational progressive landslides,progressive flow slides,and spread failures.Additionally,inclined shear bands occur within spread failures,often adopting distinctive ploughing patterns characterised by triangular shapes.The sensitive clay thickness at the base,the clay strength gradient,the sensitivity,and the softening rate significantly influence the failure mechanisms and the extent of diffused displacement.Remarkably,some of these effects mirror those observed in upward progressive landslides,underscoring the interconnectedness of these phenomena.This study contributes valuable insights into the complex dynamics of sensitive clay landslides,shedding light on the intricate interplay of factors governing their behaviour and progression.

Large-scale model testing of high-pressure grouting reinforcement for bedding slope with rapid-setting polyurethane

Bedding slope is a typical heterogeneous slope consisting of different soil/rock layers and is likely to slide along the weakest interface.Conventional slope protection methods for bedding slopes,such as retaining walls,stabilizing piles,and anchors,are time-consuming and labor-and energy-intensive.This study proposes an innovative polymer grout method to improve the bearing capacity and reduce the displacement of bedding slopes.A series of large-scale model tests were carried out to verify the effectiveness of polymer grout in protecting bedding slopes.Specifically,load-displacement relationships and failure patterns were analyzed for different testing slopes with various dosages of polymer.Results show the great potential of polymer grout in improving bearing capacity,reducing settlement,and protecting slopes from being crushed under shearing.The polymer-treated slopes remained structurally intact,while the untreated slope exhibited considerable damage when subjected to loads surpassing the bearing capacity.It is also found that polymer-cemented soils concentrate around the injection pipe,forming a fan-shaped sheet-like structure.This study proves the improvement of polymer grouting for bedding slope treatment and will contribute to the development of a fast method to protect bedding slopes from landslides.

Using negative exponential function to characterize built-up land density along slope gradient for 19 urban areas in China

In China,numerous cities are expanding into sloping land,yet the quantitative distribution patterns of urban built-up land density along the slope gradient remain unclear,limiting the understanding of sloping land urbanization.In this paper,a simple negative exponential function was presented to verify its applicability in 19 typical sloping urban areas in China.The function fits well for all case urban areas(R2≥0.951,p<0.001).The parameters of this function clearly describe two fundamental attributes:initial value a and decline rate b.Between 2000 and 2020,a tends to increase,while b tends to decrease in all urban areas,confirming the hypothesis of mutual promotion between flatland densification and sloping land expansion.Multiple regression analysis indicates that the built-up land density and the ruggedness of background land can explain 70.7%of a,while the average slope ratio of built-up land to background land,the built-up land density and the built-up land area can explain 82.1%of b.This work provides a quantitative investigative tool for distribution of urban built-up land density along slope gradient,aiding in the study of the globally increasing phenomenon of sloping land urbanization from a new perspective.

Effects of freeze-thaw cycles on sandstone in sunny and shady slopes

A growing rock engineering activity in cold regions is facing the threat of freeze-thaw(FT)weathering,especially in high mountains where the sunny-shady slope effects strongly control the difference in weathering behavior of rocks.In this paper,an investigation of the degradation of petrophysical characteristics of sandstone specimens subjected to FT cycle tests to simulate the sunny-shady slope effects is presented.To this aim,non-destructive and repeatable testing techniques including weight,ultrasonic waves,and nuclear magnetic resonance methods on standard specimens were performed.For the sunny slope specimens,accompanied by the enlargement of small pores,100 FT cycles caused a significant decrease in P-wave velocity with an average of 23%,but a consistent rise of 0.18%in mass loss,34%in porosity,67%in pore geometrical mean radius,and a remarkable 14.5-fold increase in permeability.However,slight changes with some abnormal trends in physical parameters of the shady slope specimens were observed during FT cycling,which can be attributed to superficial granular disaggregation and pore throat obstruction.Thermal shocks enhance rock weathering on sunny slopes during FT cycles,while FT weathering on shady slopes is restricted to the small pores and the superficial cover.These two factors are primarily responsible for the differences in FT weathering intensity between sunny and shady slopes.The conclusions derived from the interpretation of the experimental results may provide theoretical guidance for the design of slope-failure prevention measures and the selection of transportation routes in cold mountainous regions.

Large-scale toppling slope under water level fluctuation of reservoir:A case of Yunnan Province,China

Landslides induced by reservoir inundation are common in Southwest China,negatively influencing hydropower stations.TheWunonglong hydropower station dam was constructed in the upper reaches of the Lancang River,accordingly causing the water level at the Lajinshengu slope to increase by 30 m.A tension crack with a visible depth of 8 m was observed in the upper sector of the Lajinshengu slope after reservoir impoundment for 170 d.In the following days,numerous cracks appeared on the surface of the slope,and the maximum displacement of the slope reached 3.22 m.Then,a large-scale active deformation body within the Lajinshengu slope formed with an area of 2.62×10^(5)m^(2)and a volume of 1.65×10^(7)m^(3).Detailed field investigations,on-site monitoring,and centrifugal model tests were carried out to analyze the surface features,deformation characteristics,and failure mechanism of the Lajinshengu slope.The results show that the slope is an ancient landslide,divided into two parts(i.e.zone A and zone B)by the gully.Zone B is a traction landslide caused by the displacement of zone A.The longterm inundation weakens the soft rock at the slope foot,intensifying the toppling of bedrock and consequently triggering the sliding of the overburden in zone A.The failure mode of the Lajinshengu slope is a typical case of toppling-sliding failure,and the underlying rock toppling drives the overlying sliding.In addition,early identification methods for toppling deformation covered by overburdened soil were proposed based on monitoring data and deformation signs.

Modified Sadowski formula-based model for the slope shape amplification effect under multistage slope blasting vibration

Blasting operations,which are crucial to open-pit mine production due to their simplicity and efficiency,require precise control through accurate vibration velocity calculations.The conventional Sadowski formula mainly focuses on blast center distance but neglects the amplification effect of blasting vibration waves by terraced terrain,from which the calculated blasting vibration velocities are smaller than the actual values,affecting the safety of the project.To address this issue,our model introduces the influences of slope and time into Sadowski formula to measure safety through blast vibration displacement.In the northern section of the open-pit quartz mine in Jinchang City,Gansu Province,China,the data of a continuous blasting slope project are referred to.Our findings reveal a noticeable vibration amplification effect during blasting when a multi-stage slope platform undergoes a sudden cross-sectional change near the upper overhanging surface.The amplification vibration coefficient increases with height,while vibration waves within rocks decrease from bottom to top.Conversely,platforms without distinct crosssectional changes exhibit no pronounced amplification during blasting.In addition,the vibration intensity decreases with distance as the rock height difference change propagates.The results obtained by the proposed blast vibration displacement equation incorporating slope shape influence closely agree with real-world scenarios.According to Pearson correlation coefficient(PPMCC)analysis,the average accuracy rate of our model is 88.84%,which exceeds the conventional Sadowski formula(46.92%).

Influence of joint spacing and rock characteristics on the toppling stability of cut rock slope through a simplified limit equilibrium method

Toppling failure of rock mass/soil slope is an important geological and environmental problem.Clarifying its failure mechanism under different conditions has great significance in engineering.The toppling failure of a cutting slope occurred in a hydropower station in Kyushu,Japan illustrates that the joint characteristic played a significant role in the occurrence of rock slope tipping failure.Thus,in order to consider the mechanical properties of jointed rock mass and the influence of geometric conditions,a simplified analytical approach based on the limit equilibrium method for modeling the flexural toppling of cut rock slopes is proposed to consider the influence of the mechanical properties and geometry condition of jointed rock mass.The theoretical solution is compared with the numerical solution taking Kyushu Hydropower Station in Japan as one case,and it is found that the theoretical solution obtained by the simplified analysis method is consistent with the numerical analytical solution,thus verifying the accuracy of the simplified method.Meanwhile,the Goodman-Bray approach conventionally used in engineering practice is improved according to the analytical results.The results show that the allowable slope angle may be obtained by the improved Goodman-Bray approach considering the joint spacing,the joint frictional angle and the tensile strength of rock mass together.

Limit load and failure mechanisms of a vertical Hoek-Brown rock slope

The problem considered in this short note is the limit load determination of a vertical rock slope.The classical limit theorem is employed with the use of adaptive finite elements and nonlinear programming to determine upper and lower bound limit loads of a Hoek-Brown vertical rock slope.The objective function of the mathematical programming problem is such as to optimize a boundary load,which is known as the limit load,resembling the ultimate bearing capacity of a strip footing.While focusing on the vertical slope,parametric studies are carried out for several dimensionless ratios such as the dimensionless footing distance ratio,the dimensionless height ratio,and the dimensionless rock strength ratio.A comprehensive set of design charts is presented,and failure envelopes shown with the results explained in terms of three identified failure mechanisms,i.e.the face,the toe,and the Prandtl-type failures.These novel results can be used with great confidence in design practice,in particularly noting that the current industry-based design procedures for the presented problem are rarely found.

Stability analysis of loose accumulation slopes under rainfall:case study of a high‑speed railway in Southwest China

The high and steep slopes along a high-speed railway in the mountainous area of Southwest China are mostly composed of loose accumulations of debris with large internal pores and poor stability,which can easily induce adverse geological disasters under rainfall conditions.To ensure the smooth construction of the high-speed railway and the subsequent safe operation,it is necessary to master the stability evolution process of the loose accumulation slope under rainfall.This article simulates rainfall using the finite element analysis software’s hydromechanical coupling module.The slope stability under various rainfall situations is calculated and analysed based on the strength reduction method.To validate the simulation results,a field monitoring system is established to study the deformation characteristics of the slope under rainfall.The results show that rainfall duration is the key factor affecting slope stability.Given a constant amount of rainfall,the stability of the slope decreases with increasing duration of rainfall.Moreover,when the amount and duration of rainfall are constant,continuous rainfall has a greater impact on slope stability than intermittent rainfall.The setting of the field retaining structures has a significant role in improving slope stability.The field monitoring data show that the slope is in the initial deformation stage and has good stability,which verifies the rationality of the numerical simulation method.The research results can provide some references for understanding the influence of rainfall on the stability of loose accumulation slopes along high-speed railways and establishing a monitoring system.

Infiltration,runoff,and slope stability behaviors of infinite slope with macropores based on an improved Green–Ampt model

Infiltration–runoff–slope instability mechanism of macropore slope under heavy rainfall is unclear.This paper studied its instability mechanism with an improved Green–Ampt(GA)model considering the dual-porosity(i.e.,matrix and macropore)and ponding condition,and proposed the infiltration equations,infiltration–runoff coupled model,and safety factor calculation method.Results show that the infiltration processes of macropore slope can be divided into three stages,and the proposed model is rational by a comparative analysis.The wetting front depth of the traditional unsaturated slope is 17.2%larger than that of the macropore slope in the early rainfall stage and 27%smaller than that of the macropore slope in the late rainfall stage.Then,macropores benefit the slope stability in the early rainfall but not in the latter.Macropore flow does not occur initially but becomes pronounced with increasing rainfall duration.The equal depth of the wetting front in the two domains is regarded as the onset criteria of macropore flow.Parameter analysis shows that macropore flow is delayed by increasing proportion of macropore domain(ω_(f)),whereas promoted by increasing ratio of saturated permeability coefficients between the two domains(μ).The increasing trend of ponding depth is sharp at first and then grows slowly.Finally,when rainfall duration is less than 3 h,ωf andμhave no significant effect on the safety factor,whereas it decreases with increasingωf and increases with increasingμunder longer duration(≥3 h).With the increase ofω_(f),the slope maximum instability time advances by 10.5 h,and with the increase ofμ,the slope maximum instability time delays by 3.1 h.

44%of steep slope cropland in Europe vulnerable to drought

Steep-slope cropland plays a vital role in food production,economic development,ecosystem diversity,and Eu-ropean cultural heritage.However,these systems are susceptible to extreme weather events.The 2022 summer drought significantly impacted European agriculture,but the specific effects on steep-slope crops remain uncer-tain.Clarifying this is essential for comprehending similar future events and for implementing effective water management strategies to ensure the sustainability of steep-slope agriculture and associated ecosystem services.This study quantitatively analyzes the spatial distribution of twelve major European steep-slope(>12%)crops and assesses agricultural drought severity during the 2022 events using open-access spatial data.The satellite-based Vegetation Health Index(VHI)is utilized to identify critical hotspots.Results show that olive grove is the most widespread crop in steep slope agriculture(34%of total area),followed by wheat(24%),maize(16%),and vineyard(11%).Almost half of the steep-slope agriculture in Europe suffered drought during summer 2022.Vineyards were hardest affected at 79%,primarily in northern Portugal,northern Spain,southern France,and central Italy.Sunflowers followed at 62%,mainly in Spain,central Italy,southern France,and northern Roma-nia.Olive groves ranked third at 59%,with the most impact in northern Portugal,southern and central Spain,and southern Italy.Maize was also significantly affected at 54%.In this paper,we therefore highlight the need to increase steep-slope agriculture resilience by improving water management and promoting sustainable land practices.

Classification and rating of disintegrated dolomite strata for slope stability analysis

Although disintegrated dolomite,widely distributed across the globe,has conventionally been a focus of research in underground engineering,the issue of slope stability issues in disintegrated dolomite strata is gaining increasing prominence.This is primarily due to their unique properties,including low strength and loose structure.Current methods for evaluating slope stability,such as basic quality(BQ)and slope stability probability classification(SSPC),do not adequately account for the poor integrity and structural fragmentation characteristic of disintegrated dolomite.To address this challenge,an analysis of the applicability of the limit equilibrium method(LEM),BQ,and SSPC methods was conducted on eight disintegrated dolomite slopes located in Baoshan,Southwest China.However,conflicting results were obtained.Therefore,this paper introduces a novel method,SMRDDS,to provide rapid and accurate assessment of disintegrated dolomite slope stability.This method incorporates parameters such as disintegrated grade,joint state,groundwater conditions,and excavation methods.The findings reveal that six slopes exhibit stability,while two are considered partially unstable.Notably,the proposed method demonstrates a closer match with the actual conditions and is more time-efficient compared with the BQ and SSPC methods.However,due to the limited research on disintegrated dolomite slopes,the results of the SMRDDS method tend to be conservative as a safety precaution.In conclusion,the SMRDDS method can quickly evaluate the current situation of disintegrated dolomite slopes in the field.This contributes significantly to disaster risk reduction for disintegrated dolomite slopes.

Slope-climbing of cropland and its effects in China

In recent decades,the spatio-temporal patterns of China’s croplands have been reshaped by disturbances from anthropogenic activities,with complex changes in the topographic characteristics of croplands.Slope-climbing of croplands(SCCL)is an important issue that threatens sustainable agricultural development.While providing land with prominent location advantages,SCCL weakens the water and fertilizer retention capacity for cropland,intensifies various geological disasters,and adversely affects the ecological environment and food yield of these croplands.It is crucial to determine the spatio-temporal variation features and effects of SCCL in China to formulate more accurate cropland protection policies and to maintain food security;however,the current lack of relevant studies is detrimental for capturing trends in cropland resources and sustainable cropland use.In this study,we constructed a multi-scale slope spectrum for cropland and total terrain to explore the spatial differences and trends of SCCL from a three-dimensional view.We evaluated the natural and socioeconomic effects of SCCL in China from multiple perspectives.Results indicate that the proportion of cropland with slopes below 2°,5°,and 6°in China decreased by 0.43%,0.47%,and 0.50%from 1980 to 2020,respectively.SCCL became apparent during 1980-1990 and 2010-2020,especially over the recent decade.The cropland climbing index(CCI)and upper limited slope change(ULSC)to measure the spatio-temporal pattern of SCCL were 0.99%and 1.17°,respectively,during 2010-2020.At the agricultural regional scale,the SCCL was also concentrated in 1980-1990 and 2010-2020,and it is more pronounced in the southern areas.The proportion of provinces and prefecture-level cities with high-intensity SCCL during 1980-2020 were 87.10%and 49.73%,respectively.SCCL was comparatively more pronounced and broader from 2010 to 2020.During this period,17.84%of prefecture-level cities had no SCCL,and the average CCI for all prefecture-level cities peaked at 1.62%.In this study,we also evaluated the pros and cons of SCCL and provided targeted suggestions for decision makers and farmers to refine cropland protection policy systems and further develop the sustainable use of croplands.