Estimating the Unit Weight of Local Organic Soils from Laboratory Tests Using Artificial Neural Networks

The estimation of the unit weight of soil is carried out using laboratory methods; however, it requires high-quality research material in the form of samples with undisturbed structures, the acquisition of which, especially in the case of organic soils, is extremely difficult, time-consuming and expensive. This paper presents a proposal to use artificial neural networks to estimate the unit weight of local organic soils as leading parameters in the process of checking the load capacity of subsoil, under a direct foundation in drained conditions, in accordance with current standards guidelines. The initial recognition of the subsoil, and the locating of organic soils at the Theological and Pastoral Institute in Rzeszow, was carried out using a mechanical cone penetration test (CPTM), using various interpretation criteria, and then, material for laboratory tests was obtained. The analysis of the usefulness of the artificial intelligence method, in this case, was based on data from laboratory tests. Standard multi-layer backpropagation networks were used to predict the soil unit weight based on two leading variables: the organic content LOIT and the natural water content w. The applied neural model provided reliable prediction results, comparable to the standard regression methods.


Introduction
Soil unit weight is a geotechnical parameter that is necessary to determine, among other things, the design bearing resistance of soils in drained conditions and the vertical component of geostatic stress. Nowadays, the basic and most reliable method to estimate this parameter is by direct tests conducted in laboratory conditions. However, these tests always require the highest quality samples with undisturbed structures, collection of which is very difficult, time-consuming and expensive. Alternative methods to determine the unit weight of soil in a direct way, in situ, using various types of penetrometers, mainly based on cone penetration test (CPTM) static probes, have been sought for a long time [1][2][3][4][5][6]. Despite the widespread use of various types of penetrometers in engineering practice and the development of appropriate calculation formulas, their results have not always been strictly verified by laboratory methods; therefore, they cannot be treated as universal in relation to all types of soil, especially organic soil. Preliminary analyses in this area have shown that, at present, in relation to selected, locally-occurring, Polish organic soils, existing global relationships are of little use, and it is advisable to search for other, more accurate solutions [7].
Alternatively, from a practical point of view, these methods seem to be designed to determine soil unit weight based on the research of samples with a disturbed structure; however, they are designed for soil with natural moisture, which is cheaper and relatively easier to collect. This is particularly

Natural Water Content
In the case of organic soils, the natural water content may be extremely high and reach over 100%. The water content of the soil had the ratio of the mass of free water to the mass of dry soil, and was calculated according to Formula 1 [39]: where m w is the mass of water, m d is the mass of the dried test specimen, m 1 is the mass of the container and moist test specimen, m 2 is the mass of the container and dried test specimen and m c is the mass of the container.

Organic Content
Organic content in soils is most often determined using a mass loss on ignition test (LOI T ). The loss on ignition is normally determined on a representative sample of the soil, finer than 2.0 mm, as the mass lost by ignition on a prepared specimen at a specified temperature. The organic content is calculated on the assumption that the organic mass is totally burned by the ignition and that the mass loss is only due to the ignition of the organic matter. Unfortunately, the standard guidelines [40,41] do not specify detailed parameters of the assay, but only set very wide ranges of temperatures (400-900 • C) and roasting times (3-9 hours). This issue was analyzed in the example of local organic soils by Marut and Straż [42], and the effect was the statement that, in the considered cases, roasting for a shorter time at a given temperature gives a comparable effect to longer use time periods. The difference in this case did not exceed 1%; therefore, it was assumed that, due to similar, local origins, the samples of organic soils will be ignited at a temperature of 800 • C for 3 hours. The loss on ignition of the soil was calculated according to Formula 2 [43]: where DW 105 is the weight after oven-drying the soil to a constant weight at ca. 105 • C and DW T is the dry weight of the sample after heating to the ignition temperature (at ca. 800 • C).

Soil Unit Weight
The bulk density of the soil is the mass of soil per unit volume of the material, including any water or gas it contains. The term unit weight, γ t , is often used and is calculated by multiplying the bulk density by the acceleration due to gravity [44]: where ρ is the bulk density, m is the specimen mass, V is the volume of the specimen and g is the acceleration due to gravity.

Characteristics of the Study Area at the Rzeszow Site
The study area and data used in this paper show the effect of soil exploration on the premises where an educational building of the Theological and Pastoral University in Rzeszow is located. The site where the diagnosis was conducted, in terms of morphology, is located in the valley of the Młynówka River and reaches around 206.0 m above sea level. The area of the conducted investigation is geographically located in the Rzeszow Foothills, within the macro-region of the Sandomierz Basin, whereas geologically speaking, it is located in the south part of the Carpathian Foredeep [7]. In the study area, three locations were selected in the immediate vicinity of the CPTM soundings, where the subsoil is primarily comprised of organic soils, at a depth of about 8 m below ground level. The boreholes were drilled in the vicinity of the selected CPTM sounding points, and 135 samples with disturbed and undisturbed structures of organic soils were taken for laboratory tests. Macroscopic analysis and laboratory tests confirmed the test results with cone penetration tests, and an example of the local subsoil profile is shown in Figure 1 [45].

Artificial Neural Network
One of the most commonly used tools to describe dependencies in geotechnics was standard regression. Recently, non-standard methods, including an artificial neural network and a fuzzy logic regression, have also been increasingly used as approximation tools. Moreover, the most common neural networks [46] used for function approximation (nonlinear regression) are multi-layer perceptron (MLP), radial basis function (RBF) networks and support vector machines (SVM). The feedforward networks consist of a series of layers ( Figure 2). The first layer has a connection from the network input. Each subsequent hidden layer has a connection from the previous layer. The final output layer produces the network's output. The training network process requires a set of examples of proper network adjust (relations of parameter inputs and target outputs). The process of training a neural network involves changing the values of the weights and biases of the network to optimize network performance, as defined by the network performance function. The default performance function for feedforward networks is mean square error.
For training multilayer feedforward networks, any standard numerical optimization algorithm can be used to optimize the performance function, but there are a few that have shown excellent

Artificial Neural Network
One of the most commonly used tools to describe dependencies in geotechnics was standard regression. Recently, non-standard methods, including an artificial neural network and a fuzzy logic regression, have also been increasingly used as approximation tools. Moreover, the most common neural networks [46] used for function approximation (nonlinear regression) are multi-layer perceptron (MLP), radial basis function (RBF) networks and support vector machines (SVM). The feedforward networks consist of a series of layers ( Figure 2). The first layer has a connection from the network input. Each subsequent hidden layer has a connection from the previous layer. The final output layer produces the network's output.

Artificial Neural Network
One of the most commonly used tools to describe dependencies in geotechnics was standard regression. Recently, non-standard methods, including an artificial neural network and a fuzzy logic regression, have also been increasingly used as approximation tools. Moreover, the most common neural networks [46] used for function approximation (nonlinear regression) are multi-layer perceptron (MLP), radial basis function (RBF) networks and support vector machines (SVM). The feedforward networks consist of a series of layers ( Figure 2). The first layer has a connection from the network input. Each subsequent hidden layer has a connection from the previous layer. The final output layer produces the network's output. For training multilayer feedforward networks, any standard numerical optimization algorithm can be used to optimize the performance function, but there are a few that have shown excellent performance for neural network training. These optimization methods use either the gradient of the For training multilayer feedforward networks, any standard numerical optimization algorithm can be used to optimize the performance function, but there are a few that have shown excellent performance for neural network training. These optimization methods use either the gradient of the network performance with respect to the network weights or the Jacobian of the network errors with respect to the weights [47].
The Levenberg-Marquardt [48,49] method generally is the fastest training method. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton [50] backpropagation method is also quite fast. This algorithm requires more computation in each iteration and more storage than the conjugate gradient methods, although it generally converges in fewer iterations. Both of these methods tend to be less efficient for large networks (with thousands of weights), since they require more memory and more computation time. The Bayesian regularization [51] algorithm requires more time, but can result in good generalization for difficult, small or noisy data sets.
The architecture of the networks used can be described as I-H-O, where I is the number of inputs, H is the number of neurons in the hidden layer and O is the number in the output layer. Mostly, the definition of the network uses several types of transfer (activation) functions: log-sigmoid, hyperbolic tangent sigmoid, radial basis and some linear.
All neural network computation was performed using the Neural Network Toolbox for Matlab [52]. In all the considered examples, a shallow multilayer feedforward network with one hidden layer was applied. In input of nets, we take into account two parameters, or one of them, the organic content LOI T and the natural water content w. The output of nets was soil unit weight of organic soils γ t . In the calculations, five to eight neurons were used in the hidden layers. The Levenberg-Marquardt method was used in the training process. We used a log-sigmoid transfer function in the hidden layer and a linear function in the output layer.

Results
Based on the analysis of the exploratory research carried out in the study area at the Rzeszów site, with the CPTM probe used at a maximum depth 15.4 m below ground level, frequently occurring and varied organic soils were identified in the subsoil. The test results show the values of the cone resistance q c and sleeve friction f s , in addition to the friction ratio R f , which was used for the classification of the organic soil type using known criteria (p.2.1), as presented in the diagrams in The Levenberg-Marquardt [48,49] method generally is the fastest training method. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton [50] backpropagation method is also quite fast. This algorithm requires more computation in each iteration and more storage than the conjugate gradient methods, although it generally converges in fewer iterations. Both of these methods tend to be less efficient for large networks (with thousands of weights), since they require more memory and more computation time. The Bayesian regularization [51] algorithm requires more time, but can result in good generalization for difficult, small or noisy data sets.
The architecture of the networks used can be described as I-H-O, where I is the number of inputs, H is the number of neurons in the hidden layer and O is the number in the output layer. Mostly, the definition of the network uses several types of transfer (activation) functions: log-sigmoid, hyperbolic tangent sigmoid, radial basis and some linear.
All neural network computation was performed using the Neural Network Toolbox for Matlab [52]. In all the considered examples, a shallow multilayer feedforward network with one hidden layer was applied. In input of nets, we take into account two parameters, or one of them, the organic content LOIT and the natural water content w. The output of nets was soil unit weight of organic soils γt. In the calculations, five to eight neurons were used in the hidden layers. The Levenberg-Marquardt method was used in the training process. We used a log-sigmoid transfer function in the hidden layer and a linear function in the output layer.

Results
Based on the analysis of the exploratory research carried out in the study area at the Rzeszów site, with the CPTM probe used at a maximum depth 15.4 m below ground level, frequently occurring and varied organic soils were identified in the subsoil. The test results show the values of the cone resistance qc and sleeve friction fs, in addition to the friction ratio Rf, which was used for the classification of the organic soil type using known criteria (p.2.1), as presented in the diagrams in                Existing laboratory research has shown that, in the subsoil, a full spectrum of organic soils are present, from the low to high organic soil, which, due to ambiguous and simplified extant classification criteria [53,54], in comparison to the withdrawn national standard guidelines [55], is currently difficult to exactly define [56,57]. The organic content of the researched soils ranged from ca. 5% to ca. 85%, with a natural water content ranging from ca. 24% to ca. 418% and bulk density from 1.17 kN/m 3 to 2.25 kN/m 3 . The complete results of the index properties of organic soils at the study area at the Rzeszow site are summarized in Table 1.   Existing laboratory research has shown that, in the subsoil, a full spectrum of organic soils are present, from the low to high organic soil, which, due to ambiguous and simplified extant classification criteria [53,54], in comparison to the withdrawn national standard guidelines [55], is currently difficult to exactly define [56,57]. The organic content of the researched soils ranged from ca. 5% to ca. 85%, with a natural water content ranging from ca. 24% to ca. 418% and bulk density from 1.17 kN/m 3 to 2.25 kN/m 3 . The complete results of the index properties of organic soils at the study area at the Rzeszow site are summarized in Table 1.  Existing laboratory research has shown that, in the subsoil, a full spectrum of organic soils are present, from the low to high organic soil, which, due to ambiguous and simplified extant classification criteria [53,54], in comparison to the withdrawn national standard guidelines [55], is currently difficult to exactly define [56,57]. The organic content of the researched soils ranged from ca. 5% to ca. 85%, with a natural water content ranging from ca. 24% to ca. 418% and bulk density from 1.17 kN/m 3 to 2.25 kN/m 3 . The complete results of the index properties of organic soils at the study area at the Rzeszow site are summarized in Table 1.

Evaluation of the Soil Unit Weight Based on Laboratory Test Results
The next step after performing laboratory tests was to search for the relationships between the dependent (soil unit weight γ t ) and independent variables (natural water content w; organic matter content LOI T ).
The one-factor relationship between the soil unit weight and natural water content is presented and described in the formula in Figure 9a and the organic matter content in Figure 9b. It has been shown that there is a very strong relationship between the variables included, as evidenced by the high value of the coefficient of determination R 2 , for natural water content (0.978) and only slightly lower for organic matter content (0.946). The differences between the values of the soil unit weight, calculated separately for each of the variables, with the value determined in a laboratory, were comparable. In the case of the first variable-natural water content-the max. difference was 21.76%, and for the second variable-soil matter content-the max. difference was 20.46%, in relation to the expected value of the soil unit weight.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 17 The next step after performing laboratory tests was to search for the relationships between the dependent (soil unit weight γt) and independent variables (natural water content w; organic matter content LOIT). The one-factor relationship between the soil unit weight and natural water content is presented and described in the formula in Figure 9a and the organic matter content in Figure 9b. It has been shown that there is a very strong relationship between the variables included, as evidenced by the high value of the coefficient of determination R 2 , for natural water content (0.978) and only slightly lower for organic matter content (0.946). The differences between the values of the soil unit weight, calculated separately for each of the variables, with the value determined in a laboratory, were comparable. In the case of the first variable -natural water content-the max. difference was 21.76%, and for the second variable-soil matter content-the max. difference was 20.46%, in relation to the expected value of the soil unit weight.
(a) (b) Figure 9. The value of the soil unit weight estimated from laboratory tests was in reference to (a) natural water content and (b) organic matter content.
Additionally, based on the regression analysis, the relationship between empirically determined values of soil unit weight and leading parameters (natural water content and organic content) was determined and is described by the following formula: where γ t is the soil unit weight, LOIT is the organic content and w is the natural water content.
The differences between the values determined on the basis of laboratory tests and the model determined by the regression method, taking into account two variables, did not exceed 2.95 kN/m 3 (22.44%), with the factor R 2 = 0.797. Variable significance studies (w, LOIT) were conducted and showed that the variable representing the organic content was not significant. After eliminating the insignificant variable, it was found that the differences between the set value and the expected volume of the soil unit weight were much higher at 3.96 kN/m 3 (38.00%), in comparison to the case of the model with two variables. Due to the larger differences and the fact that the model with one explanatory variable was not very representative, the model with two variables was adopted for further analysis. The comparison of the results is presented in Figure 10. Additionally, based on the regression analysis, the relationship between empirically determined values of soil unit weight and leading parameters (natural water content and organic content) was determined and is described by the following formula: where γ t is the soil unit weight, LOI T is the organic content and w is the natural water content. The differences between the values determined on the basis of laboratory tests and the model determined by the regression method, taking into account two variables, did not exceed 2.95 kN/m 3 (22.44%), with the factor R 2 = 0.797. Variable significance studies (w, LOI T ) were conducted and showed that the variable representing the organic content was not significant. After eliminating the insignificant variable, it was found that the differences between the set value and the expected volume of the soil unit weight were much higher at 3.96 kN/m 3 (38.00%), in comparison to the case of the model with two variables. Due to the larger differences and the fact that the model with one explanatory variable was not very representative, the model with two variables was adopted for further analysis. The comparison of the results is presented in Figure 10.
For the analytical estimation of the soil unit weight, personally determined one and two-factor relationships were used; therefore, it was necessary to determine the value of the relative error in the relationships in relation to the values determined empirically in the laboratory.
The maximum and minimum values of relative error were calculated according to Formula (5): where P is number of cases, p = {1, . . . , P}, d p is the measured value and y p is the calculated value. For the analytical estimation of the soil unit weight, personally determined one and two-factor relationships were used; therefore, it was necessary to determine the value of the relative error in the relationships in relation to the values determined empirically in the laboratory.
The maximum and minimum values of relative error were calculated according to Formula (5): where P is number of cases, p={1,…, P}, dp is the measured value and yp is the calculated value. The comparison carried out for the relationship described earlier indicates a much larger value of the relative error of soil unit weight, which was noted for the one-factor relationship based on organic content (LOIT), while the lowest, comparable values were obtained for one-factor (w) and two-factor (w, LOIT) relationships. A summary of the calculated values of relative error (RE) of the soil unit weight is shown in Table 2. Table 2. Relative error (RE) of the evaluated soil unit weight of organic soil for use in single and two-factor relationships. show that evaluation of the soil unit weight, based only on organic content, is not recommended, but two-factors relationships, based on organic content and natural water content, result in better accuracy and a better way to describe the model.

Standard Regreesion
The most important task of our analyses is to compare the results of standard regression to the neural network approach in the problem of soil unit mass identification. For this reason, a division of the entire set of 135 measurements into two groups (70%, 30%) was assumed. Ninety-four of them were used as basic data to calculate the parameters of the regression fit function. The remaining 41 were used to test the predictions of the soil unit mass from the regression model with measured values. To determine the statistical evaluation during the analysis, the basic data set was randomly selected 250 times.
Four regression models of one variable were used: the two-parameter linear model (6), three-parameter polynomial model (7), two-parameter power model (8) and four-parameter power model (9). The comparison carried out for the relationship described earlier indicates a much larger value of the relative error of soil unit weight, which was noted for the one-factor relationship based on organic content (LOI T ), while the lowest, comparable values were obtained for one-factor (w) and two-factor (w, LOI T ) relationships. A summary of the calculated values of relative error (RE) of the soil unit weight is shown in Table 2. Table 2. Relative error (RE) of the evaluated soil unit weight of organic soil for use in single and two-factor relationships.

REγ (w) [%]
REγ The results show that evaluation of the soil unit weight, based only on organic content, is not recommended, but two-factors relationships, based on organic content and natural water content, result in better accuracy and a better way to describe the model.

Standard Regreesion
The most important task of our analyses is to compare the results of standard regression to the neural network approach in the problem of soil unit mass identification. For this reason, a division of the entire set of 135 measurements into two groups (70%, 30%) was assumed. Ninety-four of them were used as basic data to calculate the parameters of the regression fit function. The remaining 41 were used to test the predictions of the soil unit mass from the regression model with measured values. To determine the statistical evaluation during the analysis, the basic data set was randomly selected 250 times.
Four regression models of one variable were used: the two-parameter linear model (6), three-parameter polynomial model (7), two-parameter power model (8) and four-parameter power model (9).
F3(x) = p 1 e p 2 x (8) Additionally, the next four regression models of the two variables were used: the three-parameter surface model (10), five-parameter surface model (11,12) and eight-parameter surface model (13).
F5(x, y) = p 1 + p 2 x + p 3 y (10) F6(x, y) = p 1 + p 2 x + p 3 y + p 4 x 2 + p 5 x 2 y (11) F7(x, y) = p 1 + p 2 x + p 3 y + p 4 y 2 + p 5 xy 2 (12) F8(x, y) = p 1 + p 2 x + p 3 y + p 4 x 2 + p 5 y 2 + p 6 x 2 y + p 7 xy 2 + p 8 x 2 y 2 (13) The goodness of fit was checked using the coefficient of determination (14), mean relative error (15) and mean squared error (16): where n is number of cases, d p is the measured values, y p is the fitted values (predicted) and d p is the mean of the measured values. The comparison of the median of the goodness of fit obtained for the regression of the one variable model is presented in Table 3. The median is the most resistant statistic datum and is of central importance in robust statistics. The best model obtained in the analyses was the power model F4, using natural water content. It has 1.83% of the median average relative error for testing the model. The comparison of the MRE models, included in Table 3 for laboratory data, is presented using box-and-whisker diagrams in Figure 11. On each box, the central mark is the median, the edges of the box are the 25th and 75th percentiles, the whiskers extend to the most extreme data points, which are not considered outliers, and outliers are plotted individually.
Generally, we obtained better results using natural water content. Table 4 presents a comparison of the median of the goodness of fit obtained for the two variables regression models (F5-F8). The best result we obtained was for the F6 regression model. Generally, we obtained better results using natural water content. Table 4 presents a comparison of the median of the goodness of fit obtained for the two variables regression models (F5-F8). The best result we obtained was for the F6 regression model.  Table 4 for two laboratory parameters (LOIT, w) and are presented in Figure 12. The model F6 had the least outliers.   (LOIT, w). The left-hand plot (Figure 13a), with the blue line, is related to the base data, and the right-hand plot (Figure 13b), with the green line, is related to the test data. The coefficients of determination are very similar.  The comparison of the test results from the MRE of models are included in Table 4 for two laboratory parameters (LOI T , w) and are presented in Figure 12. The model F6 had the least outliers. Generally, we obtained better results using natural water content. Table 4 presents mparison of the median of the goodness of fit obtained for the two variables regression mode 5-F8). The best result we obtained was for the F6 regression model.  Table 4 for tw oratory parameters (LOIT, w) and are presented in Figure 12. The model F6 had the least outliers   (LOIT, w). The left-han ot (Figure 13a), with the blue line, is related to the base data, and the right-hand plot ( Figure 13b ith the green line, is related to the test data. The coefficients of determination are very similar.

Artificial Neural Networks Analysis
Next, we applied a neural regression model. In all the examples, standard multi-layer perception with one hidden layer was applied. In this case, the nets have only one element in the output vector (soil unit weight). The number of hidden neurons is obtained from a cross-correlation procedure. In the calculations, five to eight neurons were used in the hidden layers. Exactly the same pattern divided as that for the standard regression was used to learn and test the networks. The comparison of the median of the goodness of fit, obtained for a few architectures, is presented in Table 5.  20 repetitions of the network training were taken into account for each of the 250 pattern divisions. In this way, 1000 results were finally taken into account. Better prediction was obtained, like that in the standard regression for natural water content, in the input vector. A comparison of the MRE of the nets, included in Table 5, is presented in Figure 14.

Artificial Neural Networks Analysis
Next, we applied a neural regression model. In all the examples, standard multi-layer perception with one hidden layer was applied. In this case, the nets have only one element in the output vector (soil unit weight). The number of hidden neurons is obtained from a cross-correlation procedure. In the calculations, five to eight neurons were used in the hidden layers. Exactly the same pattern divided as that for the standard regression was used to learn and test the networks. The comparison of the median of the goodness of fit, obtained for a few architectures, is presented in Table 5. The presented results include the median obtained parameters for one element in the input vector. During the calculation of the values, 20 repetitions of the network training were taken into account for each of the 250 pattern divisions. In this way, 1000 results were finally taken into account. Better prediction was obtained, like that in the standard regression for natural water content, in the input vector. A comparison of the MRE of the nets, included in Table 5, is presented in Figure 14.
In this particular case, the network architecture does not affect the accuracy of approximation. Obtained mean relative errors of testing using organic matter content were in the range of 2.73-2.76%. Respectively, the MRE of testing using natural water content were in the range of 1.65-1.67%. In Table 6, the comparison of the median result, obtained for the nets with two elements in the input vector, is presented. In that approach, we also obtained a better result, like that in the standard regression with two independent variables. The mean relative errors of testing using two parameters were in the range of 1. In this particular case, the network architecture does not affect the accuracy of approximation. Obtained mean relative errors of testing using organic matter content were in the range of 2.73-2.76%. Respectively, the MRE of testing using natural water content were in the range of 1.65-1.67%. In Table 6, the comparison of the median result, obtained for the nets with two elements in the input vector, is presented. In that approach, we also obtained a better result, like that in the standard regression with two independent variables. The mean relative errors of testing using two parameters were in the range of 1.40-1.44%.  Table 6 for two laboratory parameters (LOIT, w), are presented in Figure 15. There is no clear difference in the results. The obtained results for the two inputs are a little better, which is the opposite case for the nets with one input. The smallest median average relative error test of the ANN was equal to 1.40%.     The comparison test results from the MRE of models, included in Table 6 for two laboratory parameters (LOI T , w), are presented in Figure 15. There is no clear difference in the results. The obtained results for the two inputs are a little better, which is the opposite case for the nets with one input. The smallest median average relative error test of the ANN was equal to 1.40%. In this particular case, the network architecture does not affect the accuracy of approximation. Obtained mean relative errors of testing using organic matter content were in the range of 2.73-2.76%. Respectively, the MRE of testing using natural water content were in the range of 1.65-1.67%. In Table 6, the comparison of the median result, obtained for the nets with two elements in the input vector, is presented. In that approach, we also obtained a better result, like that in the standard regression with two independent variables. The mean relative errors of testing using two parameters were in the range of 1.40-1.44%.  Table 6 for two laboratory parameters (LOIT, w), are presented in Figure 15. There is no clear difference in the results. The obtained results for the two inputs are a little better, which is the opposite case for the nets with one input. The smallest median average relative error test of the ANN was equal to 1.40%.    Figure 15. Comparison of the MRE of the nets with two elements in the input vector. Figure 16 shows detailed results for the ANN's prediction of the soil unit weight using two variables (LOI T , w) in the input vector. Figure 16a shows learning data and Figure 16b shows the test data. The coefficients of determination are very similar. Figure 16 shows detailed results for the ANN's prediction of the soil unit weight using two variables (LOIT, w) in the input vector. Figure 16a shows learning data and Figure 16b shows the test data. The coefficients of determination are very similar. Generally, the use of neural networks has allowed the soil unit weight values to be predicted based on laboratory tests with very high accuracy. The ANN regression models are slightly better than in the considered regression models. There is no clear difference in the architecture nets used. The best prediction neural networks were determined based on the lowest average relative error.

Conclusions
The presented results from the investigation of organic soils allowed for the estimation of the parameters necessary for geotechnical design and soil unit weight, without samples with undisturbed structures. The application of standard regression and neural networks to predict soil unit weight using laboratory tests was possible with acceptable prediction qualities. The maximum median values of the coefficient of determination obtained were equal, respectively, to 0.970 and 0.986. The sampling and preparation laboratory tests were very difficult, time-consuming and expensive. Therefore, the search for cheaper, faster and easier methods to determine soil unit weight in engineering practice is expected. The use of CPTM measurement results using standard regression and ANNs should be verified in the near future. Generally, the use of neural networks has allowed the soil unit weight values to be predicted based on laboratory tests with very high accuracy. The ANN regression models are slightly better than in the considered regression models. There is no clear difference in the architecture nets used. The best prediction neural networks were determined based on the lowest average relative error.

Conclusions
The presented results from the investigation of organic soils allowed for the estimation of the parameters necessary for geotechnical design and soil unit weight, without samples with undisturbed structures. The application of standard regression and neural networks to predict soil unit weight using laboratory tests was possible with acceptable prediction qualities. The maximum median values of the coefficient of determination obtained were equal, respectively, to 0.970 and 0.986. The sampling and preparation laboratory tests were very difficult, time-consuming and expensive. Therefore, the search for cheaper, faster and easier methods to determine soil unit weight in engineering practice is expected. The use of CPTM measurement results using standard regression and ANNs should be verified in the near future.