Probabilistic modelling of Escherichia coli concentration in raw milk under hot weather conditions

Climate change is one of the threats to the dairy supply chain as it may affect the microbiological quality of raw milk. In this context, a probabilistic model was developed to quantify the concentration of Escherichia coli in raw milk and explore what may happen to France under climate change conditions. It included four modules: initial contamination, packaging, retailing, and consumer refrigeration. The model was built in R using the 2nd order Monte Carlo mc2d package to propagate the uncertainty and analysed its impact independently of the variability. The initial microbial counts were obtained from a dairy farm located in Saudi Arabia to reflect the impact of hot weather conditions. This country was taken as representative of what might happen in Europe and therefore in France in the future due to climate change. A large dataset containing 622 data points was analysed. They were fitted by a Normal probability distribution using the fitdistrplus package. The microbial growth was determined across various scenarios of time and temperature storage reflecting the raw milk supply-chain in France. Existing growth rate data from literature and ComBase were analysed by the Ratkowsky secondary model. Results were interpreted using the nlstools package. The mean E. coli initial concentration in raw milk was estimated to be 1.31 [1.27; 1.35] log CFU/ mL and was found to increase at the end of the supply chain as a function of various time and temperature conditions. The estimations varied from 1.73 [1.42; 2.28] log CFU/mL after 12 h, 2.11 [1.46; 3.22] log CFU/mL after 36 h, and 2.41 [1.69;3.86] log CFU/mL after 60 h of consumer storage. The number of milk packages exceeding the 2-log French hygiene criterion for E. coli increased from 10% [8;12%] to 53% [27;77%] during consumer storage. In addition, the most significant factors contributing to the uncertainty of the model outputs were identified by running a sensitivity analysis. The results showed that the uncertainty around the Ratkowsky model parameters contributed the most to the uncertainty of E. coli concentration estimates. Overall, the model and its outputs provide an insight on the possible microbial raw milk quality in the future in France due to higher temperatures conditions driven by climate change.