Predicting and validating the load-settlement behavior of large-scale geosynthetic-reinforced soil abutments using hybrid intelligent modeling

Settlement prediction of geosynthetic-reinforced soil(GRS)abutments under service loading conditions is an arduous and challenging task for practicing geotechnical/civil engineers.Hence,in this paper,a novel hybrid artificial intelligence(AI)-based model was developed by the combination of artificial neural network(ANN)and Harris hawks’optimisation(HHO),that is,ANN-HHO,to predict the settlement of the GRS abutments.Five other robust intelligent models such as support vector regression(SVR),Gaussian process regression(GPR),relevance vector machine(RVM),sequential minimal optimisation regression(SMOR),and least-median square regression(LMSR)were constructed and compared to the ANN-HHO model.The predictive strength,relalibility and robustness of the model were evaluated based on rigorous statistical testing,ranking criteria,multi-criteria approach,uncertainity analysis and sensitivity analysis(SA).Moreover,the predictive veracity of the model was also substantiated against several large-scale independent experimental studies on GRS abutments reported in the scientific literature.The acquired findings demonstrated that the ANN-HHO model predicted the settlement of GRS abutments with reasonable accuracy and yielded superior performance in comparison to counterpart models.Therefore,it becomes one of predictive tools employed by geotechnical/civil engineers in preliminary decision-making when investigating the in-service performance of GRS abutments.Finally,the model has been converted into a simple mathematical formulation for easy hand calculations,and it is proved cost-effective and less time-consuming in comparison to experimental tests and numerical simulations.

Influence of groundwater level changes on the seismic response of geosynthetic-reinforced soil retaining walls

Geosynthetic-reinforced soil retaining walls(GSRWs)have been widely used in civil engineering projects.However,as the climate changes,extreme weather conditions and natural hazards are likely to become more frequent or intense,posing a huge threat to the stability of GSRWs.In this paper,the effect of groundwater level fluctuations on the seismic response of GSRWs is investigated.First,a dynamic numerical model was established and validated through centrifugal shaking-table test results.Using the established numerical model,the seismic response of GSRWs under four different groundwater level conditions was then investigated,i.e.,an earthquake occurring at a low groundwater level(Case LW),an earthquake occurring when the groundwater level rises(Case RW),an earthquake occurring at a high groundwater level(Case HW),and an earthquake occurring when the groundwater level drops(Case DW).The results show that the GSRW in Case DW has the worst seismic stability because of the drag forces generated by the water flowing to the outside of the GSRW.For Case RW,deformation of the GSRW under earthquake forces was prevented by the drag forces generated by the water flowing to the inside of the GSRW and the water pressure acting on the outside of the facing,giving the GSRW the best seismic stability in this case.Compared with Case LW,the seismic stability of a GSRW in Case HW is worse,because the high groundwater level will generate excess pore-water pressure during an earthquake.On this basis,we provide engineering design suggestions to be considered by practitioners.

Energy analysis of geosynthetic-reinforced slope in unsaturated soils subjected to steady flow

Soils are actually unsaturated in nature. In the present study, a stability analysis of a geosynthetic-reinforced slope in unsaturated soils subjected to various steady flow conditions is conducted based on limit analysis. Work rate by apparent cohesion due to matric suction is calculated based on the effective stress-based equation. Analytical expression of the required cohesion/stability number of slope is derived from the energy balance equation. An optimization code is programmed to capture the optimized solution of the stability number. Comparison is made to verify the present work and a parametric analysis is conducted to investigate the effects of soil type, infilitration rate, reinforcement strength and soil suction on slope stability afterwards. A set of numerical solutions is presented at the end of the paper for preliminary design purposes.

Numerical analysis of geosynthetic-reinforced embankment performance under moving loads

The performance of geosynthetic-reinforced embankments under traffic moving loads is always a hotspot in the geotechnical engineering field.A three-dimensional(3D)model of a geosynthetic-reinforced embankment without drainage consolidation was established using the finite element software ABAQUS.In this model,the traffic loads were simulated by two moving loads of rectangular pattern,and their amplitude,range,and moving speed were realized by a Fortran subroutine.The embankment fill was simulated by an equivalent linear viscoelastic model,which can reflect its viscoelasticity.The geogrid was simulated by the truss element,and the geocell was simulated by the membrane element.Infinite elements were utilized to weaken the boundary effect caused by the model geometry at the boundaries.Validation of the established numerical model was conducted by comparing the predicted deformations in the cross-section of the geosynthetic-reinforced embankment with those from the existing literature.On this basis,the dynamic stress and strain distribution in the pavement structure layer of the geosynthetic-reinforced embankment under a moving load was also analyzed.Finally,a parametric study was conducted to examine the influences of the different types of reinforcement,overload,and the moving load velocity on the geosynthetic-reinforced embankment.

Seismic effects on reinforcement load and lateral deformation of geosynthetic-reinforced soil walls

Current design methods for the internal stability of geosynthetic-reinforced soil(GRS)walls postulate seismic forces as inertial forces,leading to pseudo-static analyses based on active earth pressure theory,which yields unconservative reinforcement loads required for seismic stability.Most seismic analyses are limited to the determination of maximum reinforcement strength.This study aimed to calculate the distribution of the reinforcement load and connection strength required for each layer of the seismic GRS wall.Using the top-down procedure involves all of the possible failure surfaces for the seismic analyses of the GRS wall and then obtains the reinforcement load distribution for the limit state.The distributions are used to determine the required connection strength and to approximately assess the facing lateral deformation.For sufficient pullout resistance to be provided by each reinforcement,the maximum required tensile resistance is identical to the results based on the Mononobe-Okabe method.However,short reinforcement results in greater tensile resistances in the mid and lower layers as evinced by compound failure frequently occurring in GRS walls during an earthquake.Parametric studies involving backfill friction angle,reinforcement length,vertical seismic acceleration,and secondary reinforcement are conducted to investigate seismic impacts on the stability and lateral deformation of GRS walls.

Compaction-induced stress in geosynthetic-reinforced granular base course--A discrete element model

A discrete element method(DEM) model was used to simulate the development of compaction-induced stress in a granular base course, with and without geogrid reinforcement. The granular base course was modeled as a mixture of uniformly sized triangular particles. The geogrid was modeled as a series of equally spaced balls that interact with each other through long-range interaction contacts. The longrange interaction contact was also used to simulate a deformable subgrade. The compactor was modeled as a solid cylinder rolling at a constant speed. The DEM model shows that the geogridreinforced granular base course gains additional compaction-induced stress due to the residual tensile stress developed in the geogrid. The residual tensile stress in the geogrid increases with the number of compaction passes. Parametric analyses were also conducted to assess the effects of geogrid stiffness and subgrade modulus on the compaction-induced stress.

The competition between Bidens pilosa and Setaria viridis alters soil microbial composition and soil ecological function

Bidens pilosa is recognized as one of the major invasive plants in China.Its invasion has been associated with significant losses in agriculture,forestry,husbandry,and biodiversity.Soil ecosystems play an important role in alien plant invasion.Microorganisms within the soil act as intermediaries between plants and soil ecological functions,playing a role in regulating soil enzyme activities and nutrient dynamics.Understanding the interactions between invasive plants,soil microorganisms,and soil ecological processes is vital for managing and mitigating the impacts of invasive species on the environment.In this study,we conducted a systematic analysis focusing on B.pilosa and Setaria viridis,a common native companion plant in the invaded area.To simulate the invasion process of B.pilosa,we constructed homogeneous plots consisting of B.pilosa and S.viridis grown separately as monocultures,as well as in mixtures.The rhizosphere and bulk soils were collected from the alien plant B.pilosa and the native plant S.viridis.In order to focus on the soil ecological functional mechanisms that contribute to the successful invasion of B.pilosa,we analyzed the effects of B.pilosa on the composition of soil microbial communities and soil ecological functions.The results showed that the biomass of B.pilosa increased by 27.51% and that of S.viridis was significantly reduced by 66.56%.The organic matter contents in the bulk and rhizosphere soils of B.pilosa were approximately 1.30 times those in the native plant soils.The TN and NO_(3)^(-)contents in the rhizosphere soil of B.pilosa were 1.30 to 2.71 times those in the native plant soils.The activities of acid phosphatase,alkaline phosphatase,and urease in the rhizosphere soil of B.pilosa were 1.98-2.25 times higher than in the native plant soils.Using high-throughput sequencing of the16S rRNA gene,we found that B.pilosa altered the composition of the soil microbial community.Specifically,many genera in Actinobacteria and Proteobacteria were enriched in B.pilosa soils.Further correlation analyses verified that these genera had significantly positive relationships with soil nutrients and enzyme activities.Plant biomass,soil p H,and the contents of organic matter,TN,NO_(3)^(-),TP,AP,TK,and AK were the main factors affecting soil microbial communities.This study showed that the invasion of B.pilosa led to significant alterations in the composition of the soil microbial communities.These changes were closely linked to modifications in plant traits as well as soil physical and chemical properties.Some microbial species related to C,N and P cycling were enriched in the soil invaded by B.pilosa.These findings provide additional support for the hypothesis of soil-microbe feedback in the successful invasion of alien plants.They also offer insights into the ecological mechanism by which soil microbes contribute to the successful invasion of B.pilosa.Overall,our research contributes to a better understanding of the complex interactions between invasive plants,soil microbial communities,and ecosystem dynamics.

Biochar Decreases Soil Cadmium(Cd) Availability and Regulates Expression Levels of Cd Uptake/Transport-Related Genes to Reduce Cd Translocation in Rice

Biochar has been used as a soil amendment for heavy metal-contaminated soils, and it has the potential to mitigate Cd accumulation in plants. In this study, we used rice straw biochar(RSB) and kitchen waste biochar(KWB) to clarify the effect of biochar on Cd-contaminated neutral soil, the physiological responses to biochar application, and the gene regulatory networks in a rice genotype.

Subsoil tillage enhances wheat productivity,soil organic carbon and available nutrient status in dryland fields

Tillage practices during the fallow period benefit water storage and yield in dryland wheat crops.However,there is currently no clarity on the responses of soil organic carbon(SOC),total nitrogen(TN),and available nutrients to tillage practices within the growing season.This study evaluated the effects of three tillage practices(NT,no tillage;SS,subsoil tillage;DT,deep tillage)over five years on soil physicochemical properties.Soil samples at harvest stage from the fifth year were analyzed to determine the soil aggregate and aggregate-associated C and N fractions.The results indicated that SS and DT improved grain yield,straw biomass and straw carbon return of wheat compared with NT.In contrast to DT and NT,SS favored SOC and TN concentrations and stocks by increasing the soil organic carbon sequestration rate(SOCSR)and soil nitrogen sequestration rate(TNSR)in the 0-40 cm layer.Higher SOC levels under SS and NT were associated with greater aggregate-associated C fractions,while TN was positively associated with soluble organic nitrogen(SON).Compared with DT,the NT and SS treatments improved soil available nutrients in the 0-20 cm layer.These findings suggest that SS is an excellent practice for increasing soil carbon,nitrogen and nutrient availability in dryland wheat fields in North China.

A dynamic soil freezing characteristic curve model for frozen soil

The soil freezing characteristic curve(SFCC)plays a fundamental role in comprehending thermohydraulic behavior and numerical simulation of frozen soil.This study proposes a dynamic model to uniformly express SFCCs amidst varying total water contents throughout the freezing-thawing process.Firstly,a general model is proposed,wherein the unfrozen water content at arbitrary temperature is determined as the lesser of the current total water content and the reference value derived from saturated SFCC.The dynamic performance of this model is verified through test data.Subsequently,in accordance with electric double layer(EDL)theory,the theoretical residual and minimum temperatures in SFCC are calculated to be-14.5℃to-20℃for clay particles and-260℃,respectively.To ensure that the SFCC curve ends at minimum temperature,a correction function is introduced into the general model.Furthermore,a simplified dynamic model is proposed and investigated,necessitating only three parameters inherited from the general model.Additionally,both general and simplified models are evaluated based on a test database and proven to fit the test data exactly across the entire temperature range.Typical recommended parameter values for various types of soils are summarized.Overall,this study provides not only a theoretical basis for most empirical equations but also proposes a new and more general equation to describe the SFCC.

Straw return increases crop production by improving soil organic carbon sequestration and soil aggregation in a long-term wheat-cotton cropping system

Straw return is a promising strategy for managing soil organic carbon(SOC)and improving yield stability.However,the optimal straw return strategy for sustainable crop production in the wheat(Triticum aestivum L.)-cotton(Gossypium hirsutum L.)cropping system remains uncertain.The objective of this study was to quantify the long-term(10 years)impact of carbon(C)input on SOC sequestration,soil aggregation and crop yields in a wheat-cotton cropping system in the Yangtze River Valley,China.Five treatments were arranged with a single-factor randomized design as follows:no straw return(Control),return of wheat straw only(Wt),return of cotton straw only(Ct),return of 50%wheat and 50%cotton straw(Wh-Ch)and return of 100%wheat and 100%cotton straw(Wt-Ct).In comparison to the Control,the SOC content increased by 8.4 to 20.2%under straw return.A significant linear positive correlation between SOC sequestration and C input(1.42-7.19 Mg ha^(−1)yr^(−1))(P<0.05)was detected.The percentages of aggregates of sizes>2 and 1-2 mm at the 0-20 cm soil depth were also significantly elevated under straw return,with the greatest increase of the aggregate stability in the Wt-Ct treatment(28.1%).The average wheat yields increased by 12.4-36.0%and cotton yields increased by 29.4-73.7%,and significantly linear positive correlations were also detected between C input and the yields of wheat and cotton.The average sustainable yield index(SYI)reached a maximum value of 0.69 when the C input was 7.08 Mg ha^(−1)yr^(−1),which was close to the maximum value(SYI of 0.69,C input of 7.19 Mg ha^(−1)yr^(-1))in the Wt-Ct treatment.Overall,the return of both wheat and cotton straw was the best strategy for improving SOC sequestration,soil aggregation,yields and their sustainability in the wheat-cotton rotation system.

Assessment of soil erosion in the Irga watershed on the eastern edge of the Chota Nagpur Plateau,India

Human activities to improve the quality of life have accelerated the natural rate of soil erosion.In turn,these natural disasters have taken a great impact on humans.Human activities,particularly the conversion of vegetated land into agricultural land and built-up area,stand out as primary contributors to soil erosion.The present study investigated the risk of soil erosion in the Irga watershed located on the eastern fringe of the Chota Nagpur Plateau in Jharkhand,India,which is dominated by sandy loam and sandy clay loam soil with low soil organic carbon(SOC)content.The study used the Revised Universal Soil Loss Equation(RUSLE)and Geographical Information System(GIS)technique to determine the rate of soil erosion.The five parameters(rainfall-runoff erosivity(R)factor,soil erodibility(K)factor,slope length and steepness(LS)factor,cover-management(C)factor,and support practice(P)factor)of the RUSLE were applied to present a more accurate distribution characteristic of soil erosion in the Irga watershed.The result shows that the R factor is positively correlated with rainfall and follows the same distribution pattern as the rainfall.The K factor values in the northern part of the study area are relatively low,while they are relatively high in the southern part.The mean value of the LS factor is 2.74,which is low due to the flat terrain of the Irga watershed.There is a negative linear correlation between Normalized Difference Vegetation Index(NDVI)and the C factor,and the high values of the C factor are observed in places with low NDVI.The mean value of the P factor is 0.210,with a range from 0.000 to 1.000.After calculating all parameters,we obtained the average soil erosion rate of 1.43 t/(hm^(2)•a),with the highest rate reaching as high as 32.71 t/(hm^(2)•a).Therefore,the study area faces a low risk of soil erosion.However,preventative measures are essential to avoid future damage to productive and constructive activities caused by soil erosion.This study also identifies the spatial distribution of soil erosion rate,which will help policy-makers to implement targeted soil erosion control measures.

Manure substitution improves maize yield by promoting soil fertility and mediating the microbial community in lime concretion black soil

Synthetic nitrogen(N)fertilizer has made a great contribution to the improvement of soil fertility and productivity,but excessive application of synthetic N fertilizer may cause agroecosystem risks,such as soil acidification,groundwater contamination and biodiversity reduction.Meanwhile,organic substitution has received increasing attention for its ecologically and environmentally friendly and productivity benefits.However,the linkages between manure substitution,crop yield and the underlying microbial mechanisms remain uncertain.To bridge this gap,a three-year field experiment was conducted with five fertilization regimes:i)Control,non-fertilization;CF,conventional synthetic fertilizer application;CF_(1/2)M_(1/2),1/2 N input via synthetic fertilizer and 1/2 N input via manure;CF_(1/4)M_(3/4),1/4 N input synthetic fertilizer and 3/4 N input via manure;M,manure application.All fertilization treatments were designed to have equal N input.Our results showed that all manure substituted treatments achieved high soil fertility indexes(SFI)and productivities by increasing the soil organic carbon(SOC),total N(TN)and available phosphorus(AP)concentrations,and by altering the bacterial community diversity and composition compared with CF.SOC,AP,and the soil C:N ratio were mainly responsible for microbial community variations.The co-occurrence network revealed that SOC and AP had strong positive associations with Rhodospirillales and Burkholderiales,while TN and C:N ratio had positive and negative associations with Micromonosporaceae,respectively.These specific taxa are implicated in soil macroelement turnover.Random Forest analysis predicted that both biotic(bacterial composition and Micromonosporaceae)and abiotic(AP,SOC,SFI,and TN)factors had significant effects on crop yield.The present work strengthens our understanding of the effects of manure substitution on crop yield and provides theoretical support for optimizing fertilization strategies.

Enhancing Deep Learning Soil Moisture Forecasting Models by Integrating Physics-based Models

Accurate soil moisture(SM)prediction is critical for understanding hydrological processes.Physics-based(PB)models exhibit large uncertainties in SM predictions arising from uncertain parameterizations and insufficient representation of land-surface processes.In addition to PB models,deep learning(DL)models have been widely used in SM predictions recently.However,few pure DL models have notably high success rates due to lacking physical information.Thus,we developed hybrid models to effectively integrate the outputs of PB models into DL models to improve SM predictions.To this end,we first developed a hybrid model based on the attention mechanism to take advantage of PB models at each forecast time scale(attention model).We further built an ensemble model that combined the advantages of different hybrid schemes(ensemble model).We utilized SM forecasts from the Global Forecast System to enhance the convolutional long short-term memory(ConvLSTM)model for 1–16 days of SM predictions.The performances of the proposed hybrid models were investigated and compared with two existing hybrid models.The results showed that the attention model could leverage benefits of PB models and achieved the best predictability of drought events among the different hybrid models.Moreover,the ensemble model performed best among all hybrid models at all forecast time scales and different soil conditions.It is highlighted that the ensemble model outperformed the pure DL model over 79.5%of in situ stations for 16-day predictions.These findings suggest that our proposed hybrid models can adequately exploit the benefits of PB model outputs to aid DL models in making SM predictions.

Time-lagged Effects of the Spring Atmospheric Heat Source over the Tibetan Plateau on Summer Precipitation in Northeast China during 1961–2020:Role of Soil Moisture

The spring atmospheric heat source(AHS)over the Tibetan Plateau(TP)has been suggested to affect the Asian summer monsoon and summer precipitation over South China.However,its influence on the summer precipitation in Northeast China(NEC)remains unknown.The connection between spring TP AHS and subsequent summer precipitation over NEC from 1961 to 2020 is analyzed in this study.Results illustrate that stronger spring TP AHS can enhance subsequent summer NEC precipitation,and higher soil moisture in the Yellow River Valley-North China region(YRVNC)acts as a bridge.During spring,the strong TP AHS could strengthen the transportation of water vapor to East China and lead to excessive rainfall in the YRVNC.Thus,soil moisture increases,which regulates local thermal conditions by decreasing local surface skin temperature and sensible heat.Owing to the memory of soil moisture,the lower spring sensible heat over the YRVNC can last until mid-summer,decrease the land–sea thermal contrast,and weaken the southerly winds over the East Asia–western Pacific region and convective activities over the South China Sea and tropical western Pacific.This modulates the East Asia–Pacific teleconnection pattern,which leads to a cyclonic anomaly and excessive summer precipitation over NEC.

Semi-analytical solution for drained expansion analysis of a hollow cylinder of critical state soils

The expansion of a thick-walled hollow cylinder in soil is of non-self-similar nature that the stress/deformation paths are not the same for different soil material points.As a result,this problem cannot be solved by the common self-similar-based similarity techniques.This paper proposes a novel,exact solution for rigorous drained expansion analysis of a hollow cylinder of critical state soils.Considering stress-dependent elastic moduli of soils,new analytical stress and displacement solutions for the nonself-similar problem are developed taking the small strain assumption in the elastic zone.In the plastic zone,the cavity expansion response is formulated into a set of first-order partial differential equations(PDEs)with the combination use of Eulerian and Lagrangian descriptions,and a novel solution algorithm is developed to efficiently solve this complex boundary value problem.The solution is presented in a general form and thus can be useful for a wide range of soils.With the new solution,the non-self-similar nature induced by the finite outer boundary is clearly demonstrated and highlighted,which is found to be greatly different to the behaviour of cavity expansion in infinite soil mass.The present solution may serve as a benchmark for verifying the performance of advanced numerical techniques with critical state soil models and be used to capture the finite boundary effect for pressuremeter tests in small-sized calibration chambers.

Impact of mica on geotechnical behavior of weathered granitic soil using macro and micro investigations

The micaceous weathered granitic soil(WGS)is frequently encountered in civil engineering worldwide,unfortunately little information is available regarding how mica affects the physico-mechanical behaviors of WGS.This study prepares reconstituted WGS with different mica contents by removing natural mica in theWGS,and then mixes it with commercial mica powders.The geotechnical behavior as well as the microstructures of the mixtures are characterized.The addition of mica enables the physical indices of WGS to be specific combinations of coarser gradation and high permeability but high Atterberg limits.However,high mica content in WGS was found to be associated with undesirable mechanical properties,including increased compressibility,disintegration,and swelling potential,as well as poor compactability and low effective frictional angle.Microstructural analysis indicates that the influence of mica on the responses of mixtures originates from the intrinsic nature of mica as well as the particle packing being formed withinWGS.Mica exists in the mixture as stacks of plates that form a spongy structure with high compressibility and swelling potential.Pores among the plates give the soil high water retention and high Atterberg limits.Large pores are also generated by soil particles with bridging packing,which enhances the permeability and water-soil interactions upon immersion.This study provides a microlevel understanding of how mica dominates the behavior of WGS and provides new insights into the effective stabilization and improvement of micaceous soils.

Formation of Tianwen-1 landing crater and mechanical properties of Martian soil near the landing site

After landing in the Utopia Planitia,Tianwen-1 formed the deepest landing crater on Mars,approximately 40 cm deep,exposing precious information about the mechanical properties of Martian soil.We established numerical models for the plume-surface interaction(PSI)and the crater formation based on Computational Fluid Dynamics(CFD)methods and the erosion model modified from Roberts’Theory.Comparative studies of cases were conducted with different nozzle heights and soil mechanical properties.The increase in cohesion and internal friction angle leads to a decrease in erosion rate and maximum crater depth,with the cohesion having a greater impact.The influence of the nozzle height is not clear,as it interacts with the position of the Shock Diamond to jointly control the erosion process.Furthermore,we categorized the evolution of landing craters into the dispersive and the concentrated erosion modes based on the morphological characteristics.Finally,we estimated the upper limits of the Martian soil’s mechanical properties near Tianwen-1 landing site,with the cohesion ranging from 2612 to 2042 Pa and internal friction angle from 25°to 41°.

Physical and mechanical properties and microstructures of submarine soils in the Yellow Sea

In recent years,the exploration of seabed has been intensified,but the submarine soils of silt and sand in the Yellow Sea area have not been well investigated so far.In this study,the physical and mechanical properties of silt and sand from the Yellow Sea were measured using a direct shear apparatus and their microstructures were observed using a scanning electron microscope.The test results suggest that the shear strength of silt and sand increases linearly with the increase of normal stress.Based on the direct shear test,the scanning electron microscope was used to observe the section surface of sand.It is observed that the section surface becomes rough,with many“V”‐shaped cracks.Many particles appear on the surface of the silt structure and tend to be disintegrated.The X‐ray diffraction experiment reveals that the sand and silt have different compositions.The shear strength of sand is slightly greater than that of silt under high stress,which is related to the shape of soil particles and the mineral composition.These results can be a reference for further study of other soils in the Yellow Sea;meanwhile,they can serve as soil parameters for the stability and durability analyses of offshore infrastructure construction.

Analytical model for predicting time-dependent lateral deformation of geosynthetics-reinforced soil walls with modular block facing

To date,few models are available in the literature to consider the creep behavior of geosynthetics when predicting the lateral deformation(d)of geosynthetics-reinforced soil(GRS)retaining walls.In this study,a general hyperbolic creep model was first introduced to describe the long-term deformation of geosynthetics,which is a function of elapsed time and two empirical parameters a and b.The conventional creep tests with three different tensile loads(Pr)were conducted on two uniaxial geogrids to determine their creep behavior,as well as the a-Pr and b-Pr relationships.The test results show that increasing Pr accelerates the development of creep deformation for both geogrids.Meanwhile,a and b respectively show exponential and negatively linear relationships with Pr,which were confirmed by abundant experimental data available in other studies.Based on the above creep model and relationships,an accurate and reliable analytical model was then proposed for predicting the time-dependent d of GRS walls with modular block facing,which was further validated using a relevant numerical investigation from the previous literature.Performance evaluation and comparison of the proposed model with six available prediction models were performed.Then a parametric study was carried out to evaluate the effects of wall height,vertical spacing of geogrids,unit weight and internal friction angle of backfills,and factor of safety against pullout on d at the end of construction and 5 years afterwards.The findings show that the creep effect not only promotes d but also raises the elevation of the maximum d along the wall height.Finally,the limitations and application prospects of the proposed model were discussed and analyzed.