Structural investigations of (Ca,Sr)ZrO<Subscript>3</Subscript> and Ca(Sn,Zr)O<Subscript>3</Subscript> perovskite compounds

(Ca _x ,Sr_1?x )ZrO₃ and Ca(Sn _y ,Zr_1-y )O₃ solid solutions were synthesized by solid-state reaction at high temperature before to be studied by powder X-ray diffraction and Raman Spectroscopy. Diffraction data allow the distortion of the ABO₃ perovskite structure to be investigated according to cations substitution on A and B-sites. It is shown that distortion, characterized by Φ, the tilt angle of BO₆ octahedra, slightly increases with decreasing y content in Ca(Sn _y ,Zr_1?y )O₃ compounds and strongly decreases with decreasing x content in (Ca _x ,Sr_1?x )ZrO₃ compounds. Such results are discussed in view of the relative A and B cation sizes. Raman data show that vibrational spectra are strongly affected by the cation substitution on A-site; the frequencies of most vibrational modes increase with increasing x content in (Ca _x ,Sr_1?x )ZrO₃ compounds, i.e. with the decreasing mean size of the A-cation; the upper shift is observed for the 358 cm^?1 mode (?ν/?r = ?60.1 cm^?1/Å). On the other hand, the cation substitution on B-sites, slightly affect the spectra; it is shown that in most cases, the frequency of vibrational modes increases with increasing y content in Ca(Sn _y ,Zr_1?y )O₃ compounds, i.e. with the decreasing mean size of the B-cation, but that two modes (287 and 358 cm^?1) behave differently: their frequencies decrease with the decreasing mean size of the B-cation, with a shift respectively equal to +314 and +162 cm^?1/Å. Such results could be used to predict the location of different elements such as trivalent cations or radwaste elements on A- or B-site, in the perovskite structure.     

Permeability of granular soil employing flexible wall permeameter

Understanding the changes in permeability of soil, when soil is subjected to high confining pressure and flow pressure, which may alter the textural and geomechanical characteristics of soil, is of great importance to many geo-engineering activities such as, construction of high-rise buildings near the coast or the water bodies, earthen dams, pavement subgrades, reservoir, and shallow repositories. It is now possible to evaluate the changes in permeability of soil samples under varying conditions of confining pressure and flow pressure using flexible wall permeameter (FWP). In the present study, investigation was carried out on a cylindrical sample of granular soil employing FWP under varied conditions of confining pressure (σ₃)—50–300 kPa, which can simulate the stress conditions equivalent to depth of about 20 m under the earth’s crust, and a flow pressure (f_p)—20–120 kPa, which is mainly present near the small earthen embankment dams, landfill liners, and slurry walls near the soft granular soil with high groundwater table. The obtained results indicate a linear relationship between hydraulic conductivity (k) with effective confining pressure (σ_eff.), k, decreasing linearly with an incremental change in σ_eff.. Further, k increases significantly with an increase in f_p corresponding to each σ_eff., and q increases significantly with increase in the f_p corresponding to each (σ₃). It was also observed that corresponding to the low f_p of 20 kPa, the reduction in k is nonlinear with σ_3. The percentage reduction in k is observed to be 9, 13, and 27% corresponding to σ₃ of 50–100, 100–200, and 200-300 kPa, respectively.     

Dynamic response characteristics of a rock slope with discontinuous joints under the combined action of earthquakes and rapid water drawdown

In order to study the dynamic response characteristics of a rock slope with discontinuities under the combined action of earthquakes and rapid water drawdown, a large-scale shaking table test was performed on a rock slope with discontinuous joints. Wenchuan earthquake (WE) seismic records were performed to investigate the horizontal and vertical acceleration response and displacement response. In particular, three-dimensional optical measurement techniques was used to obtain the slope surface displacements. A comparison was made on the seismic response according to the analysis of PGD (peak ground displacement) and M _PGA (acceleration amplification coefficient) of the modeled slope. The results show that the experimental slope mainly underwent settlement and horizontal deformation when the WE records were applied in the z and x directions, respectively. The slope was first shaken by the P wave, which caused the differential settlement to occur at the surface slope; then, the slope was shaken more severely by the S wave, which led to a greater horizontal deformation. Moreover, analysis of the ΔPGD (increment of PGD) and ΔM _PGA (increment of M _PGA) under rapid drawdown suggests that the rapid water drawdown mainly impacts the deformation of the surface slope, particularly between the high and low water levels. The water infiltration through the cracks softened the material of the surface slope, and the rapid drawdown also enhanced the slope deformation. In addition, the damage evolution process of the slope can be identified, mainly including three stages: an elastic stage (<?0.168 g), a plastic stage (0.168–0.336 g), and a failure stage (>?0.336 g).     

Possible reasons for the frequency splitting of the harmonics of type II solar radio bursts

AIA/SDO data in the 193 Å channel preceding a coronal mass ejection observed at the solar limb on June 13, 2010 are used to simultaneously identify and examine two different shock fronts. The angular size of each front relative to the CME center was about 20°, and their propagation directions differed by ≈25° (≈4° in position angle). The faster front, called the blast shock, advanced the other front, called the piston shock, by R ≈ (0.02-0.03)R⊙ (R⊙ is the solar radius) and had a maximum initial speed of V_B ≈ 850 km/s (with V_P ≈ 700 km/s for the piston shock). The appearance and motion of these shocks were accompanied by a Type II radio burst observed at the fundamental frequency F and second harmonic H. Each frequency was split into two close frequencies f₁ and f₂ separated by Δf = f₂ - f₁ ? F, H. It is concluded that the observed frequency splitting Δf of the F and H components of the Type II burst could result from the simultaneous propagation of piston and blast shocks moving with different speeds in somewhat different directions displaying different coronal-plasma densities.     

Structural systematics of Ge substitution in primitive pyroxenes

Pyroxenes of general stoichiometry Mg(Ge _x Si_1?x )O₃ were encountered in attempts to synthesise Ge-substituted talcs at 0.2 GPa, 650–700 °C. Orthopyroxenes (Pbca) of compositions x = 0.21, 0.30, and 0.34 were identified, and also a P2₁/c clinopyroxene of composition x = 0.63, and C2/c clinopyroxenes of compositions x = 0.91 and 1. End-member clinoenstatite MgSiO₃-P2₁/c synthesised at 16 GPa, 1300 °C and transformed from C2/c was also included in the study. Crystal structure refinements using single-crystal XRD data showed that unit-cell parameters vary linearly with Si–Ge for the Pbca and P2₁/c pyroxenes, both of which have two symmetrically non-equivalent tetrahedral chains. Refinement of Si–Ge occupancies at tetrahedral sites showed that the two chains of all primitive pyroxenes have very different compositions, with X_Ge(TB) ? X_Ge(TA). This difference arises from the greater flexibility of the B-chain to rotate in response to tetrahedral expansion due to increasing Ge content. The TA-M2 shared polyhedral edge imposes significant constraints on the flexibility of the A-chain, which can accommodate much less Ge than the B-chain. Linear trends of cell parameters, site occupancies, and structural parameters for the primitive pyroxenes, when extrapolated to published data for MgGeO₃–Pbca, extend across the entire Si–Ge join.     

Absorption in the H(93–95)<Emphasis Type="Italic">α</Emphasis> and H(78–79)<Emphasis Type="Italic">α</Emphasis> recombination lines in the HII region of the quasar 3C 345

Absorption of the synchrotron emission of the quasar 3C 345 in the continuum and H(93–95)α and H(78–79)α radio recombination lines is studied. The upper limit for absorption in the H(93–95)α lines is T_al/T_ac < 0.7%; absorption in the H(78–79)α lines with antenna temperature T_al = 25 mK, linewidth Δf = 5.3 ± 0.08 MHz, and T_al/T_ac ≥ 0.3% has been detected. A correction to the redshift Δz = 0.00135 ± 0.00008 (z = 0.59365) has been determined.     

Characteristic scale of heterogeneity of seismically active fault and its manifestation in scaling of earthquake source spectra

Previously, similarity of source spectra of Kamchatka earthquakes with respect to the common corner frequency f_c1 and the expressed deviations from similarity for the second f_c2 and the third f_c3 corner frequencies were revealed. The value of f_c3 reflects the characteristic size L_is of fault surface; correspondingly, L_is ≈ v_rT_is, where v_r is the rupture speed and T_is ≈ 1/f_c3 is characteristic time. The estimates of f_c3 are used for normalizing f_c1 and f_c2. In this way one obtains dimensionless rupture temporal parametres τ₁ and τ₂ and can further study the dependence τ₂ (τ₁). The growth of a rupture is considered as a process of aggregation of elementary fault spots of the size L_is. The dimensionless width of the random front of aggregation is on the order of τ₂. The relationship τ₂ (τ₁) approximately follows power law with exponent β. The estimates of β derived from earthquake populations of Kamchatka, USA and Central Asia (β = 0.35–0.6) agree with values expected from the known Eden’s theory of random aggregation growth and from its generalizations.     

Enhancement of random finite element method in reliability analysis and risk assessment of soil slopes using Subset Simulation

Random finite element method (RFEM) provides a rigorous tool to incorporate spatial variability of soil properties into reliability analysis and risk assessment of slope stability. However, it suffers from a common criticism of requiring extensive computational efforts and a lack of efficiency, particularly at small probability levels (e.g., slope failure probability P _f ?<?0.001). To address this problem, this study integrates RFEM with an advanced Monte Carlo Simulation (MCS) method called “Subset Simulation (SS)” to develop an efficient RFEM (i.e., SS-based RFEM) for reliability analysis and risk assessment of soil slopes. The proposed SS-based RFEM expresses the overall risk of slope failure as a weighed aggregation of slope failure risk at different probability levels and quantifies the relative contributions of slope failure risk at different probability levels to the overall risk of slope failure. Equations are derived for integrating SS with RFEM to evaluate the probability (P _f ) and risk (R) of slope failure. These equations are illustrated using a soil slope example. It is shown that the P _f and R are evaluated properly using the proposed approach. Compared with the original RFEM with direct MCS, the SS-based RFEM improves, significantly, the computational efficiency of evaluating P _f and R. This enhances the applications of RFEM in the reliability analysis and risk assessment of slope stability. With the aid of improved computational efficiency, a sensitivity study is also performed to explore effects of vertical spatial variability of soil properties on R. It is found that the vertical spatial variability affects the slope failure risk significantly.     

Effect of Slope on <Emphasis Type="Italic">p</Emphasis>-<Emphasis Type="Italic">y</Emphasis> Curves for Laterally Loaded Piles in Soft Clay

Pile foundations are often subject to lateral loading due to various forces on a variety of structures like high rise buildings, transmission towers, power stations, offshore structures and highway and railway structures. The present investigation is to study the effect of slopes on p-y curves (where p is the static soil reaction and y is the pile deflection) due to static lateral loading in soft clay (Consistency index I_c = 0.42). A series of laboratory model tests were carried out on the instrumented model pile on sloping ground (slopes of 1V:1H, 1V:1.5H, 1V:2H, 1V:3H and 1V:5H) and with varying embedment length to diameter ratio (L/D) of 20, 25 and 30. From the experimental results, the bending moment curves along the pile shaft are double differentiated to obtain the soil resistance (p) and double integrated to obtain the deflection (y) using curve fitting method. New p-y curves for piles located on crest of soft clay with different sloping ground surface under static lateral loading are developed. Moreover, the effect of sloping angles on proposed p-y curves was studied.     

On the efficient estimation of small failure probability in slopes

The random finite element method (RFEM) combines the random field theory and finite element method in the framework of Monte Carlo simulation. It has been applied to a wide range of geotechnical problems such as slope stability, bearing capacity and the consolidation of soft soils. When the RFEM was first developed, direct Monte Carlo simulation was used. If the probability of failure (p _f ) is small, the direct Monte Carlo simulation requires a large number of simulations. Subset simulation is one of most efficient variance reduction techniques for the simulation of small p _f . It has been recently proposed to use subset simulation instead of direct Monte Carlo simulation in RFEM. It is noted, however, that subset simulation requires calculation of the factor of safety (FS), while direct Monte Carlo requires only the examination of failure or non-failure. The search for the FS in RFEM could be a tedious task. For example, the search for the FS of slope stability by the strength reduction method (SRM) usually requires much more computational time than a failure or non-failure checking. In this paper, the subset simulation is combined with RFEM, but the need for the search of FS is eliminated. The value of yield function in an elastoplastic finite element analysis is used to measure the safety margin instead of the FS. Numerical experiments show that the proposed approach gives the same level of accuracy as the traditional subset simulation based on FS, but the computational time is significantly reduced. Although only examples of slope stability are given, the proposed approach will generally work for other types of geotechnical applications.     

Laboratory parallel-beam transmission X-ray powder diffraction investigation of the thermal behavior of nitratine NaNO<Subscript>3</Subscript>: spontaneous strain and structure evolution

Present work provides in-situ structural data at a fine temperature scale from RT to the melting point of nitratine, NaNO₃. From the analysis of log e ₃₃ versus log t plots, it is possible to prove that an univocal indication on the R \( \overline{3} \) c (low temperature, LT) → R \( \overline{3} \) m (high temperature, HT) transition mechanism cannot be obtained because of the relevant role played by the arbitrary assumptions required for defining the c ₀ dependence from temperature of the HT phase. This is due to the occurrence of excess thermal expansion for the HT phase. A significantly better fit for an Ising-spin structural model over a non-Ising rigid-body one has been obtained for the LT phase. Moreover, the Ising model led to a smooth variation of the oxygen site x fractional coordinate throughout the transition. The structure of the HT polymorph has been successfully refined considering an oxygen site at x, 0, ½, with 50% occupancy. Such model was the only acceptable one from the crystal chemical point of view as the alternative model (oxygen site at x, y, z with 25% occupancy) led to unrealistically aplanar \( {\text{NO}}_{3}^{ - } \) groups.     

Quantification of photosynthetic inorganic carbon utilisation via a bidirectional stable carbon isotope tracer

The amount of bicarbonate utilised by plants is usually ignored because of limited measurement methods. Accordingly, this study quantified the photosynthetic assimilation of inorganic carbon (CO₂ and HCO₃ ^?) by plants. The net photosynthetic CO₂ assimilation (P _N), the photosynthetic assimilation of CO₂ and bicarbonate (P _N’), the proportion of increased leaf area (f _LA) and the stable carbon isotope composition (δ¹³C) of Orychophragmus violaceus (Ov) and Brassica juncea (Bj) under three bicarbonate levels (5, 10 and 15 mm NaHCO₃) were examined to determine the relationship among P _N, P _N’ and f _LA. P _N’, not P _N, changed synchronously with f _LA. Moreover, the proportions of exogenous bicarbonate and total bicarbonate (including exogenous bicarbonate and dissolved CO₂-generated bicarbonate) utilised by Ov were 2.27 % and 5.28 % at 5 mm bicarbonate, 7.06 % and 13.28 % at 10 mm bicarbonate, and 8.55 % and 17.31 % at 15 mm bicarbonate, respectively. Meanwhile, the proportions of exogenous bicarbonate and total bicarbonate utilised by Bj were 1.77 % and 3.28 % at 5 mm bicarbonate, 2.11 % and 3.10 % at 10 mm bicarbonate, and 2.36 % and 3.09 % at 15 mm bicarbonate, respectively. Therefore, the dissolved CO₂-generated bicarbonate and exogenous bicarbonate are important sources of inorganic carbon for plants.     

CO<Subscript>2</Subscript> content of andesitic melts at graphite-saturated upper mantle conditions with implications for redox state of oceanic basalt source regions and remobilization of reduced carbon from subducted eclogite

We have performed experiments to determine the effects of pressure, temperature and oxygen fugacity on the CO₂ contents in nominally anhydrous andesitic melts at graphite saturation. The andesite composition was specifically chosen to match a low-degree partial melt composition that is generated from MORB-like eclogite in the convective, oceanic upper mantle. Experiments were performed at 1–3 GPa, 1375–1550?°C, and fO₂ of FMQ ?3.2 to FMQ ?2.3 and the resulting experimental glasses were analyzed for CO₂ and H₂O contents using FTIR and SIMS. Experimental results were used to develop a thermodynamic model to predict CO₂ content of nominally anhydrous andesitic melts at graphite saturation. Fitting of experimental data returned thermodynamic parameters for dissolution of CO₂ as molecular CO₂: ln(K ⁰) = ?21.79?±?0.04, ΔV ⁰?=?32.91?±?0.65 cm³mol^?1, ΔH ⁰?=?107?±?21 kJ mol^?1, and dissolution of CO₂ as CO₃ ^2?: ln(K ⁰ ) = ?21.38?±?0.08, ΔV ⁰?=?30.66?±?1.33 cm³ mol^?1, ΔH ⁰?=?42?±?37 kJ mol^?1, where K ⁰ is the equilibrium constant at some reference pressure and temperature, ΔV ⁰ is the volume change of reaction, and ΔH ⁰ is the enthalpy change of reaction. The thermodynamic model was used along with trace element partition coefficients to calculate the CO₂ contents and CO₂/Nb ratios resulting from the mixing of a depleted MORB and the partial melt of a graphite-saturated eclogite. Comparison with natural MORB and OIB data suggests that the CO₂ contents and CO₂/Nb ratios of CO₂-enriched oceanic basalts cannot be produced by mixing with partial melts of graphite-saturated eclogite. Instead, they must be produced by melting of a source containing carbonate. This result places a lower bound on the oxygen fugacity for the source region of these CO₂-enriched basalts, and suggests that fO₂ measurements made on cratonic xenoliths may not be applicable to the convecting upper mantle. CO₂-depleted basalts, on the other hand, are consistent with mixing between depleted MORB and partial melts of a graphite-saturated eclogite. Furthermore, calculations suggest that eclogite can remain saturated in graphite in the convecting upper mantle, acting as a reservoir for C.     

Determination of the enthalpies of formation of some platinum antimonides and their phase diagrams under standard conditions

Using the method of direct synthesis calorimetry, we determined the standard enthalpy of formation of PtSb (stumpflite), Δ _f H°_298.15 (PtSb, cr) =–105.16 ± 0.84 kJ/mol and PdSb₂ (geversite), Δ _f H°_298.15 (PtSb₂,cr) =–160.92 ± 0.84 kJ/mol. Isothermal (298.15 K, p = 1 bar) phase diagrams were computed for the Pt–Sb–S and Pt–Sb–O ternary systems in the coordinates composition of the Pt–Sb binary system versus fugacity of a gaseous volatile component (O₂, S₂).     

Assessing the vegetation canopy influences on wind flow using wind tunnel experiments with artificial plants

Wind erosion causes serious problems and considerable threat in most regions of the world. Vegetation on the ground has an important role in controlling wind erosion by covering soil surface and absorbing wind momentum. A set of wind tunnel experiments was performed to quantitatively examine the effect of canopy structure on wind movement. Artificial plastic vegetations with different porosity and canopy shape were introduced as the model canopy. Normalized roughness length (Z ₀/H) and shear velocity ratio (R) were analyzed as a function of roughness density (λ). Experiments showed that Z ₀/H increases and R decreases as λ reaches a maximum value, λ _max, while the values of Z ₀/H and R showed little change with λ value beyond as λ _max.     

Water determination in iron oxyhydroxides and iron ores by Karl Fischer titration

Protohematite (Fe_2?x/3(OH) _x O_3?x 1 ≤ x < 0.5) and hydrohematite (Fe_2?x/3(OH) _x O_3?x 0.5 ≤ x < 0) are iron-defective phases containing hydroxyl groups in their structures. These species were described in prior studies mainly with the aid of X-ray diffraction and Infrared spectroscopy. The existence of these phases in soils might have influence in redox processes, and they were considered as a possible water reservoir in Martian soils. In this study, we have used for the first time the Karl Fischer titration method to determine the amount of water released after heating several synthetic samples of goethite, hematite and natural iron ores at 105, 400, 600 and 900 °C. It was found that heating at 105 °C did not remove all moisture from the samples, and higher temperatures were necessary to completely remove all the absorbed water. The water contents determined at 400, 600 and 900 °C were found to be the same within the experimental errors, suggesting the inexistence of both protohematite and hydrohematite in the investigated samples. Therefore, the above-mentioned effects of these phases in soils might have to be reevaluated.     

Behavior of 10-Å phase at low temperatures

The thermal evolution of 10-Å phase Mg₃Si₄O₁₀(OH)₂·H₂O, a phyllosilicate which may have an important role in the storage/release of water in subducting slabs, was studied by X-ray single-crystal diffraction in the temperature range 116–293 K. The lattice parameters were measured at several intervals both on cooling and heating. The structural model was refined with intensity data collected at 116 K and compared to the model refined at room temperature. As expected for a layer silicate on cooling in this temperature range, the a and b lattice parameters undergo a small linear decrease, α _a  = 1.7(4) 10^?6 K^?1 and α _b  = 1.9(4) 10^?6 K^?1, where α is the linear thermal expansion coefficient. The greater variation is along the c axis and can be modeled with the second order polynomial c _T  = c ₂₉₃(1 + 6.7(4)10^?5 K^?1ΔT + 9.5(2.5)10^?8 K^?2(ΔT)²) where ΔT = T ? 293 K; the monoclinic angle β slightly increased. The cell volume thermal expansion can be modeled with the polynomial V _T  = V ₂₉₃ (1 + 8.0 10^?5 K^?1 ΔT + 1.4 10^?7 K^?2 (ΔT)²) where ΔT = T ? 293 is in K and V in ų. These variations were similar to those expected for a pressure increase, indicating that T and P effects are approximately inverse. The least-squares refinement with intensity data measured at 116 K shows that the volume of the SiO₄ tetrahedra does not change significantly, whereas the volume of the Mg octahedra slightly decreases. To adjust for the increased misfit between the tetrahedral and octahedral sheets, the tetrahedral rotation angle α changes from 0.58° to 1.38°, increasing the ditrigonalization of the silicate sheet. This deformation has implications on the H-bonds between the water molecule and the basal oxygen atoms. Furthermore, the highly anisotropic thermal ellipsoid of the H₂O oxygen indicates positional disorder, similar to the disorder observed at room temperature. The low-temperature results support the hypothesis that the disorder is static. It can be modeled with a splitting of the interlayer oxygen site with a statistical distribution of the H₂O molecules into two positions, 0.6 Å apart. The resulting shortest O_bas–O_W distances are 2.97 Å, with a significant shortening with respect to the value at room temperature. The low-temperature behavior of the H-bond system is consistent with that hypothesized at high pressure on the basis of the Raman spectra evolution with P.     

Stokes-meter observations of the solar mean magnetic field: Probable manifestations of strong small-scale magnetic fields

For many years, information on the solar mean magnetic field (SMMF) of the Sun—an important heliophysical and astrophysical parameter—was restricted to magnetographic measurements in only one spectral line, FeI λ525.02 nm. More informative observations of the Stokes-meter parameters of the SMMF were first initiated on a regular basis at the Sayan Solar Observatory. The availability of I and V data obtained simultaneously in several spectral lines has made it possible to study fundamentally new physical problems. In this paper, based on a comparison of SMMF observations in several spectral lines, we find high correlations in the data and important systematic differences in the magnetic-field strength B, which we interpret as a manifestation of kilogauss magnetic fields in fine-structure magnetic elements. Results of theoretical modeling of the SMMF strength ratios for the FeI λ525.02 nm-FeI λ524.70 nm and FeI λ630.15 nm-FeI λ630.25 nm lines are presented. The asymmetries of the V profiles of four lines near the FeI λ525.02 nm line are examined; these lines are important diagnostics for studies of small-scale dynamical processes. The Sayan Solar Observatory SMMF measurements are in good consistency with the Wilcox Solar Observatory data for 2003: for a comparison of N = 137 pairs of points in the two data sets, the correlation coefficient ρ is 0.92 for the linear regression between the datasets B_WSO = 0.03(±0.05) + 0.93(±0.03)B_SSO.     

Integration of root systems into a GIS-based slip surface model: computational experiments in a generic hillslope environment

Root systems of trees reinforce the underlying soil in hillslope environments and therefore potentially increase slope stability. So far, the influence of root systems is disregarded in Geographic Information System (GIS) models that calculate slope stability along distinct failure plane. In this study, we analyse the impact of different root system compositions and densities on slope stability conditions computed by a GIS-based slip surface model. We apply the 2.5D slip surface model r.slope.stability to 23 root system scenarios imposed on pyramidoid-shaped elements of a generic landscape. Shallow, taproot and mixed root systems are approximated by paraboloids and different stand and patch densities are considered. The slope failure probability (P_f) is derived for each raster cell of the generic landscape, considering the reinforcement through root cohesion. Average and standard deviation of P_f are analysed for each scenario. As expected, the r.slope.stability yields the highest values of P_f for the scenario without roots. In contrast, homogeneous stands with taproot or mixed root systems yield the lowest values of P_f. P_f generally decreases with increasing stand density, whereby stand density appears to exert a more pronounced influence on P_f than patch density. For patchy stands, P_f increases with a decreasing size of the tested slip surfaces. The patterns yielded by the computational experiments are largely in line with the results of previous studies. This approach provides an innovative and simple strategy to approximate the additional cohesion supplied by root systems and thereby considers various compositions of forest stands in 2.5D slip surface models. Our findings will be useful for developing strategies towards appropriately parameterising root reinforcement in real-world slope stability modelling campaigns.