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.

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.

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.

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.

Evolution model and failure mechanisms of rainfall-induced cracked red clay slopes:insights from Xinshao County,China

Red clay landslides are widely distributed worldwide,resulting in severe loss of life and property.Although rainfall-induced red clay slopes have received extensive attention,the role of cracks in the evolutionary process of red clay slopes and their connection to failure mechanisms is still poorly understood.A comprehensive approach integrating field investigation,laboratory tests,and numerical simulations was conducted to study the 168 red clay landslides in Xinshao County,China.The results show that red clay is prone to forming cracks at high moisture content due to its low swelling and high shrinkage properties.The failure mode of red clay slopes can be summarized in three stages:crack generation,slope excavation,and slope failure.Furthermore,the retrospective analysis and numerical simulations of the typical landslide in Guanchong indicated that intense rainfall primarily impacts the shallow layer of soil within approximately 0.5 m on the intact slope.However,cracks change the pattern of rainfall infiltration in the slope.Rainwater infiltrates rapidly through the preferential channels induced by the cracks rather than uniformly and slowly from the slope surface.This results in a significant increase in both the depth of infiltration and the saturated zone area of the cracked slope,reaching 3.8 m and 36.2 m^(2),respectively.Consequently,the factor of safety of the slope decreases by 13.4%compared to the intact slope,ultimately triggering landslides.This study can provide valuable insights into understanding the failure mechanisms of red clay slopes in China and other regions with similar geological settings.

Stability assessment of tree ring growth of Pinus armandii Franch in response to climate change based on slope directions at the Lubanling in the Funiu Mountains,China

Global warming will affect growth strategies and how trees will adapt.To compare the response of tree radial growth to climate warming in different slope directions,samples of Pinus armandii Franch were collected and treering chronologies developed on northern and western slopes from the Lubanling in the Funiu Mountains.Correlation analyses showed that two chronologies were mainly limited by temperatures in the previous June–August and the combination of temperatures and moisture in the current May–July.The difference of the climate response to slopes was small but not negligible.Radial growth of the LBL01 site on the northern slope was affected by the combined maximum and minimum temperatures,while that of the LBL02 site was affected by maximum temperatures.With regards to moisture,radial growth of the trees on the north slope was influenced by the relative humidity in the current May–July,while on the western slope,it was affected by the relative humidity in the previous June–August,the current May–July and the precipitation in the current May–July.With the change in climate,the effects of the main limiting factors on growth on different slopes were visible to a certain extent,but the differences in response of trees on different slopes gradually decreased,which might be caused by factors such as different slope directions and the change in diurnal temperature range.These results may provide information for forest protection and ecological construction in this region,and a scientific reference for future climate reconstruction.

Bio-cementation for tidal erosion resistance improvement of foreshore slopes based on microbially induced magnesium and calcium precipitation

In most coastal and estuarine areas,tides easily cause surface erosion and even slope failure,resulting in severe land losses,deterioration of coastal infrastructure,and increased floods.The bio-cementation technique has been previously demonstrated to effectively improve the erosion resistance of slopes.Seawater contains magnesium ions(Mg^(2+))with a higher concentration than calcium ions(Ca^(2+));therefore,Mg^(2+)and Ca^(2+)were used together for bio-cementation in this study at various Mg^(2+)/Ca^(2+)ratios as the microbially induced magnesium and calcium precipitation(MIMCP)treatment.Slope angles,surface strengths,precipitation contents,major phases,and microscopic characteristics of precipitation were used to evaluate the treatment effects.Results showed that the MIMCP treatment markedly enhanced the erosion resistance of slopes.Decreased Mg^(2+)/Ca^(2+)ratios resulted in a smaller change in angles and fewer soil losses,especially the Mg^(2+)concentration below 0.2 M.The decreased Mg^(2+)/Ca^(2+)ratio achieved increased precipitation contents,which contributed to better erosion resistance and higher surface strengths.Additionally,the production of aragonite would benefit from elevated Mg^(2+)concentrations and a higher Ca^(2+)concentration led to more nesquehonite in magnesium precipitation crystals.The slopes with an initial angle of 53°had worse erosion resistance than the slopes with an initial angle of 35°,but the Mg^(2+)/Ca^(2+)ratios of 0.2:0.8,0.1:0.9,and 0:1.0 were effective for both slope stabilization and erosion mitigation to a great extent.The results are of great significance for the application of MIMCP to improve erosion resistance of foreshore slopes and the MIMCP technique has promising application potential in marine engineering.

An approach for determination of lateral limit angle in kinematic planar sliding analysis for rock slopes

Planar sliding is one of the frequently observed types of failure in rock slopes.Kinematic analysis is a classic and widely used method to examine the potential failure modes in rock masses.The accuracy of planar sliding kinematic analysis is significantly influenced by the value assigned to the lateral limit angleγlim.However,the assignment ofγlim is currently used generally based on an empirical criterion.This study aims to propose an approach for determining the value ofγlim in deterministic and probabilistic kinematic planar sliding analysis.A new perspective is presented to reveal thatγlim essentially influences the probability of forming a potential planar sliding block.The procedure to calculate this probability is introduced using the block theory method.It is found that the probability is correlated with the number of discontinuity sets presented in rock masses.Thus,different values ofγlim for rock masses with different sets of discontinuities are recommended in both probabilistic and deterministic planar sliding kinematic analyses;whereas a fixed value ofγlim is commonly assigned to different types of rock masses in traditional method.Finally,an engineering case was used to compare the proposed and traditional kinematic analysis methods.The error rates of the traditional method vary from 45%to 119%,while that of the proposed method ranges between 1%and 17%.Therefore,it is likely that the proposed method is superior to the traditional one.

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.

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.

Stability assessment of tree ring growth of Pinus armandi Franch in response to climate change based on slope directions at the Lubanling in the Funiu Mountains,China

Global warming will affect growth strategies and how trees will adapt.To compare the response of tree radial growth to climate warming in different slope directions,samples of Pinus armandii Franch were collected and tree-ring chronologies developed on northern and western slopes from the Lubanling in the Funiu Mountains.Correlation analyses showed that two chronologies were mainly limited by temperatures in the previous June-August and the com-bination of temperatures and moisture in the current May-July.The difference of the climate response to slopes was small but not negligible.Radial growth of the LBLO1 site on the northern slope was affected by the combined maximum and minimum temperatures,while that of the LBLO2 site was affected by maximum temperatures.With regards to moisture,radial growth of the trees on the north slope was influenced by the relative humidity in the current May-July,while on the western slope,it was affected by the relative humidity in the previous June-August,the current May-July and the precipitation in the current May-July.With the change in climate,the effects of the main limiting factors on growth on different slopes were visible to a certain extent,but the differences in response of trees on different slopes gradually decreased,which might be caused by factors such as different slope directions and the change in diurnal temperature range.These results may provide information for forest protection and ecological construction in this region,and a scientific reference for future climate reconstruction.

A semi-infinite beam theoretical model on predicting rock slope subsidence induced by underground mining

When the mining goaf is close to the cliff,rock slope subsidence induced by underground mining is significantly affected by its boundary conditions.In this study,an analytical method is proposed by considering the key strata as a semi-infinite Euler-Bernoulli beam rested on a Winkler foundation with a local subsidence area.The analytical solutions of deflection are derived by analyzing the boundary and continuity conditions of the cliff.Then,the analytical solutions are verified by the results from experimental tests,FEM and InSAR,respectively.After that,the influence of changing parameters on deflections is studied with sensitivity analysis.The results show that the distance between goaf and cliff significantly affects the deflection of semi-infinite beam.The response of semi-infinite beam is obviously determined by the length of goaf and the bending stiffness of beam.The comparisons between semi-infinite beam and infinite beam illustrate the ascendancy of the improved model in such problems.

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.

Seismic stability of expansive soil slopes reinforced by anchor cables using a modified horizontal slice method

Earthquake-induced slope failures are common occurrences in engineering practice and pre-stressed anchor cables are an effective technique in maintaining slope stability,especially in areas that are prone to earthquakes.Furthermore,the soil at typical engineering sites also exhibit unsaturated features.Explicit considerations of these factors in slope stability estimations are crucial in producing accurate results.In this study,the seismic responses of expansive soil slopes stabilized by anchor cables is studied in the realm of kinematic limit analysis.A modified horizontal slice method is proposed to semi-analytically formulate the energy balance equation.An illustrative slope is studied to demonstrate the influences of suction,seismic excitations and anchor cables on the slope stability.The results indicate that the stabilizing effect of soil suction relates strongly to the seismic excitation and presents a sine shape as the seismic wave propagates.In higher and steeper slopes,the stabilizing effect of suction is more evident.The critical slip surface tends to be much more shallow as the seismic wave approaches the peak and vice versa.

Efficient slope reliability analysis under soil spatial variability using maximum entropy distribution with fractional moments

Spatial variability of soil properties imposes a challenge for practical analysis and design in geotechnical engineering.The latter is particularly true for slope stability assessment,where the effects of uncertainty are synthesized in the so-called probability of failure.This probability quantifies the reliability of a slope and its numerical calculation is usually quite involved from a numerical viewpoint.In view of this issue,this paper proposes an approach for failure probability assessment based on Latinized partially stratified sampling and maximum entropy distribution with fractional moments.The spatial variability of geotechnical properties is represented by means of random fields and the Karhunen-Loève expansion.Then,failure probabilities are estimated employing maximum entropy distribution with fractional moments.The application of the proposed approach is examined with two examples:a case study of an undrained slope and a case study of a slope with cross-correlated random fields of strength parameters under a drained slope.The results show that the proposed approach has excellent accuracy and high efficiency,and it can be applied straightforwardly to similar geotechnical engineering problems.

Probabilistic Analysis of Slope Using Finite Element Approach and Limit Equilibrium Approach around Amalpata Landslide of West Central, Nepal

The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have different slope inclinations. The lower bench, located above the basement, consistently fails and sets others up for failure. The fluctuating water level of the slope, which travels down the slope masses, exacerbates the slide problem. The majority of these rocks are Amalpata landslide area experiences several structural disruptions. The area’s stability must be evaluated in order to prevent and control more harm from occurring to the nearby agricultural land and people living along the slope. The slopes’ failures increase the damages of house existing in nearby area and the erosion of the slope. Two modeling techniques the finite element approach and the limit equilibrium method were used to simulate the slope. The findings show that, in every case, the terrace above the basement is where the majority of the stress is concentrated, with a safety factor of near unity. Using probabilistic slope stability analysis, the failure probability was predicted to be between 98.90% and 100%.

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.

Hydrodynamic characteristics and particle tracking of 90° lateral intakes at an inclined river slope

Lateral intakes are common in rivers.The pump effciency and sediment deposition are determined by the local hydrodynamic characteristics and mainstream division width.The hydraulic characteristics of lateral withdrawal from inclined river slopes at different intake elevations should be investigated.Meanwhile,the division width exhibits significant vertical non-uniformity at an inclined river slope,which should be clarified.Hence,a three-dimensional(3-D)hydrodynamic and particle-tracking model was developed with the Open Source Field Operation and Manipulation(Open FOAM),and the model was validated with physical model tests for 90°lateral withdrawal from an inclined side bank.The flow fields,withdrawal sources,and division widths were investigated with different intake bottom elevations,withdrawal discharges,and main channel velocities.This study showed that under inclined side bank conditions,water entered the intake at an oblique angle,causing significant 3-D spiral flows in the intake rather than two-dimensional closed recirculation.A lower withdrawal discharge,a lower bottom elevation of the intake,or a higher main channel velocity could further strengthen this phenomenon.The average division width and turbulent kinetic energy were smaller under inclined side bank conditions than under vertical bank conditions.With a low intake bottom elevation,a low withdrawal discharge,or a high main channel velocity,the sources of lateral withdrawal were in similar ranges near the local inclined bank in the vertical direction.Under inclined slope conditions,sediment deposition near the intake entrance could be reduced,compared to that under vertical slope conditions.The results provide hydrodynamic and sediment references for engineering designs for natural rivers with inclined terrains.

Mechanism of Diabatic Heating on Precipitation and the Track of a Tibetan Plateau Vortex over the Eastern Slope of the Tibetan Plateau

Existing studies contend that latent heating(LH)will replace sensible heating(SH)to become the dominant factor affecting the development of the Tibetan Plateau vortex(TPV)after it moves off the Tibetan Plateau(TP).However,in the process of the TPV moving off the TP requires that the airmass traverse the eastern slope of the Tibetan Plateau(ESTP)where the topography and diabatic heating(DH)conditions rapidly change.How LH gradually replaces SH to become the dominant factor in the development of the TPV over the ESTP is still not very clear.In this paper,an analysis of a typical case of a TPV with a long life history over the ESTP is performed by using multi-sourced meteorological data and model simulations.The results show that SH from the TP surface can change the TPV-associated precipitation distribution by temperature advection after the TPV moves off the TP.The LH can then directly promote the development of the TPV and has a certain guiding effect on the track of the TPV.The SH can control the active area of LH by changing the falling area of the TPV-associated precipitation,so it still plays a key role in the development and tracking of the TPV even though it has moved out of the main body of the TP.

Uphill or downhill?Cropland use change and its drivers from the perspective of slope spectrum

The continuous decrease of low-slope cropland resources caused by construction land crowding poses huge threat to regional sustainable development and food security.Slope spectrum analysis of topographic and geomorphic features is considered as a digital terrain analysis method which reflects the macro-topographic features by using micro-topographic factors.However,pieces of studies have extended the concept of slope spectrum in the field of geoscience to construction land to explore its expansion law,while research on the slope trend of cropland from that perspective remains rare.To address the gap,in virtue of spatial analysis and geographically weighted regression(GWR)model,the cropland use change in the Yangtze River Basin(YRB)from 2000 to 2020 was analyzed and the driving factors were explored from the perspective of slope spectrum.Results showed that the slope spectrum curves of cropland area-frequency in the YRB showed a first upward then a downward trend.The change curve of the slope spectrum of cropland in each province(municipality)exhibited various distribution patterns.Quantitative analysis of morphological parameters of cropland slope spectrum revealed that the further down the YRB,the stronger the flattening characteristics,the more obvious the concentration.The province experienced the greatest downhill cropland climbing(CLC)was Shannxi,while province experienced the highest uphill CLC was Zhejiang.The most common cropland use change type in the YRB was horizontal expansion type.The factors affecting average cropland climbing index(ACCI)were quite stable in different periods,while population density(POP)changed from negative to positive during the study period.This research is of practical significance for the rational utilization of cropland at the watershed scale.