An experimental investigation of the influence of plasticity on creep degradation rate

Intrinsic soil properties, such as the Atterberg limits, are essential factors influencing the mechanical behaviour of the fine-grained soils. In this study, a series of long-term multiple-stage loading oedometer tests were performed on alluvial organic soils to investigate the creep behaviour. The plasticity ratios ranged from 0.4 to 0.63. The smaller value of the plasticity ratio Rp indicated higher soil plasticity. The results showed that the coefficient of secondary compression Cαe of alluvial organic soils was stress- and strain-rate-dependent. The coefficient of secondary compression change index m was derived using a double-logarithmic approach for a creep degradation and was related to the plasticity and clay percentage to fines. Based on the results, it was found that high plasticity soils exhibit slow creep degradation rate during one-dimensional straining under normally consolidated state. The results show that the higher soil plasticity expressed by the plasticity index, plasticity ratio and clay percentage to fines, smaller the coefficient of secondary compression change index. Moreover, the correlations among a soil plasticity properties and creep parameters for the alluvial soils have also been proposed.


Introduction
Development of road infrastructure passing through the areas of soft soils occurrence imposes more and more demanding standards in terms of assessing possible behaviour of the ground and its strengthening treatment.
An example of such an area is the _ Zuławy Fens deposit, which lies in the northern part of Poland. The deposit that covers the Vistula River delta is the area of occurrence of delta formations in the subsoil, which are soft mineralorganic soils with a significant spatial variability. This kind of soils is particularly high compressible and susceptible to the time-dependent deformations which are a serious challenge for the geotechnical design [12,13,51].
Some natural fine-grained soils exhibit time-dependent mechanical properties. Numerous studies point to the significant contribution of secondary compression in the deformation of the soil subjected to an external load [5,9,11,20,72,83]. The consolidation process is traditionally divided into two successive phases: primary compression and secondary compression. During the primary phase, soil compression is controlled by the dissipation of excess pore water pressures. This dissipation is sometimes delayed by the viscous-plastic effects (creep) [7,22,67]. At the end of the primary consolidation (EOP), the excess is dissipated, and the initial total applied stress becomes fully effective. After that, the soil continues to deform, but at a strain rate controlled by the soil viscosity. This deformation is referred to as the secondary consolidation, and a point marking the transition between the two phases is the EOP state. Some suggest that viscous behaviour occurs during and after primary consolidation in accordance with Creep Hypothesis B, formally proposed by Ladd et al. [42]. Creep deformation occurs because of the breakdown of interparticle bonds [46,48], jumping of the molecule bonds [41,49], sliding among particles [39], double porosity [63,75] and the strain rate dependence of the coefficients of friction at the interparticle contacts [34,35]. The Creep Hypothesis B is referred to as the isotache theory, in which the strain at EOP increases with the thickness of the clay and leads to unique secondary compression [14-18, 24, 36, 44, 60, 65, 71]. In contrast, there is also evidence to suggest that similar strain levels at EOP, along any stage of consolidation, can occur independently of the thickness of the consolidated layer and the duration of primary consolidation [45,47,48]. This independence is often referred to as Creep Hypothesis A, which implies the same mobilized preconsolidation pressure in the field as the preconsolidation pressure determined from laboratory tests on thin samples.
Conventionally, the time-dependent behaviour of a soil can be described using the coefficient of secondary compression C ae defined as follows: To simplify the foundation design calculations, C ae is assumed as a constant value. However, this statement leads to unrealistic soil behaviour. As indicated by Yin [69], the logarithmic function has a limitation; when time is infinite, the deformation is also infinite and may overestimate the creep settlement. To overcome the limitation of the logarithmic function, some nonlinear creep formulations have been proposed. Yin et al. [73] formulated a nonlinear expression of C ae with volumetric strain and time that considers the influence of the bulk density q and e. Wu et al. [68] also considered consecutively decreasing C ae related to the q. Zhu et al. [82] and Yin et al. [74] demonstrated a decreasing C ae as q increased. They introduced a modified C ae * , defined in a double-logarithmic plane. On the other hand, Zhu et al. [80] discussed a nonlinear variation of C ae as a function of C ae at the liquid limit, structure indicator v = r' p /r' pi -1 (v reflects interparticle bonding) and the creep coefficient-based structure indicator -= C ae /C aei -1. Further, the nonlinear creep formulation accounting for both soil density and soil structure was proposed and the relationship between the LL, LP, as well as v and C ae , was suggested.
Numerous attempts have been made to identify the stress and strain rate dependency of natural fine-grained soils. Generally, rate dependency is determined using the relationship between C ae and vertical effective stress or compression. There are two opinions about the dependence of C ae on stress. Newland and Alley [54] and Barden [4] noted that C ae is almost independent of the effective stress except near the point of yield/preconsolidation pressure, where values of C ae are slightly higher. Zhang and Wang [78] also postulated that, for some remoulded soils, C ae is independent on vertical effective stress. In contrast, there are many more opinions that say C ae is stress-dependent. This point of view is connected with the occurrence of the maximum value of C ae , which is closely related to the yield/preconsolidation stress. Several authors have reported that, as the effective stress increases, C ae increases to the maximum and then decreases as the stress continues to increase [23,43,71,74,76,79]. Considering the findings of the authors mentioned above, C ae can be either very low for the vertical stresses well below the yield/preconsolidation pressure. However, C ae increases rapidly as the yield/preconsolidation pressure is approached; it reaches its maximum at the yield/preconsolidation pressure and decreases under higher stresses. There is also evidence that C ae either increases with the effective stress until the yield/ preconsolidation pressure and remains constant thereafter [18,33] or increases almost linearly with increasing effective stress [21]. Additionally, there have been some studies on the influence of the physical properties on the C ae [1,53,58,70,77,81]. However, many geotechnical correlations can only be applied in specific soil conditions and only for the limited soil types or regions of the occurrence.
In this study, the time-dependent, viscous properties of alluvial soils were related to the plasticity characteristics. Therefore, the ratio of the plastic limit PL (%) to the liquid limit LL (%) called plasticity ratio R p (-) and the plasticity index PI (%) were considered. This study's primary objective was to investigate the decreasing trend of C ae related to the density of the consolidated soil. Therefore, a double-logarithmic approach was used and its suitability for intact alluvial soils was evaluated. The preliminary results of the long-term multiple-stage oedometer tests on different intact organic soils from the same geological origin were presented in the light of the influence of the plasticity of the soil on the changes in the creep rate. Based on the results, correlational models of the coefficient of secondary compression change index m using plasticity parameters were proposed.

Test materials
This study was conducted using intact samples retrieved from the _ Zuławy Fens located near the urban area of Gdansk, Poland. The study area is shown in Fig. 1. There, a deltaic plain was created by the accumulation processes and filled with quaternary sediments. The _ Zuławy Fens subsoil is extremely heterogeneous because of a very complex depositional history and varies site to site. The geological model for the study area [6] is generalized as shown in Fig. 2. In the model weak mineral-organic soils such as aggradate organic muds with medium compact fine-grained sands (denoted as A1 or A2 due to the depth of occurrence), weak mineral soils such as a fine-grained sands with an admixture of decayed plant debris with silt and silty clay (denoted as B1 or B2 due to the depth of occurrence) and finally alluvial sands (denoted as C) were identified.
The aggradate organic muds are yellow and grey, to black or grey, to brown and composed of admixtures and interlayers of fine sand, silty sand, loamy sand, silty clay and gyttja. Figure 3 shows a sample of the loamy soil retrieved from a borehole in the sampling area. The alluvial soils of the _ Zuławy Fens are characterized by very high variability in organic content and water content. This variability translates into different physical and mechanical properties [13,59]. The term 'alluvial mud' refers to organic fine-grained soil particles transported by water and deposited on the bottom of the riverbed or flooded area as the speed of flowing water decreased. As a formation, the mud (called also ''dy'') occupies an intermediate place between peat and gyttja. It differs from peat in its high degree of plant mass humification, in which there is no macroscopically recognizable plant debris. There is also no roughly fibrous structure typical for peats. The difference between mud and gyttja is that the mass of mud is humus, while the organic mass of gyttja is detritus.

Index properties and classification of soils
To study the influence of the physical properties on the time-dependent behaviour of the fine-grained organic soils, standard index tests were performed. The primary properties (e.g., PL, LL and SL, mineral composition, grain-size distribution) and secondary properties (e.g., water content, density, and consolidation/stress conditions for all the soils) were considered. The LL was determined according to CEN-ISO/TS 17,892-6 [30], and the PL was determined in compliance with ISO/TS 17,892-12 [31]. It should be noted that determining the PL using the thread soil rolling method may have some major limitations [25,37] and is prone to human error [2,3,19,55,64]. Based on PL and LL values, the plasticity index PI = LL -PL (%) and the reference parameter named the plasticity ratio R p = PL/LL (-) [57] were determined. The PI of the soil indicates the size of the range of water contents where the soil exhibits plastic properties. The R p also points to the plasticity of the soil and can be correlated to the specific surface area SSA (m 2 g -1 ), clayey fraction content CF (%) and the liquidity index LI (-) [62] or its equivalent water content ratio WCR (-). Based on the sizable data gathered by Spagnoli and Shimobe [57], the smaller the R p , the higher the plasticity of the soil. The particle-size distribution was identified prior to CEN-ISO- TS 17,, and the specific gravity tests were performed following the procedure in CEN-ISO-TS 17,892-3 [28]. The soil organic matter SOM (%) was determined on the basis of the total carbon content in soil after dry combustion according to ISO 10,694 [27]. In this work the soils were classified based on the Atterberg limits and the grain-size distribution consistent with the European Soil Classification System (ESCS) [38]. The results of the classification were also verified under the criteria given by Moreno-Maroto and Alonso-Azcárate [50]. Following their definition, the term clay refers to the fine-grained material whose PI is equal to or higher than LL/2.

Consolidation tests
Consolidation characteristics were investigated based on the results of the long-term (ten days) multiple-stage loading (MSL 10 ) oedometer tests. The dimension of the specimen was 20 mm high and 60 mm in diameter. The tests with uniform stress distribution and drainage conditions in the vertical direction were done using pervious top and pervious bottom (PTPB). Two loading procedures were followed. For the one-dimensional consolidation tests on samples with the highest plasticity, the loading scheme was 6.25 ? 12.5 ? 25 ? 50 ? 100 ? 200 ? 400 kPa. For the low plasticity samples, the loading scheme was 12.5 ? 25 ? 50 ? 150 ? 250 ? 500 kPa. It was assumed that the compression curves (in semilogarithmic planes) exhibit two regions with compressibility that varied depending on the stress level. In general due to the soil structure, the soil's stress-compression behaviour develops along the recompression/ reloading line (overconsolidation line). The effective vertical stress r' v at the intersection of the overconsolidation line and the normal consolidation line indicates the preconsolidation pressure r' p [10,26,52]. In this study the r' p values are determined on the basis of the bilogarithmic method [10], which utilizes the plot of the lnV and log r' v . To assess the stress and strain rate dependency, C ae was calculated from the linear slope of the e = f(logt) curve beyond the transition from the primary to secondary compression and plotted with the effective vertical stress and the void ratio.
Following the methodology given by Zhu et al. [80], the C ae for each applied stress level was plotted against the e in a double-logarithmic scale for the normally consolidated state. This approach was originally referred to reconstituted Fig. 2 Geological cross-section of the study area. Note: A1 or A2 denotes aggradate organic muds with medium compact fine-grained sands; B1 or B2 denotes fine-grained sands with an admixture of decayed plant debris with silt and silty clay; C denotes alluvial sands  fine-grained soils. However, it is interesting to check whether this approach is also valid for an intact soils. The proposed model assumes consecutive decreases of the C ae as the applied stress level increased and the compression decreased. These changes were attributed to increased density resulting from the compression. Within the projection of the Log C ae -Log e plane, there is a straight line with gradient m for the normally consolidated state (see Fig. 4). In this work parameter m is the coefficient of the secondary compression change index and is used to calculate the change in the creep rate as the void ratio decreases. Its determination is analogous to the permeability change index C k , which illustrates the change in permeability with changes in the e. The parameter m describes the rate at which creep degradation process occurs. High values of the m indicate that the creep degradation is rapidly slowing down. In turn, low values of the m illustrate slow degradation rate of C ae during compression. The adopted formulation for the creep degradation is expressed as follows [74]: In this study, the reference point (C ae ref , e ref ) was arbitrarily selected and corresponded to the soil's LL. Since the reference point is chosen, Eq. (2) can be rewritten as: Application of the given formulation (Eq. 2) implies: (i) a decreasing C ae over time for an applied stress level as the void ratio decreases during creep, (ii) the repeated decreases in C ae with applied stresses as the void ratio is consecutively lowered during loading, (iii) a positive void ratio at the end of the creep phase resulting from compression.

Results and discussion
This section presents the results of the tests performed on fine-grained alluvial soils. Particular emphasis was given to the aspects of compression behaviour, time dependency and correlations between the creep parameters and the plasticity of the soil. Table 1 presents some physical properties of the selected soils. Figure 5 shows the plasticity for each classification under the ESCS, combined with the C-Line and M-Line to further delimit the clays and other less-plastic soils. The PI values between LL/3 and LL/2 were defined as moderately or slightly clayey material. The CF of the analysed samples varied from 9 to 42%. The activity of these soils was evaluated using the relationship between CF and PI in the activity chart (see Fig. 6). Different liquid limits have an impact on determining soil activity [56] and the degree of expansiveness [61]. As a result, the investigated soils were found to be very different in terms of activity. The general expansiveness of these soils was established as low and medium, except for the sandy intermediate plasticity clay (orsaClI (C) sample) and very high plasticity silt (orSiV sample).

Compression behaviour
The void ratio-effective vertical stress and specific volume-effective vertical stress relationships for all soil samples are shown in Fig. 7. Most of the compression curves were characterized by a clearly marked breakdown associated with preconsolidation pressure or the influence of the destructuration process. The consequence of this was a lower constant slope of the curves at r' p values and a significant increase in the constant slope after exceeding this boundary values. In the tested soils, a tendency was observed to increase the slope of the curves with increasing stress, which may result from both the high presence of dispersed organic substance (e.g., sample orSiV and orsa-SiH (B)), high natural water content and the plastic state of relatively young sediments (e.g., sample orClI and orsaClI (B)). The r' p of the ten intact investigated soils ranged from 18 to 70 kPa. In addition, the values of the compression index C c pertain to the normally consolidated region were calculated from the compression curves. The C c is an important compressibility parameter with particular relevance to primary settlement calculations for normally consolidated or lightly overconsolidated soils [45]. In the geotechnical practice the C c is often related to the LL [46]. The C c of the investigated soils for the normally consolidated state is plotted against LL as shown in Fig. 8, and the equation of the straight line: C c-= 0.0096(LL) -0.249 has been obtained. The performance of the empirical model is commonly evaluated with the regression coefficient R 2 defined as follows [66]: In the analysis, the value of R 2 was 0.90 and indicates good predictive ability of the model due to the fact that the compressibility of soft soils closely associates with the mineral compositions which also control the LL of the soil. Figure 9 shows procedure for identifying parameters from the consolidation results of the sample orsaSiH (E). The relationship between C ae and e as an example illustrating the calculation procedure of the m value was also presented. For those purposes the results of three tests on the same soil material were used. As can be seen from Fig. 9a, the threshold stress r' T was found to exist and is corresponded to the maximum value of the C ae . After reaching the r T , the gradual decrease in C ae with the vertical effective stress was observed. In the discussed case r' p was determined by the bilogarithmic method and was about 20 kPa (see Fig. 9b). Thereafter, values of r' T and r' p were compared. From the comparison it was observed Fig. 5 Classification of soils based on liquid limit and plasticity index. Note that symbols Cl and Si are clay and silt, respectively; L, l, H and V are low plasticity, intermediate plasticity, high plasticity and very high plasticity, respectively that the r' T was 5 times the r' p (see Fig. 9c). Moreover, it is not surprising that the C ae decreases with decreasing void ratio e. Then, the relationship between Log C ae and Log e was used for the evaluation of the changes in C ae during compression (see Fig. 9d). As presented for the normally consolidated state the C ae -e was plotted in double log space, where linear relationship was found to exist with the gradient m = 5.48. Parameter m was taken from the equation of that linear trend line. Relatively high value of the m indicates rapid changes in C ae and the fast creep degradation rate. Figure 10 shows the variation of C ae under increasing effective vertical stress owing to the destructuration effect. As can be seen, C ae is stress-dependent and curve of the function C ae = f(logr 0 v ) has a concave downward shape with the maximum at the threshold stress r' T . In general, for each soil sample, the shape of the soil characteristic curves C ae = f(r 0 v ) was developed, and each exhibited a well-defined maximum C ae . In some cases latter part of the curve tended to be stable. With the increase in pressure, the C ae of the orSiV, orClI, orsaSiH (A), orsaSiH (E) and orsaSiH (F) sample increased faster and their C ae was significantly higher than that of other samples, except for orsaSiH (E) sample. For the above-mentioned samples, the C ae reached its peak when the pressure was 100 kPa. For other samples, the peak was reached, when the pressure was 50 kPa. The values of the r' T , for each soil sample, were significantly higher than values of the r' p . Table 2 provides the threshold values of the r' T /r' p ratio for the studied soils. As can be seen, threshold ratio ranged from 1.43 to 5 and indicates that the soil structure has a significant influence on the variation pattern of the C ae rather than yielding or preconsolidation. As a result, this variation is due to the destructuration effect with the peak value of C ae at a void ratio corresponding to the vertical stress wellabove the r' p . Increased consolidation pressure induced the gradual destructuration of the soil structure and, hence, C ae increased. In this pre-peak stage the collapse of the structure occurs with the grains readjustment and the damage to the grain cementation increases the creep [81]. Near the r' T , the soil structure broke down, and the destructuration process slowed (post-peak stage), resulting in a lower C ae value at higher stresses.

Coefficient of secondary compression change index m
The soil data were plotted in a double-logarithmic plane, and the relationship between Log(C ae ) and Log(e) was then evaluated (see Fig. 11). For all tested soils, Log(C ae ) was linearly related to Log(e) in the normally consolidated state. The resulting regression coefficient R 2 was between 0.87 and 0.99. As shown in Fig. 9, the material constant m depended strongly on the type of material (grain-size distribution, PL and LL) and initial void ratio. In general the soils classified as very high plasticity silt, high plasticity clay and sandy high plasticity clay were characterized by low values of m and corresponded to the wide range of LL and PI. However, in the case of sample orsaClI (C) is clearly an exception to this general trend and this soil behaves more like the silty soil than the clayey soil. According to the ESCS classification, the soil is classified as organic sandy intermediate plasticity clay. Following the classification principles, if the fines represent more than or equal to 50% in the sample and LL lies between 35 and 50%, the soil is classified as clay. Moreover, this soil was found to plot slightly above an A-line in the plasticity chart. However, taking into account the definition of clay given by Moreno-Maroto and Alonso-Azcárate [50] this soil does not meet their criteria, because the discussed soil was found to plot below the C-line (PI is not equal to or higher than LL/2). In that context, classification is reinterpreted and the soil may be regarded as a moderately clayey material, whose PI value is between LL/3 and LL/2, what brings its behaviour closer to the silt soil. The smaller initial void ratio and lower compressibility of this sample compared to that of samples orsaClI (A) and orsaClI show that the plasticity of the soil greatly affects a rate at which C ae change during one-dimensional compression.
Subsequently, the averaged C ae was calculated within the normally consolidated range for each soil and plotted against the LL, PI and SOM (see Fig. 12). Figure 12 shows a tendency for an increasing C ae with the increasing LL, PI and SOM. These results were consistent with those of Bjerrum et al. [8], Nakase et al. [53], Yin [70], Jesmani et al. [32] and Zhu et al. [82]. Present study and the results of the above-mentioned researchers showed that plasticity and organic matter to a great extent influence the timedependent behaviour of fine-grained soils.

The effect of the PI and R p on m
To examine the effect of the soil plasticity on m in the Log C ae -Log e relationship, two approaches were considered: (i) m as a function of plasticity index PI. (ii) m as a function of the plasticity ratio R p . Figure 13 shows the plots of m with PI and R p for the available test results on 11 intact alluvial organic soils. It was found that the accurate description of the dependency Fig. 8 Relationship between compression index and liquid limit Fig. 9 An example of the identification parameter for the analysis: a relationship between C ae and effective vertical stress; b bilogarithmic approach for the determination of the preconsolidation pressure; c relationship between C ae and ratio of r' v / r' p ; d double-logarithmic approach for the determination of the reference values of C ae . Note that: symbols e L and e P are void ratios at the liquid limit and plastic limit, respectively; C asL and C aeP are relevant coefficients of secondary compression index Acta Geotechnica (2022) 17:803-817 811 between the plasticity expressed as PI or R p and parameter m is an exponential equation. Therefore, the exponential regression analysis was executed on the data, and the m-PI and m-R p relationships were obtained. As shown in Fig. 13a, the parameter m decreased with the increase in the PI. Thus, soils with considerable size of the range of water contents where the soil exhibits plastic properties have small value of m. In this case, soils with higher plasticity exhibit slow creep degradation rate. From the plot of m with PI, the following relation could be obtained with a correlation coefficient R 2 = 0.82: Similarly, Fig. 13b shows a general trend for the m-R p relationship: with decreasing R p the m value decreases. As can be seen, the m values were lower for the soils with higher plasticity. According to the experimental evidence, the R p values of the tested soil were in the range of 0.4-0.63. This results correspond to the ordinary general sediment soils range of approximately 0.2-0.8 [57]. From the plot of m with R p , the following relation could be achieved with a correlation coefficient R 2 = 0.87: The influence of the plasticity expressed by the R p on the creep degradation rate is significant and larger than those determined with the PI. PI and R p characterize plasticity of the soil and both the high value of the PI and low value of the R p indicate high plasticity. Therefore, in this case resulting C ae degradation rate is low. The value of R 2 in Eq. (6) is relatively higher compared to its value in Eq. (5). From these results, it appears that the choice of R p as the correlating parameter for the m is preferable, because it gives the higher correlation coefficient.
Interestingly, m-PL and m-LL relationships do not give the definite trends and the data are more scattered. Yin et al. [74] correlated the PL with the slopes of Log C ae -Log e plot and observed that for a database of mainly reconstituted Finnish clays the parameter m can be approximately determined by the PL of soils by a power-type formula (R 2 = 0.8423). However, the direct comparison of the test results for soils that are in a different state  (reconstituted and intact) as well as from different geological origin is not appropriate and cannot give clear conclusive remarks. Furthermore, soils could be from the same geological origin; however, with a different R p , they can behave very differently [40]. The differences might be due to the different clay minerals, their proportions and the different clay size proportions. Experimental evidence from the literature shows the clear trend, where lower R p values correspond to higher clay percentage to fines PC (%) values [57]. Figure 14 shows the R p -PC semilogarithmic relationship for the investigated soils. As a result, the smaller is the R p , the higher the PC. Therefore, it can be assumed that the PC also affects the value of m. Hence, m was plotted with PC in the linear plane (see Fig. 15   strongly depends on PC. Considering the types of soil, m is smaller for clay (lower R p and higher PI along with PC) than for silt (higher R p and lower PI along with PC). Figure 16 compares the calculated m values using a double-logarithmic approach to that derived using developed correlations. To examine the robustness and assess the performance of the empirical models (bearing in mind that the R 2 may not give sufficient information to determine the validity of the correlation [66]), the standard error of estimate S y,x is used: S y,x is the measure of the scatter around the regression curve [62] and explains how wrong the regression model is on average, using the units of the response variable. Approximately 95% of the observations should fall within plus/minus 2S y,x from the fitted line, which is also an approximation of a 95% confidence interval. Smaller values are better because it indicates that the observations are closer to the fitted line. The scattering of the points in both studied relationships was clearly along the equality line. The correlations using PI and R p produced values for m close to the experimental values derived from the double-logarithmic approach in the case of soils with higher plasticity. In other words correlations are more suitable for fine-grained soils with the small values of m. The largest discrepancies were observed for the interval 4 \ m \ 6 for the values obtained using R p . Considering data produced with the use of PI, the scattering was high for the m [ 4. It led to the conclusion that an empirical models are more accurate for high plasticity soils; however, it must be used with great care. The invented correlations were established for soils from one area and of the same geological origin. Consequently, more research and more research data are needed to confirm their validity. Nevertheless, the trends found and the relationship between changes in soil creep parameters and plasticity parameters allow for a better understanding of the mechanism of the time-dependent visco-plastic straining.

Conclusions
The time-dependent behaviour of intact alluvial organic soils was investigated by means of the basic soil properties such as PL, LL and PC, which are commonly identified during laboratory analysis. This paper presents the results of experimental research to determine the effect of the soil plasticity on the creep behaviour of natural clayey soils. The results of experiments on 11 types of soil with the same geological origin showed that C ae was both stress and rate-dependent. It was found that soil structure has a significant influence on the pattern of the C ae -r 0 v curve with the characteristic threshold stress which corresponds to the maximum value of C ae . The pattern of the C ae -r 0 v curve has a concave downward shape for all the soils. To model  the creep degradation rate in the normally consolidated states, the double-logarithmic approach was successfully applied. The results of the analysis indicate that the smaller value of the parameter m refers to the low creep degradation rate which typically occurs in clays and clayey soils. It should also be noted that m strongly depends on the plasticity of the soil expressed by PI and R p . The following conclusions were drawn based on the findings of the study: 1. The PI and R p have a significant influence on the value of m. Soils with both higher value of PI and lower value of R p indicated higher plasticity. For the investigated soils, a definite trends were observed: while PI increased, m decreased and while R p decreased, m also decreased. 2. Clay percentage to fines affects soil's plasticity and plays an important role in the creep degradation process. Soils with higher values of PC have smaller values of m and lower creep degradation rates. 3. The relationships are valid only for soils with the same geological origin, and in the limited range of PI and R p . Further investigations should be conducted that include intact fine-grained soils with different geological origins. 4. Considering the type of the soil, low creep degradation rates occur in clays than in silts and other less-plastic soils.

Declarations
Conflicts of interest The author declares that they have no conflicts of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.