A Conceptual Model and Evaluating Experiments for Studying the Effect of Soil Deformation on Its Permeability

Soil structure and void ratio are the major factors that control the permeability changes during soil deformation. In this research, we proposed and tested a conceptual model considering these two factors based on the concept of permeability anisotropy. This model, which is expressed as k ( e ) graph, determines the total k values that soil can achieve and shows that as deformation proceeds, soil permeability passes through a speciﬁc zone in the k ( e ) graph. Thus, by deforming a soil sample, measuring its permeability during deformation, and comparing the results using the k ( e ) graph, it might be possible to predict deformation eﬀects on the permeability. To evaluate this conceptual model, we designed and built a special apparatus to carry out two sets of experiments. The ﬁrst set was performed to achieve the k ( e ) graph during static compression based on the conceptual model; and the second set was conducted to investigate the permeability changes relative to k ( e ) graph during simple shear deformation in constant volume condition. Our results show that the theoretical k ( e ) graph agrees more with the measured k ( e ) graph in medium to dense samples that might have no macropore. In addition, particles’ preferential orientation and/or anisotropic permeability were not changed during shear deformation due to three possible causes: deformation done in constant volume deformation, relatively low shear strain, and shearing along particle orientation. Void ratio and particle orientation are associated with each other, and soil shearing with constant void ratio might cause the anisotropy of permeability to be relatively constant. Thus, it is needed to design and build a new complex apparatus or use a special method for testing how permeability changes within the k ( e ) graph zone during soil deformation. structure changes on the soil permeability, we can investigate the eﬀect of soil deformation on the permeability.


Introduction
When a soil body deforms, its permeability might change. Settlement, compression, consolidation, and shearing are the common causes of deformations in soils. However, less extensive deformations such as bending, torsion, and other complicated forms may also occur. Although in some of these deformations, permeability changes are negligible; they would play a key role when fluid infiltration and seepage is a matter of concern.
us, having a good knowledge of deformation effects could be a great assistance in realistic judgment about permeability variation for designing or foresight events.
In unsaturated conditions, when soil deforms due to changes in void ratio, the degree of saturation varies leading to consequent changes in its permeability [1][2][3]. Green et al. [4] reviewed the effects of various types of soil deformation on hydraulic conductivity in order to study the agricultural management practices such as tillage and reconsolidation and wheel-track soil compaction. Yeo et al. [5], Berilgen et al. [6], Fan et al. [7], and Royal et al. [8] conducted some experiments to investigate the permeability changes of slurry materials due to loading or compression occurred in the cutoff walls, cut-off barriers, and storing contaminant soils. Bryant et al. [9] investigated the permeability changes due to the consolidation of fine-grained marine sediments in various depths and introduced empirical relations to estimate the permeability from the void ratio. Morrow et al. [10], Arch and Maltman [11], Zhang and Cox [12], Crawford et al. [13], and Kimura et al. [14,15] investigated the effect of shearing on the permeability in order to find out permeability behavior of faults and gouge materials.
e review of previous works shows that researchers often have tried to investigate the effect of one specific deformation on the permeability, and generally, there is no comprehensive study on deformation effects. Most researchers had studied the effect of soil compression or consolidation on the permeability through studying void ratio changes on the hydraulic conductivity. is group of researchers often studied one-dimensional compression type of deformation by using the oedometer apparatus or similar devices that have rigid wall cells [5][6][7][8][9][16][17][18][19][20][21][22][23][24][25][26]. In addition, some researchers studied axial compression with lateral free walls using either a triaxial compression device or compression cell with flexible walls [1,3,[27][28][29]. Another type of study focuses on soil permeability changes during exerted shear strain. e main methods employed in this regard are ring shear [14,18,30], rotary shear [31], and shear of a thin layer sample [10,12,13,32].
is study is an attempt to investigate the effect of any deformation on the permeability. In this way, we described the anisotropic permeability changes during soil static compression as a conceptual model that could help to understand how any soil deformation affects the permeability. e validity of the conceptual model was investigated through sets of laboratory tests by applying an apparatus designed specifically for this research.

Parameters Controlling the k during Soil Deformation.
Soil permeability is controlled by several parameters such as grain size, specific surface area, clay content, porosity, temperature, and pore geometry [33][34][35][36][37][38][39][40][41][42][43][44]. When a soil deforms, considering laboratory constant conditions, among these effective parameters, the only varying parameters are void ratio and soil structure (including particle arrangement, orientation, and pore geometry). In this view, soil deformations can be divided into two general groups. e first group involves deformations that occur with the changing volume and both void ratio and soil structure vary. Another group comprises deformations developed with fixed volume suggesting that the void ratio is constant and only soil structure varies. Accordingly, by knowing the influence of void ratio variations and soil structure changes on the soil permeability, we can investigate the effect of soil deformation on the permeability.

Effect of Void Ratio and Soil
Structure on the k. Void ratio plays a major role in permeability such that a decrease in this parameter causes a rapid declination in permeability values. Many researchers have studied the effect of void ratio variation on the permeability and some of them proposed empirical relationships to estimate permeability via void ratio [22,42,[45][46][47][48][49].
In order to understand the effect of soil structure changes on the permeability during deformation, the soil should be studied in constant volume deformation condition so that the effect of void ratio could be negligible. At a certain void ratio, permeability in a specific direction could have different values depending on soil structure. In other words, in a specific void ratio, the value of k depends on soil particles orientation. In a soil body, if particles have random orientation, the soil would be structurally isotropic and value of permeability in all directions is equal. However, in the same void ratio, when particles are oriented, the soil would be structurally anisotropic, leading to different permeability values in different directions. In this case, permeability in the direction parallel to soil particles orientation is more than that in the direction perpendicular to them. For these conditions, a researcher reported permeability or hydraulic anisotropy in soil, represented by the anisotropy ratio (r k ) [11, 16-18, 20, 21, 31, 32, 50, 51], which is the ratio of horizontal to vertical permeability (k h /k v ) or ratio of larger to smaller permeability value as follows:

k h(max) and k v(min)
Concept. For a given soil with the specific void ratio, if the preferential orientation is at the maximum level, the horizontal permeability (k h ) would be the maximum at this void ratio (k h � k h(max) ) while the vertical permeability (k v ) would be the minimum (k v � k v(min) ). In this situation, the anisotropy ratio is at the maximum state possible (r k(max) ) and permeability values in other directions would be between k h(max) and k v(min) .

e Conceptual Model.
Basak [52], Delage and Lefebvre [53], Al-Tabbaa and Wood [16], Chapuis et al. [54], Brown and Moore [32], and Cetin [55] showed that low compression soils have less anisotropic structures and increasing compression causes the development of the anisotropic structure in the soil. Sloane and Kell [56] and Chapuis and Gill [17] reported that static compression creates a more anisotropic structure in fine-grained and sandy soils than dynamic or other compression methods. As a result, when a saturated soil sample with a high void ratio is one-dimensionally compressed gradually, soil particles can be displaced and rotated perpendicular to the applying load axis. Since particle rotation and structural anisotropy development continues from the beginning of soil compression, the r k is expected to be maximum at each void ratio. us, it is acceptable that the measured k h and k v at each void ratio during compression are the k h(max) and k v(min) . Plotting the obtained k h(max) and k v(min) in different void ratios during one-dimensionally static compression in a log-linear scale of k(e) graph would give two k h(max) and k v(min) curves, respectively ( Figure 1). Although, the soil may be isotropic (r k ≈ 1) at the beginning of compression, with a decrease in void ratio and increase in structural anisotropy, differences of these two curves increase, leading to the increase in anisotropic permeability (i.e. r k > 1). A similar behavior can be seen in the results reported by Basak [52], Wilkinson and Shipley [25], and Al-Tabbaa and Wood [16]. e k(e) graph shown in Figure 1 represents a conceptual model according to which the soils composed of particles with unequal dimensions, produce anisotropic permeability during static compression. is model is valid until there is no fracture or canal in the soil. e k(e) graph is unique for each soil sample because particle shape and geometry in different soils vary and lead to different structural anisotropies. e inferable point from the k(e) graph ( Figure 1) is that there is a range of k value for each void ratio, where k varies between k h(max) and k v(min) depends on the soil structure as follows: In fact, when a soil body with having k h(max) and k v(min) (i.e., r k(max) ) deformed, it is expected to cause the soil structure changed; consequently, anisotropy structure decreased from the initial state caused k varied between k h(max) and k v(min) . erefore, when a soil body deforms with a constant volume, it is expected that k varies only between k h(max) and k v(min) values; however, when it deforms with a changing volume, it is needed to estimate the k in relation to its corresponding void ratio. e dashed area between the two curves in Figure 1 presents all possible values of k that soil can achieve under normal conditions. As can be seen, k is not located out of this region except in special conditions, e.g., when soil grading or particle shape changes due to some process like grain crash as a result of excessive loading or water movement along produced fractures.
us, any deformation that could occur in the soil, in either a constant or a changing volume, would lead to changes in the permeability within the dashed area in Figure 1. Accordingly, the pattern of k h and k v variation during soil deformation with respect to k h(max) and k v(min) curves in the k(e) graph can show the effect of soil deformation on its permeability.
Although, the conceptual model maybe theoretically valid for all soil types, in practice, it may not work properly for cohesionless granular soils, because particles in these soils cannot rotate freely under static compression. After a slight compression, they interlock and resist against volumetric compaction [57].
is behavior prevents the development of anisotropic structures in these soil types.

Methodology and the Procedure
We tried to study the permeability changes during soil deformation via laboratory tests based on the conceptual model ( Figure 1). For this purpose, a procedure with the following steps was defined as follows: (1) Preparing the k(e) graph (k h(max) and k v(min) curves) of soil samples (2) Determining the k h and k v variations of soil samples during progressive simple shear deformation (3) Comparison of obtained results from steps 1 and 2 in order to determine how the soil permeability changed during deformation within k(e) graph zone Step 1 includes experiments conducted to obtain k h(max) and k v(min) values in different void ratios during one-dimensional static compression in order to prepare the k(e) graph.
In step 2, tests were carried out to measure k h and k v variations during simple shear deformation. rough this deformation type, it is possible to deform the soil sample without changing its volume. As a result, we can investigate only the effect of soil structure changes on the permeability without the effect of void ratio variation.
Finally, in step 3, the results of the second step of experiments were plotted on the k h(max) and k v(min) curves were obtained from the first step of experiments to display how the permeability changed with respect to its initial state during simple shear deformation. So, it can be determined that whether measured k h and k v changed between k h(max) and k v(min) curves according to the conceptual model or not. e mechanism of simple shear deformation is schematically illustrated in Figure 2. In this deformation type, theoretically, material shear and deformation occur without volume change. Here, shear strain (c) is defined as the ratio of shear displacement (D) to the height of the sample or thickness of shear (H) as follows:

Samples Specifications.
To examine the conceptual model, we need a remolded homogeneous sample in order to eliminate the effect of cracks and fractures on the permeability. us, two fine-grained soils with no organic matter were selected from loess soil deposits which are widespread in Golestan province, Northern Iran (Table 1). Rezaee et al. [58] showed that the sphericity of the loess particles is between 0.66 and 0.8 with an average of 0.73.

Decreasing void ratio
Decreasing permeability

Test Apparatus.
In order to implement the two sets of tests (for steps 1 & 2), it was necessary to design and built a special test apparatus. is apparatus must have two main capabilities; first, it should be able to apply an incremental static force to consolidate one-dimensionally the soil sample from a high void ratio to measure k h(max) and k v(min) values at the end of each compression step. Second, it should be able to apply simple shear strain to the soil sample and measure the permeability at the end of each simple shear strain. e schematic of the apparatus is illustrated in Figure 3.
is apparatus can measure k h and k v values separately through both falling and constant head methods. e cell has a square cross-section with a dimension of 8 cm and a height of 20 cm. Here, k h can be measured with two components assigned as k h1 and k h2 . e shear zone thickness in the sample during shear deformation is 5 cm; so, the dimensions of the sample deformed during the simple shear deformation are 8 cm × 8 cm × 5 cm.
e cell is composed of three parts ( Figure 3): the upper part is always fixed during compression or shear deformation, the lower part is fixed during the compression test. But, during shear deformation, it can move on a precision rail to the left or right with a fixed strain rate of 0.05 mm.sec −1 . e middle part of the cell is a section at which the test sample bears shear deformation during the displacement of the lower part (Figures 3(a) and 3(b)).
Shear deformation in this paper was reported as a strain percentage such that each 1 mm of lateral displacement of the lower part of the cell on the rail is equal to a 2% shear strain.

Samples Preparation and Test Method.
Samples with a high void ratio can be prepared by two slurry and powder methods. Since performing desired tests on the slurry sample is difficult and involves high mistakes, we decided to prepare the samples using the powder method. For this purpose, soil samples were dried and crushed and then passed through sieve no. 10 (2 mm) such that to obtain a homogeneous soil powder. e prepared powder was put in the cell using a long pipe funnel placed at a constant 2 cm height. Afterward, when the loading piston is located in the cell, the cap would be fitted, closed, and sealed.
Before the beginning of the test, the soil is saturated from the lower to upper openings through the cell valves (Figure 3); as a result, air exits from the upper cell valve. After saturating and deairing the sample, it is ready for the test.
Horizontal permeability in all tests measured with two components: k h1 and k h2 . In permeability experiments, during shear deformation, k h1 is the permeability along the shear direction of the sample; thus, k h2 is perpendicular to the shear direction.
For the first set of tests (step 1), after each increment stage of the static load and soil compression, the sample is allowed to reach a new equilibrium condition until displacements are zero. Next, horizontal and vertical permeability values are measured separately.  For the second set of tests (step 2), the sample compressed until reaching the desired void ratio. en, to prevent the soil volume changes during the shear deformation, especially in dense samples, their vertical load was released and set to about 1 kg/cm 2 . However, the applied load was not completely released because removing the load can cause detachment between the sample body and cell wall or produce opening fractures in the sample. After compressing the soil to the desired void ratio, setting the vertical load, and reaching an equilibrium state, the k h and k v values were measured. en, soil sheared by 1 mm shear displacement of the lower part of the cell, which is equivalent to 2% shear strain. After this shearing and reaching an equilibrium state, k h and k v values were measured. is process comprises shearing, reaching an equilibrium state, and measuring k h and k v was repeated until the cell sealing is not damaged by shear displacement.
It has to be noted that all horizontal and vertical permeability were measured after water flow reached a steady state. Measurements in each step were replicated at least 5 times, and the final value was reported as their average.

First
Step Results: e k(e) Graph of Samples. Variations of k h(max) and k v(min) corresponding to the void ratio in each step of progressive compression for Ι and ΙΙ samples (the k(e) graphs) are illustrated in Figures 4 and 5, respectively.
Results show that permeability variations of samples Ι and ΙΙ during the static compression were similar, in general.
ere is a linear relationship between the void ratio and the logarithm of k v(min) , similar to previously published data [9,16,20,26,42,[59][60][61][62]. However, the decrease in k h(max) differs from k v(min) and is not linear. As illustrated in Figure 6, the anisotropic permeability ratio (r k ) of samples has a varying trend during the compression, despite what is expected from the conceptual model. Here, r k is expressed as k h(av) /k v which k h(av) denotes the average of k h1 and k h2 .
As presented in Figure 6, the sample's r k increases rapidly from the beginning of the compression and reaches a maximum level. en, the increase in compression resulted in a rapid decrease in r k until reaching near a minimum level. After this drop, at medium void ratios, r k almost has minimum values with relatively low variations during the compression but it increases again at low void ratios.
Soil structure observations during the compression show that r k changes were affected by the variation in the geometry of large pores and particle orientation (Figure 7). Macrostructure of the soil was examined using the images captured by a camera installed on the cell wall ( Figure 3). Because of the sample preparation method, the samples have a granular structure at the high void ratios and mostly are composed of large units such as ped grains and silt-size grains [63] as well as other fine units such as clays and clay size particles (Figures 7(a) and 7(b)). Additionally, pores with various sizes are seen in the samples; macropores distribute between large units while micropores distribute within these large and fine units [53,63,64].
When the soil has a high void ratio, water flows mostly through the macropores. In this regard, permeability is controlled by these pores rather than micropores that are influenced by the capillary force [22,65,66]. But at low void ratios, macropores are absent and hydraulic conductivity is controlled by micropores [67].

Shock and Vibration
Macropores are more sensitive to compaction compared to micropores [23,68]. During soil compaction or consolidation, larger pores close first [53,[69][70][71][72][73]. is evolution of soil structure affects the permeability and r k variations during soil static compression. Following these variations, generally, three stages can be recognized (Figures 4-6) as follows: Stage A. A high anisotropic permeability in low compressions due to macropores effects (pseudoanisotropic permeability) Stage B. e lowest amount of anisotropic permeability in medium compressions due to the elimination of macropores affects and undeveloped orientation of the soil particles Stage C. Gradual increase of anisotropic permeability at high compressions due to developing particle orientation In Stage A, the soil has a granular structure and water flow is controlled by macropores (Figure 7(b)). Anisotropic permeability comes from flattening the ped grains and interped macropores perpendicular to the applied load axis. Although fine particles are resorted within larger units (i.e., ped grains), they do not show a clear preferential orientation. Considering the importance of macropores in control of water flow and permeability, there are maximum anisotropy ratios in this stage, and therefore, there is pseudoanisotropy in this stage [63].
At the beginning of compression in this stage, r k is relatively low but later with a small compression, ped grains, and interped macropores started to deform, and due to the affection of interconnected interped macropores, r k reached a maximum level. In the following, as compression proceeded, peds started to merge and macropore flow paths were disconnected. is process caused permeability and r k to decline rapidly.
In B stage, ped grains disappeared so there are no large pores in the soil body but few remnant macropores may exist as trapped spaces that started to complete the closure (Figure 7(c)). So decrease in the void ratio during compression in this stage is expected to occur mostly through decreasing the micropores size and particle rotation. ere is the least anisotropic permeability in this stage because of the absence of macropores and the lack of evolution in particle orientation. Rotation of particles still is not enough to have a noticeable effect on the increasing r k .
In Stage C, soil completely is compressed and there is no macropore in the soil body (Figure 7(d)). e decreasing void ratio continues by a decrease in micropore size and particle rotation. e r k increases gently due to the development of particle orientation during the compression.
Soil structure evolution through these three stages indicates that, in stages A and B, soil particles do not have a preferential orientation and it is only in C stage that a preferential orientation is developed. So only in the C stage, a   Shock and Vibration variation of k h(max) and k v(min) during compression coincide with the conceptual model. is result shows that the sample preparation method affects the k h(max) and k v(min) curves in the k(e) graph at relatively high void ratios, which are marked by A and B stages. Preparing the sample by the powder method, which was employed in this research, caused the creation of large pores in the soil body and prevented the development of particles preferential orientation until these pores were closed completely. is behavior coincides with the results of Yong and Warkentin [63].

Second
Step Results: Variations of k during Progressive Simple Shear Deformation. Experiments of the permeability changes during simple shear deformation conducted for all defined three stages A, B, and C, so that to determine variations of k h and k v related to k h(max) and k v(min) curves in the k(e) graph. e measured k v , k h1 , and k h2, after each step of 2% incrementing simple shear strain for samples Ι and ΙΙ are presented in Figures 8 and 9, respectively. Specifications of samples at the beginning of the tests are presented in Table 2.
In stages A and B, sealing of cell was failed at low shear strain (due to the low compression of the soil sample) thus the experiment was terminated and the strain developed in the sample did not exceed 10%.

ird Step Results: Comparison of the Results from Steps 1 and 2.
e measured k values in step 2 (Figures 8 and 9) are plotted on the related k(e) graph and shown in Figures 10  and 11. e pattern of points scattering in these figures shows that except in some cases, other plotted data generally are consistent with k h(max) and k v(min) curves. Examples of these contradictory behaviors are seen in data such as k h2 in ΙΙ sample in Stage A or k h1 and k h2 in sample Ι in Stage A, which are failed data or arisen due to the closure effect of large pores and conduits (summary of the comparison results are presented in Table 2).
is result is unlike the behavior expected from the conceptual model that had expressed structure changes caused permeability changed between k h(max) and k v(min) curves. For example, in Stage C of samples in which soil particles orientation is developed, even slight increases in     void ratios due to swelling caused by shearing led to a scattering data trend similar to k h(max) and k v(min) curves ( Figures 10 and 11). is result proves that permeability mainly is affected by void ratio changes and structural changes have a negligible effect on it. is result is in line with the findings of other researchers [32,74].
Lack of consistency between the second set of experiments results and conceptual model predictions arises from the fact that in constant volume condition; the induced simple shear deformation cannot produce significant changes in particle orientation. Soil macrostructure observations during shear deformation confirm this point in all shear deformation tests. Figure 12 illustrates the angularity changes of a small fragment of a gastropod shell in sample Ι during the static compression. ese remnant shell fragments could be found in Golestan province loess deposits. It was noted that fragment angle related to the horizontal plane (perpendicular to applied load or parallel to shear direction) decreased by 23°w ith a decrease in a void ratio from 0.95 to 0.68. Also, a simple shear strain was applied to this sample at a void ratio of 0.8 (stage A of sample Ι in Table 2). e structure observed for this sample revealed that, during the shear deformation up to 8% strain, the angle of the shell did not change. us, it is obvious that the effect of simple shear deformation on particle orientation is negligible compared to the effect of compression; so the shear deformation under constant volume conditions cannot change the initial soil particles orientation i.e., soil structural anisotropy. For this reason, permeability variations during the shear deformation have the most dependency to the void ratio. As a result, plotted data in Figures 10 and 11 are consistent with the k h(max) and k v(min) curves and show a similar trend. Another two effective factors might be mentioned that could cause the anisotropic permeability relatively constant during soil shear deformation. First, exerted shear strain was inadequate to develope anisotropic permeability with comparison to severe shearing in other phenomena such as fault zones and shear plane of landslides . Secondly, exerted shear strain took place tangential to particle orientation.
is caused the particles to slide against each other with no change in their orientations or soil volume.
So, it is necessary to change the void ratio during shear deformation to rotate the soil particles and increasing anisotropic permeability. According to Arch and Maltman [11], shear deformation with the decreased void ratio produces the highest amount of orientation compared with compression or shear deformation alone. It is reported in some previous studies that the shear deformation produces much anisotropy rather than compression [11,13,18,25,32,75] but according to findings of this research, it reveals that shear deformation and compression occurred simultaneously.

Summary and Conclusions
e results of the first set of experiments revealed that when the soil is statically compressed from a high to low void ratio, generally three stages of anisotropic permeability variations could be recognized on the k(e) graph based on the soil structure evolution. Toward the end of stage B at medium to low void ratio, the variations of k h(max) and k v(min) , as well as r k , has a good consistency with the conceptual model. is is most likely due to the absence of large pore spaces in the soil body as well as soil particles orientation perpendicular to the applied load. In samples with a high void ratio (both samples in stage A), the macropores play a main role in affecting the permeability, and thus, affect the k(e) graph. So, for these samples, macropore properties, due to the sample preparation method, seem to have the main effect on the k(e) graph. It means that the k(e) graph varied with the sample preparation method for each soil type.
Results of the second set of tests show that measured k values during simple shear deformation coincide with k h(max) and k v(min) curves, suggesting that permeability is mainly affected by void ratio and permeability anisotropy not changed significantly; so the measured k h and k v are not located within the area between k h(max) and k v(min) curves. It seems that the lack of changes in the anisotropic permeability arises from the three points: first, deformation was created under constant volume conditions, second, shear strain is relatively low, and third, shear strain exerted along particle orientations. Deformation in fixed volume conditions may produce few changes in the soil structure but cannot significantly change the orientation of soil particles; because the creation of particles with preferential orientation eliminates spaces and could lead to volume reduction. On the contrary, the reduction in particles' preferential orientation caused the volume increase. erefore, the key point is that how the soil deformation might change the void ratio and soil particle preferential orientation. e constant void ratio indicates that the soil particles have no preferred orientation changes. In this case, anisotropic permeability has not changed and the values of k h and k v follow the k h(max) and k v(min) curves. However, if the soil shows any void ratio changes, it should be due to the changes in anisotropic permeability that may cause the values of k h and k v not to follow the k h(max) and k v(min) curves. e results show that measured permeability's (k h1 , k h2 and k v ) have good coincidence with the k h(max) and k v (min) curves that validate the conceptual model. e preparation method of samples affects the k h(max) and k v(min) curves at very high void ratios (approximately larger than 1), but at lower void ratios, its effects on these curves are negligible. Both conceptual model and validating tests show that, in any soil deformation, the soil particles orientation or structural anisotropy do not change in a fixed void ratio. So the void ratios have a significant role in the permeability for two reasons in soil deformation; the void ratio changes have a larger effect on the permeability relative to the soil structural anisotropy changes. Also, the soil particles' orientation and anisotropy do not change while the void ratio is relatively constant.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.