Next Article in Journal
Comparison of Immediate Neuromodulatory Effects between Focal Vibratory and Electrical Sensory Stimulations after Stroke
Previous Article in Journal
The Nitrogen Content in the Fruiting Body and Mycelium of Pleurotus Ostreatus and Its Utilization as a Medium Component in Thraustochytrid Fermentation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Multivariate Regression and Artificial Neural Network Modelling for Prediction of Physicochemical Properties of Grape-Skin Compost

by
Tea Sokač Cvetnić
1,
Korina Krog
1,
Davor Valinger
1,
Jasenka Gajdoš Kljusurić
1,*,
Maja Benković
1,
Tamara Jurina
1,
Tamara Jakovljević
2,
Ivana Radojčić Redovniković
1 and
Ana Jurinjak Tušek
1
1
Faculty of Food Technology and Biotechnology, University of Zagreb, Pierottijeva 6, 10 000 Zagreb, Croatia
2
Croatian Forest Research Institute, Cvjetno naselje 41, 10 450 Jastrebarsko, Croatia
*
Author to whom correspondence should be addressed.
Submission received: 21 February 2024 / Revised: 7 March 2024 / Accepted: 13 March 2024 / Published: 16 March 2024

Abstract

:
The reusability of by-products in the food industry is consistent with sustainable and greener production; therefore, the aim of this paper was to evaluate the applicability of multiple linear regression (MLR), piecewise linear regression (PLR) and artificial neural network models (ANN) to the prediction of grape-skin compost’s physicochemical properties (moisture, dry matter, organic matter, ash content, carbon content, nitrogen content, C/N ratio, total colour change of compost samples, pH, conductivity, total dissolved solids and total colour change of compost extract samples) during in-vessel composting based on the initial composting conditions (air-flow rate, moisture content and day of sampling). Based on the coefficient of determination for prediction, the adjusted coefficient of determination for calibration, the root-mean-square error of prediction (RMSEP), the standard error of prediction (SEP), the ratio of prediction to deviation (RPD) and the ratio of the error range (RER), it can be concluded that all developed MLR and PLR models are acceptable for process screening. Furthermore, the ANN model developed for predicting moisture and dry-matter content can be used for quality control (RER >11). The obtained results show the great potential of multivariate modelling for analysis of the physicochemical properties of compost during composting, confirming the high applicability of modelling in greener production processes.

1. Introduction

The modern way of life has led to the need to solve the problem of the large quantities of food waste produced worldwide. Looking at the amount of waste itself, the biggest problem is biodegradable waste, which consists of the leftovers from the fruits and vegetables that we consume in our daily lives [1,2]. Its improper disposal can pose a risk to the environment and human health, as decomposition of this type of waste in landfills leads to the production of methane, which consequently contributes to the creation of the greenhouse effect [3]. Despite this fact, biodegradable waste has significant potential to contribute to the circular bioeconomy because it can be used in the production of fertilizers, non-fossil fuels and soil improvers [4,5]. One of the processes developed with this goal in mind is composting [6]. Composting is considered an environmentally friendly and profitable method for processing organic waste. During the composting process, polymeric waste materials are degraded by the growth of various micro-organisms such as fungi and bacteria. This processing method is quite complex in many ways, as it involves a variety of processes ranging from microbiological to thermodynamic, which are mutually dependent [7].
Grape pomace is the most important by-product of wine production. It represents 20–25% of the total mass of processed grapes and varies according to grape variety, degree of ripeness and type of press used in processing [8]. As described by Martinez Salgado et al. [9], after pressing, grape pomace has a moisture content of around 20–30% (w/w), a C:N ratio between 40 and 45:1, a pH between 3 and 6, low electrical conductivity, a large amount of organic matter and organic forms of micro- and macro-nutrients that are mineralised over time and can therefore be considered a good substrate for composting. Grape pomace can be composted aerobically to create a soil conditioner that increases the soil’s ability to retain water and replenish its organic carbon and nutrient content [10]. The effectiveness of composting is influenced by environmental variables such as aeration, temperature, humidity and pH, as well as substrate parameters, which include the C/N ratio (carbon/nitrogen ratio) and nutrient content [11].
There are several ways to compost food waste, including in-vessel procedures, aerated static piles and windrows [12]. A number of issues, including land availability, operational complexity, capital and running costs and the possibility of annoyances, should be taken into account while choosing a composting process [13]. Windrow systems are the least expensive, but they are also noted for taking a long time to complete—typically more than 60 days. For maturing and curing, windrow systems have a high land need. In summary, in-vessel composting is a significant technique, especially in contexts where rapid and controlled decomposition, odour control and space efficiency are essential considerations. It plays a crucial role in managing organic waste in urban and industrial settings, contributing to sustainable waste-management practices [1].
In order to achieve the best possible efficiency when setting up the process, it is necessary to consider all the above listed variables and their interactions. As the specific relationships are often non-linear, a variety of effects needs to be thoroughly investigated both theoretically and practically. As a result, both the creation of a logical framework for process design and the understanding of complicated dynamic interactions can benefit from the use of mathematical modelling tools [14]. Mathematical models of the composting process are helpful to obtain information on how different process factors and conditions such as substrate composition, oxygen concentration, pollutant concentrations, composting time, temperature, etc. affect compost quality during the process [15,16,17,18,19]. In the literature, there are available examples of efficient application of multivariate techniques like cluster analysis and principle component analysis (PCA) for distinguishing compost samples of different origins [20], application of PCA for analysis of compost quality assessment [21], application of partial least squares regression (PLSR) for evaluating compost maturity [22], application of multiple linear regression for estimation of maturity of compost from food waste and agro-residues [23] and application of machine-learning methods (including random forest, extreme gradient boosting, Light Gradient Boosting Machine and Multilayer Perceptron networks) for prediction of germination index and carbon-to-nitrogen ratio [24].
Taking into account that understanding the grape-skin composting process variables’ interactions can significantly contribute to the final compost quality, in this work, the effects of the initial composting conditions (initial moisture content, air-flow rate and day of sampling) on physicochemical properties of grape-skin compost (moisture, dry matter, organic matter, ash content, carbon content, nitrogen content, C/N ratio, total colour change of compost samples, pH, conductivity, total dissolved solids, total colour change of compost extract samples) were analysed. The relationships between the listed input and output process variables were modelled using multiple linear regression (MLR), piecewise linear regression (PLR) and artificial neural network modelling (ANN). The predictive efficiencies of the mentioned models were estimated and compared with the motivation to use developed models for compost properties prediction based on different initial process conditions.

2. Materials and Methods

2.1. Materials

2.1.1. Biowaste for Composting

Biowaste used for composting was grape skins from the white grape variety Vitis vinifera cv. Graševina, harvested in 2021 (Kutjevo, Croatia). Grape pomace was stored in a freezer at −18 °C before conducting experiments. Prior to starting the composting process, seeds were separated from the skin by sieving and used in another experiment.

2.1.2. Chemicals

Sodium hydrogen carbonate was supplied by Kemika (Zagreb, Croatia). Urea was purchased from Gram-Mol (Zagreb, Croatia).

2.2. Methods

2.2.1. Grape-Skin Composting Process

The composting of m = 1.9 kg of grape skins was carried out in laboratory batch reactors with a volume of V = 5 L. The dimensions of the reactor were: diameter d = 16 cm and height L = 25 cm. The reactors were insulated with a 5 cm thick insulating layer. Initial moisture content range (IMC) and air flow (AF) range were selected according to the available literature data [25,26] to ensure optimal composting conditions. Composting experiments were performed according to central composite design (Statistica 14.0 Tibco Software Inc., Palo Alto, CA, USA) of experiment as follows: (1) IMC = 50% and AF = 0.50 L/min; (2) IMC = 50% and AF = 1.25 L/min; (3) IMC = 65% and AF = 0.50 L/min; (4) IMC = 65% and AF = 1.40 L/min; (5) IMC = 65% and AF = 0.35 L/min; (6) IMC = 50% and AF = 2.00 L/min; (7) IMC = 57.5% and AF = 1.70 L/min; (8) IMC = 57.5% and AF = 0.43 L/min; and (9) IMC = 57.5% and AF = 1.06 L/min
To ensure aerobic conditions in the reactor, continuous aeration was carried out for the entire duration of the process (30 days). Temperature was monitored with digital thermometers (ST-9263B, Sagar Scientific & Instruments, Maharashtra, India) placed in the centre of the reactor. The ambient temperature was approximately 19 °C during all experiments. During the composting process, samples were taken at 48 h intervals, resulting in 16 samples from a single reactor, and, in total, 144 samples.

2.2.2. Physicochemical Analyses of the Compost Samples

Physicochemical analyses of the compost samples included dry-matter content (DM) and moisture content (MC) measurement, pH value, total dissolved solids (TDS) and conductivity (S) measurement, organic matter (OM) and ash content measurement (AC), carbon content (CC), nitrogen content (NC) and total colour change (ΔE) measurements.

Dry-Matter Content and Moisture Content of the Compost Samples

The dry-matter and moisture content of the compost samples was determined gravimetrically by convectional drying at T = 105 °C, for t = 24 h [27].

Extraction Procedure, pH Value, Conductivity and Total Dissolved Solids of Compost Samples

A compost sample was mixed with distilled water in a ratio of 1:10 (w/v) and the extraction procedure was carried out on a magnetic stirrer at 150 rpm for 1 h. The resulting mixture was centrifuged (10,000 rpm) for 10 min and then filtered. In the filtrate, pH [28,29] was determined using a pH meter (914, Metrohm, Herisau, Switzerland) and conductivity and total dissolved solids were determined using a conductometer (SevenCompact, MettlerToledo, Greifensee, Switzerland).

Organic Matter Content and Ash Content of Compost Samples

The compost samples dried at T = 105 °C for t = 24 h were burned at T = 500 °C for t = 5 h in a muffle furnace (B410, Nabertherm, Lilienthal, Germany). The percentage loss of volatiles was expressed as a fraction of total organic matter [30], while the mass remaining after burning was expressed as the ash fraction [28].

Carbon and Nitrogen Content of the Compost Samples

Total carbon and nitrogen were determined using an elemental analyser CNS 2000 (Leco, St Joseph, MI, USA) and a spectrophotometer (LaboMed UV-VIS, Los Angeles, CA, USA) [31].

Total Colour Change of the Compost Samples and Compost Extract Samples

The colour of all compost samples during the composting process and the prepared water extracts from compost was determined using a PCE-CSM3 colorimeter (PCE Instruments, Meschede, Germany). The total colour change of the compost and the corresponding compost extracts (ΔE) was determined according to the Equation (1):
Δ E = L * L 0 * 2 + a * a 0 * 2 + b * b 0 * 2
where L0*, a0*, and b0* are the values of the Hunter coordinates of the samples/extracts of the initial substrate samples, L*, a* and b* are the values of the Hunter coordinates of the compost/compost extracts during the composting process.

2.2.3. Statistical Analyses and Mathematical Modelling

Descriptive Statistics

All measurements were repeated three times, and the results are reported as mean ± standard deviation. Descriptive statistical analyses, including minima, maxima, averages and standard deviations of the analysed variables were performed using Statistica 14.0 (Tibco Software Inc., Palo Alto, CA, USA). The normality of data distribution was tested by the Shapiro–Wilk test and the uniformity of variance with Levene’s test. Spearman’s correlation matrix was used to analyse the relationships between the process variables and all analysed properties of the compost samples because data distribution analysis showed that data were not normally distributed.

Multiple Linear Regression Modelling, Pricewise Linear Regression Modelling and Artificial Neural Network Modelling

In further calculations, it was assumed that the measured physicochemical properties of the compost samples (i = 1, …, 12: moisture content (Y1), dry-matter content (Y2), organic matter content (Y3), ash content (Y4), carbon content (Y5), nitrogen content (Y6), C/N ratio (Y7), total colour change of compost samples (Y8), pH (Y9), conductivity (Y10), total dissolved solids (Y11), total and colour change of compost extract samples (Y12)) gathered through nine independent experiments can be described as a function of initial moisture content (X1), air flow (X2) and day of sampling (X3) according to Equation (2):
Y i = f X 1 , X 2 , X 3
Multiple linear regression (MLR) (Equation (3)), piecewise linear regression (PLR) (Equation (4)) and artificial neural network (ANN) models were used to evaluate the relationship between input and output variables.
Y i = b 0 + b 1 · X 1 + b 2 · X 2 + b 3 · X 3
Y i = b 01 + b 11 · X 1 + b 21 · X 2 + b 31 · X 3 for   Y i b n + b 02 + b 12 · X 1 + b 22 · X 2 + b 32 · X 3 for   Y i > b n
The MLR model parameters (Equation (3)) and PLR model parameters (Equation (4)) were calculated using the Levenberg-Marquardt algorithm implemented in Statistica 14.0 (Tibco Software Inc., Palo Alto, CA, USA). Using the least squares method, the programme searches for optimal solutions in the parameter space of the function. The calculations were repeated 50 times with a convergence level of 10−6 and a confidence interval of 95% [32]. The data set (432 data points for each output variable) for MLR and PLR modelling was randomly split 70:30 into a calibration and a prediction data set. The applicability of the developed calibration models was estimated using the coefficient of determination for calibration (Rcal2), the adjusted coefficient of determination for calibration (Rcal2adj), the root-mean-square error for calibration (RMSE) and the F-value of the model. Predictive performance of the models was estimated using the coefficient of determination for prediction (Rpred2), the adjusted coefficient of determination for calibration (Rpred2adj), the root-mean-square error of prediction (RMSEP), the standard error of prediction (SEP), the ratio of prediction to deviation (RPD) and the ratio of the error range (RER) [33].
In addition, multilayer perceptron (MLP) ANNs were used to predict the physicochemical properties of the compost samples. ANN models were developed separately for each analysed output-process variable. The ANN models contained an input layer, a hidden layer and an output layer. The input layer had three neurons representing the conditions of the composting process (moisture content, air flow and day of sampling), the output layer had one neuron (moisture content, dry-matter content, organic-matter content, ash content, carbon content, nitrogen content, C/N ratio, total colour change of compost samples, pH value, conductivity, total amount of dissolved solids or total colour change of compost extract samples) and the number of neurons in the hidden layer varied between four and 13 and was selected by the algorithm. For the activation functions of the hidden layer and the output layer, the identity, logistic, hyperbolic tangent and exponential activation functions were randomly selected. For ANN modelling, the data set was split 70:30 into a calibration and a prediction data set. In addition, the calibration dataset was split into 70% for network training, 15% for network testing and 15% for model validation. The backpropagation algorithm was used for model training. The applicability of the developed calibration models was estimated using the coefficient of determination for calibration (Rcal2), the adjusted coefficient of determination for calibration (Rcal2adj), and the root-mean-square error for calibration (RMSE). Prediction performance of the models was estimated based on the coefficient of determination for prediction (Rpred2), the adjusted coefficient of determination for calibration (Rpred2adj), the root-mean-square error for prediction (RMSEP), the standard error of prediction (SEP), the ratio of prediction to deviation (RPD) and the ratio of the error range (RER) [33].

3. Results and Discussion

3.1. Physicochemical Properties of Compost Samples

Composting promotes the recycling of organic waste, yielding a final product with substantial bioenergy potential and considerable nutritional benefits for the soil. The shift towards more effective compost production and management necessitates a comprehensive grasp of the entire process, the materials utilized, and the physical attributes of those materials [34]. Therefore, in this work, the effects of different initial moisture contents, air-flow rates and days of sampling on the physicochemical properties of the compost samples (moisture content, dry-matter content, organic-matter content, ash content, carbon content, nitrogen content, C/N ratio, total colour change of the compost samples, pH, conductivity, total dissolved solids and total colour change of the compost-extract samples) during the grape-skin composting process were analysed. The minimum, maximum and standard deviations of the analysed physicochemical properties are listed in Table 1.
As shown in Table 1, the moisture content in the experiments ranged between 53.612 and 65.540%, which corresponds to the optimal range for the composting process [35]. The highest average value for moisture content was obtained in experiment 9, which had the lowest initial value for moisture content at the beginning of the composting process and the highest value for air-flow rate. Also, the highest values for total dissolved solids (TDS) and conductivity were achieved in experiment 9. Similar results were obtained by Peng et al. [36]. Higher aeration rate accelerates degradation reactions that result in salts and, consequently, the total dissolved solids and conductivity values are higher. While the release of mineral salts like phosphate and ammonia salts through the breakdown of organic materials can cause an increase in conductivity values, the decline in conductivity values during the composting process is directly related to increased concentrations of nutrients like nitrates and nitrites [37]. Moreover, experiment 2 was conducted at the lowest initial moisture content and higher air-flow rate, which resulted in a greater colour change of the compost compared to the other trials, confirming substrate degradation [38]. In addition, the highest average pH value was achieved in the experiment 3. The pH value is related to the microbial degradation during the composting process. In the first phase, it decreases due microbial degradation of organic matter and formation of organic acids. Later, the pH increases due to acid consumption by micro-organisms [28,35]. Experiment 7 had the lowest values for dry-matter and organic-matter content and, consequently, lower carbon content and C/N ratio, but this experiment was characterised by the highest value for ash content. According to the literature [28], organic-matter content and ash content have a reciprocal relationship, with a higher organic-matter content leading to a lower ash content.
Spearman’s correlation matrix was used to determine the relationships between the conditions of the composting process and physicochemical properties of the compost during the process (Table 2). Significant correlations are marked in bold. Results showed that initial moisture content was positively correlated with compost moisture (r = 0.6831), ash content (r = 0.2220) and compost pH (r = 0.1827). The results also showed that the initial moisture content was negatively correlated with the dry-matter content of the compost (r = −0.6831) and with the organic-matter content (r = −0.2204). The results obtained are in agreement with the results of Makan et al. [39] who showed that the initial moisture content had a significant effect on aerobic composting and that higher moisture contents are better for composting organic waste in vascular bioreactors. Similarly, Yeh et al. [40] showed that an initial moisture content of 55–70% is optimal for effective composting of food waste. The correlation matrix also showed that air-flow rate was negatively correlated with organic-matter content (r = −0.1297), with carbon content (r = −0.2236), with C/N ratio (r = −0.1727) and with total compost colour change (r = −0.2782) and positively correlated with ash content (r = 0.1306), nitrogen content (r = 0.1533), pH (r = 0.1057), TDS (r = 0.2835), conductivity (r = 0.2537) and overall colour change of the compost extracts (r = 0.2127). The results presented are consistent with the literature which states that aeration rate is one of the most important variables in the composting process, and that variations in aeration rate affect temperature, moisture content and oxygen-supply rate, among other variables [41]. The compost material will dry out and cool down if the aeration rate is higher than the optimal rate; on the other hand, an aeration rate that is too low will result in an oxygen deficit that prevents the microbes from receiving enough oxygen to support their activity [42]. The correlations between sampling day and the variables of the composting process were also analysed using the correlation matrix. It can be seen that sampling day was negatively correlated with the dry-matter content (r = −0.3170), organic matter content (r = −0.5370), carbon content (r = −0.2548) and the C/N ratio (r = −0.3527), while it was positively correlated with moisture content (r = 0.3170), ash content (r = 0.5362), nitrogen content (r = 0.3032), pH value (r = 0.5751), TDS value (r = 0.5503), conductivity (r = 0.5539) and the total colour change of the compost extracts (r = 0.4024). In general, it can be concluded that a higher initial moisture content, a higher air-flow rate and a longer duration of the composting process lead to greater decomposition of organic matter.

3.2. Multiple Linear Regression, Piecewise Linear Regression and Artificial Neural Network Models for Prediction of Physicochemical Properties of Compost during the Composting Process

The effects of initial moisture content, air flow and day of sampling on physicochemical properties of compost during the composting process were analysed. Multiple linear regression, partial linear regression and artificial neural network models were developed and their performance for predicting the physicochemical properties of grape-skin compost was evaluated. The applicability of the calibration models to describe the physicochemical properties of compost was estimated using R2, R2adj and RMSE. The applicability of the model prediction was estimated based on Rpred2, Rpred2adj and RMSEP, SEP, RPD and RER.

3.2.1. Multiple Linear Regression Models

Parameters of the MLR models and the PLR models were estimated using the Levenberg–Marquardt algorithm implemented in Statistica 14.0 with a confidence interval of 95% and the values are given in Supplementary Table S1. For the MLR models, the best agreement was obtained between the experimental data and the data predicted by the model for moisture content (Figure 1a) and dry-matter content (Figure 1b) (Rcal2 = 0.779, Rcal2adj = 0.777, RMSE = 2.772%, Rpred2 = 0.738, Rpred2adj = 0.738, RMSEP = 2.781%, SEP = 0.245, RPD = 1.948, RER = 7.274), followed by pH (Figure 1i) (Rcal2 = 0.447, Rcal2adj = 0.441, RMSE = 0.883, Rpred2 = 0.391, Rpred2adj = 0.376, RMSEP = 0.915, SEP = 0.081, RPD = 1.278, RER = 5.167) and the total dissolved solids content in the compost (Figure 1j) (Rcal2 = 0.422, Rcal2adj = 416, RMSE = 318.47 mg/L, Rpred2 = 0.396, Rpred2adj = 0.382, RMSEP = 342.322 mg/L, SEP = 30.139 mg/L, RPD = 1.275, RER = 5.308). The highest scatter between model-predicted data and experimental data was found for nitrogen content (Figure 1f) (Rcal2 = 0.158, Rcal2adj = 0.149, RMSE = 0.424%, Rpred2 = 0.123, Rpred2adj = 0.102, RMSEP = 0.583%, SEP = 0.051%, RPD = 0.806, RER = 3.789). The MLR calibration model’s quality was also estimated using residual analysis. MLR model residual-analysis results are presented in Supplementary Figure S1 for the model with the highest Rcal2 (moisture content model) and for the model with the lowest Rcal2 (total dissolved solids). This approach is based on examining the residuals’ patterns [43]. It became apparent that the residuals for MLR models were normally distributed (Supplementary Figure S1) because the normality condition was fulfilled because the residual plots were distributed approximately in a straight line. The normal distribution of the residuals was also confirmed by the bell-shaped histograms that show the measurement distribution [43]. The plots of the predicted values vs. residuals the show that there is no pattern in the residuals, suggesting that the models adequately describe the experimental data. Furthermore, it was discovered that the residuals ranged around the central value without any clear outliers, indicating that the degree of randomisation was suitable and that the order of testing had no effect on the results [44].
Statistical analysis of the MLR models showed that the model parameters b1 (coefficient with the initial moisture content), b2 (coefficient with the air flow rate) and b3 (coefficient with the sampling day) had a significant influence (p < 0.05) on all 12 analysed model outputs (Table S2). It is important to mention that Palechor-Trochez et al. [45] reported strong correlations between the overall change in organic carbon and the change in colour coordinates of the compost during the composting process. Furthermore, the colour changes during the composting process could be due to the presence of dissolved and particulate organic matter [46].
A large F-value (greater that F-critical = 2.0838 for analysed data set), with a small p-value (p < 0.001) for developed MLP models implies that there is a general relationship between the response and the predictors [47]. The relationship between observed and model predicted data was also estimated based on the R2 value. According to Henseler et al. [48], Hair et al. [49] and Hussain et al. [50] an R2 value of 0.75 is considered substantial, an R2 value of 0.50 is considered moderate and an R2 value of 0.26 is considered as weak. The suitability of the developed MLR models for predicting the physicochemical properties of grape skins during composting was also estimated using the ratio of prediction to deviation (RPD) and the ratio of the error range (RER). Models with RPD < 1.4 are considered non-reliable, those with RPD in range from 1.4 to 2 are considered fair, while models with RPD > 2 are described as excellent models [51]. Furthermore, models with RER > 4 are acceptable for data screening, models with RER > 10 can be used for quality control, and models with RER > 15 can be used for quantification [52]. Therefore, based on Rpred2, only the MLP models developed for the prediction of moisture content and dry-matter content can be considered moderately suitable and, based on the RPD values, only these two models can be considered reliable. On the other hand, based on RER values, all MLP models can be accepted for screening except the model predicting nitrogen content. Therefore, the proposed models need to be improved. Similar results were presented by Chikae et al. [23] where multiple linear regression modelling was applied for the prediction of germination index based on pH value, ammonium concentration, acid phosphate activity and esterase activity and R2 of 0.791 was obtained. On the other hand, Fouguira et al. [53] showed very efficient application (R2 > 0.990) of nonlinear multiple regression modelling for prediction of pH, carbon-to-nitrogen ratio and organic-matter content based on initial waste composition.

3.2.2. Piecewise Linear Regression Models

The basic idea behind PLR is that data should be modelled with the regression function piecewise when following different linear trends in different regions of the data [54]. The results obtained in this paper show that the developed PLR models describe the experimental results with higher accuracy than the proposed MLP models (Supplementary Table S2). For all PLR models, developed Rpred2 value was greater than 0.63.
The best agreement between the experimental data and the data predicted by the model was again obtained for the moisture content (Figure 2a) and the dry-matter content (Figure 2b) (Rcal2 = 0.837, Rcal2adj = 0.835, RMSE = 2.501%, Rpred2 = 0.834, Rpred2adj = 0.831, RMSEP = 2.144%, SEP = 0.189%, RPD = 2.526, RER = 9.433), followed by the total colour change of the compost (Figure 2h) (Rcal2 = 0.781, Rcal2adj = 0.779, RMSE = 1.764, Rpred2 = 0.777, Rpred2adj = 0.772, RMSEP = 1.925, SEP = 0.169, RPD = 2.159, RER = 8.802) and pH (Figure 2i) (Rcal2 = 0.784, Rcal2adj = 0.782, RMSE = 0.548, Rpred2 = 0.766, Rpred2adj = 0.745, RMSEP = 0.692, SEP = 0.043, RPD = 2.116, RER = 8.611). The highest scatter between model and experimental data was found for the ash content (Figure 3d) (Rcal2 = 0.674, Rcal2adj = 0.670, RMSE = 3.322%, Rpred2 = 0.635, Rpred2adj = 0.624, RMSEP = 3.434%, SEP = 0.302%, RPD = 1.708, RER = 7.111). The improvement in the applicability of the PLR models is particularly evident in the PLR model, which was developed to describe the nitrogen content in the compost. The MLP model developed to describe the nitrogen content was described by Rpred2 = 0.128, RMSEP = 6.738%, RPD = 1.065 and RER = 4.575 while the PLR model developed to describe the nitrogen content was described by Rpred2 = 0.647, RMSEP = 0.261%, RPD = 1.947 and RER = 7.316. It should be noted that the PLR model guarantees an Rpred2 value that is about five times higher, an RMSEP value that is about 25 times lower and RPD and RER values that are about two times higher than those for the MPL models. The residual analysis was also performed for the PLR calibration models for moisture content and ash content (Figure S1). PLR models’ goodness of fit was confirmed through a normal probability plot of the residuals, the predicted values versus residuals versus plot, histogram of the residuals and residuals versus the order of the data plot (Figure S1). The statistical analysis of the PLR models showed the same trend as for the MLR models (Table S1); in particular, the parameters b1 (coefficient connected to the initial moisture content), b2 (coefficient connected to the air flow rate) and b3 (coefficient connected to the sampling day), which had a significant influence (p < 0.05) on all 12 analysed model outputs. Based on Rpred2 (Table S2) only the PLR models developed for the prediction of moisture content and dry-matter content can be considered substantial, while the other proposed PLR models can be considered moderate (Rpred2 > 0.5). Furthermore, based on the RPD values the PLR models developed for the prediction of moisture content, dry-matter content, total colour change of the compost samples and prediction of pH can be considered excellent (RPD > 2), while the other models can be considered fair (1.4 < RPD < 2). As for the MLR models, based on RER values, all PLR models can be acceptable for screening (4 < RER < 10). Based on the obtained results, it can be concluded that PLR models can be used for the description and prediction of physiochemical properties of compost during the composting process. Similarly, Costello et al. [55] developed an efficient pricewise model (R2 > 0.75) for the prediction of pH and electrical conductivity based on shoot dry weight. The mentioned authors compared MLR and PLR modelling performance and showed that PLR models are more efficient.

3.2.3. Artificial Neural Network Models

Multilayer perceptron (MLP) neural networks were developed to further improve the prediction of physiochemical properties of grape-skin compost during the composting process. ANNs were compared to PLR modelling as nonlinear models, and it was expected that the ANNs could better and more accurately describe the experimental data [56]. ANN models were developed individually for each of the selected physicochemical properties of the compost. As shown in Table 3 and Figure 3, the developed ANN models show better performance compared to PLR models and, especially, MLR models. The developed ANN models provided good agreement between the experimental values and the values predicted by the model at learning, testing and validation levels, as the coefficients of determination were above 0.70 in all three phases and the errors of the model were low (RMSE < 0.014).
Table 3 shows the ANNs that were selected as optimal for predicting specific composting properties based on R2 and RMSE for the training, test and validation datasets, as well as for considering the number of neurons in the hidden layer. A lower number of neurons in the hidden layer was considered advantageous as it implies a simpler network structure [57]. Results showed that the best agreement between the experimental data and the data predicted by the ANN model was obtained for moisture content (Rpred2 = 0.9050, Rpred2adj = 0.9028, RMSEP = 1.7078%, SEP = 0.1504%, RPD = 3.1786, RER = 11.8455), followed by dry-matter content (Rpred2 = 0.9038, Rpred2adj = 0.9015, RMSEP = 1.7173%, SEP = 0.1512%, RPD = 3.1550, RER = 11.7799) and pH value of the compost (Rpred2 = 0.8322, Rpred2adj = 0.8282, RMSEP = 0.4896, SEP = 0.0431, RPD = 2.3890, RER = 9.6603). The highest dissipation between model and experimental data was obtained for the C/N ratio (Rpred2 = 0.6542, Rpred2adj = 0.66483, RMSEP = 3.6098, SEP = 0.3178, RPD = 1.9895, RER = 8.5399). Based on Rpred2, only the ANN models developed for the prediction of moisture content, dry-matter content, organic-matter content, pH and conductivity can be considered substantial (Rpred2 > 0.75), while the other ANN models can be considered moderate (Rpred2 > 0.65). The goodness of fit of the ANN calibration models for moisture content and C/N content (Figure S1) was confirmed through a normal probability plot of the residuals, the predicted values versus residuals versus plot, histogram of the residuals and residuals versus the order of the data plot (Figure S1). Moreover, the ANN model developed for the prediction of moisture content, dry-matter content, organic-matter content and the pH of grape-skin compost samples can be considered excellent (RPD > 2) based on the RPD values, while the other models can be considered fair (RPD > 1.6). The applicability improvement of the ANN models is also clear from the analysis of the RER values. Based on the RER values, the ANN models developed for the prediction of moisture and dry-matter content can be used for quality control (RER > 11), while all other models are acceptable for screening (RER > 7.8). Based on the results obtained, it can be concluded that the ANN models can be used for the description and prediction of the physiochemical properties of compost during the composting process with higher precision than the MLP or PLR models. Usage of the developed ANN models can contribute to reduction of number of the necessary experiments and provide an efficient prediction of the properties of the compost depending on the initial conditions of the composting process. The results obtained are in agreement with the results presented by Hosseinzadeh et al. [58] in which the ANN models provided a better prediction for the recovery of total nitrogen and total phosphorus from waste by vermicomposting than the MLR models. Furthermore, the superior prediction performance of ANN modelling compared to multiple linear regression was also shown by: (i) Dumenci et al. [59] for evaluation of olive mill waste compost based on composting mixture composition, (ii) Singh et al. [60] for modeming of compost production under different climate conditions, (iii) by Shi et al. [61], for prediction of humic acid content in the final compost based on the carbon-to-nitrogen content, initial moisture content, type of inoculant and composting day and (iv) Abdi et al. [56] for prediction of electrical conductivity, pH, carbon-to-nitrogen ratio and germination index, based on inlet-air rates, initial carbon-to-nitrogen ratios of 18 and the addition of coco peat biochar.
Comparing the results obtained using all three modelling approaches, it can be concluded that grape-skin compost’s physicochemical properties can be predicted based on the initial composting conditions. For both the MLR and PLR models, results showed that all estimated regression coefficients are significant, indicating the importance of selected input variables on the composting efficiency. Furthermore, it can also be noticed that the ANN modelling approach was superior in prediction performance, indicating the highly nonlinear nature of the interactions between input and output compositing variables.
To evaluate the importance of the input variables for the results of the ANN model, global sensitivity analysis was performed. As shown in Figure 4, moisture content was found to be the most important variable for the composting process of grape skins. According to Gurusamy et al. [62], maintaining an efficient composting process requires optimisation of moisture content; while high moisture content limits the mass transfer of air, low moisture content inhibits microbial activity, leading to the development of unstable and immature compost. The global sensitivity coefficient for initial moisture content was above 50% for all output variables except total dissolved solids and conductivity values. For carbon content and nitrogen content, the global sensitivity coefficient for initial moisture content was above 90%. This is in line with the results presented by Ghanney et al. [63] where the importance of moisture content for the change in carbon content during the composting process was also confirmed. These authors found that the organic-matter content of composting straw and cow manure decreased significantly at a moisture content of 65%, while the loss was lower at a moisture content of 45%. The global sensitivity analysis of the ANN models for the total amount of dissolved solids and the conductivity values of the compost showed that all three process variables have approximately the same effect on the analysed output variables. This can be explained by the fact that the change in electrical conductivity during the composting process correlates with the decomposition of organic matter [64], which is influenced by the compost moisture and oxygen-supply rate [65].

4. Conclusions

The obtained results indicate the importance of all these analysed process variables on the compost’s physical and chemical properties. Furthermore, based on RER values, all the MLP models can be accepted for screening except the model predicting nitrogen content, while the PLR and ANN models can be used for the description and prediction of the physiochemical properties of compost during the composting process. Application of the modelling approach in composting process analysis can contribute to a reduction in the number of necessary experiments and provide an efficient prediction of the properties of the compost depending on the initial conditions of the composting process.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/bioengineering11030285/s1, Table S1. Coefficients and statistics of multiple linear regression models for prediction of physicochemical properties of compost during the composting process based on initial moisture content, air-flow rate and sampling day. Significant coefficients and standard errors (SE) are marked bold. Figure S1. Analysis of the residuals for the: (a) MLR model describing moisture content, (b) MLR model describing total dissolved solids values, (c) PLR model describing moisture content, (d) PLR model describing ash content, (e) ANN model describing moisture content and (f) ANN model describing C/N ratio. (1) normal probability plot, (2) predicted value plots vs. residual plot, (3) histogram of residuals, (4) case vs residual plot. Table S2. Coefficients and statistics of piecewise linear regression models for prediction of physicochemical properties of compost during the composting process based on initial moisture content, air flow and sampling day. Significant coefficients and standard errors (SE) are marked bold.

Author Contributions

Conceptualization, A.J.T. and T.S.C.; methodology, A.J.T., T.S.C. and T.J. (Tamara Jakovljević); software, A.J.T.; validation, T.S.C., K.K. and A.J.T.; formal analysis, T.S.C. and K.K.; investigation, T.S.C. and K.K.; resources, I.R.R.; data curation, T.S.C., K.K. and A.J.T.; writing—original draft preparation, T.S.C., K.K. and A.J.T.; writing—review and editing, T.S.C., K.K., D.V., J.G.K., M.B., T.J. (Tamara Jurina), T.J. (Tamara Jakovljević), I.R.R. and A.J.T.; visualization, T.S.C. and A.J.T.; supervision, A.J.T.; project administration, I.R.R.; funding acquisition, I.R.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the European Union through the European Regional Development Fund, Competitiveness and Cohesion 2014–2020 (KK.01.1.1.07.0007).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The authors confirm that the data supporting the findings of this study are available within the article.

Acknowledgments

The authors warmly acknowledge Kutjevo d.d. (Kutjevo, Croatia) for providing the Graševina grape pomace biomass.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Waqas, M.; Hashim, S.; Humphries, U.W.; Ahmad, S.; Noor, R.; Shoaib, M.; Naseem, A.; Hlaing, P.T.; Lin, H.A. Composting Processes for Agricultural Waste Management: A Comprehensive Review. Processes 2023, 11, 731. [Google Scholar] [CrossRef]
  2. Debrah, J.K.; Vidal, D.G.; Dinis, M.A.P. Raising Awareness on Solid Waste Management through Formal Education for Sustainability: A Developing Countries Evidence Review. Recycling 2021, 6, 6. [Google Scholar] [CrossRef]
  3. Sokač, T.; Valinger, D.; Benković, M.; Jurina, T.; Kljusurić, J.G.; Redovniković, I.R.; Tušek, A.J. Application of Optimization and Modeling for the Enhancement of Composting Processes. Processes 2022, 10, 229. [Google Scholar] [CrossRef]
  4. Zohoori, M.; Ghani, A. Municipal Solid Waste Management Challenges and Problems for Cities in Low-Income and Developing Countries. Int. J. Sci. Eng. Appl. 2017, 6, 039–048. [Google Scholar] [CrossRef]
  5. Spoann, V.; Fujiwara, T.; Seng, B.; Lay, C. Municipal Solid Waste Management: Constraints and Opportunities to Improve Capacity of Local Government Authorities of Phnom Penh Capital. Waste Manag. Res. 2018, 36, 985–992. [Google Scholar] [CrossRef] [PubMed]
  6. Sayara, T.; Basheer-Salimia, R.; Hawamde, F.; Sánchez, A. Recycling of Organic Wastes through Composting: Process Performance and Compost Application in Agriculture. Agronomy 2020, 10, 1838. [Google Scholar] [CrossRef]
  7. Fan, S.; Li, A.; ter Heijne, A.; Buisman, C.J.N.; Chen, W.S. Heat Potential, Generation, Recovery and Utilization from Composting: A Review. Resour. Conserv. Recycl. 2021, 175, 105850. [Google Scholar] [CrossRef]
  8. Spinei, M.; Oroian, M. The Potential of Grape Pomace Varieties as a Dietary Source of Pectic Substances. Foods 2021, 10, 867. [Google Scholar] [CrossRef]
  9. Martínez Salgado, M.M.; Ortega Blu, R.; Janssens, M.; Fincheira, P. Grape Pomace Compost as a Source of Organic Matter: Evolution of Quality Parameters to Evaluate Maturity and Stability. J. Clean. Prod. 2019, 216, 56–63. [Google Scholar] [CrossRef]
  10. Muhlack, R.A.; Potumarthi, R.; Jeffery, D.W. Sustainable Wineries through Waste Valorisation: A Review of Grape Marc Utilisation for Value-Added Products. Waste Manag. 2018, 72, 99–118. [Google Scholar] [CrossRef] [PubMed]
  11. Azis, F.A.; Rijal, M.; Suhaimi, H.; Abas, P.E. Patent Landscape of Composting Technology: A Review. Inventions 2022, 7, 38. [Google Scholar] [CrossRef]
  12. Malakahmad, A.; Idrus, N.B.; Abualqumboz, M.S.; Yavari, S.; Kutty, S.R.M. In-Vessel Co-Composting of Yard Waste and Food Waste: An Approach for Sustainable Waste Management in Cameron Highlands, Malaysia. Int. J. Recycl. Org. Waste Agric. 2017, 6, 149–157. [Google Scholar] [CrossRef]
  13. Walker, L.; Charles, W.; Cord-Ruwisch, R. Comparison of Static, in-Vessel Composting of MSW with Thermophilic Anaerobic Digestion and Combinations of the Two Processes. Bioresour. Technol. 2009, 100, 3799–3807. [Google Scholar] [CrossRef]
  14. Papraćanin, E.A.; Petric, I. Mathematical Modeling and Simulation of the Composting Process in a Pilot Reactor. Bull. Chem. Technol. Bosnia Herzeg. 2017, 47, 39–48. [Google Scholar]
  15. Ajmal, M.; Aiping, S.; Uddin, S.; Awais, M.; Faheem, M.; Ye, L.; Rehman, K.U.; Ullah, M.S.; Shi, Y. A Review on Mathematical Modeling of In-Vessel Composting Process and Energy Balance. Biomass Convers. Biorefinery 2022, 12, 4201–4213. [Google Scholar] [CrossRef]
  16. He, X.; Han, L.; Ge, J.; Huang, G. Modelling for Reactor-Style Aerobic Composting Based on Coupling Theory of Mass-Heat-Momentum Transport and Contois Equation. Bioresour. Technol. 2018, 253, 165–174. [Google Scholar] [CrossRef] [PubMed]
  17. Loan, L.T.T.; Takahashi, Y.; Nomura, H.; Yabe, M. Modeling Home Composting Behavior toward Sustainable Municipal Organic Waste Management at the Source in Developing Countries. Resour. Conserv. Recycl. 2019, 140, 65–71. [Google Scholar] [CrossRef]
  18. Koolivand, A.; Godini, K.; Saeedi, R.; Abtahi, H.; Ghamari, F. Oily Sludge Biodegradation Using a New Two-Phase Composting Method: Kinetics Studies and Effect of Aeration Rate and Mode. Process Biochem. 2019, 79, 127–134. [Google Scholar] [CrossRef]
  19. Jain, M.S.; Paul, S.; Kalamdhad, A.S. Kinetics and Physics during Composting of Various Organic Wastes: Statistical Approach to Interpret Compost Application Feasibility. J. Clean. Prod. 2020, 255, 120324. [Google Scholar] [CrossRef]
  20. Zbytniewski, R.; Buszewski, B. Characterization of Natural Organic Matter (NOM) Derived from Sewage Sludge Compost. Part 2: Multivariate Techniques in the Study of Compost Maturation. Bioresour. Technol. 2005, 96, 479–484. [Google Scholar] [CrossRef]
  21. Karak, T.; Bhattacharyya, P.; Paul, R.K. Assessment of Co-Compost Quality by Physico-Chemical and Exploratory Data Analysis. Clean—Soil Air Water 2014, 42, 836–848. [Google Scholar] [CrossRef]
  22. Shen, G.; Chen, Y.; Zhang, J.; Wu, Y.; Yi, Y.; Li, S.; Yin, S. Quantitative Analysis of Index Factors in Agricultural Compost by Infrared Spectroscopy. Heliyon 2023, 9, e14010. [Google Scholar] [CrossRef]
  23. Chikae, M.; Ikeda, R.; Kerman, K.; Morita, Y.; Tamiya, E. Estimation of Maturity of Compost from Food Wastes and Agro-Residues by Multiple Regression Analysis. Bioresour. Technol. 2006, 97, 1979–1985. [Google Scholar] [CrossRef]
  24. Wan, X.; Li, J.; Xie, L.; Wei, Z.; Wu, J.; Wah Tong, Y.; Wang, X.; He, Y.; Zhang, J. Machine Learning Framework for Intelligent Prediction of Compost Maturity towards Automation of Food Waste Composting System. Bioresour. Technol. 2022, 365, 128107. [Google Scholar] [CrossRef]
  25. Bertran, E.; Sort, X.; Soliva, M.; Trillas, I. Composting Winery Waste: Sludges and Grape Stalks. Bioresour. Technol. 2004, 95, 203–208. [Google Scholar] [CrossRef]
  26. Pinto, R.; Correia, C.; Mourão, I.; Moura, L.; Brito, L.M. Composting Waste from the White Wine Industry. Sustainability 2023, 15, 3454. [Google Scholar] [CrossRef]
  27. Chan, M.T.; Selvam, A.; Wong, J.W.C. Reducing Nitrogen Loss and Salinity during “struvite” Food Waste Composting by Zeolite Amendment. Bioresour. Technol. 2016, 200, 838–844. [Google Scholar] [CrossRef] [PubMed]
  28. Waqas, M.; Nizami, A.S.; Aburiazaiza, A.S.; Barakat, M.A.; Ismail, I.M.I.; Rashid, M.I. Optimization of Food Waste Compost with the Use of Biochar. J. Environ. Manag. 2018, 216, 70–81. [Google Scholar] [CrossRef] [PubMed]
  29. Nasreen, Z.; Qazi, J.I. Lab Scale Composting of Fruits and Vegetable Waste at Elevated Temperature and Forced Aeration. Pak. J. Zool. 2012, 44, 1285–1290. [Google Scholar]
  30. Diaz, M.J.; Madejon, E.; Lopez, F.; Lopez, R.; Cabrera, F. Optimization of the Rate Vinasse / Grape Marc for Co-Composting Process. Process Biochem. 2002, 37, 1143–1150. [Google Scholar] [CrossRef]
  31. Lovreškov, L.; Radojčić Redovniković, I.; Limić, I.; Potočić, N.; Seletković, I.; Marušić, M.; Butorac, L.; Jurinjak Tušek, A.; Jakovljević, T. Are Foliar Nutrition Status and Indicators of Oxidative Stress Associated with Tree Defoliation of Four Mediterranean. Plants 2022, 11, 3484. [Google Scholar] [CrossRef] [PubMed]
  32. Jurinjak Tušek, A.; Jurina, T.; Benković, M.; Valinger, D.; Belščak-Cvitanović, A.; Gajdoš Kljusurić, J. Application of Multivariate Regression and Artificial Neural Network Modelling for Prediction of Physical and Chemical Properties of Medicinal Plants Aqueous Extracts. J. Appl. Res. Med. Aromat. Plants 2020, 16, 100229. [Google Scholar] [CrossRef]
  33. Fearn, T. Assessing Calibrations: SEP, RPD, RER and R2. NIR News 2002, 13, 12. [Google Scholar] [CrossRef]
  34. Friday Oyewusi, T.; Akinfemi Osunbitan, J.; Abidemi Ogunwande, G.; Akinola Omotosho, O. Investigation into Physico-Chemical Properties of Compost Extract as Affected by Processing Parameters Investigation into Physico-Chemical Properties of Compost Extract as Affected by Processing Parameters. Environ. Chall. 2021, 5, 100370. [Google Scholar] [CrossRef]
  35. Onwosi, C.O.; Igbokwe, V.C.; Odimba, J.N.; Eke, I.E.; Nwankwoala, M.O.; Iroh, I.N.; Ezeogu, L.I. Composting Technology in Waste Stabilization: On the Methods, Challenges and Future Prospects. J. Environ. Manag. 2017, 190, 140–157. [Google Scholar] [CrossRef]
  36. Peng, L.; Tang, R.; Wang, G.; Ma, R.; Li, Y.; Li, G.; Yuan, J. Effect of Aeration Rate, Aeration Pattern, and Turning Frequency on Maturity and Gaseous Emissions during Kitchen Waste Composting. Environ. Technol. Innov. 2023, 29, 102997. [Google Scholar] [CrossRef]
  37. Irvan; Husaini, T.; Trisakti, B.; Batubara, F.; Daimon, H. Composting of Empty Fruit Bunches in the Tower Composter—Effect of Air Intake Holes Composting of Empty Fruit Bunches in the Tower Composter—Effect of Air Intake Holes. IOP Conf. Ser. Mater. Sci. Eng. 2018, 309, 012066. [Google Scholar] [CrossRef]
  38. Karanja, A.W.; Njeru, E.M.; Maingi, J.M. Assessment of Physicochemical Changes during Composting Rice Straw with Chicken and Donkey Manure. Int. J. Recycl. Org. Waste Agric. 2019, 8, 65–72. [Google Scholar] [CrossRef]
  39. Makan, A.; Assobhei, O.; Mountadar, M. Effect of Initial Moisture Content on the In-Vessel Composting under Air Pressure of Organic Fraction of Municipal Solid Waste in Morocco. J. Environ. Health Sci. Eng. 2013, 10, 3. [Google Scholar] [CrossRef]
  40. Yeh, C.K.; Lin, C.; Shen, H.C.; Cheruiyot, N.K.; Camarillo, M.E.; Wang, C.L. Optimizing Food Waste Composting Parameters and Evaluating Heat Generation. Appl. Sci. 2020, 10, 2284. [Google Scholar] [CrossRef]
  41. Nguyen, T.P.; Koyama, M.; Nakasaki, K. Effects of Oxygen Supply Rate on Organic Matter Decomposition and Microbial Communities during Composting in a Controlled Lab-Scale Composting System. Waste Manag. 2022, 153, 275–282. [Google Scholar] [CrossRef]
  42. Mejias, L.; Komilis, D.; Gea, T.; Sánchez, A. The Effect of Airflow Rates and Aeration Mode on the Respiration Activity of Four Organic Wastes: Implications on the Composting Process. Waste Manag. 2017, 65, 22–28. [Google Scholar] [CrossRef]
  43. Topp, R.; Gomez, G. Residual Analysis in Linear Regression Models with an Interval-Censored Covariate. Stat. Med. 2004, 23, 3377–3391. [Google Scholar] [CrossRef] [PubMed]
  44. Panić, M.; Radović, M.; Cvjetko Bubalo, M.; Radošević, K.; Rogošić, M.; Coutinho, J.A.P.; Radojčić Redovniković, I.; Jurinjak Tušek, A. Prediction of PH Value of Aqueous Acidic and Basic Deep Eutectic Solvent Using COSMO-RS σ Profiles’ Molecular Descriptors. Molecules 2022, 27, 4489. [Google Scholar] [CrossRef] [PubMed]
  45. Palechor-Tróchez, J.J.; Ordoñez Santos, L.E.; Villada-Castillo, H.S. Relationship between Color CIELab and Total Organic Carbon in Compost. Adv. Mater. Sci. Eng. 2018, 2018, 5474239. [Google Scholar] [CrossRef]
  46. Zahrim, A.Y.; Leong, P.S.; Ayisah, S.R.; Janaun, J.; Chong, K.P.; Cooke, F.M.; Haywood, S.K. Composting Paper and Grass Clippings with Anaerobically Treated Palm Oil Mill Effluent. Int. J. Recycl. Org. Waste Agric. 2016, 5, 221–230. [Google Scholar] [CrossRef]
  47. Tredennick, A.T.; Hooker, G.; Ellner, S.P.; Adler, P.B. A Practical Guide to Selecting Models for Exploration, Inference, and Prediction in Ecology. Ecology 2021, 102, e03336. [Google Scholar] [CrossRef]
  48. Henseler, J.; Ringle, C.M.; Sinkovics, R.R. The Use of Partial Least Squares Path Modeling in International Marketing. Adv. Int. Mark. 2009, 20, 277–319. [Google Scholar] [CrossRef]
  49. Hair, J.F.; Ringle, C.M.; Sarstedt, M. Partial Least Squares Structural Equation Modeling: Rigorous Applications, Better Results and Higher Acceptance. Long. Range Plann 2013, 46, 1–12. [Google Scholar] [CrossRef]
  50. Hussain, S.; Fangwei, Z.; Siddiqi, A.F.; Ali, Z.; Shabbir, M.S. Structural Equation Model for Evaluating Factors Affecting Quality of Social Infrastructure Projects. Sustainability 2018, 10, 1415. [Google Scholar] [CrossRef]
  51. Chang, C.-W.; Laird, D.A.; Mausbach, M.J.; Hurburgh, C.R. Near—Infrared Reflectance Spectroscopy–Principal Components Regression Analyses of Soil Properties. Soil. Sci. Soc. Am. J. 2001, 65, 480–490. [Google Scholar] [CrossRef]
  52. Sim, J.; McGoverin, C.; Oey, I.; Frew, R.; Kebede, B. Near-Infrared Reflectance Spectroscopy Accurately Predicted Isotope and Elemental Compositions for Origin Traceability of Coffee. Food Chem. 2023, 427, 136695. [Google Scholar] [CrossRef] [PubMed]
  53. Fouguira, S.; El Haji, M.; Benhra, J.; Ammar, E. Optimization of Olive Oil Extraction Wastes Co Composting Procedure Based on Bioprocessing Parameters. Heliyon 2023, 9, e19645. [Google Scholar] [CrossRef]
  54. Yang, L.; Liu, S.; Tsoka, S.; Papageorgiou, L.G. Mathematical Programming for Piecewise Linear Regression Analysis. Expert. Syst. Appl. 2016, 44, 156–167. [Google Scholar] [CrossRef]
  55. Costello, R.C.; Sullivan, D.M.; Bryla, D.R.; Strik, B.C.; Owen, J.S. Compost Feedstock and Compost Acidification Affect Growth and Mineral Nutrition in Northern Highbush Blueberry. HortScience 2019, 54, 1067–1076. [Google Scholar] [CrossRef]
  56. Abdi, R.; Shahgholi, G.; Sharabiani, V.R.; Fanaei, A.R.; Szymanek, M. Prediction Compost Criteria of Organic Wastes with Biochar Additive in In-Vessel Composting Machine Using ANFIS and ANN Methods. Energy Rep. 2023, 9, 1684–1695. [Google Scholar] [CrossRef]
  57. Chen, X.; Zheng, H.; Wang, H.; Yan, T. Can Machine Learning Algorithms Perform Better than Multiple Linear Regression in Predicting Nitrogen Excretion from Lactating Dairy Cows. Sci. Rep. 2022, 12, 12478. [Google Scholar] [CrossRef]
  58. Hosseinzadeh, A.; Baziar, M.; Alidadi, H.; Zhou, J.L.; Altaee, A.; Asghar, A.; Jafarpour, S. Bioresource Technology Application of Artificial Neural Network and Multiple Linear Regression in Modeling Nutrient Recovery in Vermicompost under Different Conditions. Bioresour. Technol. 2020, 303, 122926. [Google Scholar] [CrossRef]
  59. Aycan Dümenci, N.; Cagcag Yolcu, O.; Aydin Temel, F.; Turan, N.G. Identifying the Maturity of Co-Compost of Olive Mill Waste and Natural Mineral Materials: Modelling via ANN and Multi-Objective Optimization. Bioresour. Technol. 2021, 338, 125516. [Google Scholar] [CrossRef]
  60. Singh, T.; Uppaluri, R.V.S. Application of ANN and Traditional ML Algorithms in Modelling Compost Production under Different Climatic Conditions. Neural Comput. Appl. 2023, 35, 13465–13484. [Google Scholar] [CrossRef]
  61. Shi, C.F.; Yang, H.T.; Chen, T.T.; Guo, L.P.; Leng, X.Y.; Deng, P.B.; Bi, J.; Pan, J.G.; Wang, Y.M. Artificial Neural Network-Genetic Algorithm-Based Optimization of Aerobic Composting Process Parameters of Ganoderma Lucidum Residue. Bioresour. Technol. 2022, 357, 127248. [Google Scholar] [CrossRef]
  62. Gurusamy, N.N.; Puffer, N.; Jongh, C.D.; Rodriguez, C.; Aspray, T.J. Effect of Initial Moisture Content and Sample Storage Duration on Compost Stability Using the ORG0020 Dynamic Respiration Test. Waste Manag. 2021, 125, 215–219. [Google Scholar] [CrossRef]
  63. Ghanney, P.; Qiu, H.; Anning, D.K.; Yang, H.; Wang, Y.; Kugbe, J.X. Applied Sciences Moisture-Induced Pattern of Gases and Physicochemical Indices in Corn Straw and Cow Manure Composting. Appl. Sci. 2021, 11, 8493. [Google Scholar] [CrossRef]
  64. Abdallah, S.E.; Mazrou, Y.S.A.; Elsakhawy, T.; Elgarhy, R.; Elmetwalli, A.H.; Elsayed, S.; Elmessery, W.M. Optimizing the In-Vessel Composting Process of Sugarbeet Dry-Cleaning Residue. Agriculture 2022, 12, 427. [Google Scholar] [CrossRef]
  65. Antızar-Ladislao, B.; Lopez-Real, J.; Beck, A.J. Investigation of Organic Matter Dynamics during In-Vessel Composting of an Aged Coal—Tar Contaminated Soil Using Fluorescence Excitation—Emission Spectroscopy. Chemosphere 2006, 64, 839–847. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Comparisons between experimental data and MLR models predicted data of physicochemical properties of compost during the composting process. (a) moisture content, (b) dry-matter content, (c) organic-matter content, (d) ash content, (e) carbon content, (f) nitrogen content, (g) carbon/nitrogen ratio, (h) compost total colour change, (i) pH, (j) total dissolved solids, (k) conductivity, (l) compost extracts total colour change. (○) calibration data set and () prediction data sets. Black line represents MLP model.
Figure 1. Comparisons between experimental data and MLR models predicted data of physicochemical properties of compost during the composting process. (a) moisture content, (b) dry-matter content, (c) organic-matter content, (d) ash content, (e) carbon content, (f) nitrogen content, (g) carbon/nitrogen ratio, (h) compost total colour change, (i) pH, (j) total dissolved solids, (k) conductivity, (l) compost extracts total colour change. (○) calibration data set and () prediction data sets. Black line represents MLP model.
Bioengineering 11 00285 g001
Figure 2. Comparisons between experimental data and PLR models predicted data of physicochemical properties of compost during the composting process. (a) Moisture content, (b) dry-matter content, (c) organic-matter content, (d) ash content, (e) carbon content, (f) nitrogen content, (g) carbon/nitrogen ratio, (h) compost total colour change, (i) pH, (j) total dissolved solids, (k) conductivity, (l) compost extracts total colour change. (○) calibration data set and () prediction data sets.
Figure 2. Comparisons between experimental data and PLR models predicted data of physicochemical properties of compost during the composting process. (a) Moisture content, (b) dry-matter content, (c) organic-matter content, (d) ash content, (e) carbon content, (f) nitrogen content, (g) carbon/nitrogen ratio, (h) compost total colour change, (i) pH, (j) total dissolved solids, (k) conductivity, (l) compost extracts total colour change. (○) calibration data set and () prediction data sets.
Bioengineering 11 00285 g002
Figure 3. Comparisons between experimental data and ANN models predicted data of physicochemical properties of compost during the composting process. (a) Moisture content, (b) dry-matter content, (c) organic-matter content, (d) ash content, (e) carbon content, (f) nitrogen content, (g) carbon/nitrogen ratio, (h) compost total colour change, (i) pH, (j) total dissolved solids, (k) conductivity, (l) compost extracts total colour change. (○) training, () test, () validation and () prediction data sets for ANN models.
Figure 3. Comparisons between experimental data and ANN models predicted data of physicochemical properties of compost during the composting process. (a) Moisture content, (b) dry-matter content, (c) organic-matter content, (d) ash content, (e) carbon content, (f) nitrogen content, (g) carbon/nitrogen ratio, (h) compost total colour change, (i) pH, (j) total dissolved solids, (k) conductivity, (l) compost extracts total colour change. (○) training, () test, () validation and () prediction data sets for ANN models.
Bioengineering 11 00285 g003
Figure 4. Global sensitivity analysis of important variables of the composting process (MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon/nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples).
Figure 4. Global sensitivity analysis of important variables of the composting process (MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon/nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples).
Bioengineering 11 00285 g004
Table 1. Minimum, maximum and standard deviations of analysed physicochemical properties (MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon/nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples).
Table 1. Minimum, maximum and standard deviations of analysed physicochemical properties (MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon/nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples).
MC (%)DM (%)OM (%)AC (%)CC (%)NC (%)C/N∆EcpHTDS (mg/L)S (µS/cm)∆Ee
Exp1Min50.14441.57566.54816.52547.5001.17021.5492.0744.520491.0001002.0000.374
Max58.42549.85683.47533.45252.4002.33040.94018.3088.0801551.0003120.0004.994
St.dev.2.1432.1434.6704.6701.3040.3455.6374.6451.064389.335723.9740.875
Exp2Min48.85539.18766.78714.18147.3001.16021.1201.0194.490286.000568.0000.242
Max60.81351.14585.81933.21352.8002.50041.72419.8997.7901700.0003400.0002.922
St.dev.2.3322.3324.7534.7531.4580.3545.6734.5901.039477.961933.3590.734
Exp3Min55.44029.20157.18314.64246.4001.07019.8760.8555.130387.000799.0000.179
Max70.79944.56085.35842.81752.6002.46045.98118.5189.2501442.0002940.0004.505
St.dev.3.8933.8938.4308.4301.7180.3696.0464.7151.167375.912757.7411.262
Exp4Min56.42330.49963.18614.31547.0001.07020.6641.4664.770337.000674.0000.192
Max69.50143.57785.68536.81451.8002.41045.98114.0028.2201752.0003520.0002.378
St.dev.2.5912.5915.6945.4801.0900.3616.2533.3961.058464.734914.8620.710
Exp5Min56.42330.19258.80614.31546.7001.07020.3771.4664.660313.000627.0000.192
Max69.80841.20382.37241.19451.7002.39045.98115.0188.2701413.0002880.0002.588
St.dev.2.4762.4765.5105.3891.3060.3505.7633.5531.124417.009820.5010.706
Exp6Min49.75534.02762.07624.54143.9001.08014.8732.4445.7001044.0001991.0000.554
Max65.97350.24575.45937.92448.9003.16044.1679.3038.6002180.0004470.0005.294
St.dev.4.1354.1353.6713.6711.0440.5197.5741.4451.029380.752822.3071.390
Exp7Min56.98926.92657.97023.52644.2001.03014.6981.5275.360835.0001666.0000.167
Max73.07443.01176.47442.03049.1003.24045.53410.4629.2601867.0003930.0004.373
St.dev.3.65534.0275.2515.2511.1030.5848.1081.9671.382318.789896.3931.405
Exp8Min52.5663.65560.91220.47945.9001.08019.3501.9205.620926.0001833.0000.244
Max70.05829.94279.52139.08848.9002.47043.36411.4188.5201598.0003180.0004.022
St.dev.3.02236.8233.9033.9030.7730.3856.5341.7960.752208.822424.2831.217
Exp9Min58.6783.02261.23824.26245.6001.10015.5051.3195.230802.0001567.0000.379
Max70.74029.26075.73838.76249.3003.07043.3649.5828.6802111.0004100.0004.165
St.dev.2.56434.4603.6413.6410.8160.5177.1951.5521.166374.269720.6541.055
Table 2. Spearman’s correlation matrix for the determination of the relationships between conditions of the composting process and physicochemical properties of the compost during the process (IMC—initial moisture content, AFR—air flow rate, SD—sampling day, MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon-to-nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples). Significant correlations are marked bold.
Table 2. Spearman’s correlation matrix for the determination of the relationships between conditions of the composting process and physicochemical properties of the compost during the process (IMC—initial moisture content, AFR—air flow rate, SD—sampling day, MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon-to-nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples). Significant correlations are marked bold.
IMCAFRSDMCDMOMACCCNCC/NΔEcpHTDSSΔEe
IMC1.000
AFR−0.2341.000
SD0.0000.0001.000
MC0.683−0.0760.3171.00
DM−0.6830.076−0.317−1.001.00
OM−0.220−0.130−0.537−0.660.661.000
AC0.2220.1310.5360.66−0.66−1.0001.000
CC−0.006−0.224−0.255−0.280.280.451−0.4501.000
NC−0.0290.1530.3030.22−0.22−0.2700.2700.0451.000
C/N0.023−0.173−0.353−0.250.250.331−0.3310.071−0.9901.000
ΔEc0.150−0.2780.4320.07−0.07−0.1530.151−0.019−0.1670.1461.000
pH0.1830.1060.5850.49−0.49−0.5950.594−0.4620.326−0.4000.2031.000
TDS−0.1670.2840.5500.14−0.14−0.4460.445−0.3590.217−0.2750.1460.5121.000
S−0.1510.2530.5540.15−0.15−0.4380.437−0.3420.194−0.2520.1790.5030.9851.000
ΔEe−0.0590.2130.4020.23−0.23−0.4070.405−0.4070.024−0.0790.0580.3370.4090.4151.000
Table 3. The artificial neural network models for prediction of physicochemical properties of grape-skin compost (MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon to nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples).
Table 3. The artificial neural network models for prediction of physicochemical properties of grape-skin compost (MC—moisture content, DM—dry-matter content, OM—organic-matter content, AC—ash content, CC—carbon content, NC—nitrogen content, C/N—carbon to nitrogen ratio, ∆Ec—total colour change of compost samples, TDS—total dissolved solids, S—conductivity, ∆Ee—total colour change of compost extract samples).
CalibrationPrediction
Output Network Train. Perf. Train. ErrorTest Perf. Test ErrorValid. Perf. Valid. ErorHidden ActivationOutput ActivationRpred2Rpred2adjRMSEPSEPRPDRER
MCMLP 3-6-10.9262
1.1201
0.9248
1.1335
0.9137
1.5458
TanhLogistic0.90500.90281.70780.15043.172611.8455
DMMLP 3-10-10.9248
1.1222
0.9117
1.1258
0.9107
1.4933
TanhLogistic0.90380.90151.71730.15123.155011.7799
OMMLP 3-9-10.7978
2.2311
0.7677
2.4293
0.7594
2.8152
TanhTanh0.75310.74723.10220.27312.00478.9766
ACMLP 3-10-10.8359
1.6794
0.8344
1.8354
0.8057
2.0836
TanhTanh0.72330.71673.30300.29081.85878.4310
CCMLP 3-8-10.8623
0.5664
0.8613
0.6265
0.8522
0.7596
LogisticIdentity0.66580.65710.97700.08601.63808.3933
NCMLP 3-4-10.7556
0.0605
0.7554
0.0625
0.7236
0.0639
LogisticExponential0.65160.64330.27930.02461.68197.9125
C/NMLP 3-10-10.7367
0.2016
0.6826
0.2683
0.6672
0.3913
TanhExponential0.65420.64833.60980.31781.98958.5399
∆EcMLP 3-5-10.9219
1.3085
0.9177
1.3384
0.9056
1.3718
LogisticExponential0.73440.72812.08420.18351.99438.1286
pHMLP 3-10-10.9137
0.1148
0.8961
0.1569
0.8674
0.1640
TanhIdentity0.83220.82820.48960.04312.38909.6603
TDSMLP 3-8-10.9231
83.3771
0.8961
134.5269
0.8334
163.2544
TanhTanh0.71510.7084236.951420.86241.84187.6682
SMLP 3-10-10.8670
87.4412
0.8965
151.1360
0.8076
182.1221
TanhLogistic0.75290.7447522.993446.04701.97638.5246
∆EeMLP 3-7-10.8055
0.3691
0.7936
0.6189
0.7853
0.6371
LogisticLogistic0.69520.68560.87460.07701.95317.8557
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sokač Cvetnić, T.; Krog, K.; Valinger, D.; Gajdoš Kljusurić, J.; Benković, M.; Jurina, T.; Jakovljević, T.; Radojčić Redovniković, I.; Jurinjak Tušek, A. Application of Multivariate Regression and Artificial Neural Network Modelling for Prediction of Physicochemical Properties of Grape-Skin Compost. Bioengineering 2024, 11, 285. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering11030285

AMA Style

Sokač Cvetnić T, Krog K, Valinger D, Gajdoš Kljusurić J, Benković M, Jurina T, Jakovljević T, Radojčić Redovniković I, Jurinjak Tušek A. Application of Multivariate Regression and Artificial Neural Network Modelling for Prediction of Physicochemical Properties of Grape-Skin Compost. Bioengineering. 2024; 11(3):285. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering11030285

Chicago/Turabian Style

Sokač Cvetnić, Tea, Korina Krog, Davor Valinger, Jasenka Gajdoš Kljusurić, Maja Benković, Tamara Jurina, Tamara Jakovljević, Ivana Radojčić Redovniković, and Ana Jurinjak Tušek. 2024. "Application of Multivariate Regression and Artificial Neural Network Modelling for Prediction of Physicochemical Properties of Grape-Skin Compost" Bioengineering 11, no. 3: 285. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering11030285

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop