Next Article in Journal
Nasopharyngeal Carcinoma Burden and Its Attributable Risk Factors in China: Estimates and Forecasts from 1990 to 2050
Previous Article in Journal
Prevalence of Physical Violence in the Medical-Forensic Approach in the Years 2015–2020 in City and Neighboring Municipalities: Perspectives from Poland—Poznań Study
Previous Article in Special Issue
Assessment of Radiation Risk Perception and Interest in Tritiated Water among Returnees to and Evacuees from Tomioka Town within 20 km of the Fukushima Daiichi Nuclear Power Plant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Survey of Domestic Refrigerator Storage Temperatures in Poland for Use as a QMRA Tool for Exposure Assessment

by
Constantine-Richard Stefanou
*,†,
Anna Szosland-Fałtyn
and
Beata Bartodziejska
Institute of Agricultural and Food Biotechnology State Research Institute, 02-532 Warszawa, Poland
*
Author to whom correspondence should be addressed.
EU-FORA Fellowship Programme hosted at IBPRS-PIB.
Int. J. Environ. Res. Public Health 2023, 20(4), 2924; https://doi.org/10.3390/ijerph20042924
Submission received: 30 December 2022 / Revised: 22 January 2023 / Accepted: 5 February 2023 / Published: 8 February 2023

Abstract

:
In the framework of Quantitative Microbiological Risk Assessment, the estimation of the ingested dose of a hazard by the consumer is of paramount importance. This may be calculated by means of predictive modeling of growth/inactivation of the pathogen studied. For products that spend the majority of their shelf life in the domestic refrigerator, storage temperature will significantly impact the microbial population dynamics. To describe the variability of domestic storage temperatures in Poland, a survey including 77 participants, was carried out in Lodz, Poland. Participants were provided with temperature data loggers, which measured their refrigerator temperature for 24 h in 5-min intervals. The temperature-time profiles were used to calculate the mean working temperature, standard deviation, minimum and maximum values, and the data were statistically analyzed to find the best fitting probability distribution using R programming language. Out of the tested refrigerators, 49.35% had a mean working temperature of over 5 °C and 3.9% exceeded 10 °C. Distribution fitting scenarios were tested for goodness of fit, and the final selected distribution was a truncated normal distribution. This study can prove useful in Monte Carlo simulation analysis for stochastic quantitative food risk assessment in Poland.

1. Introduction

Perishable foods often require refrigeration to hinder microbial spoilage and/or pathogen growth. The controlled temperature of a refrigerator is meant to inhibit or delay the growth of microbes and particularly pathogens that could pose a threat to human health. Improper storing temperatures may allow for growth and increase the risk of foodborne disease. Domestic refrigerated storage is a crucial part of the cold chain, and proper storing temperatures contribute to food safety [1]. Control measures and good practice guides are in place for all other parties in the farm-to-fork path, including producers, processors, and distributors. Consumers, on the other hand, can be uninformed about proper food storage temperatures and practices [2]. Thus, temperature abuse or improper storage can occur during domestic refrigerated storage [3,4]. The importance of proper consumer storage and handling in terms of food safety is reflected in the reported source of outbreaks in the EU Zoonosis annual reports. For the period of 2016–2021, a total of 1411 household outbreaks occurred in Poland with the most commonly reported foodborne pathogen being Salmonella spp. [5]. A compounding factor for this is that most consumers have an optimistic bias when it comes to the domestic environment, viewing it as a safe place, hence, their risk perception at home is reduced. In addition to consumer behavior at home, the dietary shift towards ready-to-eat foods (RTE) draws attention to domestic storage as such products can be mildly processed with little or no preservatives and rely on respecting the cold chain to ensure safety and quality throughout the duration of their shelf-life [1].
For psychorotrophic bacteria in particular, temperature during storage can be the decisive factor for growth inhibition and control. Storage at increased temperatures could allow growth in contaminated products and lead to an increased risk of foodborne disease. Quantitative microbiological risk assessment (QMRA) requires the estimation of the ingested dose of a hazard by the consumer. This is calculated during the exposure assessment step by modeling the fate of the hazard in the food prior to consumption. Furthermore, reliable data on consumer storage are of paramount importance for establishing a shelf-life for perishable refrigerated foods [1].
Refrigerators can vary considerably, and the temperature fluctuates constantly during the cooling cycle. This variability in storage temperatures is of paramount importance, especially in stochastic Quantitative Microbiological Risk Assessment (QMRA) [6]. In deterministic QMRA, when simulating the consumer phase, the storage temperature is defined as a single point estimate, usually the mean or mode, instead of a probability distribution. This approach does not take into account the inherent variability of refrigerator working temperatures and may lead to erroneous final estimations. Therefore, the aim of this study was to collect data on domestic storage temperatures in Polish fridges and describe the variability of the storage temperature for use as input in stochastic risk assessment exercises and shelf-life determination.

2. Materials and Methods

Participants were provided with a temperature data logger (LogTag model TRIX-8, LogTag North America Inc., Auckland, New Zealand) and instructed to place it on the middle shelf of their refrigerator for 24 h. Participants were told not to change their refrigerator settings and to perform the measurement under normal usage conditions. This aimed at obtaining a realistic temperature data set that reflected everyday use and conditions. The survey carried out included 77 participants, inhabitants of Lodz in Poland, and took place from January to July of 2022. The data logger took temperature readings at 5-minute intervals, and these data points were used to calculate the mean temperature, standard deviation, minimum and maximum values for each refrigerator.
The recorded mean temperature values were statistically analyzed to find the best-fitting probability distribution to describe the variability of Polish household refrigerator working temperatures. The probability distribution fitting was performed in R with the fitdistrplus package. Fitting was performed for the Weibull, triangular, lognormal, pert, and normal distributions, with parameters being estimated by using the maximum likelihood estimation (mle) method. All distribution fitting scenarios were statistically tested for goodness-of-fit. To select the distribution which best described the data, the Chi-squared, Kolmogorov–Smirnov (KS), Cramér–von Mises (CVM), and Anderson–Darling (AD) statistics were calculated alongside the Akaike’s Information (AIC), and Bayesian Information Criteria (BIC). These values were compared for all fitted distributions to confirm the best fit.

3. Results

3.1. Domestic Refrigerator Temperature Results

Indicative temperature profiles of Polish domestic refrigerators are presented in Figure 1. The data from all collected temperature profiles were used to calculate the mean working temperature of each individual refrigerator and construct a histogram for the distribution and the empirical cumulative density plot of the mean temperatures (Figure 2). A detailed summary of the recorded data is provided in Appendix A.
Table 1 presents the frequency of observed mean temperatures for the surveyed domestic refrigerators. Out of the tested refrigerators, 49.35% had a mean working temperature of over 5 °C and 3.9% exceeded 10 °C. Out of all the collected data points of all temperature profiles, 49.42% of the temperature readings were above 5 °C and 6.08% were above 10 °C. The mean working temperature was found to be 5.1 °C with a standard deviation of 2.8 °C, and minimum and maximum mean values were −1.2 °C and 15.2 °C, respectively.
A Box and Whisker plot was used to display the central tendency and dispersion of the mean working temperatures of the surveyed refrigerators (Figure 3). The 25th percentile (Q1) was 3.49 °C, and the 75th percentile (Q3) was 6.92 °C with an interquartile range (IQR) of 3.43 °C. The median was calculated as 5 °C and the ends of the whiskers were at −1.66 °C and 12.07 °C, respectively. The dot corresponds to the outlier value of 15.2 °C.
The individual temperature profiles (Figure 1) displayed significant variation and hence, aiming to evaluate this variability for each profile, in addition to the mean working temperature, the standard deviation was also estimated from each temperature profile. The histogram of the calculated standard deviations for the 77 measured temperature profiles is shown in Figure 4. Relying on a single average value can prove inaccurate as the variability can range vastly from 0.3 °C to 3.7 °C, as seen in the figure.

3.2. Probability Distribution Fitting to Temperature Data

The fitting of probability distributions to the average temperature data was performed in R using the fitdist function. The data were fitted to the Weibull, triangular, lognormal, pert, and normal distributions. For each fitting, maximum likelihood estimation was used to estimate the probability distribution parameters. The goodness-of-fit statistics and criteria for each fitted distribution are presented in Table 2. Additionally, the bias factor for each fitted distribution was calculated using the Monte Carlo simulation method for 10,000 iterations, and the calculated mean bias factor is also presented in Table 2.
All the calculated statistics and criteria, with the exception of the Kolmogorov–Smirnov statistic pointed to the normal distribution as the best fitting distribution to describe the mean temperature of domestic refrigerators. The fitted normal would require truncation at the minimum and maximum recorded values of −1.2 °C and 15.2 °C. In contrast, the Kolmogorov–Smirnov statistic had the lowest value for the pert distribution. The suitability of the normal distribution as the preferred probability density function to describe these data was also evident in the histogram of theoretical densities and the P–P plot for all fitted probability distributions (Figure 5).

4. Discussion

Domestic storage is a key and often overlooked link in the food cold chain, and, to our knowledge, this survey is the first to focus on domestic storage temperatures in Poland. The survey conducted revealed that 88.31%, the majority of domestic refrigerators, function at over 4 °C, with an additional 3.9% having a mean temperature of above 10 °C. These numbers are in stark contrast with the European Food Safety Agency’s suggestion of below 5 °C for household refrigerators [7]. The suggested refrigerator temperature can vary among countries, but all agree on setting the limit at below 7 °C, with many countries lowering this to below 5 °C [8]. In our study, the temperature of the room in which the refrigerator was located was not considered, as Laguerre et al. noted no impact on the refrigerator temperature. Laguerre et al. also concluded that the temperature setting options of fridges did not influence the measured mean temperature and, as such, this factor was also excluded from the current study [9]. Drawing from the work of Jofré et al., door opening frequency was assumed to have little effect on the data [1]. The periodic fluctuations of temperature observed for each refrigerator, visible in Figure 1, are expected to be the result of the refrigeration cycle. Fluctuations during the course of the measurements for each refrigerator were observed, with the individual standard deviation estimated for each temperature profile ranging from 0.3 °C to 3.7 °C. James et al. note that temperature fluctuations in refrigerators can vary even for refrigerators with the same mean working temperature [8].
The surveyed households had fridges working at a mean temperature ranging from as low as −1.2 °C up to 15.2 °C. This noted variability in storage temperatures is of paramount importance for both food business operators (FBOs) and risk assessors aiming to guarantee the legally set appropriate level of protection (ALOP) for refrigerated foods. According to article 3 of Regulation (EC) No 2073/2005, FBOs are obliged to ensure that the food safety criteria applicable throughout the shelf-life of foods can be met under reasonably foreseeable conditions of distribution, storage, and use, including the consumer phase. Furthermore, the setting of use-by-dates for refrigerated pre-packed and ready-to-eat products requires applying stochastic methods for realistic estimations [10]. This is especially true in the case of psychotropic pathogens, such as Listeria monocytogenes, which are often detected in products meant for direct consumption [11].
The switch from single-point estimates for microbial growth parameters, such as storage temperature, to distributions requires a shift in perspective towards the full range of possible and probable values. Probability distributions reflect the variability of parameters, such as temperature, and aid in supporting more informed risk management decisions.
Based on similar studies in the literature, the Probability Density Function expected to best fit the data was expected to be the Normal distribution [1,10,12]. An adequate number of temperature data points was required for better distribution description and more precise fitting selection. For the statistical tests and criteria used to rank the fitted distributions, the lower the value the better fit of the distribution to the data. The Chi-squared, KS, and AD tests are often used to evaluate the acceptability of a theoretical distribution for the description of observed values [13]. The values of these statistics were used to compare the goodness-of-fit of several probability distributions to observed data. Additionally, the Akaike and Bayesian information criteria were used to estimate the prediction error and, therefore, relative quality of the fit. Applying these goodness-of-fit tests, the final domestic storage temperature distribution for Poland was described as a truncated normal distribution with a mean = 5.1 °C, sd = 2.8 °C and truncated at −1.2 °C and 15.2 °C. Roccato et al. in their analysis of domestic refrigerator temperatures in Europe, opted for describing the variability of European domestic refrigerators by a normal distribution with mean = 6.1 °C and sd = 2.8 °C for the northern European countries. Despite data from Poland not being included in the aforementioned study, our results are in relative agreement [10]. Future similar surveys could further validate our findings by including a larger number of participants and spanning a longer period of time. For all tested distributions, the bias factor estimated from a 10,000 iterations Monte Carlo simulation was within the acceptable range of 0.75–1.25.

5. Conclusions

The domestic environment can act as the setting of outbreaks of foodborne illnesses and storage temperature is a key factor for food safety at the consumer level, since it is a major controlling factor of microbial behavior in foods. A domestic refrigerator temperature survey was carried out for the first time in Poland, aiming to quantify the variability of mean storage temperatures. This study has provided a probability density distribution describing domestic refrigerator temperatures in Poland. The described truncated normal probability distribution (mean = 5.1 °C, sd = 2.8 °C, min = −1.2 °C, max = 15.2 °C) can be incorporated in Monte Carlo simulations for stochastic quantitative risk assessment and shelf-life determination. Efforts to quantify consumer phase storage variability increase the robustness and realism of risk assessments and, therefore, crucial risk management decisions critical for ensuring food safety.

Author Contributions

Conceptualization, C.-R.S.; data curation, C.-R.S. and A.S.-F.; writing—original draft, C.-R.S.; writing—review & editing, C.-R.S.; supervision, B.B. and A.S.-F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Raw data of the surveyed domestic refrigerators.
Table A1. Raw data of the surveyed domestic refrigerators.
IDRecorded Mean Temperature (°C)Standard Deviation (°C)Minimum Value (°C)Maximum Value (°C)Range
17.01.25.214.29
26.20.85.113.68.5
39.40.78.311.73.4
48.90.77.9168.1
52.71.6−6.36.112.4
63.61.32.215.713.5
76.90.56.310.54.2
87.60.76.613.16.5
94.21.12.211.79.5
105.30.83.711.37.6
112.11.50.31312.7
12−0.41.5−3.210.213.4
135.51.03.110.67.5
145.50.54.914.49.5
157.50.96.312.66.3
167.80.95.912.56.6
173.90.73.213.310.1
186.61.14.812.98.1
193.21.30.310.910.6
205.31.03.410.97.5
21−0.21.0−0.812.513.3
227.40.46.811.34.5
232.61.11.59.68.1
243.41.3−1.38.29.5
255.51.13.111.48.3
263.50.82.49.26.8
27−1.20.8−2.18.210.3
285.30.63.98.74.8
293.30.52.385.7
305.00.54.413.18.7
317.60.76.413.26.8
327.40.66.311.85.5
334.81.11.412.811.4
343.40.91.711.710
355.30.94.213.29
364.11.51.613.511.9
374.72.40.58.98.4
387.00.75.812.56.7
396.13.7−0.811.812.6
402.41.01.59.27.7
417.60.85.513.98.4
424.81.03.1117.9
43−0.10.7−1.19.911
444.91.81.211.910.7
457.00.75.511.66.1
465.50.84.8138.2
473.50.91.811.69.8
486.91.33.715.411.7
493.30.81.18.57.4
504.10.6311.38.3
510.70.9−0.511.111.6
524.41.8−5.67.413
534.60.345.81.8
544.70.93.410.47
554.71.23.311.78.4
564.10.92.885.2
575.50.64.210.46.2
5810.30.58.312.84.5
591.21.2−0.49.59.9
601.90.90.47.77.3
610.60.8−0.210.911.1
628.21.14.915.710.8
636.10.65.211.96.7
645.10.74.410.66.2
6515.20.514.421.57.1
663.71.6−0.210.911.1
677.20.56.311.14.8
6810.50.410.115.75.6
696.33.3−1.613.415
708.70.38.414.56.1
716.31.34.613.28.6
725.51.12.913.610.7
734.90.93.715.311.6
748.01.55.817.211.4
755.01.13.11511.9
762.40.91.39.68.3
77−0.60.8−1.39.510.8

References

  1. Jofré, A.; Latorre-Moratalla, M.L.; Garriga, M.; Bover-Cid, S. Domestic refrigerator temperatures in Spain: Assessment of its impact on the safety and shelf-life of cooked meat products. Food Res. Int. 2019, 126, 108578. [Google Scholar] [CrossRef]
  2. Ovca, A.; Škufca, T.; Jevšnik, M. Temperatures and storage conditions in domestic refrigerators—Slovenian scenario. Food Control 2021, 123, 107715. [Google Scholar] [CrossRef]
  3. Janjić, J.; Katić, V.; Ivanović, J.; Bošković, M.; Starčević, M.; Glamočlija, N.; Baltić, M.Ž. Temperatures, cleanliness and food storage practises in domestic refrigerators in Serbia, Belgrade: Temperatures, cleanliness and food storage practises in Serbia. Int. J. Consum. Stud. 2016, 40, 276–282. [Google Scholar] [CrossRef]
  4. Taoukis, P.S.; Giannakourou, M.C.; Koutsoumanis, K.; Bakalis, S. Modelling the Effect OF House Hold Chilled Storage Conditions on the Risk Distribution of Meat Products. Acta Hortic. 2005, 674, 435–439. [Google Scholar] [CrossRef]
  5. European Food Safety Authority; European Centre for Disease Prevention and Control. The European Union One Health 2021 Zoonoses Report. EFSA J. 2022, 20, 7666. [Google Scholar] [CrossRef]
  6. Koutsoumanis, K.; Pavlis, A.; Nychas, G.-J.E.; Xanthiakos, K. Probabilistic Model for Listeria monocytogenes Growth during Distribution, Retail Storage, and Domestic Storage of Pasteurized Milk. Appl. Environ. Microbiol. 2010, 76, 2181–2191. [Google Scholar] [CrossRef]
  7. Ricci, A.; Allende, A.; Bolton, D.; Chemaly, M.; Davies, R.; Fernández Escámez, P.S.; Girones, R.; Herman, L.; Koutsoumanis, K.; Nørrung, B.; et al. Listeria monocytogenes contamination of ready-to-eat foods and the risk for human health in the EU. EFSA J. 2018, 16, 5134. [Google Scholar] [CrossRef]
  8. James, C.; Onarinde, B.A.; James, S.J. The Use and Performance of Household Refrigerators: A Review. Compr. Rev. Food Sci. Food Saf. 2017, 16, 160–179. [Google Scholar] [CrossRef] [PubMed]
  9. Laguerre, O.; Derens, E.; Palagos, B. Study of domestic refrigerator temperature and analysis of factors affecting temperature: A French survey. Int. J. Refrig. 2002, 25, 653–659. [Google Scholar] [CrossRef]
  10. Roccato, A.; Uyttendaele, M.; Membré, J.-M. Analysis of domestic refrigerator temperatures and home storage time distributions for shelf-life studies and food safety risk assessment. Food Res. Int. 2017, 96, 171–181. [Google Scholar] [CrossRef] [PubMed]
  11. Chlebicz, A.; Śliżewska, K. Campylobacteriosis, Salmonellosis, Yersiniosis, and Listeriosis as Zoonotic Foodborne Diseases: A Review. Int. J. Environ. Res. Public Health 2018, 15, 863. [Google Scholar] [CrossRef] [PubMed]
  12. Andritsos, N.D.; Stasinou, V.; Tserolas, D.; Giaouris, E. Temperature distribution and hygienic status of domestic refrigerators in Lemnos island, Greece. Food Control 2021, 127, 108121. [Google Scholar] [CrossRef]
  13. Hassan, H.F.; Dimassi, H.; El Amin, R. Survey and analysis of internal temperatures of Lebanese domestic refrigerators. Int. J. Refrig. 2015, 50, 165–171. [Google Scholar] [CrossRef]
Figure 1. LogTag model TRIX-8 temperature profile of Polish domestic refrigerators.
Figure 1. LogTag model TRIX-8 temperature profile of Polish domestic refrigerators.
Ijerph 20 02924 g001
Figure 2. (a) Distribution of the mean working temperatures of Polish household refrigerators, (b) Empirical Cumulative Density function of the mean working temperatures of Polish household refrigerators.
Figure 2. (a) Distribution of the mean working temperatures of Polish household refrigerators, (b) Empirical Cumulative Density function of the mean working temperatures of Polish household refrigerators.
Ijerph 20 02924 g002
Figure 3. Box and whiskers plot of domestic storage temperature data.
Figure 3. Box and whiskers plot of domestic storage temperature data.
Ijerph 20 02924 g003
Figure 4. Histogram of the standard deviation of recorded temperature profiles.
Figure 4. Histogram of the standard deviation of recorded temperature profiles.
Ijerph 20 02924 g004
Figure 5. (a) Histogram of theoretical densities and (b) P–P plot for the fitted probability distributions.
Figure 5. (a) Histogram of theoretical densities and (b) P–P plot for the fitted probability distributions.
Ijerph 20 02924 g005
Table 1. Frequency table of observed mean temperatures in Polish household refrigerators.
Table 1. Frequency table of observed mean temperatures in Polish household refrigerators.
Mean Temperature [°C]T < 00 ≤ T < 22 ≤ T < 44 ≤ T < 66 ≤ T < 88 ≤ T < 1010 ≤ T < 12T > 12
Frequency54152620421
Percentiles6.49%5.19%19.48%33.77%25.97%5.19%2.60%1.30%
Table 2. Calculated goodness-of-fit statistics and criteria used for best fit selection.
Table 2. Calculated goodness-of-fit statistics and criteria used for best fit selection.
Goodness-of-Fit StatisticsWeibull DistributionTriangular DistributionLognormal DistributionPert DistributionNormal Distribution
Chi-squared15.630.531.413.912.9
Kolmogorov–Smirnov0.010.210.160.080.08
Cramer–von Mises0.140.820.550.090.08
Anderson–Darling0.874.413.230.600.51
Goodness of fit criteria
AIC379.99396.28402.76382.46378.88
BIC384.68403.30407.45391.83383.57
Bias Factor (Bf)
Mean Bf1.00101.10711.00981.00061.0093
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

Stefanou, C.-R.; Szosland-Fałtyn, A.; Bartodziejska, B. Survey of Domestic Refrigerator Storage Temperatures in Poland for Use as a QMRA Tool for Exposure Assessment. Int. J. Environ. Res. Public Health 2023, 20, 2924. https://doi.org/10.3390/ijerph20042924

AMA Style

Stefanou C-R, Szosland-Fałtyn A, Bartodziejska B. Survey of Domestic Refrigerator Storage Temperatures in Poland for Use as a QMRA Tool for Exposure Assessment. International Journal of Environmental Research and Public Health. 2023; 20(4):2924. https://doi.org/10.3390/ijerph20042924

Chicago/Turabian Style

Stefanou, Constantine-Richard, Anna Szosland-Fałtyn, and Beata Bartodziejska. 2023. "Survey of Domestic Refrigerator Storage Temperatures in Poland for Use as a QMRA Tool for Exposure Assessment" International Journal of Environmental Research and Public Health 20, no. 4: 2924. https://doi.org/10.3390/ijerph20042924

APA Style

Stefanou, C. -R., Szosland-Fałtyn, A., & Bartodziejska, B. (2023). Survey of Domestic Refrigerator Storage Temperatures in Poland for Use as a QMRA Tool for Exposure Assessment. International Journal of Environmental Research and Public Health, 20(4), 2924. https://doi.org/10.3390/ijerph20042924

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