ARTICLE

Qualitative Characteristics and Determining Shelf-Life of Milk Beverage Product Supplemented with Coffee Extracts

Ji-Woo Yoon†, Sung-Il Ahn†, Ha-Na Kim, Jun-Hong Park, Sun-Young Park, Jae-Hoon Kim1, Duk-Geun Oh1, Jin-Woo Jhoo, Gur-Yoo Kim*
Author Information & Copyright
Animal Products and Food Science Program, Division of Applied Animal Science, College of Animal Life Sciences, Kangwon National University, Chuncheon 24341, Korea
1Seoul F&B Co., Ltd., Korea
These two authors contribute equally to this work.
*Corresponding author Gur-Yoo Kim Animal Products and Food Science Program, Division of Applied Animal Science, College of Animal Life Sciences, Kangwon National University, Chuncheon 24341, Korea Tel: +82-33-250-8647 Fax: +82-33-251-7719 E-mail: gykim@kangwon.ac.kr

Copyright © 2017, Korean Society for Food Science of Animal Resources. This is an open access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Jan 31, 2017 ; Revised: Mar 22, 2017 ; Accepted: Mar 22, 2017

Published Online: Apr 30, 2017

Abstract

This study was conducted to establish the shelf-life of a milk beverage product supplemented with coffee extracts. Qualitative changes including peroxide value (PV), microorganism content, caffeine content, and sensory evaluation were measured periodically in beverages kept at 10, 20, and 30°C for 8 wk. Lipid oxidation of the product was measured by peroxide value analysis, and apparent changes were observed during a 4 wk storage period. Caffeine analysis revealed that the changes in caffeine content were negligible during the storage period. Total aerobic bacteria, Escherichia coli, yeast, and mold were not detected in the products during an 8 wk storage period. Sensory evaluation revealed that after 4 wk of storage overall acceptance was less than 3 points on a 5-point scale. In this study, PV was used as an indicator of the shelf-life of the milk beverage product. PV analysis revealed that a value of 20 meq/kg was the end of the shelf-life using the Arrhenius equation and the accelerated shelf-life test (ASLT). Assuming that the beverages are kept at 4°C during distribution, calculation of when the PV reached the quality limit point (20 meq/kg) was done with the equation ln(PV) = 0.3644X − 2.21834 and, using that equation, PV = e0.3644X-2.21834 was calculated. Therefore, 14.3086 wk was determined to be the shelf-life of the milk beverage supplemented with coffee when stored at 4°C.

Keywords: milk beverage; coffee; peroxide value; shelf-life; Arrhenius equation

Introduction

The shelf-life of foods is the maximum possible period from the date of manufacture to sale. Physicochemical and microbiological quality changes in the food during the storage period affect the value of a product, so it is very important to manufacturers and consumers to be able to predict how long quality will be preserved under specific conditions (Cho et al., 1993).

There are two ways to determine the expiration date of a food product, the accelerated shelf-life test (ASLT) and an actual experiment. Actual experimentation is the principal method used to set an accurate expiration date, but it requires a long experimental time and relatively high costs in cases where products have shelf-lives of 3 mon or more (Park et al., 2014). The accelerated shelf-life test (ASLT) makes it possible to measure shelf-life in a shorter period of time using harsher conditions, so it is used much more frequently than actual experiments to predict the expiration dates of food economically (Kim et al., 2008). There are several methods associated with the ASLT, such as temperature multiplication methods, oxygen pressure multiplication methods, moisture multiplication methods, and mixing methods (Lee et al., 1995; Lee et al., 1999; Han et al., 2005).

The ASLT with temperature is carried out using the effect of temperature on the decay rate of foods, storing samples at 3 or more temperatures, measuring indicators of quality periodically, calculating an equation by plugging the results of experiments at different temperatures into the Arrhenius equation, and then predicting the shelf-life of samples. This method is efficient in economic aspects, such as time and financial cost, especially in cases of products that have shelf-lives greater than 3 mon. It is known to be very accurate when temperature is used as a single variable (Kim et al., 2008).

As the economy has improved and people spend more time enjoying cultural life nowadays, food service industries have become more upscale and diverse. The coffee beverage sector is growing very fast, and domestic coffee beverage consumption in present-day Korea has been rising sharply (Kim et al., 2013). The palatability of coffee has been increasing, and a number of milk beverages with coffee are being produced and sold on the market. In this study, the shelf-life of a milk beverage supplemented with coffee was predicted using the accelerated test with the Arrhenius equation, which is commonly used to predict shelf-life.

Materials and Methods

Samples

A milk beverage supplemented with coffee that is currently on the market was provided by Seoul F&B Co. (Korea). It is composed of 30% raw milk, 28% coffee extract, 29% purified water, 6% sucrose, and 0.03% emulsifier. Each sample was stored in an incubator set to a different temperature (10, 20 and 30°C) to perform the accelerated test.

Peroxide value (PV)

To measure the peroxide value of each sample, lipid extraction was performed according to ISO Standard 14156: 2001 (International Organization for Standardization, 2001). One gram of extracted lipid was melted in 25 mL of acetic acid (Daejung Chemicals & Metals Co., Ltd.) and chloroform (Daejung Chemicals & Metals Co., Ltd.) (3:2). Then, the samples were left in a dark room after being mixed slightly with 1 mL of a saturated solution of potassium iodide (Junsei Chemical Co., Ltd, Japan). After admixing strongly with 30 mL of water, 0.01 N sodium thiosulfate (Daejung Chemicals & Metals Co., Ltd.) was titrated using 1 mL of a starch solution (Yakuri Pure Chemicals Co., Ltd, Kyoto, Japan) as an indicator.

Peroxide value  ( meq/kg )  =  a b × f amount of sample ( g ) × 10

Where a is the amount of 0.01 N sodium thiosulfate (mL) used in the sample test, b is the amount of 0.01 N sodium thiosulfate used in the blank test, and f is the titer of 0.01 N sodium thiosulfate.

Microbial analysis

To analyze the microbial content of samples stored at 10, 20, and 30°C for 8 wk, 1 mL of each sample was diluted 10-fold with 0.1% peptone solution (Becton, Dickinson and Company, USA). Total aerobic bacteria were counted on 3M PetrifilmTM aerobic count plates (3M Health Care, USA) that were incubated for 48 h at 35°C. E. coli and coliform were counted on 3M PetrifilmTME. coli and coliform count plates (3M Health Care, USA) that were incubated for 24 h at 35°C. Yeast and mold were counted on 3M PetrifilmTM yeast and mold count plates (3M Health Care, USA) that were incubated for 5 d at 25°C. To enumerate thermophiles, the samples were left in an 80°C water bath for 15 min, and assessed for total aerobic bacteria, E. coli, yeast, and mold. Then, red colonies created after incubation were counted and multiplied by the dilution rate.

Caffeine content

The amount of caffeine in samples was measured using HPLC (LC-20A, Shimadzu Co., Japan) with a Luna C18 column (10 μm, 250 × 4.6 mm, Phenomenex, USA). Methanol (Fisher Scientific Korea LTD., Korea), acetic acid (Daejung Chemicals & Metals Co., Ltd.), and water (Fisher Scientific Korea Ltd.) were used as mobile phases (20:1:79) and compounds obtained from the column were analyzed at a wavelength of 280 nm. The column was kept in a column oven at 35°C and a flow rate of 1.0 mL per min to achieve optimum resolution of caffeine. One gram of each sample diluted in 100 mL of water was filtered with a 0.45 μm syringe filter (Adavantec, Japan), and then the solution containing the sample was used to measure the caffeine content. The caffeine content in each sample was quantified by diluting authentic caffeine standard (Sigma, USA) in distilled water at a variety of concentrations to make a standard curve.

Sensory evaluation

A trained sensory panel of 10 researchers at the Program of Animal products and Food science, Kangwon National University evaluated the milk beverage supplemented with coffee once a week over an 8 wk period. The sensory panel consisted of 5 women and 5 men educated in how to evaluate and describe the quality and taste of the sample beverage. Sensory analysis attributes were color, body of coffee, aroma flavor, aftertaste, sweetness, sourness, bitterness, and overall acceptance; the qualities were scored using a scale ranging from 1 (very slight perception) to 5 (very intense perception).

Determination of shelf-life

To predict the shelf-life of the milk beverage supplemented with coffee, a slope, an intercept, and a correlation coefficient were calculated based on the Arrhenius equation through linear regression analysis of PV by accelerating temperature (10, 20, and 30°C).

In most cases, deterioration of the quality of foods follows a zero-order of a first-order reaction formula, indicated with the following numerical expression:

d A / d t = k ( A ) n
(1)

Where A is the quality measurement at time t, A0 is the measurement at time t0, n is the order of the reaction, t is the reaction time (wk), and k is a constant; a zero-order reaction Eq. (2) was derived from Eq. (1), then Eq. (3) was obtained by integrating the zero-order reaction Eq. (2).

d A / dt = k ( A ) n ( n = 0 )
(2)

A t = A 0 k t
(3)

Also, the first-order reaction was analyzed with Eq. (4), below, from Eq. (1), then Eq. (5) was obtained by integrating the first-order reaction Eq. (4).

d A / d t = k ( A ) n ( n = 1 )
(4)

A = A 0 exp ( k t )
(5)

Eq. (6) was derived by taking the natural logarithm (ln) of Eq. (5).

l n A t =   l n A 0 k t
(6)

The Arrhenius equation, which is universally used to determine activation energy, was used to measure the temperature dependence of quality deterioration in the samples. Eq. (7) was revealed by the Arrhenius equation.

K = k 0 exp ( E a / R T )
(7)

Eq. (8) was obtained by taking the natural logarithm of Eq. (7).

ln K = ( E a / R ) × ( 1 / T ) + ln A = S / T + I
(8)

Where T is the absolute temperature, k is the reaction constant at temperature T, R is the gas constant (1.987 Kcal/mol), Ea is the activation energy, S is the slope, and lnA (I) is the intercept.

Meanwhile, the rate of deterioration of product quality was revealed with the below equations using the prior equations and formulas.

k t = t × exp ( ln ( 100 ) exp ( S / T + I ) )
(9)

A t = exp ( ln ( A 0 k t ) )
(10)

Deterioration rate  ( % ) = 100 A t
(11)

Here, t is the storage period, A0 is the quality index at the time of initial storage, At is the quality index at a specific point in time (t), and kt is the reaction constant at a specific point in time.

Statistical analysis

All results are expressed as mean ± SD (standard deviation) of three replicates. Data were analyzed by SPSS statistics 22.0 (SPSS Inc., USA; IBM, Korea) followed by Turkey’s b test for multiple comparisons.

Results and Discussion

Peroxide value (PV)

Rancidity of edible oils and fatty foods caused by lipid oxidation is a serious problem in some sectors of the industry. Because of the unfortunate consequences of lipid oxidation in foods, it is critical that information about the oxidative stability of susceptible food items be obtained before they are marketed (Ragnarsson & Labuza, 1977).

An increase in temperature is the most common and effective means of accelerating oxidation. The rate of the reaction is exponentially related to the temperature; thus, shelf-life decreases logarithmically with increasing temperature. In single-component lipid systems, oxidation can be represented by the following equation:

  • RH + O2 → ROOH → Secondary Products = rancidity or end of shelf-life

Where RH = polyunsaturated lipid substrate; O2 = oxygen; ROOH = lipid hydroperoxide (primary product), and K = overall rate constant for ROOH production.

According to the result of lipid oxidation, it is likely that the PV of 100% lipids extracted from milk beverages stored at each temperature (10, 20, and 30°C) will increase gradually with increasing temperature and storage period (data not shown). However, PV was converted into the concentration of lipid oxidation in 2.3% of lipids in the samples (Table 1). The time point at which 20-40 meq/kg of the PV is reached is generally considered as the starting point of fat rancidity in the case of animal lipids; however, in fresh lipids, it is the time point at which 10 meq/kg of the PV is reached.

Table 1. Peroxide value of the milk beverage supplemented with coffee during an 8 wk storage period at 10, 20, and 30°C (meq/kg)
Storage Temperature (°C) Storage period (wk)

0 1 2 4 6 8
10 0.10879 0.03519 0.09637 0.39629 0.89539 2.37797
20 0.10879 0.05451 0.09591 0.61548 1.23671 2.5806
30 0.10879 0.12075 0.31602 0.91862 1.95707 3.02381
Download Excel Table

Microbial analysis

Foods have a high possibility of supporting microorganisms as time goes by; pathogenic microbes might be especially harmful to humans. Although sterilization and pasteurization are performed on food products, populations of microorganisms can be revived and grow from spores over time, so a bacterial indicator is very important in determining the end of shelf-life (Korea Food and Drug Administration, 2011).

We measured the general and thermoduric microbial count (total aerobic bacteria, E. coli, yeast, and mold) during the storage period (8 wk) and found that none of the microorganisms were detected over the course of 8 wk, which means that sterilization of the beverages was performed properly when they were first manufactured (Table 2, 3).

Table 2. Microbial analysis of the milk beverage supplemented with coffee during storage at 10, 20, and 30°C
Microbe °C Storage period (wk)

0 1 2 4 6 8
E. coli / coliform 10 N.D1) N.D N.D N.D N.D N.D
20 N.D N.D N.D N.D N.D N.D
30 N.D N.D N.D N.D N.D N.D
Aerobic microbe 10 N.D N.D N.D N.D N.D N.D
20 N.D N.D N.D N.D N.D N.D
30 N.D N.D N.D N.D N.D N.D
Yeast 10 N.D N.D N.D N.D N.D N.D
20 N.D N.D N.D N.D N.D N.D
30 N.D N.D N.D N.D N.D N.D

1N.D: Not detected.

Download Excel Table

Table 3. Thermoduric microbial analysis of the milk beverage supplemented with coffee during storage at 10, 20, and 30°C
Microbe °C Storage period (wk)

0 1 2 4 6 8
E. coli / coliform 10 N.D1) N.D N.D N.D N.D N.D
20 N.D N.D N.D N.D N.D N.D
30 N.D N.D N.D N.D N.D N.D
Aerobic microbe 10 N.D N.D N.D N.D N.D N.D
20 N.D N.D N.D N.D N.D N.D
30 N.D N.D N.D N.D N.D N.D
Yeast 10 N.D N.D N.D N.D N.D N.D
20 N.D N.D N.D N.D N.D N.D
30 N.D N.D N.D N.D N.D N.D

1N.D: Not detected.

Download Excel Table

Caffeine content

Caffeine stimulates the central nervous system and can improve alertness and reaction time and mitigate tiredness (Griffiths and Grifiths, 2005). However, higher doses of caffeine (over 200 mg/d) can also cause health problems, including anxiety, increase in blood pressure, headache, nervousness, panic attacks, tremors, sleeplessness, and so on (Griffiths and Grifiths, 2005; Hicks et al., 1996; Stephenson, 1977). Coffee is the most common beverage that people use to take in caffeine, so caffeine content and changes in caffeine content during the storage period need to be measured in beverages supplemented with coffee.

We measured the caffeine content of the milk beverage with coffee. The caffeine content of the samples tended to decrease over the course of 8 wk (Table 4). The caffeine content of samples before storage was the highest (0.43±0.01 mg/mL) among all samples. Some reports mentioned that the caffeine content of coffee or tea is likely to decrease when beverages are stored for a long time (Park and Kim, 2006; So et al., 2014).

Table 4. Caffeine content of the milk beverage supplemented with coffee during storage at 10, 20, and 30°C (mg/mL)
Storage Temperature (°C) Storage period (wk)

0 1 2 4 6 8
10 0.43±0.01A 0.41±0.01AB 0.40±0.00B 0.40±0.00Ba 0.40±0.01B 0.41±0.01ABa
20 0.43±0.01 0.39±0.05 0.41±0.00 0.38±0.01b 0.40±0.01 0.39±0.01b
30 0.43±0.01A 0.43±0.01A 0.40±0.01BC 0.39±0.01Cab 0.40±0.01BC 0.41±0.00Ba

A-CDifferent letters in the same row indicate significant differences at p<0.05.

a,bDifferent letters in the same column indicate significant differences at p<0.05.

Download Excel Table

Sensory evaluation

The results of sensory evaluation are shown in Table 5. There were no significant difference in flavor, aroma, coffee body, and aftertaste during the storage period, but color, sweetness, bitterness, sourness, and overall acceptability changed. Sweetness decreased gradually at 30°C and reached 2.1±1.2, the lowest level of all. In contrast to sweetness, bitterness increased at all temperatures, and sourness also increased during storage, especially significantly at 30°C, and reached 2.9±1.2 at 8 wk, the highest value of all. Overall acceptability likewise decreased, and was lowest at 30°C (1.7±0.8).

Table 5. Sensory evaluation of the milk beverage supplemented with coffee during storage at 10, 20, and 30°C
°C Storage period (wk)

0 1 2 4 6 8
Color 10 2.7±0.8AB 3.3±0.7A 2.7±0.7AB 2.9±0.7AB 2.6±0.7AB 2.2±0.6B
20 2.7±0.8AB 3.2±0.6A 2.5±0.9AB 3.1±0.6A 2.9±0.7AB 2.0±0.7B
30 2.7±0.8 3.3±0.7 2.9±1.1 3.4±0.8 2.9±0.7 2.3±0.9
Flavor 10 3.2±1.0 3.4±1.1 3.6±1.0 3.6±1.0 2.7±1.2 2.9±0.9
20 3.2±1.0 3.0±0.9 3.3±1.1 3.3±0.7 2.5±1.3 2.6±1.2
30 3.2±1.0 3.1±1.1 3.3±1.1 2.8±0.8 2.2±1.3 2.1±1.0
Aroma 10 3.7±1.3 3.9±1.1 3.5±1.1 4.0±0.9 3.0±1.1 3.2±0.8
20 3.7±1.3 3.6±0.8 3.2±1.2 3.6±0.8 2.8±0.9 2.8±1.0
30 3.7±1.3 3.6±1.0 3.2±1.2 3.1±1.2 2.3±1.4 2.3±1.2
Coffee body 10 2.3±1.3 2.6±1.0 2.7±1.1 2.6±0.8 1.9±0.7 2.0±0.8
20 2.3±1.3 2.5±0.9 2.6±1.1 2.5±0.8 1.7±0.7 1.7±0.5
30 2.3±1.3 2.4±0.8 2.5±0.9 2.4±1.1 1.7±0.7 1.7±0.9
After taste 10 3.1±1.4 3.3±1.1 2.9±1.1 3.6±1.1 2.8±1.1 2.6±1.0
20 3.1±1.4 3.2±1.1 3.0±0.9 3.5±1.0 2.4±1.3 2.7±1.4
30 3.1±1.4 3.3±1.3 3.0±0.9 2.8±1.2 2.4±1.3 2.2±1.0
Sweetness 10 3.6±1.0 3.9±1.0 3.9±0.7 3.9±1.0 3.2±1.1 2.8±0.8
20 3.6±1.0 3.8±0.6 3.8±0.9 3.8±1.0 2.7±1.3 2.6±0.8
30 3.6±1.0A 3.9±0.9A 3.5±0.9AB 2.9±1.0ABC 2.3±1.3BC 2.1±1.2C
Bitterness 10 1.3±0.7B 1.2±0.4B 1.0±0.0B 1.3±0.5B 1.4±0.7B 2.1±1.0A
20 1.3±0.7B 1.1±0.3B 1.0±0.0B 1.3±0.5B 1.4±0.7B 2.2±1.0A
30 1.3±0.7AB 1.1±0.3B 1.1±0.3B 1.3±0.5AB 1.5±0.8AB 2.0±1.1A
Sourness 10 1.2±0.6 1.3±0.7 1.3±0.7 1.5±0.8 1.6±1.0 2.1±1.0
20 1.2±0.6 1.2±0.6 1.4±0.7 1.7±1.1 2.2±1.2 2.3±0.8
30 1.2±0.6B 1.5±0.7B 1.8±1.0AB 2.4±1.2AB 2.3±1.2AB 2.9±1.2A
Overall acceptability 10 3.0±1.2AB 3.6±1.0A 3.8±0.6Aa 3.6±1.0Aa 2.9±1.0AB 2.1±1.0B
20 3.0±1.2A 3.3±0.8A 3.4±0.8Aab 3.1±0.7Aa 2.6±1.2AB 1.8±0.6B
30 3.0±1.2A 3.0±1.1A 2.6±1.0ABb 2.2±0.6ABb 2.0±1.2AB 1.7±0.8B

A-CDifferent letters in the same row indicate significant differences at p<0.05.

a,bDifferent letters in the same column indicate significant differences at p<0.05.

Download Excel Table

The results of sensory evaluation of the milk beverage supplemented with coffee showed that all categories tended to receive lower scores gradually during 8 wk of storage. Deterioration of quality was especially pronounced at 30°C as evaluated by taste and overall acceptability.

Determination of shelf-life

When evaluating the experiments on the milk beverage supplemented with coffee during the 8-wk storage period, PV was selected as the final indicator of shelf-life. Fig. 1(a) shows changes in PV at each storage temperature; the graph is exponential. Therefore, all PVs were converted into natural logarithm values, which produced a linear equation (Fig. 1(b)).

kosfa-37-2-305-f001
Fig. 1. Peroxide value of the milk beverage supplemented with coffee at 10, 20, and 30°C.
Download Original Figure

With the peroxide values from samples stored at 10, 20, and 30°C during storage, regression and coefficient (R2) values of each sample were calculated, and are given in Table 6. The coefficient (R2) of the first-order reaction formula was much higher than that of the zero-order reaction formula, so the first-order reaction was used to estimate the shelf-life of the milk beverage with coffee. According to Labuza and Bergquist (1983), the kinetics of PV follow a first-order reaction, and our results corresponded to those findings.

Table 6. Regressions determined by ASLT1) of the milk beverage supplemented with coffee
Regression order Temperature (°C) Regression Coefficient (R2)
0 10 Y = 0.2565X − 0.2097 0.7536
20 Y = 0.2883X − 0.1619 0.8497
30 Y = 0.3591X − 0.1141 0.9397
1 10 Y = 0.3803X − 2.3085 0.9899
20 Y = 0.3961X − 2.1681 0.9973
30 Y = 0.4259X − 2.0477 0.9714

1)ASLT: accelerated shelf-life test.

Download Excel Table

Activation energy and regression were obtained using the Arrhenius equation with the first-order reaction formula (Table 7). The reaction regression of the milk beverage was lnK = −483.79X + 0.737 and the coefficient of determination (R2) was 0.9677, which means that the regression of samples determined with the Arrhenius equation is highly accurate for quality deterioration by PV of samples kept at different temperatures during the storage period.

Table 7. Calculation of activation energy and kinetic regression from zero-order kinetics of the milk beverage supplemented with coffee
Storage temp. (°C) Absolute temp. (T) 1/T K 1) lnK 2) Regression (ln K = −(Ea/R3))(1/T)+ln A) Activation energy (Kcal/mol)
10 283 0.0035 0.3803 −0.96679 Y=−483.79X+0.737 (R2=0.9677) 961.291
20 293 0.0034 0.3961 −0.92609
30 303 0.0033 0.4259 −0.85355

1)K: kinetic constant; 2)lnK: logarithmically transformed kinetic constant; 3)R: gas constant (1.987 Kcal/mol); 4)X: 1/T.

Download Excel Table

We assumed that milk beverages can be kept at 4°C during distribution and calculated when the PV reached the quality limit point (20 meq/kg) using the equation ln(PV) = 0.3644X − 2.21834, and with this equation, PV = e0.3644X-2.21834 was calculated and indicated in Fig. 2. Therefore, 14.3086 weeks was determined to be the shelf-life of the milk beverage supplemented with coffee during storage at 4°C.

kosfa-37-2-305-f002
Fig. 2. Peroxide value of the milk beverage supplemented with coffee at 4°C.
Download Original Figure

Conclusion

This study was carried out to establish the shelf-life of a milk beverage product supplemented with coffee. During the storage period, total aerobic bacteria, E. coli, yeast, and mold were not detected and the results of overall acceptance in sensory evaluation showed less than 3 points on a 5-point scale after a 4 wk storage period. This study revealed that the PV is a suitable indicator of the shelf-life of the milk beverage product. Given the results of PV analysis, a point of 20 meq/kg was determined to be the end of shelf-life through the Arrhenius equation with an accelerated shelf-life test (ASLT). Finally, the regression for qualitative degradation by an increase of PV in the milk beverage supplemented with coffee was ln(PV) = 0.3644X − 2.21834 by the Arrhenius equation, and the shelf-life of the beverage was estimated to be 14.3086 wk at 4°C.

Acknowledgements

This research was financially supported by the Ministry of Education (MOE) and the National Research Foundation of Korea (NRF) through the Human Resource Training Project for Regional Innovation (No. 2015H1C1A103 5886).

References

1.

Cho H. Y., Kwon Y. J., Kim I. K., Pyun Y. R. Estimation of kinetic parameters of nonenzymatic browning reaction using equivalent time at reference temperature with linearly increasing temperature profile. Korean J. Food Sci. Technol. 1993; 25:178-184.

2.

Griffiths R. R., Griffiths J. W. Caffeine, Grolier multimedia encyclopedia. Scholastic Library Publishing Inc.Danbury, CT 2005http://go.grolier.com

3.

Han K. S., Lee E. J., Hong S. P. Prediction of shelf-life of commercially sterilized Korean soups using accelerated experiment. Korean J. Food Cookery Sci. 2005; 21:149-154.

4.

Hicks M. B., Hsieh Y. P., Bell L. N. Tea preparation and its influence on methylxanthine concentration. Food Res. Int. 1996; 29:325-330.

5.

International Organization for Standardization (ISO) Lait et produits laitiers - Méthodes d’extraction des lipides et des composés liposolubles, ISO 14156:2001, FIL-IDF 172: 2001. ISOGeneva, Switzerland 2001

6.

Kim T. H., Chae S. J., Kim C. W. A Study on the coffee consumption behavior by lifestyle. Korean J. Hotel Admin. 2013; 22:93-112.

7.

Kim W. S., Chung M. S., Choi S. H., Kim Y. S., Kim W. S., Lee M. S., Park S. O., Park I. S., Kim A. Y., Park Y. H. Development of experiment model for establishment of shelf-life of foods. Korea Health Industry Development Institute. 2008.

8.

Korea Food and Drug Administration Maximum residue limits of pesticide. Ohsong, Republic of Korea. 2011.

9.

Labuza T. P., Bergquist S. Kinetics of oxidation of potato chips under constant temperature and sine wave temperature conditions. J. Food Sci. 1983; 48:712-715.

10.

Lee N. G., Yoon J. Y., Lee S. R. Computation of Q10 values and shelf-life for canned and bottled orange juices. Korean Soc. Food Sci. Technol. 1995; 27:748-752.

11.

Lee S. G., Lee G. B., Son G. Y. Prediction of shelf-life and quality changes of dried noodle during storage period. Korean J. Food Cookery Sci. 1999; 15:127-132.

12.

Park J. H., An D. S., Lee D. S., Park E. J. Prediction of shelf-life of sea tangle drink. J. Korean Soc. Food Sci. Nutr. 2014; 43:784-790.

13.

Park J. H., Kim J. K. Change in chemical components of green powder tea during storage period at −5°C storage temperature. Korean Soc. Food Preserv. 2006; 13:681-685.

14.

Ragnarsson J. O., Labuza T. P. Accelerated shelf-life testing for oxidative rancidity in foods - A review. Food Chem. 1977; 2:291-308.

15.

So Y. J., Lee M. Y., Yoo K. M., Kang H. J., Hwang I. K. Physicochemical characteristics and antioxidant activity of Dutch coffee depending on different extraction conditions and storage. Korean J. Food Sci. and Technol. 2014; 46:671-676.

16.

Stephenson P. E. Physiologic and psychotropic effects of caffeine on man - A review. J. Am. Diet. Assoc. 1977; 71:240-247.