681
Views
7
CrossRef citations to date
0
Altmetric
Original Articles

Modeling the Equilibrium Moisture Sorption Data for Some Nigerian Foods

&
Pages 1-13 | Received 27 Feb 2003, Accepted 23 Mar 2004, Published online: 06 Feb 2007

Abstract

The water sorption isotherms of thirty different Nigerian food products at various temperatures were obtained from literature. The collected data were divided into three groups namely cassava-based products, grains, and others. The data were used to test the applicability of nine sorption models using three statistical parameters as the criteria. The GAB model fits most of the data well when either a mean percentage deviation or standard error estimate was used as the criteria. With the coefficient of determination as the criterion, the Exponential, Bradley, and Smith models were the best for cassava-based products, grains, and others respectively.

Introduction

Nigeria is a nation with many agricultural products. Such products include rice, sorghum, wheat, plantains, yams, okra, and cassava. Most of these products are available on a seasonal basis and normally require a short time to mature. They are, therefore, abundant during harvest with the supply tapering off as the next harvest approaches. Consequently, a good percentage of these farm products spoil during the abundant period. It should be noted that cassava alone provides the major source of dietary calories; it is used as a cheap source of carbohydrate food for man and livestock.Citation[1] For human consumption, it can either be cooked or eaten raw as a vegetable or it can be grated and roasted into flour and eaten as gari. A fundamental property of biological material, which influences dehydration, shelf life prediction, and storage stability, is the water sorption characteristics.Citation[2]

Many authors have presented the sorption isotherms of foods that can be found in Nigeria. While some of these authors have modeled their data, others have not. The authors have also used different criteria to assess the suitability of the model to fit their sorption data. The main objectives of this work are to use a set of nine models to predict the sorption isotherm for some Nigerian foods and to use three different criteria to assess their suitability. The sorption data have been obtained from literature.

Data Sources and Analysis

Water Sorption Data

A total of ninety-one sorption isotherm data sets for different Nigerian foods have been obtained from literature. The sources are summarized in Table . The collected data are divided into three groups: cassava based, grains, and others. With the exception of Onayemi and OluwamukomiCitation[3] who presented the sorption data in tabular form, the data are in graphical form. The actual values used for this work are extracted from these graphs. A minimum of six pairs of equilibrium moisture content (m) and water activity (a w ) were extracted from each graph. The ranges for m and a w for each product are also included in Table ; the actual values are shown in Table .

Table 1 Sources and range of the various sorption dataFootnote*

Table 2 Raw sorption data for the different products

Sorption Models

Available models are either dependent or independent of temperature. Those that are independent of temperature, in general, give a better fit for the sorption data at any given temperature. They are either linear or can easily be linearized. For this work, we have used nine such models to fit the extracted sorption data. The linear forms of these models are shown in Table where m o is the monolayer moisture content; c and k are related to the monolayer and multilayer properties respectively; A and B are empirical constants; m o , c, and k must be positive to have any physical meaning.Citation[4] Also, k and c are both temperature dependent and are related to the interaction energy between water and food and the interaction energy between the multiple layers of water, respectively.Citation[15]

Table 3 Linear forms of selected sorption models

As can be seen from Table , the GAB and Iglesias-Chirife models are three-parameter equations while the others are two-parameter equations. The GAB and BET models have theoretical basesCitation[2] while the others are empirical. When the parameter k is 1, the GAB equation reduces to the BET equation. For any particular product, the numerical values of m o and c obtained using the BET and GAB models are usually different. It has also been observed that most models can only predict the sorption for a particular range of activities. For example, the BET model is known to give a good fit for a w up to 0.43.Citation[15] However, our use of the nine models is to ascertain whether any of them would predict the sorption data for the whole range of interest. It should be noted that apart from Sanni et al.Citation[2] who used five of the models presented, other authors have used at most three such models.

Assessing Goodness of Fit

Different statistical parameters have been used in the literature to assess the suitability of the models for predicting the sorption data. Some of these parameters areCitation[6] Citation[28] Citation[29]: coefficient of determination (r 2), standard error of the estimate for water activity (SE), mean relative percentage deviation (P %), Root mean square deviation (B), χ 2—test, and plots of residuals.

From Table , it is apparent that the linear forms of the models are transformed sorption data variables. Hence, a value of r 2 close to 1 does not necessarily imply that the model fits the actual data very well. From a practical point of view, P permits a direct visualization of the fitting ability of a given model. The residual plot is useful for detecting inconsistencies between experimental and calculated values. A random distribution of the residuals around zero indicates that the model correctly represents the particular set of data. For comparison of the various models listed in Table , we have used only r 2, SE, and P.

Results and Discussions

For each product and at a given temperature, the model that gave the best fit using r 2, SE, and P as criteria are shown in Table . It can be seen that, except for few cases, the three criteria do not give the same best fit model. Using P as the criterion, it can be seen that the GAB model fits about 60% of the data used for this study while the Caurie model was not suitable for any of the data studied. A closer look at Table will also reveal that the GAB model fits 76% of the cassava products, 58% of the grain data, and 48% of the others. When SE is used as the criteria, a similar trend as that for P is obtained but with lower percentages. With r 2 as the criteria, the Exponential model fits 38% of the cassava-based data; the Bradley model fits 35% of the grains data, and Smith model fits 35% of the others data. The GAB, BET, and Caurie models gave poor fits if r 2 is used as the fitness criteria. The r 2 for the GAB model varies from 0.9291 to 0.9971 for cassava-based products, 0.4821 to 0.9998 for grains, and 0.018 to 0.9931 for others. This shows that although the GAB model was not the best for cassava-based products, it was also not the worst when r 2 is used as the criterion. The BET model is known to be a good fit for a w up to 0.43Citation[15]; this explains why it was a poor fit for our data since the a w s are higher than 0.50 (see Table ). It would appear that the Caurie model is not suitable for the data sets used for this work.

Table 4 Best model using the different criteria and the GAB parameters

The GAB parameters are shown in Table . For any given product, each of the parameters varies with temperature. This confirms the earlier resultsCitation[2] Citation[16] which showed that the influence of temperature on GAB parameters can be calculated with an Arrhenius form of equation. Table also shows that m o is negative for ripe plantain at 40 and 50°C and for yam flour at 25°C despite their respective P and SE values being the lowest. Since a negative m o has no physical meaning, it follows that the GAB model cannot be used for these data.

The 95% confidence interval for the average values of mo, c, and k (ignoring the negative values) is computed for each type of food products, shown in Table . The average values for k are less than 1 and are approximately constant (two decimal place accuracy) for each of the food product types. Our data, therefore, validate the fact that, theoretically, values of k should be less than 1. The confidence intervals for m o and c are much larger than that for k. It would, therefore, appear that, within the 95% confidence limit, only m o and c vary with temperature.

Table 5 95% confidence intervals for the GAB parameters

A two-way analysis of variance was also used to compare the values of P, r 2, and SE obtained for the nine models when applied to the various sorption data. The results are summarized in Table . Note that F(P), F(r 2), and F(SE) are the calculated F values for P, r 2, and SE data respectively while F crit represents the critical F-Statistic values at the specified probabilities (p = 0.05 and p = 0.01). Table shows that the computed F values are significantly greater than the critical F values at the 0.05 or 0.01 levels. We can, therefore, conclude that there is significant difference among the nine models considered, irrespective of the criteria used. The variability is highest when r 2 is used as the criterion, followed by P while SE is the lowest.

Table 6 Summary of the analysis of variance for the P, r 2 and SE data

The ability of one or more criteria to predict the same best model is shown in Table . This table shows that none of the three criteria gave the same best model for 7% of the data sets used while the three criteria gave the same best model for 15% of the data sets. It would also appear that using either P or SE would give the same best model for a majority of the data sets used. This is probably because the linear form is a transformation of m and/or a w . Thus, the r 2 is an assessment of how linear the transformed variable is rather than the actual variable. P and SE are dependent on the actual variables, and these probably explain why these two criteria gave similar results for most of the foods. The results for the worst model (not shown here) also follow similar trend.

Table 7 Ability of the different criteria to predict the same best model

Conclusion

The sorption isotherm data for thirty Nigerian foods at various temperatures have been modeled by using nine different equations. The goodness of fit for each model was assessed by using three statistical parameters. It was found that one parameter is not good enough to assess the predictive ability of the models. Using either the mean percentage deviation or standard error of the estimate as the criterion as the parameter, the GAB model fits roughly 60% of the data. For cassava-based products, the GAB model gives the best fit for 75% of the data sets. With the coefficient of determination as the criterion, the GAB model is not as good as other models such as the Exponential, Bradley, or Smith models. Except for very few cases, the three criteria did not produce the same best model for any particular data set.

References

  • Kemdirim , O.C. , Chukwu , O.A. and Achinewu , S.C. 1995 . Effect of traditional processing of cassava on the cyanide content of gari and cassava flour . Plat Foods for Human Nutrition , 48 : 335 – 339 .
  • Sanni , L.O. , Atere , C. and Kuye , A. 1997 . Moisture sorption isotherms of fufu and tapioca at different temperatures . Journal of Food Engineering , 34 ( 2 ) : 203 – 212 .
  • Onayemi , O. and Oluwamukomi , M.O. 1987 . Moisture equilibria of some dehydrated cassava and yam products . J. Food Process Eng. , 9 : 191 – 200 .
  • Salgado Cervantes , M. , Garcia Alvarado , M.A. and Waliszewski Kubiak , K.N. 1994 . Modeling of water activity and enthalpy of water sorption in cassava chips . Drying Technology , 12 ( 7 ) : 1743 – 1752 .
  • Sanni , L.O. , Kolawole , A.G. , Akingbala , J.O. and Kuye , A. 1999 . Effect of drying methods on moisture sorption isotherms of fufu at three temperature . Drying Technology , 17 ( 1 & 2 ) : 285 – 297 .
  • Kuye , A. and Sanni , L.O. 2002 . Modelling sorption isotherms of lafun and soyflour using a spreadsheet . International J. of Food Properties , 5 ( 3 ) : 599 – 610 .
  • Sanni , L.O. and Kuye , A. Effect of temperature on water sorption isotherms of some fortified cassava products . Proceeding of the 12th International Drying Symposium, IDS2000 . Armsterdam : Elsevier Science . Paper no. 23
  • Sopade , P.A. and Ajisegiri , E.S. 1994 . Moisture sorption study on Nigerian foods: maize and sorghum . J. Food Process Engineering , 17 : 33 – 56 .
  • Jideani , J.A. , Durojaiye , O.O. and Jideani , I.A. 1995 . Moisture sorption isotherms of intermediate moisture millet dough . Nigerian Food Journal , 13 : 94 – 102 .
  • Aguerre , R.J. , Suarez , C. and Viollaz , P.E. 1983 . Moisture desorption isotherms of rough rice . J. of Food Technology , 18 : 345 – 351 .
  • Andrieu , J. , Stamatopoulous , A. and Zafiropoulous , M. 1985 . Equations for fitting desorption isotherms of durum wheat pasta . J. of Food Technology , 20 : 651 – 657 .
  • Ajisegiri , E.S. , Sopade , P.A. and Abass , A.B. 1994 . Moisture sorption on Nigerian foods Kuka . J. of Stored Prod. Res. , 30 : 331 – 338 .
  • Ihekoronye , A.I. and Ngoddy , P.O. 1985 . Integrated Food Science and Technology for the Tropics N.Y. : Macmillan Pub Ltd. .
  • Ajibola , O.O. 1986 . Desorption isotherms of gari from 40°C to 70°C . J. Agric. Eng. Res. , 35 : 207 – 210 .
  • Kim , K.H. , Song , Y. and Yam , K.T. 1991 . Water sorption characteristics of dried red peppers (Capsicum annum L.) . Int J. of Food Science & Technology , 29 : 339 – 345 .
  • Kim , S.S. and Bhowmik , S.R. 1994 . Moisture sorption isotherms of concentrated yogurt and microwave vacuum dried yogurt powder . J. Food Eng. , 21 : 157 – 175 .
  • Brunauer , S. , Emmelt , P.H. and Teller , E. 1938 . Adsorption of gases in multimolecular layers . J. Am. Chem. Soc. , 60 : 309
  • Bradley , R.S. 1936 . Polymer adsorbed films Part I . The adsorption of argon on salt crystals at low temperature and the determinations of surface fields. J. Chem. Soc. , 58 : 1467 – 1474 .
  • Caurie , M. 1970 . Minimum moisture for stored food . J. Food Technol. , 5 : 301 – 307 .
  • Lazarides , H.N. 1990 . Sorption isotherm characteristics of an intermediate meat product . Lebensm-Wiss u—Techn. , 23 : 418 – 421 .
  • Guggenheim , E.A. 1966 . Applications of Statistical Mechanics Oxford Clarendon Press .
  • Anderson , R.H. 1946 . Modification of the BET equation . J. Am. Chem. Soc. , 68 : 689 – 691 .
  • de Boer , J.H. 1953 . The Dynamic Character of Adsorption, , 2nd Ed. Oxford Clarendom Press .
  • Iglesias , H.A. and Chirife , J. 1978 . An empirical equation for fitting water sorption isotherms of fruits and related products . Can Inst. Food Sci. Technol. J. , 11 : 12 – 18 .
  • Boquet , R. , Chirife , J. and Iglesias , H.A. 1978 . Equations for fitting water sorption isotherms of food. Part II: Evaluation of various two-parameter models . J. Food Technology , 13 : 319 – 328 .
  • Oswin , C.R. 1946 . The kinetics of packaging life III: the isotherm . J. Chem. Ind. London , 65 : 419
  • Smith , S.E. 1947 . The sorption of water vapour by high polymers . J. Am. Chem. Sac. , 69 : 646
  • Sopade , P.A. , Ajisegiri , E.S. and Badau , M.H. 1992 . The use of Peleg's equation to model water absorption in some cereal grains during soaking . J. of Food Engineering , 15 : 269 – 283 .
  • Sopade , P.A. 2001 . Criteria for an appropriate sorption model based on statistical analysis . Int. J. of Food Properties , 4 ( 3 ) : 405 – 418 .

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.