Research Article | Open Access

Ping Li, Luanhua Dong, Xiaowen Gao, Tonglu Li, Xiaokun Hou, "An Extension of Taylor’s *φ*-Circle Method and Some Stability Charts for Submerged Slopes", *Advances in Civil Engineering*, vol. 2020, Article ID 8252838, 10 pages, 2020. https://doi.org/10.1155/2020/8252838

# An Extension of Taylor’s *φ*-Circle Method and Some Stability Charts for Submerged Slopes

**Academic Editor:**Chao Wu

#### Abstract

Taylor’s *φ*-circle method is a classical method for slope stability calculation, which has analytical solutions. Taylor derived equations in two cases separately, namely, (i) the outlet of the critical failure surface is at the slope toe and (ii) the outlet of the failure surfaces is not at the slope toe. The method is only appropriate for two conditions (without underground water table in slopes or totally submerged slopes). In this study, a general equation that unifies the equations of the two cases is proposed and partially submerged condition is introduced. The critical failure surfaces corresponding to the minimum factor of safety are determined using the computer program proposed by the authors. The general expression of the safety factor of slopes under the following four conditions is derived, namely, (i) partly submerged, (ii) completely submerged, (iii) water sudden drawdown, and (iv) water slow drawdown. The corresponding charts for practical use are available.

#### 1. Introduction

The *φ*-circle method proposed by Taylor [1] is one of the classical methods of slope stability calculation, which is suitable for homogeneous soil slopes (i.e., without underground water table in slopes or totally submerged slopes). It is assumed that the potential failure surface is a circular arc and the resultant force of the friction and normal stress on the entire failure surface is tangent to the *φ*-circle. The *φ*-circle is a circle which is concentric with the failure surface. The radius of the *φ*-circle is the product of the radius of failure surface arc *R* and sin*φ*. Thus, this method is termed the *φ*-circle method. This *φ*-circle method regards the whole mass above the failure surface as a research object and has strict mathematical derivation and analytical solution. In addition, this method follows both the force and moment equilibriums. In the *φ*-circle method, the dimensionless term *N*, which is related to the cohesion, unit weight, the factor of safety, and the slope height, is termed the “stability number.” The expression for stability number proposed by Taylor [1] is simple and intuitive. The results are easily expressed in the chart. The method is generally used in the homogeneous and isotropic simple slope analysis and has been extended to some special cases, such as spatial variable strength parameters [2, 3], linearly increasing cohesion with depth [4–7], quasistatic seismic effect [8–10], and slope couple foundation stability research [11]. Taylor’s charts have been re-examined and developed using finite-element limit analysis and elastic-plastic finite-element method [12–15]. The advantage of the numerical methods is that they have no assumptions about the shape and the location of the failure surface. It has been found that Taylor’s original stability number is “reasonably accurate” when the failure surface is circular, but differs more than 10% when the failure surface is not circular. Some numerical solutions have been developed to extend the slope stability chart in three dimensional slopes [16–18] and heterogeneous slopes [19].

Cojean and Fleurisson [20] and Cojean and Cai [21] have well documented the effect of water level fluctuation on the factor of safety of slopes. They analyzed the factor of safety under the conditions of static water level and water sudden drawdown in graphical form. Li et al. [22] presented a landslide case on reservoir bank which shows that the factor of safety decreases firstly and then increases with the rise of water level, which agrees with the results of Cojean and Fleurisson [20]. The analytical solution for water drawdown slopes still is not available.

In this study, the original Taylor’s *φ*-circle method is modified, which unifies the original equations in two cases into one general equation. Then, analytical solutions for the factor of safety of slopes subjected to the water drawdown are derived. Four cases are considered, namely, (i) partially submerged slopes, (ii) completely submerged slopes, (iii) water sudden drawdown, and (iv) water slow drawdown. Searching for the critical failure surface according to the minimum factor of safety was carried out using the computer program that the authors proposed. The charts of the factor of safety versus slope inclinations and water levels are provided to illustrate the effect of water level on the slope stability.

#### 2. Modified Taylor’s *φ*-Circle Method

Figure 1 shows a typical slope model. In this model, AB_{1}C_{1} represents the potential sliding mass; *H*_{1} is the slope height; *β* is the slope angle; 2*α*_{1} is the central angle of the failure surface arc AB_{1}; *λ*_{1} is the inclination angle of the chord AB_{1}; and nH_{1} is the distance between the slope toe and the outlet of potential failure surface. When *n* = 0, the outlet of the potential failure surface is just at the slope toe. The radius *r* of the *φ*-circle equals *R*sin* φ*. *φ* is the internal friction of the failure surface.

From the geometry of Figure 1, the radius *R* of the failure surface arc is

According to Taylor’s derivation, the gravity of the potential sliding mass AB_{1}*C*_{1} iswhere *γ*_{1} is the unit weight of the soil above the water level.

Let

Equation (2) is simplified as

The gravity moment of potential sliding mass AB_{1}*C*_{1} iswhere *d*_{1} is the arm of gravity to the center of the *φ*-circle.

Taylor has derived the equations in two cases, namely, (i) the outlet of the failure surface is at the slope toe (i.e., *n* = 0) and (ii) the outlet of the failure surfaces is not at the slope toe (i.e., *n* > 0). When *n* > 0, the failure arc has the minimum factor of safety when the right vertical tangent of the *φ*-circle passes through the midpoint of the slope *C*_{1}*E*. When the center of the failure arc is limited to the vertical line on the left of the midpoint of *C*_{1}*E*, *n* can be calculated with the condition. It will reduce one degree of freedom in searching the failure surface, significantly alleviating the computation work. However, with the development of the computation ability of the computer, it is not difficult to search the failure surface with total degrees of freedom nowadays. In this study, a general equation of moment is derived here which unifies the equations in the two cases, namely, *n* = 0 and *n* > 0. Thus, *n* is one degree of freedom for searching the failure surface, which is convenient for programming.

In Figure 1, AEC_{1}*F* is a parallelogram, whose center of gravity lies in the midpoint of the line AC_{1}. The arm of AEC_{1}*F* to the center of the *φ*-circle can be generally expressed as

The gravity of parallelogram AEC_{1}F is *γ*_{1}, and substituting equation (1) into equation (6), the gravity moment is

Replacing the last term of equation (5) with the equation (7), the general expression for the gravity moment of the sliding mass AB_{1}*C*_{1} is

Let

Then, equation (8) can be simplified as

With modified equation (11), the stability number can be derived in exactly the same way as Taylor’s derivation.

#### 3. The Factor of Safety in Water Drawdown Slopes

The submerged slope is shown in Figure 2. When water level declines from a height of *H*_{2} to *H*_{3} as shown in Figure 3, the corresponding depth of water increases from *D*_{2} to *D*_{3}. Supposing there is no water recharge at the inner side of the slope, the pore water pressure on the failure surface can be expressed by the following equation:where *D* is the depth at the failure surface between the water level before drawdown (*H*_{2}) and water level after drawdown (*H*_{3}), *u* is the pore water pressure at the depth of *D*, *r*_{u} is the pore water pressure coefficient, and is the unit weight of water.

Two extreme situations are considered. One is that the sliding mass is completely impermeable, so the soil between the water level before drawdown (*H*_{2}) and after drawdown (*H*_{3}) should be saturated. The pore water pressure on the failure surface can be calculated by the difference of *H*_{2} and *H*_{3}. The other case is that the sliding mass is completely permeable, and the pore water pressure on the failure surface is approximately zero. The general case is between the two extreme conditions, so the average unit weight of the soil between the water levels before and after drawdown is estimated by interpolating between the saturated unit weight and the unsaturated unit weight aswhere *γ*_{a} is the average unit weight of soil between the water levels before and after drawdown, is the unit weight of water, *e* is the void ratio, is water content of the soil, and *r*_{u} is the pore water pressure coefficient.

When the pore water pressure has not yet dissipated, *γ*_{a} equals saturated unit weight. When the pore water pressure has completely dissipated, *γ*_{a} equals the unsaturated unit weight, which means that the water level in the sliding mass is consistent with the reservoir.

The factor of safety can be obtained by superimposing three independent slope models with the same slope angle *β*. Model 1 is the simple homogenous slope AB_{1}C_{1} with unsaturated unit weight as the original Taylor’s model shown in Figure 1. Model 2 is the slope AB_{2}C_{2}, where *H*_{2} is the slope height, 2*α*_{2} is the central angle of the failure surface arc AB_{2}, and *λ*_{2} is the inclination angle of the chord AB_{2}. The unit weight of AB_{2}C_{2} equals the difference between *γ*_{a} and *γ*_{1}, as shown in Figure 2. Model 3 is the slope AB_{3}C_{3}, where *H*_{3} is the slope height, 2*α*_{3} is the central angle of the failure surface arc AB_{3}, and *λ*_{3} is the inclination angle of the chord AB_{3}. The unit weight of AB_{3}C_{3} equals the difference of *γ*′ (the effective unit weight) and *γ*_{a}, as shown in Figure 3. It is necessary to take the pore water pressure between the water levels before and after drawdown into consideration in model 3. Because the pore water pressure passes through the center of the *φ*-circle, its moment is zero. Superimposing the gravity and gravity moment of the three models, we can obtain the practical gravity W and its moment *W*_{d}. Based on the derivation of Taylor, the analytical expression for the factor of safety in water drawdown slopes can be established.

Deduction of the gravity and gravity moment for the latter two models is listed in Appendix A. The polygon of force equilibrium in the case of water declining is shown in Figure 4, in which there is an added pore water pressure on the failure surface compared with Taylor’s original model. There are four forces acting on the sliding mass, the weight *W*, the resultant cohesion *C*, the resultant normal force *P*′, and the pore water pressure . The magnitude and direction of force *C* are the same as Taylor’s original model. The magnitude of *W* is obtained by equation (A.11) in Appendix A. *P*′ is tangent to the *φ*-circle and has an angle with *W*. The expression for angle is the same as Taylor’s model. passes through the center of *φ*-circle. The resultant of *P*′ and is assumed to be *P*. The angle *µ* can be calculated using equation (14). The detailed deduction of equation (14) and the related terms in equation (14) is listed in Appendix B.

The general expression for the factor of safety in the case of water drawdown is shown in equation (15). The detailed reduction is listed in Appendix C.where

Based on the general form, the factor of safety for six cases can be simplified as below.

For a simple homogenous slope as Taylor’s model, under the conditions of *H*_{2} = *H*_{3} = 0, *K*_{2} = *K*_{3} = 0, = 0, and = , equation (15) can be simplified as

For the completely submerged slope, under the conditions of *H*_{2} = *H*_{3} = *H*_{1}, *K*_{2} = *K*_{3}=1, *r*_{u} = 0, *γ*_{2} = 0, *γ*_{3} = *γ*′ − *γ*_{1}, = 0, and = , equation (15) can be simplified as the same as equation (17) in which *γ*_{1} is replaced with *γ*′.

For the partially submerged slope, under the conditions of *H*_{2} = *H*_{3}, *K*_{2} = *K*_{3}, *r*_{u} = 0, *γ*_{2} = 0, *γ*_{3} = *γ*′−*γ*_{1}, = 0, and = , equation (15) can be simplified as

For the slope where water level declines very slowly or the soil is highly permeable, the water level in the slope is synchronous with the reservoir water level. In this case, the expression for the factor of safety is the same as equation (18).

For the slope that experiences water drawdown with a decrease of water level, the soil in the zone between the water levels before and after the water drawdown is regarded to be saturated due to capillarity. In this case, *γ*_{2} = *γ*_{sat}−*γ*_{1}, *γ*3 = −, = 0, and = , and equation (15) can be simplified as

For the slope having water sudden drawdown, the soil between the water levels before and after drawdown is saturated and the pore water pressure on the failure arc has not dissipated. In this case, *r*_{u} = /*γ*, *γ*_{2} = *γ*_{sat} − *γ*_{1}, and *γ*_{3} = −, and equation (15) can be simplified as

#### 4. Charts for the Factor of Safety

The failure surface arc has three degrees of freedom, namely, *n*, *α*_{1}, and *λ*_{1}. Given reasonable domains for the three degrees of freedom, the minimum factor of safety corresponding to the critical failure surface can be searched with any optimum method such as the golden-section method. In this study, the failure surface is determined using the computer program accomplished by the authors.

Considering the simple homogenous slope first, the factor of safety has a linear relationship with *c*/*γ* H_{1} as shown in Equation (15). Given *c*/*γ H*_{1} = 0.01 and 0.10, respectively, the minimum factors of safety against a range of slope inclinations and internal friction angles are shown as charts in Figure 5.

**(a)**

**(b)**

An example is described to explain how to use the charts. The slope height H1 is 50 m; the slope inclination angle *β* is 30°; the cohesion *c* is 30 kPa; the internal friction *φ* is 25°; and the unit weight *γ* is 17 kN/m^{3}. With the given *β* and *φ* values, *F* is 1.04 when *c*/*γH*_{1} = 0.01 and is 1.89 when *c*/*γH*_{1} = 0.10 based on Figure 5(a) and Figure 5(b). Afterward, *c*/*γH*_{1} is 0.034 based on the given *c*, *γ*, and *H*_{1} values. *F* of the slope is 1.27 by interpolation between 1.04 and 1.89.

In Taylor’s chart, the factor of safety is defined as the ratio of actual cohesion to critical cohesion. To find the reliable factor of safety, it needs to adjust *F* for getting equal values to fit both the friction angle and the stability number, which is a tedious process in using the chart. Figure 5 provides a direct and simple way to determine the factor of safety.

For a submerged slope with the static water level, *F* has a linear relationship with *c*/*H*_{1} as shown in equation (19). Given *c*/*H*_{1} = 1.0 and 2.0, *γ* = 17. 5 kN/m³, *e* = 0.67, = 28%, = 9.81 kN/m³, and *K*_{2} is a series of values between 0 −1.0, the minimum *F* with a range of *β* and *φ* is calculated. The results are presented as charts in Figure 6. It suggests that for the steep slope (i.e., *β* ≥ 60°), the factor of safety increases monotonously with the rise of water level; for the gentle slopes, the factor of safety decreases at the beginning and then increases with the rise of water level.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

For the case of water sudden drawdown, assume *H*_{1} = 50 m, *γ* = 16 kN/m³, *e* = 0.43, = 28%, = 9.81 kN/m³, *r*_{u} = 0, = 15%, *β* = 30°, *φ* = 20°, and a range of original water levels *H*_{2} and the terminal water levels *H*_{3} between 10 m and 45 m. Calculations are carried out with the orthogonal values of *K*_{2} and *K*_{3}. The results are presented in Figure 7. The top curve shows the factor of safety against static water where *K*_{2} = *K*_{3}. The curves that intersect with the top curve show the factor of safety against the water drawdown from *K*_{2} to *K*_{3}. It suggests in this example that the factor of safety increases with the rise of the static water level. However, it decreases with the water sudden drawdown for all the original water levels. The larger the water level drawdown is, the lower the factor of safety is. The results are in agreement with the results by Cojean and Fleurisson [20] and Cojean and Caï [21].

The results in Figures 6 and 7 are calculated under given parameters. It is difficult to make general charts for the factor of safety for all slopes. However, it is easy and practical to figure out using the computer program.

#### 5. Conclusions

Based on Taylor’s *φ*-circle method, a general expression for the factor of safety of the slope in the water drawdown case is proposed. Some typical cases can be simplified with the general expression, such as homogenous, partly submerged, and completely submerged slopes, as well as water sudden drawdown and slow drawdown slopes. A computer program written by the authors is used to search the failure surface and figure out the minimum factor of safety.

For simple homogenous slopes, the charts are designed for determining the factor of safety with aid of interpolation between two values of *c*/*γH*_{1} that has a linear relation with *F*. It avoids the tedious process of iteration with Taylor’s stability number chart.

For the slopes with static water level or water drawdown, the factor of safety is a linear function with *c*/*H*_{1}. It is difficult to give a set of general charts due to the additional free variables such as unit weights above and below the water level. However, it is easy to obtain using the proposed computer programs. Some typical cases have been calculated and shown as charts in the paper.

#### Appendix

#### A. Gravity and Gravity Moment Calculation for the Slope of Water Drawdown

For the simple slope model shown in Figure 1, the gravity and gravity moment of the sliding mass have been derived by Taylor [1] and given in the text as equations (4) and (11).

For the model shown in Figure 2, the unit weight of the sliding mass is

Let

From Figure 1, it can be seen that

Then,

For the model shown in Figure 3, the unit weight is

Similarly, we have

Let

Similar to equations (4) and (11) in the text, we have

Referring to equations (9) and (10) in the text, we have

Superimposing the gravity and gravity moment of the three models, we can obtain the practical gravity and gravity moment.

Substituting equation (A.11) into equation (A.12), we get

#### B. Deduction of the Angle *µ* between *P*′ and *P*

Based on equation (12) in the text, the resultant pore water pressure on the failure surface between water levels before and after drawdown can be estimated bywhere *l*_{2} (*l*_{3}) is the arc length of the failure surface from point *B*_{1} to the corresponding water depth *D*_{2} (*D*_{3}), 2*θ* is the center angle corresponding to the depth *D*, and 2*θ*_{p} is the center angle corresponding to the point of the resultant force .

From Figure 3, it results that

D2 and D3 can be calculated with equation (B.1) as *θ* is given the value as *θ*_{2} or *θ*_{3}.

The acting position of is approximately on (*H*_{2}−*H*_{3})/3 above the water level after drawdown or (2*D*_{3}−*D*_{2})/3 below the top ground of the slope as shown in Figure 3. (2*D*_{3}−*D*_{2})/3 can be expressed as

So,

Substituting equations (B.2), (B.3), and (B.5) into equation (B.1), can be deduced by integration and the result is

Referring to Figure 3, there is

So, the angle *δ* between and *C* is

By the cosine law in triangle ABD in Figure 4,

By the sinusoidal law in triangle ADE in Figure 4,

Substituting equation (B 9) into equation (B.10), we get

By the cosine law and sinusoidal law in triangle BDE in Figure 4, there are

So, equation (B 13) becomes

Substituting equations (B 13) and (B 16) into equation (B 14), we get

#### C. Deduction of the Factor of Safety

A common definition of the factor of safety is the ratio of the shear strength of the soil to the shear strength mobilized to retain the limit equilibrium:in which *F* is the factor of safety, *c* is the cohesion, *c*_{d} is the mobilized soil cohesion, *φ* is the internal friction angle, and *φ*_{d} is the mobilized internal friction angle.

The mobilized internal angle is

Referring to Taylor [1], we havewhere

Using sinusoidal law and referring to Figure 4, there iswhereIn which

Substituting equation (C 4) into (C 3) and then substituting equations (C 3) and (A 11) into (C 1), the factor of safety can be deduced as

The angle *µ* in equation (C 2) is deduced in Appendix B.

#### Data Availability

The data used to support the findings of this study are included in the article.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The authors acknowledge the funding received from the National Natural Science Foundation of China (program nos. 41877242 and 41772278) and the National Key R&D Program of China (2017YFC1501302) which supported this study.

#### References

- D. W. Taylor, “Stability of earth slopes,”
*Journal of the Boston Society of Civil Engineers*, vol. 24, no. 3, pp. 197–246, 1937. View at: Google Scholar - D. V. Griffiths, J. Huang, and G. A. Fenton, “Influence of spatial variability on slope reliability using 2-D random fields,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 135, no. 10, pp. 1367–1378, 2009. View at: Publisher Site | Google Scholar - H. Zhu, D. V. Griffiths, G. A. Fenton, and L. M. Zhang, “Undrained failure mechanisms of slopes in random soil,”
*Engineering Geology*, vol. 191, pp. 31–35, 2015. View at: Publisher Site | Google Scholar - R. E. Gibson and N. Morgenstern, “A note on the stability of cuttings in normally consolidated clays,”
*Géotechnique*, vol. 12, no. 3, pp. 212–216, 1962. View at: Publisher Site | Google Scholar - J. H. Hunter and R. L. Schuster, “Stability of simple cuttings in normally consolidated clays,”
*Géotechnique*, vol. 18, no. 3, pp. 372–378, 1968. View at: Publisher Site | Google Scholar - D. V. Griffiths and X. Yu, “Another look at the stability of slopes with linearly increasing undrained strength,”
*Géotechnique*, vol. 65, no. 10, pp. 824–830, 2015. View at: Publisher Site | Google Scholar - P. Li, Q. Zhang, X. Hou et al., “Stability analysis of an undrained soil slope against the linearly increasing shear strength,”
*Hydrogeology & Engineering Geology*, vol. 43, no. 3, pp. 59–65, 2016, in Chinese. View at: Google Scholar - R. Baker, R. Shukha, V. Operstein, and S. Frydman, “Stability charts for pseudo-static slope stability analysis,”
*Soil Dynamics and Earthquake Engineering*, vol. 26, no. 9, pp. 813–823, 2006. View at: Publisher Site | Google Scholar - P. P. Sahoo and S. K. Shukla, “Taylor’s slope stability chart for combined effects of horizontal and vertical seismic coefficients,”
*Géotechnique*, vol. 69, no. 4, pp. 344–354, 2019. View at: Publisher Site | Google Scholar - C. Sun, J. Chai, B. Ma, T. Luo, Y. Gao, and H. Qiu, “Stability charts for pseudostatic stability analysis of 3D homogeneous soil slopes using strength reduction finite element method,”
*Advances in Civil Engineering*, vol. 2019, Article ID 6025698, 18 pages, 2019. View at: Publisher Site | Google Scholar - W. Liu, T. Li, and P. Li, “Numerical analysis for failure mode and ultimate bearing capacity of striping footings on the crest of undrained slopes,”
*Journal of Engineering Geology*, vol. 24, no. 2, pp. 197–203, 2016, (in Chinese). View at: Google Scholar - R. Baker and Y. Tanaka, “A convenient alternative representation of Taylor’s stability chart,”
*Slope Stability Engineering*, vol. 1, pp. 253–257, 1999. View at: Google Scholar - R. Baker, “A second look at taylor’s stability chart,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 129, no. 12, pp. 1102–1108, 2003. View at: Publisher Site | Google Scholar - T. Steward, N. Sivakugan, S. K. Shukla, and B. M. Das, “Taylor’s slope stability charts revisited,”
*International Journal of Geomechanics*, vol. 11, no. 4, pp. 348–352, 2011. View at: Publisher Site | Google Scholar - D. V. Griffiths and C. M. Martin, “Critical failure mechanisms in relatively flat undrained slopes,”
*Géotechnique Letters*, pp. 1–5, 2019. View at: Google Scholar - R. L. Michalowski and T. Martel, “Stability charts for 3D failures of steep slopes subjected to seismic excitation,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 137, no. 2, pp. 183–189, 2011. View at: Publisher Site | Google Scholar - K. Lim, A. V. Lyamin, M. J. Cassidy et al., “Three-dimensional slope stability charts for frictional fill materials placed on purely cohesive clay,”
*International Journal of Geomechanics*, vol. 16, no. 2, Article ID 04015042, 2016. View at: Publisher Site | Google Scholar - C. Sun, J. Chai, Z. Xu et al., “3D stability charts for convex and concave slopes in plan view with homogeneous soil based on the strength-reduction method,”
*International Journal of Geomechanics*, vol. 17, no. 5, Article ID 06016034, 2017. View at: Publisher Site | Google Scholar - R. L. Michalowski, “Stability charts for uniform slopes,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 128, no. 4, pp. 351–355, 2002. View at: Publisher Site | Google Scholar - R. Cojean and J. A. Fleurisson, ““Influence de la structure géologique sur la stabilité de versants progressivement submergés par la montée d’un plan d’eau,” in
*Proceedings of the 6th International IAEG Congress*, pp. 2113–2132, Amsterdam, Netherlands, August 1990. View at: Google Scholar - R. Cojean and Y. J. Caï, “Analysis and modeling of slope stability in the Three-Gorges Dam reservoir (China) - the case of Huangtupo landslide,”
*Journal of Mountain Science*, vol. 8, no. 2, pp. 166–175, 2011. View at: Publisher Site | Google Scholar - T. Li, C. Zhang, P. Xu et al., “Stability ssessment and tabilizing pproaches for the Majiagou Landslide, ndergoing the ffects of ater evel luctuation in the Three Gorges Reservoir Area,”
*Landslide Disaster Mitigation in the Three Gorges Reservoir, China*, Springer, Berlin, Germany, 2009. View at: Google Scholar

#### Copyright

Copyright © 2020 Ping Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.