Methods Development for the Constrained Elastic Modulus Investigation of Organic Material in Natural Soil Conditions

Compressibility is one of the most important mechanical properties of soil. The parameter that characterizes compressibility is the constrained modulus of elasticity. Knowledge of this is important to calculate the settlement of a structure foundation on peat material. According to soil classification by EN ISO 14688-2, peat is an organic soil that contains min. 20% organic matter. It is a highly organic type of soil. Peat material has large compressibility. The value of the constrained elasticity modulus for peat is ca. 400 kPa, while it may be ca 1.0–1.6 MPa for consolidated peat. Due to the extensive range of the modulus, experimental research in this field is proposed. It is suggested to load the peat material layer with an embankment and to determine its total settlement. Based on this, a program was developed to determine the settlement–strain relationship. The authors propose an approach according to two models: the first is based on constant stress distribution in the soil with an oedometer test. The second considers the variability of stresses in the soil and the influence of the loaded area. Both methods were tested based on numerical simulations, and then an experimental field in Szczecin was used. The formulae for the constrained modulus of elasticity measurement were derived; in practical conditions, a uniaxial deformation state can be used with the combination of the total settlement.


Introduction
Peat material is a highly organic soil. It comprises decomposed remains of plants in a wet environment [1][2][3]. Peat has a dark brown to black colour, spongy consistency, and characteristic organic odour [2]. Many behavioural parameters and properties are similar to clay soils, but they are radically different from strength characteristics. Therefore, it is taken for an immensely soft soil [3,4]. Peat has a high water content. In this soil, water can comprise 90% of the mass; thus, it is a poor substrate for buildings [1,4]. The natural moisture content can reach up to 1500%. This type of soil has low density and organic matter contents. Unit weight for this type of soil is 8 to 14.3 kN m −3 [2,5]. Huge water content (>200%) and porosity are also reasons for huge compressibility of peat (including significant secondary and tertiary compression), low stiffness and low shear strength (ca. 5-20 kPa) [2,[6][7][8]. Peat also shows a tendency to creep [7]. Cohesion for this type of soil is ca. 1 to 14 kPa, and the angle of internal friction is from ca. 3 to 24 [2,9]. Other characteristic features can include high spatial variability and permeability [2]. Due to the high level of groundwater and peat's low density, the effective tensions are relatively low [7], which shows peat has weak geotechnical properties [8].
It is crucial to identify an organic soil parameter. In comparison with mineral soil, these parameters vary. Therefore, traditional methods of soil mechanical parameters testing may not be suitable for organic material such as peat.
It is possible to choose the parameters of the compressibility of peat such that they relate to the technical parameters. However, there is an example of engineering projects to It is possible to choose the parameters of the compressibility of peat such that they relate to the technical parameters. However, there is an example of engineering projects to accept peat as a soil in roads, bridges and harbours. One such example, as described by Tan [10], is part of the US Route 44 relocation project. This route covers areas of cranberry bogs with deep layers of peat. Peat in the planned roadside was dug entirely and was covered with granular soils, using sheet piling as retaining structures. The leftover peat was the area of cooperation between the organic substrate and the structure. In addition, peat can be used as steel slag-dredged soil mixtures for construction materials at harbours [11]. For this reason, it is important to know the mechanical parameter of organic soil, such as the constrained modulus of elasticity.
Problems of soil materials are statically indeterminate due to stress-strain features, which are not linear. The behaviour of soil depends on the environment, time, pressure and temperature. Soil is different in every location. Physical and mechanical properties are evaluated based on small samples received from subsoil in nearly all cases. The structure of the soil is susceptible and can be disturbed by sampling. This may be a reason for the differences in measured laboratory tests versus in situ tests [12].
One of the main features of the settlement calculation is the relationship between stress and strain [7]. The stress-strain behaviour is complex, and depends on the composition, void ratio, and method of stress applied. In analyses, this behaviour uses formulae and concepts from the theory of elasticity, which is simplified because the soil's non-linear curve of stress-strain relation is replaced by a straight line [12]. Sometimes, in stress-strain relationships, it is convenient to use a soil modulus. This modulus will be different for each different type of soil. This parameter is not constant for soil but is a parameter that describes the behaviour of soil for a given set of stresses [12]. The constrained modulus of elasticity describes soil deformation under the influence of a given load in a uniaxial state without the possibility of lateral deformation. This type of modulus is shown in Figure 1. The constrained modulus of elasticity E is one of the major compressibility parameters. It is used in the one-dimensional analysis of consolidation [7]. This parameter is an inverse of the coefficient of volume compressibility (mv). The modulus is the parameter that describes uniaxial stress to uniaxial strain for limited compression [12]. This parameter is highly dependent on pressure level [13].
The large compressibility of peat characterizes organic soil. Therefore, deformation characteristics are significant. Final settlement under a given load may take a long time due to the mechanical properties of organic soil. Therefore, long-term consolidation processes result in considerable deformation [2]. Compared to other soils, peat compression takes much longer: the primary consolidation is swift; however, secondary compression takes a long time. Factors influencing peat compressibility are fibre content, The constrained modulus of elasticity E is one of the major compressibility parameters. It is used in the one-dimensional analysis of consolidation [7]. This parameter is an inverse of the coefficient of volume compressibility (m v ). The modulus is the parameter that describes uniaxial stress to uniaxial strain for limited compression [12]. This parameter is highly dependent on pressure level [13].
The large compressibility of peat characterizes organic soil. Therefore, deformation characteristics are significant. Final settlement under a given load may take a long time due to the mechanical properties of organic soil. Therefore, long-term consolidation processes result in considerable deformation [2]. Compared to other soils, peat compression takes much longer: the primary consolidation is swift; however, secondary compression takes a long time. Factors influencing peat compressibility are fibre content, natural water content, initial permeability, void ratio, and nature and arrangement of soil particles [2].
These parameters are difficult to maintain during an extended peat test in the laboratory. Engineers to calculate the settlement have to know the constrained modulus of elasticity, which is generally defined from the oedometer test, named oedometer modulus M 0 (kPa) [2,14]. The settlement parameters of peat can be determined from standard oedometer (one-dimensional compression) tests similar to in mineral soils [13], but there may be some differences in findings. In the oedometer test, the soil sample is placed into a metal ring. The height of the metal ring depends on the apparatus, which is generally a few centimetres. Then, the sample is placed in an oedometer, and it is loaded by force. The value of the odometer modulus depends on the type of soil. For example, mineral soil such as sand has the biggest value on the oedometer modulus, i.e., ca. from 40,000 kPa for fine sand to 20,000 kPa [15] for coarse sand. For soils such as clay and silt, the oedometer modulus can be from 10,000 to 60,000 kPa [15]. The smallest value is for organic soil. In an oedometer, compression is a significant source of strain [12]. Stress is adhibited to the soil sample along the vertical axis without the possibility of strain in the horizontal direction. In this case, the axial strain is equal to the volumetric strain, e.g., in Figure 1 [12]. It influences the value of settling the sample in the oedometer and then the value of the oedometer modulus.
Oedometric modulus M 0 can be calculated from the oedometer test with Equation (1): where M 0 -oedometer constrained modulus of elasticity in the range of loads from σ i to σ i+1 (kPa); ∆σ i -stress increment (kPa); ε-unit deformation of the sample [12]. The constrained modulus of elasticity for peat and other organic soils is low. In Mazury lake district (Poland, Central Europe), the oedometric modulus is between 498-846 kPa for a load range from 0 to 400 kPa [16]. Near Konstancin (Poland, Central Europe), it is from 1660 to 1726 kPa [16], and near Warsaw (Poland, Central Europe), it is from 100 to 770 kPa for a load range from 0 to 100 kPa [16]. In West Pomerania (Poland, Central Europe), it is from 199 to 980 kPa for a load range from 25 to 50 kPa [17,18] and from 180 to 600 kPa for a load range from 50 to 100 kPa [17,18], and from 600 to 1570 kPa for a load range from 100 to 200 kPa [17,18]. Compared with mineral soil, these values are incredibly small.
As the load increases, organic soils show a non-linear nature of the stress-strain relationship. Even for low loads, the settlements are considerable. The constrained modulus of elasticity of loaded peat tends to increase with increasing load over time. The compressibility of peat is considerable [7,19,20]. In research by the Warsaw University of Life Sciences and the Swedish Geotechnical Institute [7,19], the variability of the compressibility modulus for the stress range was given for peat as M 0 = 30σ' p , where M 0 is the oedometric modulus and σ' p is the initial pre-consolidation pressure.
Soil can also be the basis for agriculture activity, for both construction and environmental purposes [21]. In each of these cases, vertical displacement of the crop yield may be significant, followed by the formation of porosity and by settlement of the soil due to the reduction of porosity [21]. With this in mind, the authors decided that an essential element to determine the susceptibility to organic soil load is the peat-constrained modulus of elasticity.
Most of the peat in Poland (Central Europe) is near Szczecin, where the experiment was carried out. This results in the need to strengthen such a substrate, e.g., by preloading it with an overload embankment [12]. This type of soil improvement improves peat properties such as compressibility, stiffness, or permeability [2]. In that method, it is important to check the settlement of embankment for all times of consolidation. Buildings can be placed on the site after the consolidation is complete, which is why preloading of embankment is also a method to know the real value of the constrained modulus of elasticity. With these values, engineers can predict settlements closer to real settlements on these sites.
Research on the deformation of organic soils loaded with experimental embankments show large vertical and horizontal displacements in the zone of the loaded subsoil. Based on the observations, it can be seen that a large part of the horizontal deformation occurs during loading and shortly after its completion. Horizontal displacements do not play a significant role in the process of organic soil consolidation, except in the early stages of construction [7,22]. On this basis, calculations of the course of settlements and dispersion of the excess pore pressure in the soil under the embankment can be carried out per the one-dimensional consolidation theory, considering the appearance of plastic strains when loading the substrate [7,19,22].
The constrained modulus of elasticity E is poorly described in the literature, and the intention is to focus on its determination in the oedometric study described previously. However, the settlement of the embankment on organic soils has been described earlier, based on the cases of embankments on organic soils at research sites in Mielimaka, Białośliwie, and Antoniny (Poland, Central Europe) [7,19,22]. Test loads with embankments of organic soils were described by Huat, Ibrahim and his team, and Indraratna [2,23,24]. An interesting approach to the load on organic soils was demonstrated by the team of Waruwu and Rifa, who loaded the peat material with an embankment under laboratory conditions [25][26][27]. In these works, a constant constrained modulus of elasticity E was adopted, focusing on the consolidation of organic soils. The authors used a new approach to determine this peat parameter based on the presented models. Taking into account changes in vertical stresses in the soil and dimensions of the loaded area, designers can use the new approach to determine the modulus, while strengthening it. When determining the peat compressibility parameter, the authors analysed the relationship between the material and the load in conditions where there is no longer excess pore pressure and consolidation is complete.
This paper presents a method for the practical determination of the compressibility modulus of an organic soil layer in natural conditions, which will refine the calculation of settlement in engineering practice. The research was based on the total settlement of the embankment, i.e., after its consolidation.
The authors present the experimental studies performed on the embankment settlement on a natural scale. The oedometric modulus of organic soil compressibility to those obtained from research on-site was compared. The authors have created an original method for the practical description of changes in the constrained modulus of elasticity E of organic soils preloaded with an embankment, based on experimental research in the field. This method may eliminate the differences between the analytical calculations of soil settlement and the actual conditions found in the field.

Materials and Methods
The authors presented two methods to determine the constrained modulus of elasticity of peat material. In each method, the organic soil layer should be loaded with an overloading embankment, and its final settlement should be determined. The first method is based on the constant stress distribution as in the oedometer test. The second method includes the influence of vertical stresses from the entire surface of the embankment on the analysed place.
The uniform loading of the organic soil with a preloaded embankment with dimensions B × L was used to determine the constrained modulus of elasticity, shown in Figure 2. The experiment does not include the inclination of slopes because their effect by such an extensive area of the embankment is not significant [6,28,29].
In the analysis of the stress-strain relationship, the embankment was divided into smaller rectangular calculation areas. Under each of them, there is a column of organic soil (peat) loaded from the top with an embankment layer, and the bottom is limited by a non-deformed compacted medium sand (Figure 3). In the analysis of the stress-strain relationship, the embankment was divided into smaller rectangular calculation areas. Under each of them, there is a column of organic soil (peat) loaded from the top with an embankment layer, and the bottom is limited by a non-deformed compacted medium sand ( Figure 3).   Compression of the organic soil layer (peat material layer) by the embankment can be made of non-cohesive soil, e.g., sand. The body of the embankment can be built of fine, medium, coarse sands or coarse-grained compacted soils. One should strive to obtain soil moisture close to the optimal value, which will enable the achievement of the maximum volume density of the soil skeleton [6,12,30].
The layer of organic soil (peat) of known thickness HT has the properties described earlier, systematized in field and laboratory tests. The underlying layer can be, for example, sand, which is not compressible. The load applied to the peat layer is assumed to be value σ0, and the resulting stresses act on the calculation column as a vertical force. As a result, the column is compressed. In the paper, the authors assume that this compression is determined by a representative modulus for uniaxial compression [6,29].
Calculations should be performed at points resulting from the division of the embankment. The force concentrated in its centre can be determined, as shown in Figure  4. According to the principle of superposition, one analysed place will be influenced by all other calculation fields [6,12,31]. Compression of the organic soil layer (peat material layer) by the embankment can be made of non-cohesive soil, e.g., sand. The body of the embankment can be built of fine, medium, coarse sands or coarse-grained compacted soils. One should strive to obtain soil moisture close to the optimal value, which will enable the achievement of the maximum volume density of the soil skeleton [6,12,30].
The layer of organic soil (peat) of known thickness H T has the properties described earlier, systematized in field and laboratory tests. The underlying layer can be, for example, sand, which is not compressible. The load applied to the peat layer is assumed to be value σ 0 , and the resulting stresses act on the calculation column as a vertical force. As a result, the column is compressed. In the paper, the authors assume that this compression is determined by a representative modulus for uniaxial compression [6,29].
Calculations should be performed at points resulting from the division of the embankment. The force concentrated in its centre can be determined, as shown in Figure 4. According to the principle of superposition, one analysed place will be influenced by all other calculation fields [6,12,31]. When an embankment acts on the subsoil, each soil element (peat material) under the embankment at any depth is bounded by adjacent elements. These elements are subjected to a state of stress [32,33]. It was assumed that there are no horizontal deformations of the soil outside the boundaries of the loaded area. Compared to the thickness of the compressible soil under the embankments, there is no ideal onedimensional compression condition. On-site, there is some settlement caused by lateral displacement of the soil from under the load area. However, inclinometer measurements of lateral soil deformation at the embankment boundaries indicate that settlements due to lateral deformations are generally small compared to compression settlements [32]. Therefore, the general method of calculating the settlement is based on the uniaxial state of strain. Consolidation is not a point of this paper, and the value of settlement is taken after the consolidation.
The embankment layer settled most in the middle and the least at the corners. Stresses from external load in the soil depend on its value, place of its application, and the area of action [6,31]. The stress is determined based on the Boussinesq theory and the principle of superposition [6,12,31], according to Equation (2): where P -force acting on the subsoil in the analysed calculation area (kN); z -depth consideration space (m); l -distance of the force from the analysed point in the embankment determined according to Equation (3) (m): where x, y -coordinates of the force applied in the embankment plane (m). When an embankment acts on the subsoil, each soil element (peat material) under the embankment at any depth is bounded by adjacent elements. These elements are subjected to a state of stress [32,33]. It was assumed that there are no horizontal deformations of the soil outside the boundaries of the loaded area. Compared to the thickness of the compressible soil under the embankments, there is no ideal one-dimensional compression condition. On-site, there is some settlement caused by lateral displacement of the soil from under the load area. However, inclinometer measurements of lateral soil deformation at the embankment boundaries indicate that settlements due to lateral deformations are generally small compared to compression settlements [32]. Therefore, the general method of calculating the settlement is based on the uniaxial state of strain. Consolidation is not a point of this paper, and the value of settlement is taken after the consolidation.
The embankment layer settled most in the middle and the least at the corners. Stresses from external load in the soil depend on its value, place of its application, and the area of action [6,31]. The stress is determined based on the Boussinesq theory and the principle of superposition [6,12,31], according to Equation (2): where P-force acting on the subsoil in the analysed calculation area (kN); z-depth consideration space (m); l-distance of the force from the analysed point in the embankment determined according to Equation (3) (m): where x, y-coordinates of the force applied in the embankment plane (m). The diagram of the stress distribution in the centre according to the Boussinesq theory is shown in Figure 5 [6,29,31].
The diagram of the stress distribution in the centre according to the Boussinesq theory is shown in Figure 5 [6,29,31]. Stresses in the column at any depth "z" and the effect of load on the state of vertical stresses in the remaining columns can be written by Equation (4): According to the superposition principle, in the analysed column of organic soil, the stresses will consider the influence of the load exerted on the analysed area and all the others ( Figure 6) [6,29,31]. Stresses in the column at any depth "z" and the effect of load on the state of vertical stresses in the remaining columns can be written by Equation (4): According to the superposition principle, in the analysed column of organic soil, the stresses will consider the influence of the load exerted on the analysed area and all the others ( Figure 6) [6,29,31]. When analysing the selected calculation area/column, it is reasonable to u uniaxial model that shows the stresses from other calculation fields. Only vertical str in the soil will be used to perform the calculations, and the analysed settlements w complete, i.e., after the subsoil consolidation process is completed.
The constrained modulus of elasticity E can be determined by a model that is b on the constant stress distribution, as in the oedometer test (first model), and reconciling the influence of the entire load area and the thickness of the organic soil u it (second model). Based on this, a map of the constrained modulus of elasticity variab under the preload embankment can be created. The least-squares method sh approximate its value based on the values at the selected points to obtain the value o When analysing the selected calculation area/column, it is reasonable to use a uniaxial model that shows the stresses from other calculation fields. Only vertical stresses in the soil will be used to perform the calculations, and the analysed settlements will be complete, i.e., after the subsoil consolidation process is completed. The constrained modulus of elasticity E can be determined by a model that is based on the constant stress distribution, as in the oedometer test (first model), and by reconciling the influence of the entire load area and the thickness of the organic soil under it (second model). Based on this, a map of the constrained modulus of elasticity variability under the preload embankment can be created. The least-squares method should approximate its value based on the values at the selected points to obtain the value of the constrained modulus of elasticity for the entire layer. This value of the constrained modulus of elasticity E can determine peat material settlement at the design stage.
The formulae for the constrained modulus of elasticity were derived, which in practical conditions, allows the use of a uniaxial deformation state with the combination of the total settlement. The authors describe all load and deformation components and transposed them into a uniaxial state of deformation, which corresponds to total settlement. Based on this, a calculation program was developed to determine the settlement-stress (load) relationship.

The First Model to Determine the Constrained Modulus of Elasticity in Peat Material
In the first model, calculations are made at selected points based on the discretization of the embankment. When using this model, it should be assumed that the external load stresses in the peat layer are equal σ z = σ 0 = const. The influence of the load from other design areas is not taken into account.
Organic soil constrained modulus of elasticity (E m_1 ) is constant throughout the column (E m_1 = const and ∂E ∂z = 0). The assumptions for the calculations are shown in Figure 7. The calculated settlement of the peat material layer in Model 1 was defined as: where s -settlement of the organic soil layer (m); σ0 -load caused by the preload embankment (kPa); HT -thickness of the organic soil layer (m); Em_1 -constrained modulus of elasticity in the considered column determined by the first model (kPa).
Knowing the settlement in the analysed place of the embankment and using the above dependence, it is possible to determine the organic soil constrained modulus of The calculated settlement of the peat material layer in Model 1 was defined as: where s -settlement of the organic soil layer (m); σ 0 -load caused by the preload embankment (kPa); H T -thickness of the organic soil layer (m); E m_1 -constrained modulus of elasticity in the considered column determined by the first model (kPa).
Knowing the settlement in the analysed place of the embankment and using the above dependence, it is possible to determine the organic soil constrained modulus of elasticity by the first model with Equation (6).

The Second Model to Determine the Constrained Modulus of Elasticity in Peat Material
The second model contains information about the influence of load acting at any point in the load area of the column and the formation of vertical stress in the soil. The calculations are performed in points determined by the division of the embankment into calculation areas.
The Boussinesq theory and the principle of superposition are applied. The distribution in the organic soil layer is not constant: σz = f(x,y,z), ∂σ z ∂z = 0. In the second model, the influence of the remaining calculation areas is taken into account, following the principle of superposition (Figure 8b).
The compressibility modulus determined by the second model is constant throughout the peat layer: E m_2 = const and in the column: ∂E ∂z = 0. The assumptions for the calculations are outlined in Figure 8a. The second model is analysed at any point A with the coordinates (x A , y A ), according to the adopted coordinate system (e.g., Figure 6), where x A ∈ 0; B and y A ∈ 0; L . The stresses in the organic soil column will be calculated according to the Bouss theory:  The stresses in the organic soil column will be calculated according to the Boussinesq theory: where P-force acting on the subsoil in the analysed calculation area [kN]; l-distance of the force from the analysed point in the embankment determined according to Equation (8) (m): x, y -coordinates of the application of the force in the plane of the embankment (m); x A , y A -coordinates of the analysed point (m). The vertical force P on any elementary area of the embankment is determined by: where A -calculated area field (m 2 ); σ 0 -embankment load (kPa). The total calculation for settlement of the peat layer in the second model is presented as: Concerning the model presented in the earlier work [6], the authors consider the influence of embankment stresses at a given point by integrating the entire embankment area, not summing up the previously calculated stress values at each assumed point.
The stresses in the organic soil layer at a point are determined based on Equation (11) [6]: Taking into account the influence of the remaining calculation areas on point A and the change of stresses in the soil, it can be concluded that: Assuming (13) should be solved as follows: where (14) to (19): The total computational settlement of organic soil containing the relationship (10) according to the second model can be determined: Knowing the load on the embankment and the settlement of organic soil, the constrained modulus of elasticity according to the second model can be represented by Equation (21):

Results
The embankment erected on organic soils was numerically analysed. Then, the peat compressibility modulus was determined under natural conditions to verify the proprietary solution.   The embankment was divided into smaller calculation areas. Points and forces acting inside each field were determined with the use of a proprietary calculation program. The division of the embankment is shown in Figure 10. The embankment was divided into smaller calculation areas. Points and forces acting inside each field were determined with the use of a proprietary calculation program. The division of the embankment is shown in Figure 10. It was assumed that the total settlement (st) of the embankment was most significant in the middle and most minor in the corners. Its maximum value is 0.48 m and the minimum at 0.12 m. Final settlement values are shown in Table 1.  It was assumed that the total settlement (s t ) of the embankment was most significant in the middle and most minor in the corners. Its maximum value is 0.48 m and the minimum at 0.12 m. Final settlement values are shown in Table 1. For the case described above, the constrained modulus of elasticity of the peat materials at the designated points was determined, as described in the first method. The constrained moduli of elasticity of peat loaded with a 5 × 10 m embankment are included in Table 2 and are shown in Figure 11.  From the obtained results, the largest constrained modulus of elasticity Em_1 is at the corners at 2485.33 kPa, and the smallest was at the centre of the embankment at 621.33 kPa. The average is 1164.27 kPa ( Figure 11 and Table 2). The difference between the maximum and minimum value is 1864.00 kPa.

Numerical Simulations for the First and Second Method
The shape of the diagram shows the inverse variability to the settlement diagram. The lowest settlement occurred where the module is most significant, and the greatest is the lowest. This opposite is the result of the assumptions used in the first model.

Determining the Constrained Modulus of Elasticity with the Second Model in Peat Material for Numerical Simulation
The numerical analysis was performed for the assumed case. Then, the constrained modulus of elasticity of the peat material was determined using the second method. A parameter Ji (kN/m) was first established to apply this method. The values of the Ji (kN/m) parameter can be found in Figure 12 and Table 3.  From the obtained results, the largest constrained modulus of elasticity E m_1 is at the corners at 2485.33 kPa, and the smallest was at the centre of the embankment at 621.33 kPa. The average is 1164.27 kPa ( Figure 11 and Table 2). The difference between the maximum and minimum value is 1864.00 kPa.
The shape of the diagram shows the inverse variability to the settlement diagram. The lowest settlement occurred where the module is most significant, and the greatest is the lowest. This opposite is the result of the assumptions used in the first model.

Determining the Constrained Modulus of Elasticity with the Second Model in Peat Material for Numerical Simulation
The numerical analysis was performed for the assumed case. Then, the constrained modulus of elasticity of the peat material was determined using the second method. A parameter J i (kN/m) was first established to apply this method. The values of the J i (kN/m) parameter can be found in Figure 12 and Table 3.  The solution presented in Figure 12 shows that the most significant values of Ji are below the embankment's centre (Ji_max = 97.97 kN/m), and the smallest are at the corners (Ji_min = 38.82 kN/m). The distribution of Ji confirms the correctness of these findings.
Then, the constrained modulus of elasticity value was determined using the Ji parameter and the final settlement value at individual points. The results are shown in Table 4 and Figure 13.   The solution presented in Figure 12 shows that the most significant values of Ji are below the embankment's centre (J i_max = 97.97 kN/m), and the smallest are at the corners (J i_min = 38.82 kN/m). The distribution of Ji confirms the correctness of these findings.
Then, the constrained modulus of elasticity value was determined using the Ji parameter and the final settlement value at individual points. The results are shown in Table 4 and Figure 13.

Full-Scale Experiment
After checking the models during numerical simulations, a full-scale experiment was performed.
The experimental tests were carried out at the test site on the Ostrów Brdowski Island in Szczecin, at the construction site of a factory for windmills. An experimental plot with a preload embankment was used. The embankments were founded on a layer of low peat, under which there are sands of valley bottoms and floodplain terraces. The peat material thickness is 9 m [6,31,34].
The subject of the experiment was the settlement process of the preload embankment. The data on the stress-strain relationship of peat material was obtained by reading the settlement of the disc benchmarks at fixed surface points at specific measurement points. The embankment is shown in the photo (Figure 14).

Full-Scale Experiment
After checking the models during numerical simulations, a full-scale experiment was performed. The experimental tests were carried out at the test site on the Ostrów Brdowski Island in Szczecin, at the construction site of a factory for windmills. An experimental plot with a preload embankment was used. The embankments were founded on a layer of low peat, under which there are sands of valley bottoms and floodplain terraces. The peat material thickness is 9 m [6,31,34].
The subject of the experiment was the settlement process of the preload embankment. The data on the stress-strain relationship of peat material was obtained by reading the settlement of the disc benchmarks at fixed surface points at specific measurement points. The embankment is shown in the photo (Figure 14). The calculations took into account the variability of the constrained modulus of elasticity and vertical stress. The result obtained in this way reflects the features found in natural conditions. Organic soil samples were taken to determine geotechnical parameters. For comparison, the authors also performed traditional tests using an oedometer.
The embankment size is B = 63.58 m, L = 110.2 m, and it has two sections: A2, which contains measuring points from A2-P01 to A2-P16, and A2 ext, with the A2-P17 measuring point. The scheme of the embankment with the locations of the measurement points is shown in Figure 15. The calculations took into account the variability of the constrained modulus of elasticity and vertical stress. The result obtained in this way reflects the features found in natural conditions. Organic soil samples were taken to determine geotechnical parameters. For comparison, the authors also performed traditional tests using an oedometer.
The embankment size is B = 63.58 m, L = 110.2 m, and it has two sections: A2, which contains measuring points from A2-P01 to A2-P16, and A2 ext, with the A2-P17 measuring point. The scheme of the embankment with the locations of the measurement points is shown in Figure 15. The embankment loaded the subsoil. The load value was equal to σ0 = 35 kPa. The embankments were observed for 84 days. A vertical drainage system was used to complete the total settlement faster. The values of total settlement at the indicated points are summarized in Table 5. The total settlement of the embankment is shown in Figure 16. The embankment loaded the subsoil. The load value was equal to σ 0 = 35 kPa. The embankments were observed for 84 days. A vertical drainage system was used to complete the total settlement faster. The values of total settlement at the indicated points are summarized in Table 5. The total settlement of the embankment is shown in Figure 16.   In the embankment, the highest settlement value reached 1.175 m at the A2-P10 point and reached the lowest of 0.092 m at A2-P17, with a mean value of 0.410 m.
The authors used the above values to determine the constrained modulus of elasticity of the peat material layer based on its overload with an embankment by the described methods.
One of the essential elements when changing the volume under load is rheological properties. They make it possible to consider changes in the behaviour of the sample, even if they have only structural changes in a fragmented medium [35]. An example of this is self-propelled colloidal populations that induce time-dependent three-dimensional fluid flows, which change the rheological properties of the centres [35]. The authors did not consider this issue in this paper.

Determining the Constrained Modulus of Elasticity with the First Model in Peat Material for a Full-Scale Experiment
Based on the previously described method, the constrained modulus of elasticity was determined for the presented case of a full-scale experiment.
The constrained modulus results for the actual appendage in Szczecin are summarized in Table 6 and Figure 17. The authors used the above values to determine the constrained modulus of elasticity of the peat material layer based on its overload with an embankment by the described methods.
One of the essential elements when changing the volume under load is rheological properties. They make it possible to consider changes in the behaviour of the sample, even if they have only structural changes in a fragmented medium [35]. An example of this is self-propelled colloidal populations that induce time-dependent three-dimensional fluid flows, which change the rheological properties of the centres [35]. The authors did not consider this issue in this paper.

Determining the Constrained Modulus of Elasticity with the First Model in Peat Material for a Full-Scale Experiment
Based on the previously described method, the constrained modulus of elasticity was determined for the presented case of a full-scale experiment.
The constrained modulus results for the actual appendage in Szczecin are summarized in Table 6 and Figure 17.   For the full-scale experiment, the largest modulus is at the A2-P17 point and amounts to 3423.91 kPa, while the smallest is at the A2-P10 point (Em_1 = 268.09 kPa). The average of all points is 1107.02 kPa.

Determining the Constrained Modulus of Elasticity with the Second Model in Peat Material for a Full-Scale Experiment
For the case described above and following the previously adopted assumptions, the Ji parameter was calculated in the calculation points of the embankment. The result is shown in Table 7 and Figure 18.  For the full-scale experiment, the largest modulus is at the A2-P17 point and amounts to 3423.91 kPa, while the smallest is at the A2-P10 point (E m_1 = 268.09 kPa). The average of all points is 1107.02 kPa.

Determining the Constrained Modulus of Elasticity with the Second Model in Peat Material for a Full-Scale Experiment
For the case described above and following the previously adopted assumptions, the J i parameter was calculated in the calculation points of the embankment. The result is shown in Table 7 and Figure 18.   Figure 18 and Table 7 show that the highest value Ji, Ji_max = 278.34 kN/m, is located at A2-P12 below the centre of the embankment, and the smallest, Ji_min = 204.01 kN/m, in A2-P01 (Ji_min = 204.01 kN/m). The span between these results is 74.32 kN/m.
Based on the known total settlement of the embankment at the points and the value of the Ji parameter, the constrained modulus of elasticity of the peat material was determined using the second method. The results are shown in Table 8 and Figure 19.   Figure 18 and Table 7 show that the highest value J i , J i_max = 278.34 kN/m, is located at A2-P12 below the centre of the embankment, and the smallest, J i_min = 204.01 kN/m, in A2-P01 (J i_min = 204.01 kN/m). The span between these results is 74.32 kN/m.
Based on the known total settlement of the embankment at the points and the value of the J i parameter, the constrained modulus of elasticity of the peat material was determined using the second method. The results are shown in Table 8 and Figure 19.  According to Figure 19, the largest modulus, Em_2 = 2743.17 kPa, appears at the point A2-P17, the place of the lowest embankment settlement. The smallest, Em_2 = 232.16 kPa, in A2-P10, where settlement is the most intense. The mean value obtained from the calculations presented above is 877.62 kPa According to Figure 19, the largest modulus, E m_2 = 2743.17 kPa, appears at the point A2-P17, the place of the lowest embankment settlement. The smallest, E m_2 = 232.16 kPa, in A2-P10, where settlement is the most intense. The mean value obtained from the calculations presented above is 877.62 kPa

Discussion
The constrained modulus of elasticity of the peat material was determined via both methods presented for the numerical simulation and the actual case study in Szczecin.
In numerical simulations, the constrained modulus of elasticity was obtained from the first model from 621. 33  For the simulation, the span between the maximum and minimum values for the first method is 1864 kPa, and the second is 141.65 kPa. It may prove that the influence of the loaded surface is of great importance. It is the result of made assumptions. The modulus is 4.43 times greater on the mean of the first method than on the second.
An embankment with B = 63.58 m, L = 110.2 m, and with 17 measuring points was made. The settlement was measured at all these points until the process stabilized and consolidation was completed. The total settlement was the basis for selecting the peatconstrained modulus of elasticity, which will consider the triaxial state of stresses and the uniaxial state of strains. Figure   For the simulation, the span between the maximum and minimum values for the first method is 1864 kPa, and the second is 141.65 kPa. It may prove that the influence of the loaded surface is of great importance. It is the result of made assumptions. The modulus is 4.43 times greater on the mean of the first method than on the second.
An embankment with B = 63.58 m, L = 110.2 m, and with 17 measuring points was made. The settlement was measured at all these points until the process stabilized and consolidation was completed. The total settlement was the basis for selecting the peatconstrained modulus of elasticity, which will consider the triaxial state of stresses and the uniaxial state of strains. Figure 21 presents the minimum, maximum and average values for both methods and the results of the oedometric tests for the full-scale experiment on the Brdowski Island in Szczecin (Poland).
For a full-scale experimental test, the spread between the maximum and minimum value is, respectively, 3155.83 kPa, and the second method is 2546.28 kPa. The constrained modulus of elasticity is 1.26 times greater on the mean of the first method than on the second. The average modulus tested in the oedometer is bigger than both the first (29% greater) and the second (63% greater). For a full-scale experimental test, the spread between the maximum and minimum value is, respectively, 3155.83 kPa, and the second method is 2546.28 kPa. The constrained modulus of elasticity is 1.26 times greater on the mean of the first method than on the second. The average modulus tested in the oedometer is bigger than both the first (29% greater) and the second (63% greater).
By performing experimental tests and analyses, we know the dimensions of the embankment (B = 63.58 m and L = 110.2 m) and its height (hn = 1.9 m). The unit weight of the embankment material was 18.42 kN m −3 based on the research. The load on the peat was 35 kPa. The final settlement was known in the specified ones after completing the consolidation process. Subsequently, the strictness modules were determined by the methods described. On this basis, the constrained modulus of elasticity in the peat layer was optimized using the least-squares method. In engineering practice, those representative values of the constrained modulus of elasticity are determined for the entire soil layer to determine the predicted settlement of the object. For this purpose, the optimized value of the peat material constrained modulus of elasticity was determined using the least-squares method, referring to the calculated value and the measured settlement. The calculated settlement from each point was compared to the measurements. The results are shown for the first method in Figure 22, the second method in Figure 23, and the oedometer test in Figure 24. By performing experimental tests and analyses, we know the dimensions of the embankment (B = 63.58 m and L = 110.2 m) and its height (h n = 1.9 m). The unit weight of the embankment material was 18.42 kN m −3 based on the research. The load on the peat was 35 kPa. The final settlement was known in the specified ones after completing the consolidation process. Subsequently, the strictness modules were determined by the methods described. On this basis, the constrained modulus of elasticity in the peat layer was optimized using the least-squares method. In engineering practice, those representative values of the constrained modulus of elasticity are determined for the entire soil layer to determine the predicted settlement of the object. For this purpose, the optimized value of the peat material constrained modulus of elasticity was determined using the least-squares method, referring to the calculated value and the measured settlement. The calculated settlement from each point was compared to the measurements. The results are shown for the first method in Figure 22, the second method in Figure 23, and the oedometer test in Figure 24.
For the first method, the representative value is 708.49 kPa and 64% of the average value. In the second method, this value is 561.68 kPa and constitutes 64% of the average value. In an oedometer test, the representative value is 715.33 kPa and is 50% of the mean value. A summary of these values is shown in Figure 25.    For the first method, the representative value is 708.49 kPa and 64% of the average value. In the second method, this value is 561.68 kPa and constitutes 64% of the average value. In an oedometer test, the representative value is 715.33 kPa and is 50% of the mean value. A summary of these values is shown in Figure 25. The first method is based on the constant stress distribution. It was based on the oedometric model of the consolidation description relationship, assuming the effects of soil deformation and vertical filtration in the column [32]. In this method, the stress distribution is constant for the entire layer, σz = const., in the peat column. Due to the constant stress distribution, the constrained modulus of elasticity results can be compared to those obtained in the oedometer. It can be seen in Figure 25 where the value of the representative constrained modulus of elasticity for the layer of peat material considered in the first method in the oedometer test has a similar value.
The second method makes the expected results more likely. Its determinants are: The first method is based on the constant stress distribution. It was based on the oedometric model of the consolidation description relationship, assuming the effects of soil deformation and vertical filtration in the column [32]. In this method, the stress distribution is constant for the entire layer, σ z = const., in the peat column. Due to the constant stress distribution, the constrained modulus of elasticity results can be compared to those obtained in the oedometer. It can be seen in Figure 25 where the value of the representative constrained modulus of elasticity for the layer of peat material considered in the first method in the oedometer test has a similar value.
The second method makes the expected results more likely. Its determinants are: defining the compressibility modulus as a constant value over the entire thickness of the organic soil column, considering the variable vertical stress in the peat, and the influence of the load on the entire embankment. It is in line with Boussinesq' s theory and is used for the practical calculation of settlement. In the case of peat material, this is of great importance. The second method considers more details than the first.
For the optimized tightness modulus determined by Models 1 and 2, the settlement was calculated and then compared with the actual settlement of the embankment at each measured point (from A2-P01 to A2-P17). The comparison is shown in Figure 26. From Figure 26, it can be seen that there are differences between the measured settlement and that calculated from the optimized value of modulus. One of the significant differences is that the settlement calculated from the compressibility modulus of Model 1 is a constant value at each point; however, this is not so. Model 1 assumes constant stresses in the soil, which allows the value to be related to the oedometer conditions. These conditions have an impact on the calculated settlement value. Nevertheless, the most similar results were obtained in points A2-P02, A2-P05, A2-P07, A2-P14 and A2-P15. The pair of points A2-P02 and A2-P07, as well as A2-P14 and A2-P15, are located next to each other. The greatest plains are at points A2-P10 and A2-P11 under the centre of the embankment. In Model 2, the stress values in the peat layer varied, and the influence of the entire loaded area was taken into account. Adding only this effect made the calculated settlement values variable, bringing the results closer to the real case. In Model 2, the stress values in the peat layer varied, and the influence of the entire loaded area was taken into account. Adding only this effect made the calculated settlement values variable. This value brings the results closer to the real case. Similar settlement results are in points A2-P14 and A2-P15, and lesser at points A2-P02, A2-P03, A2-P04, A2-P05, A2-P06 and A2-P07. However, the lowest compliance is at points A2-P14 and A2-P15. Despite the significant difference, the measured settlement and the calculated settlement (the modulus determined by Method 2) have the highest values near points A2-P10 and A2-P11. The discrepancy between the first and the second is a consequence of the simplifications used. Those simplifications are: the external load stresses in the peat layer are equal σz=σ0=const; the influence of the load from other design areas is not taken into account (first model); the Boussinesq theory and the principle of superposition are applied; the distribution in the organic soil layer is not constant (second model); organic soil constrained modulus of elasticity (Em_1 and Em_2) is constant throughout the layer. Obtaining such high results in the first method causes the creation of a column of organic soil over the entire thickness. These differences result from considering uniaxial deformations of the organic soil layer and the triaxial stress state in the parameter based on the embankment settlement. The proprietary method allows to obtain the results of From Figure 26, it can be seen that there are differences between the measured settlement and that calculated from the optimized value of modulus. One of the significant differences is that the settlement calculated from the compressibility modulus of Model 1 is a constant value at each point; however, this is not so. Model 1 assumes constant stresses in the soil, which allows the value to be related to the oedometer conditions. These conditions have an impact on the calculated settlement value. Nevertheless, the most similar results were obtained in points A2-P02, A2-P05, A2-P07, A2-P14 and A2-P15. The pair of points A2-P02 and A2-P07, as well as A2-P14 and A2-P15, are located next to each other. The greatest plains are at points A2-P10 and A2-P11 under the centre of the embankment. In Model 2, the stress values in the peat layer varied, and the influence of the entire loaded area was taken into account. Adding only this effect made the calculated settlement values variable, bringing the results closer to the real case. In Model 2, the stress values in the peat layer varied, and the influence of the entire loaded area was taken into account. Adding only this effect made the calculated settlement values variable. This value brings the results closer to the real case. Similar settlement results are in points A2-P14 and A2-P15, and lesser at points A2-P02, A2-P03, A2-P04, A2-P05, A2-P06 and A2-P07. However, the lowest compliance is at points A2-P14 and A2-P15. Despite the significant difference, the measured settlement and the calculated settlement (the modulus determined by Method 2) have the highest values near points A2-P10 and A2-P11. The discrepancy between the first and the second is a consequence of the simplifications used. Those simplifications are: the external load stresses in the peat layer are equal σ z =σ 0 =const; the influence of the load from other design areas is not taken into account (first model); the Boussinesq theory and the principle of superposition are applied; the distribution in the organic soil layer is not constant (second model); organic soil constrained modulus of elasticity (E m_1 and E m_2 ) is constant throughout the layer. Obtaining such high results in the first method causes the creation of a column of organic soil over the entire thickness. These differences result from considering uniaxial deformations of the organic soil layer and the triaxial stress state in the parameter based on the embankment settlement. The proprietary method allows to obtain the results of calculations that are consistent with those observed on site.
Neglecting the effect of horizontal pressures is based on laboratory experiments, which were performed in an oedometer. The research procedure in the oedometer is adapted to involve the horizontal pressure when comparing the results of the proposed model with that of which was performed in the oedometer, where horizontal pressure is not included. This justifies the assumption that Model 1 gives results of E values close to that which derives from the oedometer. This was the reason that the authors, as a first step, consider the module while neglecting the effect of the horizontal pressure.
The lower the constrained modulus of elasticity, the greater the settlement. Determining the appropriate modulus values affects the designed amount of settlement. Overestimated values of the module may result in underestimating the amount of settlement, which might be more significant. It can be seen that when comparing the full-scale experimental studies to the oedometer studies, the oedometric modulus is greater than that of the modulus determined by both methods.

Conclusions
Organic soil, as in peat, is a complex medium. Therefore, the methods of determining its parameters were examined to find which best reflect the soil' s natural state, susceptibility, and variables. Full-scale experimental research was conducted in the field (in practice) and simulations were performed to prove that the proposed models can be used in the practical calculations of peat-constrained modulus of elasticity.
The authors proposed a method to optimize the constrained modulus of elasticity of peat material by presenting two methods. The first is based on a constant stress distribution, such as in the oedometer test, and the second considers the influence of the loaded area. Both methods were verified numerically based on field experiments on a natural scale and then used to determine the constrained modulus of elasticity of peat material for a real case. The modulus for both methods was determined for peat from the Szczecin area (Poland) based on the described method.
The methods recommended for the settlement calculation are based on the value of the constrained modulus of elasticity. These methods are based on the assumption of uniaxial settlement. The experimental methods use the uniaxial compression method, which considers the 3-axis stress state. When determining the constrained modulus of elasticity, the presented method takes this fact into account. Using this allows engineers to determine the settlements that occur on site.
The field-determined constrained modulus of elasticity is lower than the oedometric modulus coming from laboratory tests, which may indicate the subsequent oversizing of the settlement when using the oedometer-based method. If the constrained modulus of elasticity is lower, the settlement is more significant.
Based on the comparison of the measured and calculated settlement values, it can be seen that there is a need to expand Model 2 with additional variables. It is currently the subject of the authors' research on the issue of the peat compressibility modulus.

Data Availability Statement:
The data used for the research are included in the report [34]. The company ST3 Offshore consented to the use of the report for research purposes.