Issue 
Int. J. Metrol. Qual. Eng.
Volume 12, 2021



Article Number  12  
Number of page(s)  11  
DOI  https://doi.org/10.1051/ijmqe/2021010  
Published online  28 May 2021 
Research Article
Behaviour desorption study of the essential oil of Cedrus atlantica in a porous clay versus insecticidal activity against Sitophilus granarius: explanation of the phenomenon by statistical studies
^{1}
Superior School of Technology, University of Sultan Moulay Slimane, PB 170, Khenifra 54000, Morocco
^{2}
Medicinal Research Institute, Center for Research and Study of Djibouti, BP 486, Djibouti
^{*} Corresponding author: t.ainane@usms.ma
Received:
14
February
2021
Accepted:
15
May
2021
Substances of natural origin and more particularly essential oils currently represent an alternative solution in the fight for the protection of stored foodstuffs. The way to use these essential oils and the storage conditions remain the main handicap in dealing with this subject. This article develops the use of porous supports as a medium for the study of the behavior of essential oils with respect to insecticidal activities. The process relates to the fixing and/or grafting of essential oils in porous clay media in a welldefined geometric form included in the storage methods. The study of transfers of the essential oil of Cedrus atlantica in a porous clay medium in the case of desorption was made by analytical and numerical models of diffusion process, to know the behavior of the oil and to determine some physical parameters (diffusivity D_{z}, activation energy E_{a}, evaporation rate F and constant of evaporation K) which explain the mechanisms involved and to try to exploit them in parallel with the insecticidal activities against of Sitophilus granarius the main cereal pest by statistical approaches such as: design of experiments and principal component analysis.
Key words: Clay / essential oil / behaviour desorption / diffusion / insecticidal activity
© A. Ainane et al., Published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Cereal storage is designed to keep grains for later use from many external factors, such as humidity, temperature, and light, that can lead to their degradation [1]. In addition, good storage also protects its content from other influences such as: microorganisms, shocks, odours, vibrations, dust and compressive forces [2]. The Technological development has brought the possibility to access to countless foods from every part of the globe, so the importance of protecting their organoleptic properties is food quality conformity [3–5].
More precisely, stored cereals are prone to postharvest loss in quality and quantity due to infestation by different groups of insects, which is considered the principal problem of storage [6]. Among the stored product beetles, Sitophilus granaries (Coleoptera: Dryophthorinae) are a primary pest in storage of grainbased products, principally wheat, but also attack others cereals: oats, rye, barley, corn, as well as derived products (flour) and other seeds, especially Fabaceae [7,8]. Control of these insects relies heavily on the treatment by synthetic insecticides such as pyrethroids, organochlorines, organophosphates, carbamates and fumigants (mainly methyl bromide and phosphine) [9]. However, due to toxicity to consumers, insecticide resistance and resurgence of pests associated with synthetic insecticides, alternative solutions are requiring sought [10]. Pesticides from natural sources, relatively cheaper, biodegradable, which are means accessible and available, less toxic to nontarget organisms and less prone to resistance by insect species are considered potential resources [11]. Among these resources, Essential oils produced from aromatic and medicinal plants have received much work due to their broad spectrum of activities, they are volatile, complex compounds characterized by a strong odour and are formed by aromatic plants as secondary metabolites (monoterpenes, sesquiterpenes and phenylpropanoid compounds) which show different insecticide properties (larvicidal, antifeedant, repellent, fumigant and ovicidal activities) [12]. Due forward described properties, essential oils are potentially suitable for integrated pest management programmes [13].
Cedar extracts considered an effective natural pesticide, especially their essential oil. This last obtained from different geographically regions, such as Morocco, Algeria, Cyprus, Lebanon, Syria, Turkey, Afghanistan and Himalayas, have been reported with their use in entomology to obtain many treatments traditionally ascribed to insecticidal activities [14]. In terms of previous reports on the chemical composition of this plant, the oil from the wood of this species has been extensively studied, sesquiterpenoids of the himachalane compounds being the major components [15].
The objective of this work is to carry out a desorption transfer study of the essential oil of Cedrus atlantica after fixation in a porous clay, by the use of experimental designs (complete factorial design) and analytical models of diffusion process (2nd Fick's Law), to know the process of essential oil and to determine the relation of certain physical parameters between them (diffusivity, constant of evaporation and rate of evaporation) which can explain the mechanism in play with the insecticidal activity against Sitophilus granaries.
2 Materials and methods
2.1 Materials
2.1.1 Clay
The studied natural porous clay RC was selected from Bejaad, Morocco (32°47'01.7”N 6°13'52.3”W). The chemical composition of the clay raw material (Tab. 1) shows that the major elements are SiO_{2}, Fe_{2}O_{3} and Al_{2}O_{3}. The amounts of K_{2}O, Na_{2}O and MgO are significant. From XRD analyses (Fig. 1), it is found that this clay is mainly constituted of quartz, calcite, dolomite, kaolinite, illite and hematite [16].
Chemical composition (mass %) of the clay.
Fig. 1 XRD diffractograms of the clay RC. 
2.1.2 Essential oil of cedrus atlantica
The essential oil of Cedrus atlantica (EOCa), which was obtained by steam distillation from the sawdust of Cedrus atlantica Man. from Morocco, was analysed by chromatography–mass spectrometry (GC–MS). The GC analysis was performed on SHIMADZU GC14B equipped with and FID detector and a LM5 (30 m × 0.25 mm × 0.3 mm) capillary column. The components in the oil were identified by comparison of the mass spectra obtained by CG–MS and literature data. The results obtained by the analysis are demonstrated in Table 2. In total were identified 22 constituents (92.10%).
The data obtained with the chemical analysis are in accord with previous works reported in literature of the compounds analysed, the sesquiterpenes hydrocarbons αhimachalene (15.63%), βhimachalene (31.24%) and γhimachalene (14.46%) are the majority compounds, in which the component βhimachalene is the most abundant and represent practically the twothirds of a percentage of the oil composition [17].
Chemical constitution of the EOCa.
2.1.3 Preparation of the mixture clay/essential oil
The RC clay was homogenized, finely ground (<63 µm) and heated for 2 h at 200 °C for the organic complete combustion. a welldefined volume of EOCa has been mixed with a quantity of clay powder, prepared between 10 and 40% (v/w) in order to have a total fixation of the oil on the powder. Then, the prepared materials (RC+EOCa) are moved to metal cylinders (diameters 1 cm and 2 cm), for the continuation of the experimental works.
2.2 Theoretical considerations
2.2.1 Assumptions
Material transfer of EOCa is considered with onedimensional diffusion in the cylindrical coordinate system.
The diffusion takes place under transient conditions with a constant diffusivity and the rate of evaporation.
During sorption, the concentration of EOCa on clay RC face reaches the equilibrium value as soon as the EOCa is mixed with RC.
At the beginning of desorption, the concentration of EOCa throughout the clay is constant.
2.2.2 Model of transfer of matter by diffusion
Fick's second law for diffusion in terms of cylindrical coordinates (r, θ, z), can be written as [18,19]:$$\frac{\partial C}{\partial t}=\frac{1}{r}\left\{\frac{\partial}{\partial r}\left(r{D}_{r}\frac{\partial C}{\partial r}\right)+\frac{\partial}{\partial \theta}\left(\frac{{D}_{\theta}}{r}\frac{\partial C}{\partial \theta}\right)+\frac{\partial}{\partial z}\left(r{D}_{z}\frac{\partial C}{\partial z}\right)\right\}.$$(1)
Since, the cylinder height is smaller in front of the diameter and the lateral cylinder surface, the diffusion transfer will be directed in the direction of the height, which allows to consider the transfer model is onedimensional. The problem is made easier and the diffusion equation becomes Fick's 2nd Law:$$\frac{\partial C}{\partial t}={D}_{z}\frac{{\partial}^{2}C}{\partial {z}^{2}}.$$(2)
Considering the following boundary conditions:$$\{\begin{array}{ccc}\hfill t=0\hfill & \hfill C={C}_{0}\hfill & \hfill 0\le z\le h\hfill \\ \hfill t>0\hfill & \hfill C={C}_{t}\hfill & \hfill z=h.\hfill \end{array}$$(3)
The analytical solution of equation (2) has been demonstrated by Crank (1979) [18] according to the conditions displayed in equation (3), it can be written as follows:$$C\left(z,t\right)=\frac{4}{\pi}{\displaystyle \sum _{n=0}^{\infty}}\mathrm{e}\mathrm{x}\mathrm{p}\left({D}_{z}\cdot \frac{{\left(2n+1\right)}^{2}}{{h}^{2}}\cdot {\pi}^{2}\cdot t\right)\cdot \mathrm{c}\mathrm{o}\mathrm{s}\left(\frac{\left(2n+1\right)\pi}{2h}\cdot z\right).$$(4)
The total mass M_{t} of EOCa in the clay RC at an instant t is obtained by integration of the variable C (concentration) on the thickness of the material, and on a surface S subjected to the concentration flow:$${M}_{t}={\displaystyle \underset{0}{\overset{h}{\int}}}S.C\left(z,t\right)dz.$$(5)
The analytical solution of these equations taking into account the boundary conditions and after an infinite desorption time tends towards equilibrium is given by:$$\frac{{M}_{t}{M}_{\infty}}{{M}_{0}{M}_{\infty}}=\frac{8}{{\pi}^{2}}{\displaystyle \sum _{n=0}^{\infty}}\frac{1}{{\left(2n+1\right)}^{2}}\mathrm{e}\mathrm{x}\mathrm{p}\left({D}_{z}\cdot \frac{{\left(2n+1\right)}^{2}}{{h}^{2}}\cdot {\pi}^{2}\cdot t\right)$$(6)
with:
M_{0}: The initial mass of EOCa absorbed.
M_{t}: The mass of EOCa desorbed after a given time of aging.
M_{∞}: the mass of liquid at equilibrium.
The simplifying the equation (6) by taking the first term of the series solution and assuming M_{∞} = 0:$$\frac{{M}_{t}}{{M}_{0}}=\frac{8}{{\pi}^{2}}\mathrm{e}\mathrm{x}\mathrm{p}\left(\frac{{D}_{z}\cdot {\pi}^{2}\cdot t}{{h}^{2}}\right).$$(7)
The term (8/π^{2}) is assumed to be equal to 1. The last equation becomes:$$\frac{{M}_{t}}{{M}_{0}}=\mathrm{e}\mathrm{x}\mathrm{p}\left(\frac{{D}_{z}\cdot {\pi}^{2}\cdot t}{{h}^{2}}\right).$$(8)
During our study of the modelling and in consideration of the theoretical approximations already cited, we will focus on the three parameters which explain the desorption process: Diffusivity D_{z}, Evaporation constant K and evaporation rate F:
The diffusivity does practically obtaining by the linear, logarithmic extrapolation of equation (8).
the constant of evaporation K approximated to the diffusivity by the following equation:
The evaporation rate F was determined using the initial value of the desorption flux, it is the value of the initial gradient of the desorption curve as a function of time, as indicated in the equation (10):
2.3 Experimentation and modelling
2.3.1 Insecticidal activity test
An insecticidesusceptible population of the adult weevil species (Sitophilus granarius) were used in this study. The population of Sitophilus granarius was obtained from the Department of environmental engineering of ESTKhenifra (University of Sultan Moulay Slimane, Khenifra, Morocco). The strain was reared on wheat grains free of insecticide residue in glass containers (1 L) within growth chambers at 27 ± 2 °C, 70 ± 10% of relative humidity (RH), and 12:12 h photoperiod (D: L) [20].
EOCa or EOCa + RC were placed in steel cylinders (h = 0.5 cm; Φ = 1 cm or 2 cm or 3 cm), then the cylinders are included in glass petri dishes containing 10 insects. The assembly with a negative control (without product) are introduced into a fumigation chamber included in the experimental enclosure (temperature and relative humidity controlled). The mortality of insects is recorded as a function of time for 24 h. Repetitions were performed in triplicate for each trial to minimize errors.
Corrected mortality in treating insects is expressed by equation (11):$$M\%=\left(\frac{{M}_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}}{M}_{\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{t}\mathrm{r}\mathrm{o}\mathrm{l}}}{100{M}_{\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{t}\mathrm{r}\mathrm{o}\mathrm{l}}}\right)\times 100$$(11)
with:
M%: Mortality corrected;
M_{test}: Mortality observed during the test;
M_{control}: mortality observed in the control.
The determination of the lethal dose of 50% LD50 is determined by linear interpolation on curves giving the percentage of mortality as a function of the logarithm of the concentration tested.
2.3.2 Optimization of transfer conditions
In order to minimize the experimental conditions of diffusion transfer of EOCa essential oil in RC porous clay media, we performed a statistical approach based on experimental designs.
The transfer was presented as the evaporation rates F [21]. The design matrix (Full Factorial Design) was performed according to 4 factors X_{i} such as:
Factor 1 = C: Essential oil concentration (0.01 μL/cm^{3} and 0.02 μL/cm^{3});
Factor 2 = T: Temperature (25 and 30 °C);
Factor 3 = D: Cylinder diameter (1 cm and 2 cm);
Factor 4 = M: Mass of clay rock (0.05 g and 0.10 g).
The number of tests carried out is calculated by the formula:$$\text{number}\text{\hspace{0.17em}}\text{of}\text{\hspace{0.17em}}\text{tests}={2}^{k}.$$(12)
The number 2 corresponds to the number of two levels −1 and +1 and k stands for the number of factors studied.
In our study, the number of factors equals 4:
number of tests = 2^{4} = 16
This corresponds to 16 tests.
Table 3 presents the matrix of the experiments and the evaporation rates obtained, from which the rows present the tests coded for the factors studied according to the high (+1) or low (−1) level of each factor.
The polynomial model presents:$$F={b}_{0}+{\displaystyle \sum _{i=1}^{n}}{b}_{i}{X}_{i}+{\displaystyle \sum _{i=1}^{n}}{\displaystyle \sum _{j=1}^{n1}}{b}_{ij}{X}_{i}{X}_{j}+{\displaystyle \sum _{i=1}^{n}}{\displaystyle \sum _{j=1}^{n1}}{\displaystyle \sum _{k=1}^{n2}}{b}_{ijk}{X}_{i}{X}_{j}{X}_{k}+{b}_{ijkl}{X}_{i}{X}_{j}{X}_{k}{X}_{l}$$(13)
Or:
Mean: b_{0};
4 main effects: b_{i};
6 effects present the 2nd order interactions: b_{ij};
4 effects present the 3rd order interactions: b_{ijk};
1 interaction of the 4th order: b_{ijkl}.
Matrix of experiments and evaporations rate (flux) F obtained for each test.
2.3.3 Study of insecticidal activity behaviour of essential oil mixture in clay porous media
The interaction between the behaviour of EOCa in porous clay media RC and insecticidal activity was achieved by processing the matter transfer by diffusion parameters vs. mortality of insecticidal activity, according to the statistical tool of Principal Component Analysis (PCA) which is an extremely powerful tool for compressing and synthesizing information, useful when there is a sum of quantitative data to be processed, to be interpreted and to know the correlations between them [22].
2.3.4 Computer and statistical processing
All computer processing was done using MATLAB software. The numerical modelling results obtained will be compared with the experimental results, in order to find the parameters to be searched for and to make the maximum interpretation. Other computer tools have been used in parallel to perform certain calculations or graphical representation, such as: EXCEL, ORIGIN LAB, DESIGN PLAN and XLSTAT.
3 Results and analysis
3.1 Insecticidal activities in porous clay media
All the results of the insecticidal tests of the EOCa alone or with the porous clay media RC and GC, against Sitophilus granarius are displayed in Table 4 in the form of lethal doses of 50% (LD_{50}) depending on the parameters studied such as: the diameter of cylinder D and the temperature of the incubation T. From this table, we deduce that the LD_{50} decrease in the presence of porous clayey RC and GC media, which proves that the insecticidal power has changed when the essential oil has been fixed on the porous media. In addition, the studied factors such as cylinder diameter and temperature keep the same properties of insecticidal activities, hence increasing cylinder diameter and increasing temperature promotes good activity. Also, it can be concluded that the insecticidal activity of the EOCa fixed on porous media and their persistence depends on the nature of the interactions between the constituent compounds of essential oils and the media studied [23].
LD_{50} of the insecticidal activities of the EOCa fixed on porous media.
3.2 Optimization of essential oil transfer conditions in porous media
As we described previously, the experimental design chosen for the optimization of the flux parameter of the evaporation of the EOCa fixed in the porous medium is a full factorial design of four (4) factors: the concentration of l essential oil C, temperature T, cylinder diameter D and mass of clay medium M. Each factor requires 2 levels: high level is coded +1 and low level is coded −1, which implies a total number of 16 simulations. Therefore, the mathematical model associated with the 1st degree polynomial is presented as a function of the factors plus the second, third and fourth order interactions of these factors.
The coefficients of the flux model of the essential oil fixed on the RC medium are displayed in Table 5. Following the results obtained by the optimization route, the mean flux is 2.6061.10^{−4}. If we take into account only the factors studied on the medium, we deduce that:
The increasing concentration of the essential oil increases the flux;
The increase in temperature increases the flow;
The increase in cylinder diameter decreases the flow;
The growth of the mass of porous medium decreases the flux.
The set of graphical representations of flux as a function of two factors in 3D are shown in Figure 2. The interaction of two of the four factors also influences the flux values. The flux increases with respect to the interaction of two concentrationtemperature parameters, it decreases with respect to a parameter with the interactions of concentrationdiameter, concentrationmass, masstemperature and temperaturediameter and it decreases with respect to the interaction from two parameters to diametermass interactions. All the graphical representations of flow as a function of two factors in 3D are shown in Figure 2.
In general, and according to the results of Tables 6 and 7, the coefficients of determination R^{2} for the two media show the correct explanation of the optimized models. The Fisher test value is F_{test} = 2.42 which explains why the model is not significant with respect to noise (uncertainties). The predictive coefficient of determination R_{pred}^{2} gives a negative value which implies that the total mean of flux presented can be a better predictor of the response of the model studied [24]. The variation ratio of the response measures the flux/uncertainty ratio, is equal to 6.68 which confirms that the response is adequate (ratio greater than 4 which is desirable).
Parameters of the presented polynomial flux model for desorption.
Fig. 2 3D representation of flow as a function of two factors. 
Statistical study of the presented polynomial model of flux during desorption.
Descriptive and predictive qualities of polynomial models.
3.3 Modelling of the essential oil diffusion process in porous media
To fully understand the behavior and the diffusion process of EOCa in the RC medium, two intermediate paths were made: the study of the desorption kinetics and the simulation of the concentration profiles inside the cylinder. This study makes it possible to explain the mechanism involved, as well as a means of precisely determining the parameters: the diffusion coefficient D_{z}, the evaporation flux F, the evaporation constant K and the activation energy E_{a}.
Several kinetic models cited in the literature can study the desorption mechanism [25–27], in our case the model tested is the analytical solution of the Fickian diffusion model, which is described previously. The methodology adopted consists in controlling the mass percentage of the oil as a function of the time of the experimental results with simulations of analytical treatments during desorption. To make the simulation according to this model, it is necessary to seek the coefficient of diffusion D which is the principal parameter. The kinetic study was carried out under the optimal conditions corresponding to the tests carried out in the previous paragraph (C = 0.01 μL/cm^{3}, T = 30 °C, M = 0.10 g and D = 2 cm). The study of the linearity of equation (6) of the experimental results obtained by the gravimetric technique, gives us approximate values of the diffusion coefficients D_{z} according to the optimal conditions, it is worth 9.81.10^{−6} cm^{2}/h. The experimental results of the simulation of the desorption kinetics of the EOCa in the RC medium are shown in Figure 3. The curves obtained show a very clear branch which reflects a decrease in the mass percentage of the essential oil, which explains the desorption mechanism, we see that the amount desorbed is 5.77% over 24 h.
The simulation made of the model is closer to the experimental results, which proves the validity of the analytical model developed and on the other hand the ability to give a reliable estimate of the other parameters to be sought. In general, the simulation has important interests, it can deduce the quantities of oil that remained during the time of biological applications (insecticidal activities). Figure 4 shows the simulation of the profiles of the concentration of essential oil in the medium. These concentration profiles are used to give good information about the transfer of EOCa within porous media, so that the concentrations can be determined for each altitude and for a more precise time.
Fig. 3 Analytical simulation of the desorption kinetics of EOCa in the RC medium. (•): Experimental results; (…): Simulation, C = 0.01 µ L/cm^{3}; T = 30 ° C; M = 0.10 g; D = 2 cm. 
Fig. 4 Simulation of concentration profiles − Case of desorption of the EOCa in the RC medium. 
3.4 Interaction of the behavior of essential oils in porous media and insecticidal power
As we described previously, the behavior of the EOCa in the RC medium is identified on the one hand by the diffusion parameters: the evaporation flux F, the diffusivity D_{z}, the evaporation constant K and the activation energy E_{a}, and insecticidal activity is presented by M% mortality. The calculated values of all parameters are mentioned in Table 8, the interpretation of these results remains complex. On the other hand, the values of the energy of the activation are constant and low, which confirms that the two clay porous media have the same energetic behavior, moreover, the negative value proves the opposite direction of desorption (adsorption) and it shows the great capacity of porous media to fix essential oils [28].
Principal Component Analysis (PCA) is a statistical tool to study the behaviour of EOCa in RC media with respect to insecticidal activity and to determine hidden information. The case that we are going to study crosses 16 subjects (the matrices produced during the study of the optimization of the flow of evaporation in the two porous media) and the four (4) variables: the flow of evaporation F, diffusivity D_{z}, evaporation constant K and mortality M% (Tab. 8).
To fully explain the results, we must start with the choice of the axes that are interesting for our analysis. Table 9 present the eigenvalues of the four (4) presentative axes obtained. According to the Kaiser criterion [29], the axes F1 and F2 were retained whose eigenvalues are greater than 1 (λ_{1} = 1.89 and λ_{2} = 1.38), with variabilities of 47.31% and 34.70% respectively, and corresponding to the total of 82% of the information. These eigenvalues of the axes selected return a “good proportion” of the analysis, this means that the sum of the inertia explained by each of the axes represents a significant part of the total inertia.
In the mapping displayed in Figure 5, we can see the following interpretations:
The parameters: the diffusion coefficient D_{z} and the evaporation constant K are very close to the correlation circle towards the F1 axis and therefore very well represented on the mapping. The rather closed angle (starting from the origin) formed by the points (D_{z}: cos^{2}(θ) = 0.78 and K: cos^{2}(θ) = 0.72), indicates that these two variables are well correlated with each other. These results come back to the dependence of these two parameters according to equation (9);
The parameters: the flow of evaporation F and the percentage of mortality M% are close to the circle of correlation towards the axis F2, so they are well represented on the mapping. The rather closed angle formed by the points (F: cos^{2}(θ) = 0.24 and M%: cos^{2}(θ) = 0.15), indicates that these two variables are well correlated between them. These results show that the observed M% mortality of insecticide activity depends on the flow of evaporation F;
On the other hand, the almost right angle formed by the four parameters indicates that the two variables of the diffusion coefficient D_{z} and of the evaporation constant K are independent of the two other variables which are the evaporation flux F and the percentage mortality M%.
Parameters of diffusion and of insecticidal activity of EOCa desorption in the RC.
Eigenvalues of the presentative axes of the PCA.
Fig. 5 Mapping of correlations between parameters of EOCa desorption and insecticidal activity. 
4 General discussion
The essential oil of Cedrus atlantica presented a good insecticidal activity against the stored commodity pest of wheat Sitophilus granarius, moreover, this activity was evolved when it was fixed in the clay porous medium. The mechanism of action of this activity is fumigation, hence the appropriate model of the process is a transfer of the material.
To explain the mechanism, we have carried out a study of the relationship between the behavior of the essential oil in the porous medium with the insecticidal activity, hence the study of transfer of essential oil was done by using mathematical analytic models of Fickian diffusion processes under assumptions and theoretical approximations that facilitate numerical calculations. As a result, some physical parameters: the flux F, the diffusivity D_{z}, the evaporation constant K were determined in parallel with the biological parameter of a mortality M% of insecticidal activities against Sitophilus granarius. The use of the statistical tool of the principal component analysis (PCA) which consists in the determination of the correlation between all the obtained parameters allows to show that the biological activity is dependent on the flux F.
In parallel, the optimization of the diffusion transfer of the essential oil was done according to the design of experiments method with a Full factorial design. The results show that the flux F increases with the concentration of the essential oil and with the temperature, but decreases with the cylinder diameter (contact surface) and the mass of the porous medium.
In general, to have a good insecticidal activity against Sitophilus granarius it is necessary to play on the direct increase of the flow of the essential oil fixed on the porous medium by the increase of the concentration of the essential oil and the temperature and by the decrease of the cylinder diameter and the mass of the porous medium.
Finally, the limitations of this study can be cited as follows:
The constituents of essential oils do not have the same affinity with respect to the porous medium;
The insecticidal activity of the essential oil fixed in the porous medium and their persistence depends on the nature of interactions between the constituent compounds of the essential oil and the medium studied.
5 Conclusion
The study of the desorption behavior of Cedrus atlantica essential oil fixed on porous clay medium in relation to its insecticidal activities against Sitophilus granarius was carried out according to mathematical studies of Fick's law and welldefined statistical approaches such as principal component analysis and design of experiments. Several data and information were exploited to explain the mechanism of observation of the remarkable biological activity in the case of oil fixation in the porous medium.
The observed mortality of this insecticidal activity depends on the evaporation flow, and this last parameter also depends on the concentration of essential oil, the temperature, the diameter of the cylinder used and the mass of the porous medium. All the results of the data collected explain that the surface of the porous medium delays the affinity of the constituents of the oil (bioactive compounds) during the desorption.
References
 S.A. Mir, B.N. Dar, M.A. Shah, S.A., Sofi, A.M. Hamdani, C.A. Oliveira, A.S. Sant'Ana, Application of new technologies in decontamination of mycotoxins in cereal grains: challenges, and perspectives, Food Chem. Toxicol. 111976 (2021) [Google Scholar]
 Â. Luís, F. Domingues, A. Ramos, Production of hydrophobic zeinbased films bioinspired by the lotus leaf surface: characterization and bioactive properties, Microorganisms 7, 267 (2019) [Google Scholar]
 R.V. George, H.O. Harsh, P. Ray, A.K. Babu, Food quality traceability prototype for restaurants using blockchain and food quality data index, J. Clean. Prod. 240, 118021 (2019) [Google Scholar]
 R.P. Guiné, S.G. Florença, M.J. Barroca, O. Anjos, The link between the consumer and the innovations in food product development, Foods 9, 1317 (2020) [Google Scholar]
 S. Perera, S. Nanayakkara, M.N.N. Rodrigo, S. Senaratne, R. Weinand, Blockchain technology: Is it hype or real in the construction industry? J. Ind. Inf. Integr. 17, 100125 (2020) [Google Scholar]
 J. Bauchet, S. Prieto, J. RickerGilbert, Improved drying and storage practices that reduce aflatoxins in stored maize: experimental evidence from smallholders in Senegal, Am. J. Agric. Econ. 103, 296–316 (2021) [Google Scholar]
 W.R. Morrison III, A. Lanba, B. Hall, A. Bruce, Novel implementation of laser ablation tomography as an alternative technique to assess grain quality and internal insect development in stored products, J. Stored Prod. Res. 86, 101552 (2020) [Google Scholar]
 T. Vendl, M. Frankova, R. Aulicky, V. Stejskal, First record of the development of Sitophilus oryzae on two rodent bait formulations and literature overview of stored product arthropods infestations in rodent baits, J. Stored Prod. Res. 86, 101557 (2020) [Google Scholar]
 J. Rösner, B. Wellmeyer, H. Merzendorfer, Tribolium castaneum: a model for investigating the mode of action of insecticides and mechanisms of resistance, Curr. Pharm. Des. 26, 3554–3568 (2020) [Google Scholar]
 A.A. Warra, M.N.V. Prasad, African perspective of chemical usage in agriculture and horticulture—their impact on human health and environment, in: Agrochemicals Detection, Treatment and Remediation. ButterworthHeinemann, 2020, pp. 401–436 [Google Scholar]
 H. Musarurwa, N.T. Tavengwa, Emerging green solvents and their applications during pesticide analysis in food and environmental samples, Talanta 121507 (2020) [Google Scholar]
 P. Pavez, R. Figueroa, M. Medina, D. Millán, R.D. Falcone, R.A. Tapia, Choline [amino acid] ionic liquid/water mixtures: a triple effect for the degradation of an organophosphorus pesticide, ACS Omega 5, 26562–26572 (2020) [Google Scholar]
 M.B. Isman, Commercial development of plant essential oils and their constituents as active ingredients in bioinsecticides, Phytochem. Rev. 19, 235–241 (2020) [Google Scholar]
 N. Aoussar, N. Rhallabi, R.A. Mhand, R. Manzali, M. Bouksaim, A. Douira, F. Mellouki, Seasonal variation of antioxidant activity and phenolic content of Pseudevernia furfuracea, Evernia prunastri and Ramalina farinacea from Morocco, J. Saudi Soc. Agric. Sci. 19, 1–6 (2020) [Google Scholar]
 S. Zrira, M. Ghanmi, Chemical composition and antibacterial activity of the essential of Cedrus atlantica (Cedarwood oil), J. Essent. Oil Bear. Plants 19, 1267–1272 (2016) [Google Scholar]
 A. Ainane, M. Taleb, F. ElHajjaji, B. Hammouti, A. Chetouani, T. Ainane, dependence between two types of most abundant natural clays in Bejaad province (Central Morocco) using a statistical approach, Mor. J. Chem. 9, 210–220 (2021) [Google Scholar]
 A. Ainane, F. Khammour, S. Charaf, M. Elabboubi, M. Elkouali, M. Talbi, T. Ainane, Chemical composition and insecticidal activity of five essential oils: Cedrus atlantica, Citrus limonum, Rosmarinus officinalis, Syzygium aromaticum and Eucalyptus globules, Mater. Today Proc. 13, 474–485 (2019) [Google Scholar]
 J. Crank, The mathematics of diffusion (Oxford University Press, 1979), pp. 105–254 [Google Scholar]
 M.B. Teixeira, S. Tobinaga, A diffusion model for describing water transport in round squid mantle during drying with a moisturedependent effective diffusivity, J. Food Eng. 36, 169–181 (1998) [Google Scholar]
 A. Ainane, F. Khammour, M. Elkouali, M. Talbi, A. Oussaid, A. Lemhidi, T. Ainane, Evaluation of the toxicity of the essential oils of certain mints grown in the region of Settat (Morocco): Mentha piperita, Mentha pulegium and Mentha spicata against, Sitophilus Granarius, Sitophilus Oryzae and Sitophilus Zeamais, J. Anal. Sci. Appl. Biotechnol. 1, 1–10 (2019) [Google Scholar]
 L.Á. Oliveira, J.C. Santos, T.H. Panzera, R.T. Freire, L.M. Vieira, F. Scarpa, Evaluation of hybridshortcoirfibrereinforced composites via full factorial design, Composite Structures 202, 313–323 (2018) [Google Scholar]
 D. Granato, J.S. Santos, G.B. Escher, B.L. Ferreira, R.M. Maggio, Use of principal component analysis (PCA) and hierarchical cluster analysis (HCA) for multivariate association between bioactive compounds and functional properties in foods: A critical perspective, Trends Food Sci. Technol. 72, 83–90 (2018) [Google Scholar]
 S. Gharby, A. Asdadi, M. Ibourki, A. Hamdouch, T. Ainane, L.A.I. Hassani, Chemical characterization of the essential oil from aerial parts of Lavandula rejdalii Upson & Jury, a medicinal plant endemic to Morocco, J. Essent. Oil Bear. Plants 23, 1422–1427 (2020) [Google Scholar]
 F. Yang, C. Long, Z. Wei, L. Long, Optimization of medium using response surface methodology to enhance the growth of Effrenium voratum (Symbiodiniaceae, Dinophyceae), J. Phycol. 56, 1208–1215 (2020) [Google Scholar]
 A. Tursi, F. Chidichimo, R. Bagetta, A. Beneduci, BTX removal from open aqueous systems by modified cellulose fibers and evaluation of competitive evaporation kinetics, Water 12, 3154 (2020) [Google Scholar]
 G.C. Leandro, C. Capello, B.L. Koop, J. Garcez, A.R. Monteiro, G.A. Valencia, Adsorptiondesorption of anthocyanins from jambolan (Syzygium cumini) fruit in laponite® platelets: Kinetic models, physicochemical characterization, and functional properties of biohybrids, Food Res. Int. 109903 (2020) [Google Scholar]
 F.A. Houle, R.E. Miles, C.J. Pollak, J.P. Reid, A purely kinetic description of the evaporation of water droplets, J. Chem. Phys. 154, 054501 (2021) [Google Scholar]
 J. Foroozesh, S. Kumar, Nanoparticles behaviors in porous media: application to enhanced oil recovery, J. Mol. Liq. 113876 (2020) [Google Scholar]
 S. AsanteOkyere, C. Shen, Y.Y. Ziggah, M.M. Rulegeya, X. Zhu, Principal component analysis (PCA) based hybrid models for the accurate estimation of reservoir water saturation, Comput. Geosci. 145, 104555 (2020) [Google Scholar]
Cite this article as: Ayoub Ainane, Fatouma Mohamed AbdoulLatif, Jalludin Mohamed, Wissal Attahar, Manar Ouassil, Zine Laabidine Shybat, Adnane El Yaacoubi, Tarik Ainane, Behaviour desorption study of the essential oil of Cedrus atlantica in a porous clay versus insecticidal activity against Sitophilus granarius: explanation of the phenomenon by statistical studies, Int. J. Metrol. Qual. Eng. 12, 12 (2021)
All Tables
Parameters of diffusion and of insecticidal activity of EOCa desorption in the RC.
All Figures
Fig. 1 XRD diffractograms of the clay RC. 

In the text 
Fig. 2 3D representation of flow as a function of two factors. 

In the text 
Fig. 3 Analytical simulation of the desorption kinetics of EOCa in the RC medium. (•): Experimental results; (…): Simulation, C = 0.01 µ L/cm^{3}; T = 30 ° C; M = 0.10 g; D = 2 cm. 

In the text 
Fig. 4 Simulation of concentration profiles − Case of desorption of the EOCa in the RC medium. 

In the text 
Fig. 5 Mapping of correlations between parameters of EOCa desorption and insecticidal activity. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.