Next Article in Journal
Diversity of a Lactic Acid Bacterial Community during Fermentation of Gajami-Sikhae, a Traditional Korean Fermented Fish, as Determined by Matrix-Assisted Laser Desorption/Ionization Time-of-Flight Mass Spectrometry
Previous Article in Journal
Improvement of Intramuscular Fat in longissimus Muscle of Finishing Thai Crossbred Black Pigs by Perilla Cake Supplementation in a Low-Lysine Diet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cold Chain Food and COVID-19 Transmission Risk: From the Perspective of Consumption and Trade

1
Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
3
Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
4
Key Laboratory of Carrying Capacity Assessment for Resource and Environment, Ministry of Natural Resources, Beijing 100101, China
*
Author to whom correspondence should be addressed.
Foods 2022, 11(7), 908; https://doi.org/10.3390/foods11070908
Submission received: 2 December 2021 / Revised: 8 February 2022 / Accepted: 15 March 2022 / Published: 22 March 2022
(This article belongs to the Section Food Security and Sustainability)

Abstract

:
Since the outbreak of the coronavirus disease 2019 (COVID-19), political and academic circles have focused significant attention on stopping the chain of COVID-19 transmission. In particular outbreaks related to cold chain food (CCF) have been reported, and there remains a possibility that CCF can be a carrier. Based on CCF consumption and trade matrix data, here, the “source” of COVID-19 transmission through CCF was analyzed using a complex network analysis method, informing the construction of a risk assessment model reflecting internal and external transmission dynamics. The model included the COVID-19 risk index, CCF consumption level, urbanization level, CCF trade quantity, and others. The risk level of COVID-19 transmission by CCF and the dominant risk types were analyzed at national and global scales as well as at the community level. The results were as follows. (1) The global CCF trade network is typically dominated by six core countries in six main communities, such as Indonesia, Argentina, Ukraine, Netherlands, and the USA. These locations are one of the highest sources of risk for COVID-19 transmission. (2) The risk of COVID-19 transmission by CCF in specific trade communities is higher than the global average, with the Netherlands–Germany community being at the highest level. There are eight European countries (i.e., Netherlands, Germany, Belgium, France, Spain, Britain, Italy, and Poland) and three American countries (namely the USA, Mexico, and Brazil) facing a very high level of COVID-19 transmission risk by CCF. (3) Of the countries, 62% are dominated by internal diffusion and 23% by external input risk. The countries with high comprehensive transmission risk mainly experience risks from external inputs. This study provides methods for tracing the source of virus transmission and provides a policy reference for preventing the chain of COVID-19 transmission by CCF and maintaining the security of the global food supply chain.

1. Introduction

The recent COVID-19 event related to frozen foods has caused public concern about the cold food chain (CCF) [1]. Many studies have focused on the impact of the pandemic on the food supply chain [2] and the external supply of food [3]. Since the outbreak of COVID-19, there have been continuous studies of the transmission path [4]. Air-borne transmission through breathing is a widely accepted route [5]. To our knowledge, the study of the seafood market in Wuhan, China initially considered seafood to be a potential carrier of viruses, and SARS-CoV-2 has been detected in imported CCF [6,7]. As a result, the questions of whether SARS-CoV-2 is foodborne and whether the food supply chain is a transmission path for the virus have aroused public concern and scholarly attention [8].
Testing for COVID-19 is a challenge because there is a lack of reliable, inexpensive, and easy-to-use detection tools, specifically for food [9]. Laboratory studies have shown that COVID-19 can remain highly stable under refrigeration at 4 °C, and in freezing conditions from −10 to −80 °C, on fish, meat, poultry, and swine skin, for 14–21 days [10,11,12]. From July to August 2020, at least 9 cases of live SARS-CoV-2 were reported being detected on the surface of frozen food in China, triggering a new round of COVID-19 spread in Beijing and Qingdao [11]. Particularly, in the Beijing Xinfadi epidemic (June, 2020), 45 people in close contact with the CCF environment were tested positive in Xinfadi Market [1]. In the Qingdao epidemic (September 2020), two positive patients are the stevedores engaged in CCF, both carried out loading and unloading of frozen cod in bulk. Gene sequencing results show that: the virus genome sequence of the stevedores is highly homologous to the original piece from the frozen cod outer package, with two distinct nucleotides. Once the destination is reached, the COVID-19 may infect stevedores or other staff working in the port through direct contact with the goods according to Liu and his colleague [13]. Lately, in the latest round of the COVID-19 in Dalian, China (in early November 2021), 16 of the 45 positive patients were workers in designated cold storage for imported CCF.
These findings indicate that frozen food is an important and potential source of the outbreak of COVID-19 [14]. In China, many batches of CCF have been detected for surviving COVID-19 (Table S1). At least eight cases of the COVID-19 spread are related to the imports of CCF in 2020. All the zero patients had contact with the contaminated cold chain food and they were all cold chain practitioners. All the SARS-CoV-2 sources can be traced from the outer packages or container samples of the cold chain (Table 1). In some countries with significant COVID-19 transmission rates, such as in the USA, Canada, Brazil, Germany, and Ireland, many production and processing workers in the CCF environment have been infected and have transmitted COVID-19 [15,16,17,18]. If there is COVID-19 infection in a certain part of a processing plant among processing personnel, cross-contamination in the CCF supply chain environment may occur, increasing the probability of COVID-19 transmission [18,19,20] (Table S2).
Meanwhile, there is also evidence of a correlation between the survival environment of COVID-19 and the cold chain transportation environment in international trade [10,21], providing environmental conditions that could allow SARS-CoV-2 to be transported through the international cold chain [17]. COVID-19 may be transported and exported on frozen food packaging over a long period of time. Cold chain transportation has the environmental attribute of low temperatures throughout the process, allowing COVID-19 to remain stable, prolonging its survival time, and increasing the difficulty of elimination [22]. There is research evidence that cold chain transportation in the frozen food industry may have caused a recurrence of COVID-19 cases in the destination (Table S2) [13,23,24]. In high-risk areas of the COVID-19 pandemic, all links of the CCF production and marketing chain could be contaminated by people with latent and asymptomatic infections [25].
Based on the existing research and our analysis, we hypothesize that meat and meat products, dairy products, fruits, vegetables, ready-to-eat food, and other cold chain transport food may be potential carriers and transmitters of COVID-19 [24,26,27]. Different parts of the CCF supply chain, such as harvesting, processing, transportation, and retail, have physical environments that may spread infections. Contaminated food in the cold chain has become a focus of COVID-19 transmission research in food chain circulation [11] and reflects the possibility that COVID-19 can be transmitted from farm to table, highlighting a source of virus transmission between countries.
The continuous increase in globalization has increased the intensity of food trade links between countries, leading the global food system to become increasingly integrated. There has been a gradual increase in the dependence of countries with weak food production capacity on external food supplies, with an increased demand for transnational long-distance food transportation [28]. Urbanization leads to the formation of large-scale urban agglomerations. As a result, local agricultural systems cannot provide enough food to meet the demand, and cities increasingly rely on food distribution systems outside the region. Cross-regional food allocation has become important in supporting the food system of big cities.
The acceleration of globalization and urbanization has also changed food consumption demands, including the rise of fast-food culture (e.g., McDonald’s and KFC), and an increase in the demand for fresh food (e.g., fruits, vegetables, and aquatic products). Food consumption in different places and the increased demand for the convenient preservation of food had led to increased requirements for food transportation, advancing the development of cold chain transportation [29]. Cold chain transportation, with its properties of low temperature and preservation, plays an important role in local-to-local food flow and global food trade and plays an increasingly important role in daily life [22]. It has become an indispensable part of the food supply chain, particularly in urban areas. As a result, the spread of the virus through the cold chain between and within countries could heighten local food safety risks around the world across the trade network [30].
In countries with weak COVID-19 detection capability, the transmission of COVID-19 through CCF may not be found in time and may not receive adequate attention [20]. Some scholars believe that cold chain food is a potential and non-negligible transmission route in the occurrence and development of COVID-19 [31,32], but few studies on the risk of viral transmissions through the CCF trade [33]. Despite this, the scientific identification of global CCF COVID-19 transnational infection sources, and the risk of transmissions in the global CCF chain, is important for cutting off COVID-19 transmission through CCF and maintaining the security of the food supply chain [26,34]. Therefore, by applying the perspective of CCF consumption levels and the global CCF trade network (CCFTN), this study explores the risk level of COVID-19 transmission by CCF. The study combines the COVID-19 risk index (CRI), CCF consumption level, and trade level of different countries to address the following questions: (1) What are the trade network patterns and core node characteristics of the global CCFTN? (2) What is the “source” of the intensity pattern of COVID-19 transmission by CCF in the trade network? (3) What is the risk level of SARS-CoV-2 transmission by CCF in different countries? What are the differences in dominant risk factors? What are the regional differences?
Based on the global CCF consumption quantity and trade flow under general conditions, as well as the current COVID-19 risk of different countries, here, CCF consumption quantity and trade flow pattern and the risk of COVID-19 transmission by CCF were studied. We used consumption and trade matrix data from the Food and Agriculture Organization of the United Nations (FAO), specifically on meat, eggs, milk, vegetables, and fruits. Nutritive conversion factors (NCF) were used to convert the trade quantity into energy (kcal) equivalent. Based on an analysis of the CCF consumption level, the complex network analysis (CNA) method was used to analyze the CCF trade network features and to identify the core nodes. To identify the risk “source” of COVID-19 transmission by CCF, a COVID-19 transmission model for CCF was constructed, including the risk level, consumption level of CCF, urbanization rate (UR), and trade quantity of CCF in different countries. The risk level of CCF transmission was evaluated across the globe—trade communities—country dimensions, and the dominant risk types of different countries were identified. This study provides a scientific basis for preventing COVID-19 transmission through the cold chain, scientifically preventing and controlling the risk of transmission, and protecting the security of the global food supply chain.

2. Methods and Data Processing

Considering the quantity of food consumption in cold chain and the flow of international trade network, this study builds on the following research framework (Figure 1), including analysis of CCFTN characteristics—, evaluation of COVID-19 diffusion source intensity—, Assessment of COVID-19 transmission risk level—, classification of dominant risk types, and—Policy recommendations”. The study starts with two perspectives, i.e., internal diffusion and migration diffusion. First, the internal diffusion risk (IDR) was evaluated based on the UR, CCF consumption level [35], and domestic CRI. Second, the CNA method was used to analyze the characteristics and community structure pattern of global CCFTN; this method has been widely applied to study trade characteristics and virus spread [36,37]. The global migration diffusion source intensity was then analyzed, based on the out-strength degree and CRI level. Third, the external input risk (EIR) was evaluated based on CCF imports, the number of partners, and the CRI of partner countries. Based on the IDR and EIR, the comprehensive risk assessment and dominant risk type identification were determined using different dimensions: globe, community, and country. Finally, proposes adaptive strategies to prevent the risk of COVID-19 transmission through CCF and maintain the security of the global food supply chain.

2.1. Methods

2.1.1. Cold Chain Food Trade Network Analysis

(1)
Network properties
Based on the CCF trade network, we analyzed the characteristics of the CCFTN using three indicators that are closely related to COVID-19 transmission: node degree, node strength, and betweenness centrality.
Node degree ( K i ) measures the number of contacts maintained by each node. Here, it represents the number of trade partners a country has. The node in-degree counts the number of links incoming to a node; it is measured by
K i i n = j = 1 N a j i
The node out-degree counts the number of links emanating from a node; it is measured by
K i o u t = j = 1 N a i j
where a is an element of CCFTN. Generally speaking, the higher the K i i n and S i i n of a country, the higher the risk of migration diffusion risk. The higher the K i o u t and S i o u t of a country, the higher the intensity of transnational diffusion risk.
Node strength ( S i ) is the weighted corollary to the node degree. It measures the sum of the weights for nodal links. Node in-strength sums the value of links incoming to a node; it is measured by
S i i n = j = 1 N w j i
The node out-strength sums the value of links emanating from a node; it is measured by
S i o u t = j = 1 N w i j
where w is an element of weighted CCFTN. If a country has a higher value of K i o u t and S i o u t , it indicates the country has a larger number of export partners and volumes. By extension, the country has a higher level of influence in the CCFTN. In contrast, a country with a high value of K i i n and S i i n indicates that the country has a larger number of import partners and volumes [38].
The Betweenness Centrality (BC) of a node is defined as a number of shortest paths through the node. BC is measured as follows:
BC i = 1 N 2 i , t n s i i g s i
In the equation, g s i is the number of shortest paths from node s to node i; and n s i i is the number of shortest paths through node i in the g s i shortest paths from node s to node i. A higher BC value for a node indicates that the node is more important in the CCFTN, and that there is a higher level of COVID-19 transmission efficiency by CCF [37].
(2)
Community detection
Several independent and highly connected modules (i.e., communities) generally form in a complex network; the connection strength between countries in the same community is higher than the connection strength with different communities. To better analyze the source of EIR, we identified the community structure of the CCFTN. The specific calculation formula of modularity is as follows:
Q = 1 2 m i j ( A i j k i k j 2 m ) ( c i c j )
In the equation, c i and c j represent the communities in which the i and j nodes are located, respectively. The variable represents the binary functions for estimating whether two points are in the same community. If c i = c j , then the value is 1; otherwise, the value is 0. The variable A i j is the weight of the connection between nodes i and j. The expression
k i = j A i j
is the sum of all connection weights that contain node i. The expression
m = 1 2 i j A i j
is the total contact weight of the entire network, where the modularity is a standardized index, and the value interval is (−1, 1). In this study, the Modularity index is used to evaluate the degree of community separation. Q represents the Modularity index, and its value ranges from (0, 1). The larger the value, the more obvious the differentiation of network community and the higher the quality of community division [39].

2.1.2. Analysis of the Intensity of COVID-19 Transmission Sources

We constructed the intensity of the COVID-19 transnational diffusion index (CTDI) to evaluate the probability that COVID-19 occurs as an output. This is defined as follows:
C T D I i = S i o u t × C R I i  
If a country has a high value of C T D I , it indicates that the country has a higher probability and range of COVID-19 transmission through the CCFTN to other countries.

2.1.3. Risk Assessment Model of COVID-19 Migration Diffusion

As noted above, COVID-19 transmission risk was analyzed from the perspective of internal diffusion and external input. The internal diffusion risk index (IDRI) is defined as:
I D R I i = C L C C F i × U R i × C R I i
where C L C C F i represents the CCF consumption level of country i; U R i represents the UR of country i; and C R I i represents the COVID-19 risk index of country i. Generally speaking, the higher the level of urbanization, the higher the demand for CCF transfer and transportation. The greater the consumption of CCF, the higher the probability that residents will be exposed to CCF. Therefore, the higher the level of urbanization, CCF consumption, and the risk of the COVID-19, the higher the internal diffusion risk of COVID-19 transmission by CCF.
Next, the external input risk index (EIRI) is defined as:
E I R I i = j I Q j × C R I j
where I Q j represents the quantity of CCF imported from country j to country i; and C R I j represents the COVID-19 risk index of country j. Correspondingly, the greater the import volume, the higher the risk of the COVID-19 risk in the country of the import pattern country, and the higher external input risk. Next, we build a comprehensive transmission risk index (CTRI) using IDRI plus EIRI, as follows:
C T R I i = I D R I i + E I R I i
The range in values of IDRI and EIRI strongly and equally influence CTRI; therefore, we applied the Jenks breaks in ArcGIS 10.5 to divide IDRI and EIRI into 10 levels, assigning the numbers 1–10 according to the level from low to high.
Finally, the contribution rate of IDRI and EIRI to CTRI was used to classify the dominant risk types. This was calculated based on the proportion of IDRI and EIRI scores to CTRI. If the contribution rate of IDRI was more than 50%, the country was classified as having an internal diffusion risk type (IDRT). If the contribution rate of EIRI was more than 50%, the country was classified as having an external input risk type (EIRT). If the contribution rate of IDRI and EIRI is equal, the country is classified as having a double risks type (DRT). The criteria for classification of all different elements (i.e., IDRI, EIRI, CTRI, CFF Consumption, and UR) above list in Table 2.

2.2. Data Processing

2.2.1. Cold Chain Food Data

Our study selected 175 items in these five foods (see Table S3), namely meat, eggs, milk, vegetables, and fruits. In order to make the statistical data of consumption quantity and trade flow of different types of CCF comparable as well as better reflect the spatial differences of consumption level and trade flow of CCF, the nutritive factors were used. Through nutritive factors, the consumption and trade volume of CCF are uniformly converted into energy (kcal) to reflect its quantity [40]. Since the global cold chain food consumption-trade data have not yet been released during the COVID-19 pandemic, we use the latest data that can be found, which are the data from 2018, to present the global CCF consumption and trade situation under normal circumstances. The metadata of food trade and consumption were collected from the Trade and Food Balance database of FAOSTAT (http://www.fao.org/faostat/en/#data, accessed on 20 May 2021). In particular, missing data were replaced and supplemented by the conversion of import and export data, under the premise that one country’s export is another country’s import. For example, when querying import data, if the import data for a reporting country were missing for the partner country, the missing data were collected by querying the export data of other partner countries exporting to that reporting country (e.g., the FAOSTAT “Partner” and “Reporter” countries). After this data supplementation step (needed for about 30 countries), trade data for 198 countries were finally obtained. Network analysis and visualization were conducted using Gephi (version 0.9.2, https://gephi.org/, accessed on 1 June 2021).

2.2.2. COVID-19 Risk Index Data

The COVID-19 risk index data were collected from the website at the http://covid19-risk-index.com/ (accessed on 15 May 2021); the data were compiled by teams at the Huashan Hospital of Fudan University in China and from other institutions. To mitigate the impact of data fluctuations, the average values from 1 January to 30 April 2021, were used If data were missing for the countries in the study sample, we extrapolated the data based on the average of their sub-regions (using United Nations classification standards). Finally, the data of UR (urban population of the total population) in 2018 were collected from the World Bank.

3. Results

3.1. Consumption Level and International Trade Network Characteristics of Cold Chain Food

3.1.1. Consumption Level and Urbanization Rate

From the perspective of CCF consumption patterns, Figure 2a shows that 31 countries belong to very high consumption levels, such as European countries, the USA, Argentina, Australia and some European countries. A total of 51 countries are classified as having a high consumption level, mainly distributed in Americas and Asia. Moreover, a total of 116 countries are classified as having a medium consumption level and below, mainly in Africa, West Asia, Southeast Asia, and Southeast Asia. From the perspective of the urbanization level and distribution pattern, the countries in Europe and America are classified as very high and high; the countries in Southeast Asia, West and Central Africa, and the Southeast Asian islands are generally classified as having a medium level; and the land countries in South Asia and Southeast Asia are generally classified as having a low and very low level (Figure 2b).

3.1.2. International Trade Network and Community Structure

(1)
International trade network characteristics
The CCFTN has 198 nodes and 11,924 lines, with an annual trade quantity of 80,249.05 × 1010 kcal. The probability distribution fitting results illustrating different attributes of the nodes in the cold chain transportation network show that the out-strength distribution of nodes in CCFTN has typical scale-free distribution characteristics with strong heterogeneity (Figure 3a). Export core nodes play a leading role in the structure and normal operation of the trade network.
In terms of trade flow (Table 3 and Figure 4), the top 20 countries in the CCFTN account for 66.06% of the export quantities; with the exception of the out-degree of Papua New Guinea, Bolivia, and Belarus, which are lower than 70, the out-degree of other top 20 countries exceed 100. According to Table 3, the top 20 countries account for 77.40% of the world’s total out-strength, with the top 10 countries accounting for 56.86% of the globe’s totality. With respect to the trade quantity and out-strength, the top 10 countries affect the CCFTN pattern, given their large export quantities and relatively extensive external links. These countries also form the core node of the trade network and the transition point of COVID-19 risk. The average BC value for the top 20 countries is 902.68 (Table 3 and Figure 3b), which is 6.59 times the global average; and the average value of the out-degree is 164, which is 2.72 times the global average. The countries with high and middle BC levels have an important position in a trade hub, and profoundly impact the global cold chain transport trade network structure. Therefore, such countries are more likely to become potential nodes for transmitting COVID-19.
(2)
Community structure
The results of community structure analysis show that the CCFTN includes 6 trade communities, dominated by big exporting countries (Figure 4). The total quantity of intra-community trade is 48,477.35 × 1010 kcal, accounting for 60.37% of the trade across the total network (Table 4). Community division results and internal energy flow characteristics are as follows (Figures S1 and S2). (1) The quantity of trade for the NLD-DEU community accounts for 28.43% of the total CCFTN. The Netherlands, Germany, Spain, Belgium, and France are the main exporting countries; CCF mainly flows to Germany, the Netherlands, Italy, France, Belgium, and other European countries, and also radiate to parts of Africa. (2) The quantity of trade for the IND-UKR community accounts for 15.07% of the total CCFTN. CCF mainly flows from Indonesia, Ukraine, Malaysia, and New Zealand to India, China, and Malaysia, radiating to Southeast Asia, Australia, and Africa. (3) The quantity of trade for the USA-CAN community accounts for 8.54% of the total CCFTN. The USA and Canada are the main exporting countries. CCF mainly flows into each other between the USA and Mexico as well as Canada. South Korea, Japan, and other American countries are important CCF export destinations. (4) The quantity of trade for the RUS-ARG community accounts for 6.33% of the total CCFTN. Russia, Argentina, Brazil, and Belarus are the main CCF outflow sources. CCF flows mainly to Russia, Algeria, Egypt, Iran, Norway, Venezuela, and some South American, Central Asia, and West Asia countries. (5) The quantity of trade for the TUR-ARE community accounts for 1.85% of the total CCFTN. CCF mainly flows from the United Arab Emirates, Turkey, and Saudi Arabia to Iraq, Yemen, the United Arab Emirates, and other Western Asian countries. (6) The quantity of trade for the SRB-MKD community accounts for 0.15% of the total CCFTN, with generally singular flows of CCF from Serbia to the other four countries.

3.2. Risk Assessment of COVID-19 Transmission by Cold Chain Food

The results of transnational diffusion risk intensity evaluation showed that 5 countries are associated with very high-risk levels (Figure 5), namely, Indonesia, Argentina, Ukraine, Netherlands, and the USA. There are 8 countries classified at high risk level, and 11 countries are at a medium risk level; 27 and 147 countries are classified as having low and very low risk levels, respectively, accounting for 13.64% and 74.24% of the total (see Figure 5 and Figure S3).

3.3. COVID-19 Transmitted by Cold Chain Simulation

3.3.1. Spatial Pattern of Internal Diffusion Risk

The results of the IDR assessment show that the average IDRI is 24.63 globally, with significant differences among trade communities and countries. For specific trade communities (Table 4), the averages of the SRB-MKD and NLD-DEU IDR indices are about 1.5 and 2 times the global average, respectively; the IDR within the communities is high. The averages of the USA-CAN, TUR-ARE, and ARG-RUS community IDR values are 25.93, 25.16, and 23.32 respectively, which are close to the global average; the IDR within the communities is not high. The IDR value in the IDN-UKR community is only 11.12, which is lower than the global level.
From the perspective of national scale (Figure 6a and Table 5), 19 countries are classified as having a very high IDR level, mainly in Europe, but also including the USA, Argentina, Brazil, Oman, Kuwait, and other countries. A total of 37 countries are classified as having a high IDR level, mainly in Europe. There are 40 countries with medium IDR levels and 49 with low IDR levels scattered in Asia and Africa, respectively. There are 53 countries with very low levels of IDR, mainly distributed in Africa, but also including India, Pakistan, and some Southeast Asian countries.

3.3.2. Spatial Pattern of External Input Risk

The results of the EIR assessment show that the average global EIRI is 23.01. The EIRI of the NLD-DEU and USA-CAN communities are higher than the global average level, at 36.57 and 23.81, respectively. The EIR indices of the IDN-UKR, ARG-RUS, and TUR-ARE communities range between 15–20, lower than the global average. The EIRI of the SRB-MKD community is only 4.63, reflecting a low level of external risk (Table 4).
There are 5 countries (i.e., Indonesia, China, the Netherlands, the USA, and Germany) with high EIRI; this also includes the major CCF importing countries and BC core countries. There are 7 countries with high-level EIRI, including Italy, Britain, France, Belgium, Spain, Mexico, and Bangladesh. Countries with medium and low EIR levels account for 12.12% and 16.67% of the total number of countries, respectively. The countries with medium risk are mainly distributed in West Asia and South America. The countries with low-level risk mainly include European countries. Approximately 65% of the countries in the world have a very low EIR level; these are mainly distributed in Africa (Figure 6b and Table 5).

3.3.3. Spatial Pattern of Comprehensive Risk

The CTR assessment showed that the average value of global CTR index is 9.24. For each community, except for the IDN-UKR the average of the other five communities (i.e., the NLD-DEU, the USA-CAN, the ARG-RUS, the TUR-ARE, the SRB-MKD) is higher than the global average. Of these 5 communities, the NLD-DEU group has the highest CTR index, followed by SRB-MKD; the CTR index of the two communities was about 11.5. These are followed by TUR-ARE, USA-CAN, and ARG-RUS, with CTR indices ranging between 9.5 and 9.9 (See Table 4).
From a country-level perspective, 11, 45, 51, 51, and 40 countries have CTR indices from very high level, to high, to middle, to low, to very low level, respectively. The few countries with a very high CTR level include the Netherlands, the USA, Germany, Belgium, France, Spain, Britain, Italy, Mexico, Poland, and Brazil. There is little difference in the indices across the other CTR levels. In terms of spatial distribution, European and American countries generally have very high and high CTR levels (Figure 7a). The CTR index of countries in Africa is low, with the exception of Algeria and South Africa. In Central Asia, South Asia, and Southeast Asia, the CTR level is medium or lower. Some countries in Western Asia, China, the Philippines, Malaysia, and Australia are at a high CTR level (See Figure 7a and Table 5).

3.4. Types of Dominant Risk

The identification of dominant risk types revealed that 61.62% of countries are dominated by IDR, with 22.73% of countries are experiencing primarily external input risks. Also, some a few countries face double risks, and accounting for 15.66% of the total country count. From a spatial distribution perspective (Figure 7b), the IDR types are mainly distributed in Africa, Southeast Asia, Central Asia. In contrast, North America, South America, and most parts of Europe are dominated by external risk. A few countries in Eastern Europe, North Africa, and South Asia are dominated by double risks.

4. Discussion

There have been instances of COVID-19 spread related to CCF events, and there is a chance that CCF is a carrier that facilitates the transmission [14]. This study uses 2018 global matrix data of CCF consumption and trade, NCF, and a CAN method to identify the core nodes and communities of the global CCFTN. The study also explored the potential source intensity for COVID-19 transmission by CCF. The risk assessment model evaluated the extended diffusion and migration diffusion risk of COVID-19; the variables included the CRI, UR level, CCT consumption level, CCF import volumes, and CRI of exporting countries. The risks associated with COVID-19 transmission based on CCF trade were evaluated at three different dimensions i.e., national, and global scales, as well as community level. This led to the identification of the dominant risk types.
The main contributions of this study are as follows. (1) CNA method was applied to analyze global CCFTN patterns and core nodes. (2) The risk of COVID-19 through CCF was analyzed from two perspectives: local-to-local extended diffusion, and country-to-country migration diffusion. The dominant risk type was also identified for different countries. Our research could provide a scientific basis for preventing the transmission COVID-19 by CCF and maintaining the global food supply chain and food safety. Also, another contribution of this study is that it focused on the unique context of the outbreak events related to CCF. CCF has special requirements related to the storage environment, making it a unique research object. Within this unique context, factors such as UR, the consumption level of CCF, the import quantities of CCF, and the CRI of importing partner countries were all considered to comprehensively evaluate the risk of COVID-19 transmission in CCF. In contrast to previous studies on cereal or commodity trade, the study describes global CCF trade pattern characteristics from the perspective of physical quantities of food, using the NCF and CNA method.
Our study found that Indonesia, Argentina, Ukraine, Netherlands, the USA, and other countries have a high level of intensity with respect to CCF migration diffusion sources. The NLD-DEU community, USA-CAN community, countries importing CCF, and core countries with a high BC level (e.g., Indonesia, China, Netherlands, the USA, and Germany) face a high EIR. There is a high level of CCF consumption and urbanization in Europe and North America. The SRB-MKD and NLD-DEU communities, European countries, the USA, Argentina, Brazil, Oman, and Kuwait, face high IDR.
In terms of the CTR, the CTR level of most trade communities is higher than the global average, with the NLD-DEU community having the highest levels. The CTR levels of 11 countries are at a very high level, including 8 European countries (namely, Netherlands, Germany, Belgium, France, Spain, Britain, Italy, and Poland) and 3 American countries (i.e., the USA, Mexico, and Brazil). The European region faces a high CTR level, with most countries dominated by ETI. Globally, IDR dominates the transmission risk of COVID-19 by CCF in many countries (accounting for about 62% of all countries, mainly in African countries), followed by the countries with high levels of EIR (accounting for about 23% of the countries); fewer countries faced both risks (internal and external) (approximately 15% of the countries).
There are three key results highlighted here. Firstly, the CCFTN has typical scale-free network characteristics, and the stability of the trade network is significantly affected by the large trading countries. Given this, stable CCF output from large exporting countries significantly impacts the CCFTN. The COVID-19 pandemic negatively affected global food production, processing, and transportation. Stable CCF exports also significantly impact the maintenance of the global CCF supply chain; while large CCF exports also increase the probability of COVID-19 transmission [7]. Therefore, the large CCF exporting countries have become the key nodes for maintaining global food trade security and transnational export controls during the COVID-19 pandemic. At the same time, due to the functions served by transfer hubs, the countries with high and middle BC levels may become transfer hubs for COVID-19 during food importation. This occurs through the import and export of contaminated processed food. This is risk associated with countries having a BC value more than 1000, which includes the USA, the Netherlands, Britain, Malaysia, Canada, France, Belgium, and China. These countries are key nodes of food transfer in the global cold chain and could be core hubs of COVID-19 transnational spread—and for cutting the viral transmittal off.
Secondly, after years of development, CCF has become an important part of the global food supply chain, and the security of CCF impacts global food security. Once the pandemic began, it was inevitable that COVID-19 would impact the global food supply chain [7]. The world now faces increased risks from farm and food supply chain interruption, with increased risks of unemployment and food insecurity in food-related fields. The interruption of the food supply chain also increases global food waste. Cold chain technology reduces food loss and improves the efficiency of the food supply chain and food quality standards [35]. During the COVID-19 pandemic, managing CCF inventories has become an important way to address food shortages [41,42,43].
However, cold chain technology also increases the scope of food trade circulation, improves food accessibility and availability. During the COVID-19 pandemic, quarantine protocols have increased people’s demand for CCF in their homes [44]. Compared with non-cold food chain transportation, CCF faces the dual pressure of controlling COVID-19 transmission and maintaining a smooth food supply chain [11]. Therefore, strengthening monitoring of the virus and maintaining the normal circulation of CCF is of great significance for maintaining food security. This deserves increased attention from policymakers [45].
Thirdly, this study found that internal diffusion risk is the dominant risk for many countries with respect to possible COVID-19 spread by CCF. It is vital to strengthen virus monitoring and daily protection of CCF workers to prevent and control COVID-19. The virus can be transmitted through respiration, contact, and other routes; as such, food contamination may occur during CCF production, processing, storage, transportation, and retail. As such, virus infection may occur at any point in the CCF life cycle, from farm to fork [11]. This highlights the possibility that the transmission chain may extend from human infection to food pollution back to human infection.
In terms of environmental factors (cold and humidity) and operational design (workers standing side-by-side on crowded and noisy processing lines), the situation in meat processing plants has been shown to facilitate the spread of the virus (Table S2). The crowded living conditions of temporary foreign workers in fruit and vegetable production operations can also support virus transmission and spread [46,47]. The risk associated with CCF is higher compared to non-CCF because of these unique operational variables. The COVID-19 pandemic has increased the risk of global unemployment, and cold chain workers face higher risks of infection while also facing unemployment risks [48]. They may also not report medical concerns, exacerbating the risk of virus transmission [49]. This also deserves more attention and again highlights the need to strengthen monitoring tools for detecting COVID-19 transmission in the cold chain environment, and to introduce daily prevention and control practices for workers involved in the cold chain [4]. These steps would help maintain the normal operation of the cold chain food system and prevent and control future transmissions.
Countries with a strong source intensity value with respect to COVID-19 migration diffusion should work to strengthen their load-level monitoring of exported CCF, to prevent transmission through contaminated food. Countries with a high EIRI that import CCF should also strengthen their monitoring of imports to prevent COVID-19 transmission. We acknowledge that this testing recommendation is challenging, due to the lack of a reliable, inexpensive, and easy-to-use COVID-19 detection tool specifically for food [9]. We recommend, therefore, that international organizations develop guidelines for COVID-19 detection in CCF and apply inter-disciplinary methods to develop new tools of prevention [5]. It is also important to assist countries that have an insufficient detection capacity to gradually eliminate COVID-19 transmission.
There are five countries (i.e., Indonesia, Argentina, Ukraine, Netherlands, and the USA) with a very high level of diffusion risk intensity. Those five countries are the countries of global concern in preventing COVID-19 transmission by CCF trade networks. The country with high import quantity and in-degree and high BC (e.g., India, China, Netherlands, the USA, Germany) facing high-level external input risk. If lack of detection ability of COVID-19, those countries can easily become a transmission hub for global COVID-19 diffusion by CCF trade network. So, those countries are the countries of special concern in COVID-19 diffusion in the world.
Most countries with the high consumption level of CCF are faced with high comprehensive transmission risk, and most of them belong to external input risk types (mainly European countries). On the one hand, these countries need to strengthen the COVID-19 prevention and control in the processing, transportation, and preparation of CCF, popularize the knowledge of prevention in CCF, and reduce the risk of infection through contact with CCF. On the other hand, it is necessary to strengthen the port transshipment of CCF and the storage management of imported CCF, cut off the chain of external input risk in time, and maintain the security of CCF and the food supply chain.
For the countries that belong to the internal diffusion risk type (mainly in Africa), their CCF COVID-19 detection ability is insufficient. We suggest that international organizations should give priority to strengthening assistance to enhance the COVID-19 detection intensity of their CCF and reduce the risk of internal diffusion.
Also, some relevant studies showed that the lack of reliable, cheap, and easy-to-use COVID-19 detection tools for food was still one of the challenges [9,50,51], and we then suggest that global medical institutions and enterprises strengthen the research and development of cold chain food virus monitoring. We also believed that policymakers formulate technical guidelines and disinfection for SARS-CoV-2 prevention in processing-packaging-loading and uploading(stevedore)-retailing-cooking of CCF, to guide residents to scientifically handle imported CCF in daily life and reduce the risk of COVID-19 epidemic caused by exposure to CCF.
It is worth nothing that our study had two main methodological limitations. (1) The data were limited to the data in the FAO Database. As such, we included only 5 kinds of food that need cold chain transportation. These foods do account for most CCF and are generally representative of that food category. However, we were not able to obtain good data about aquatic product trade quantities. As such, future studies could focus on aquatic products as the research object. Future research could mine aquatic data and physical conversion parameters to study the risk associated with SARS-CoV-2 transmission in aquatic products.
It should be noted that our study did not take the trade data after the COVID-19 outbreak into consideration due to accessibility. It is pointed out that after the outbreak, the international trade will be affected in a short time, and the trade volume will decrease, while the domestic CCF will increase with the implementation of the isolation policy [2,46,52,53], which may affect the accuracy of the results of this study. This study is based on the general CCF consumption and trade volume. Considering that the COVID-19 pandemic is continuing, a comparative study of CCF and the risk of SARS-CoV-2 transmission needs to be carried out around the COVID-19 Outbreak in the near future.
(2) Different countries are in different development stages of cold chain transportation, and the prevention and control measures for COVID-19 differ across the CCF flow process. Although the comprehensive COVID-19 risk index (e.g., CRI) covers the factor of prevention and control measures, the research did not more granularly detect COVID-19 in the CCF process flow. This may affect the accuracy of the assessment results. Additional industry data about cold chain industry development could improve the accuracy of the risk index and improve the risk assessment of COVID-19 transmission by CCF. In addition, future studies would be helpful in validating that the incidences of COVID-19 transmission have occurred across the hypothesized communities, to validate or refute the model in real-life contexts. The lack of actual testing to validate the model is an acknowledged research gap that deserves further exploration [9,50,51]. As we all know, the difference of food matrix will also affect the survival time of virus [26]. Future research can focus on the spread of CCF viruses from the perspective of food matrix.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/foods11070908/s1, Figure S1: Flow of the emery within the global CCFTN communities. A total of 198 countries are ranked according to the total trade quantity and plotted clockwise in descending order. The size of the out bar indicates the total trade quantity (unit: 106 kcal). Export quantity is indicated with links emanating from the out bar of the same color. (a) NLD-DEU community, (b) IDN-UKR community, (c) USA-CAN community, (d) ARG-RUS community, (e) TUR-ARE community, (f) SRB-MKD community, Figure S2: Spatial pattern of trade community of cold chain foods trade networks (CCFTN), Figure S3: Spatial pattern of the COVID-19 risk index (CRI), Table S1: Imported cold-chain food tested positive for SARS-CoV-2 in China in 2020, Table S2. Cold-chain food tested positive for SARS-CoV-2 in countries in 2020, Table S3: The cold chain food products and the calorie factors, Table S4: The country name and ISO3 code.

Author Contributions

C.Z. wrote the manuscript text and contributed to the graphics; Y.Y., Z.F., C.X., T.L., Y.L. and X.S. contributed to the revision of the methods, data processing, results, and discussion of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA20010201, XDA20010203) and the Chinese Academy of Sciences (ZDRW-XH-2021-3).

Institutional Review Board Statement

Not available.

Informed Consent Statement

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

Data Availability Statement

The data presented in this study are available on request to authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pang, X.H.; Ren, L.L.; Wu, S.S.; Ma, W.T.; Yang, J.; Di, L.; Li, J.; Xiao, Y.; Kang, L.; Du, S.C.; et al. Cold-chain food contamination as the possible origin of COVID-19 resurgence in Beijing. Natl. Sci. Rev. 2020, 7, 1861–1864. [Google Scholar] [CrossRef] [PubMed]
  2. Laborde, D.; Martin, W.; Swinnen, J.; Vos, R. COVID-19 risks to global food security. Science 2020, 369, 500–502. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, C.; Yang, Y.; Feng, Z.; Xiao, C.; Lang, T.; Du, W.; Liu, Y. Risk of Global External Cereals Supply under the Background of the COVID-19 Pandemic: Based on the Perspective of Trade Network. Foods 2021, 10, 1168. [Google Scholar] [CrossRef] [PubMed]
  4. Sun, C.; Cheng, C.; Zhao, T.; Chen, Y.; Ayaz Ahmed, M. Frozen food: Is it safe to eat during COVID-19 pandemic? Public Health 2021, 190, e26. [Google Scholar] [CrossRef]
  5. Ceniti, C.; Tilocca, B.; Britti, D.; Santoro, A.; Costanzo, N. Food Safety Concerns in “COVID-19 Era”. Microbiol. Res. 2021, 12, 53–68. [Google Scholar] [CrossRef]
  6. Li, Q.; Guan, X.; Wu, P.; Wang, X.; Zhou, L.; Tong, Y.; Ren, R.; Leung, K.S.; Lau, E.H.; Wong, J.Y. Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. N. Engl. J. Med. 2020, 382, 1199–1207. [Google Scholar] [CrossRef]
  7. Galanakis, C.M. The Food Systems in the Era of the Coronavirus (COVID-19) Pandemic Crisis. Foods 2020, 9, 523. [Google Scholar] [CrossRef]
  8. Duda-Chodak, A.; Lukasiewicz, M.; Ziec, G.; Florkiewicz, A.; Filipiak-Florkiewicz, A. COVID-19 pandemic and food: Present knowledge, risks, consumers fears and safety. Trends Food Sci. Technol. 2020, 105, 145–160. [Google Scholar] [CrossRef]
  9. Adelodun, B.; Ajibade, F.O.; Tiamiyu, A.O.; Nwogwu, N.A.; Ibrahim, R.G.; Kumar, P.; Kumar, V.; Odey, G.; Yadav, K.K.; Khan, A.H.; et al. Monitoring the presence and persistence of SARS-CoV-2 in water-food-environmental compartments: State of the knowledge and research needs. Environ. Res. 2021, 200, 111373. [Google Scholar] [CrossRef]
  10. Fisher, D.; Reilly, A.; Kang Eng Zheng, A.; Cook, A.; Anderson, D. Seeding of outbreaks of COVID-19 by contaminated fresh and frozen food. BioRxiv 2021. [Google Scholar] [CrossRef]
  11. Han, J.; Zhang, X.; He, S.; Jia, P. Can the coronavirus disease be transmitted from food? A review of evidence, risks, policies and knowledge gaps. Environ. Chem. Lett. 2020, 19, 5–16. [Google Scholar] [CrossRef] [PubMed]
  12. Kiss, I.Z.; Berthouze, L.; Taylor, T.J.; Simon, P.L. Modelling approaches for simple dynamic networks and applications to disease transmission models. Proc. R. Soc. A-Math. Phys. Eng. Sci. 2012, 468, 1332–1355. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, P.; Yang, M.; Zhao, X.; Guo, Y.; Wang, L.; Zhang, J.; Lei, W.; Han, W.; Jiang, F.; Liu, W.J. Cold-chain transportation in the frozen food industry may have caused a recurrence of COVID-19 cases in destination: Successful isolation of SARS-CoV-2 virus from the imported frozen cod package surface. Biosaf. Health 2020, 2, 199–201. [Google Scholar] [CrossRef] [PubMed]
  14. Han, S.; Liu, X. Can imported cold food cause COVID-19 recurrent outbreaks? A review. Environ. Chem. Lett. 2021, 20, 119–129. [Google Scholar] [CrossRef] [PubMed]
  15. Porter, K.A.; Ramaswamy, M.; Koloski, T.; Castrodale, L.; McLaughlin, J. COVID-19 Among Workers in the Seafood Processing Industry: Implications for Prevention Measures—Alaska, March-October 2020. MMWR Morb. Mortal. Weekly. Rep. 2021, 70, 622–626. [Google Scholar] [CrossRef] [PubMed]
  16. Pedreira, A.; Taskin, Y.; Garcia, M.R. A Critical Review of Disinfection Processes to Control SARS-CoV-2 Transmission in the Food Industry. Foods 2021, 10, 283. [Google Scholar] [CrossRef] [PubMed]
  17. Hu, L.; Gao, J.; Yao, L.; Zeng, L.; Liu, Q.; Zhou, Q.; Zhang, H.; Lu, D.; Fu, J.; Liu, Q.S.; et al. Evidence of Foodborne Transmission of the Coronavirus (COVID-19) through the Animal Products Food Supply Chain. Environ. Sci. Technol. 2021, 55, 2713–2716. [Google Scholar] [CrossRef]
  18. Dyal, J.W.; Grant, M.P.; Broadwater, K.; Bjork, A.; Waltenburg, M.A.; Gibbins, J.D.; Hale, C.; Silver, M.; Fischer, M.; Steinberg, J.; et al. COVID-19 Among Workers in Meat and Poultry Processing Facilities-19 States, April 2020. Mmwr-Morb. Mortal. Wkly. Rep. 2020, 69, 557–561. [Google Scholar] [CrossRef]
  19. Mardones, F.O.; Rich, K.M.; Boden, L.A.; Moreno-Switt, A.I.; Caipo, M.L.; Zimin-Veselkoff, N.; Alateeqi, A.M.; Baltenweck, I. The COVID-19 Pandemic and Global Food Security. Front Vet. Sci. 2020, 7, 578508. [Google Scholar] [CrossRef]
  20. Benkeblia, N. In the landscape of SARS-CoV-2 and fresh fruits and vegetables: The fake and hidden transmission risks. J. Food Saf. 2021, 41, e12898. [Google Scholar] [CrossRef]
  21. Dai, M.; Li, H.; Yan, N.; Huang, J.; Zhao, L.; Xu, S.; Wu, J.; Jiang, S.; Pan, C.; Liao, M. Long-term Survival of SARS-CoV-2 on Salmon as a Source for International Transmission. J. Infect. Dis. 2020, 223, 537–539. [Google Scholar] [CrossRef]
  22. Zhang, X.R. COVID-19 transmission in cold chain: A safe and green new-generation cold chain is demanded. Int. J. Energy Res. 2021, 45, 6483–6488. [Google Scholar] [CrossRef]
  23. Gonzalez, N.; Marques, M.; Domingo, J.L. Respiratory viruses in foods and their potential transmission through the diet: A review of the literature. Environ. Res. 2021, 195, 110826. [Google Scholar] [CrossRef] [PubMed]
  24. Lewis, D. Can COVID spread from frozen wildlife? Scientists probe pandemic origins. Nature 2021, 591, 18–19. [Google Scholar] [CrossRef] [PubMed]
  25. Chi, Y.; Zheng, S.; Liu, C.; Wang, Q. Transmission of SARS-CoV-2 on cold-chain food overpacks: A new challenge. J. Glob. Health 2021, 11, 03071. [Google Scholar] [CrossRef]
  26. Yekta, R.; Vahid-Dastjerdi, L.; Norouzbeigi, S.; Mortazavian, A.M. Food products as potential carriers of SARS-CoV-2. Food Control 2021, 123, 107754. [Google Scholar] [CrossRef]
  27. Ceylan, Z.; Meral, R.; Cetinkaya, T. Relevance of SARS-CoV-2 in food safety and food hygiene: Potential preventive measures, suggestions and nanotechnological approaches. Virusdisease 2020, 31, 154–160. [Google Scholar] [CrossRef]
  28. Velebit, B.; Djordjevic, V.; Milojevic, L.; Babic, M.; Grkovic, N.; Jankovic, V.; Yushina, Y. The common foodborne viruses: A review. In Proceedings of the 60th International Meat Industry Conference Meatcon2019; 2019; Volume 333, p. 012110. [Google Scholar] [CrossRef]
  29. Rodrigue, J. The Geography of Transport Systems, 5th ed.; Taylor & Francis: London, UK, 2020. [Google Scholar]
  30. Galimberti, A.; Cena, H.; Campone, L.; Ferri, E.; Dell’Agli, M.; Sangiovanni, E.; Belingheri, M.; Riva, M.A.; Casiraghi, M.; Labra, M. Rethinking Urban and Food Policies to Improve Citizens Safety After COVID-19 Pandemic. Front Nutr. 2020, 7, 569542. [Google Scholar] [CrossRef]
  31. Ji, W.; Li, X.; Chen, S.; Ren, L. Transmission of SARS-CoV-2 via fomite, especially cold chain, should not be ignored. Proc. Natl. Acad. Sci. USA 2021, 118, e2026093118. [Google Scholar] [CrossRef]
  32. Zhou, P.; Shi, Z.L. SARS-CoV-2 spillover events. Science 2021, 371, 120–122. [Google Scholar] [CrossRef]
  33. Han, X.; Xu, Y.; Fan, L.; Huang, Y.; Xu, M.; Gao, S. Quantifying COVID-19 importation risk in a dynamic network of domestic cities and international countries. Proc. Natl. Acad. Sci. USA 2021, 118, e2100201118. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, F.; Liu, Y.; Ya, B.; He, J.; Leng, T.; Sun, Y.; Xin, Q. Latest Effective Measures to Combat COVID-19: A Review. Iran. J. Public Health 2021, 50, 640–648. [Google Scholar] [CrossRef] [PubMed]
  35. Reimagining the Cold Chain in a Post-Pandemic World; June 2020. Available online: https://www.corporate.carrier.com/Images/Reimagining-the-Cold-Chain-White-Paper_tcm558-82726.pdf (accessed on 10 September 2021).
  36. Torreggiani, S.; Mangioni, G.; Puma, M.J.; Fagiolo, G. Identifying the community structure of the food-trade international multi-network. Environ. Res. Lett. 2018, 13, 054026. [Google Scholar] [CrossRef]
  37. Ercsey-Ravasz, M.; Toroczkai, Z.; Lakner, Z.; Baranyi, J. Complexity of the international agro-food trade network and its impact on food safety. PLoS ONE 2012, 7, e37810. [Google Scholar] [CrossRef]
  38. Newman, M. Network; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  39. Newman, M.E.; Girvan, M. Finding and evaluating community structure in networks. Phys. Rev. E Stat. Nonlin Soft Matter Phys. 2004, 69, 026113. [Google Scholar] [CrossRef] [Green Version]
  40. D’Odorico, P.; Carr, J.A.; Laio, F.; Ridolfi, L.; Vandoni, S. Feeding humanity through global food trade. Earth’s Future 2014, 2, 458–469. [Google Scholar] [CrossRef]
  41. Gustavsson, J.; Cederberg, C.; Sonesson, U. Global Food Losses and Food Waste; FAO: Rome, Italy, 2011. [Google Scholar]
  42. Yahia, E.M.; Smolak, J. Developing the Cold Chain for Agriculture in the Near East and North Africa (NENA). Policy Brief; FAO: Cairo, Egypt, 2014. [Google Scholar]
  43. Lange, B.; Priesemann, C.; Geiss, M.; Lambrecht, A. Promoting Food Security and Safety via Cold Chains; Internationale Zusammenarbeit (GIZ) GmbH: Bonn, Germany, 2016. [Google Scholar]
  44. COVID-19 Spurs Demand in the Frozen Food Industry. Available online: https://foodtank.com/news/2021/02/covid-spurs-demand-in-the-frozen-food-industry/ (accessed on 1 September 2021).
  45. de Paulo Farias, D.; Dos Santos Gomes, M.G. COVID-19 outbreak: What should be done to avoid food shortages? Trends Food Sci. Technol. 2020, 102, 291–292. [Google Scholar] [CrossRef]
  46. Aday, S.; Aday, M.S. Impact of COVID-19 on the food supply chain. Food Qual. Saf. 2020, 4, 167–180. [Google Scholar] [CrossRef]
  47. Hobbs, J.E. The COVID-19 pandemic and meat supply chains. Meat Sci. 2021, 181, 108459. [Google Scholar] [CrossRef]
  48. Dang, H.H.; Viet Nguyen, C. Gender inequality during the COVID-19 pandemic: Income, expenditure, savings, and job loss. World Dev. 2021, 140, 105296. [Google Scholar] [CrossRef]
  49. Ursachi, C.S.; Munteanu, F.D.; Cioca, G. The Safety of Slaughterhouse Workers during the Pandemic Crisis. Int J. Environ. Res. Public Health 2021, 18, 2633. [Google Scholar] [CrossRef] [PubMed]
  50. Bosworth, A.; Whalley, C.; Poxon, C.; Wanigasooriya, K.; Pickles, O.; Aldera, E.L.; Papakonstantinou, D.; Morley, G.L.; Walker, E.M.; Zielinska, A.E.; et al. Rapid implementation and validation of a cold-chain free SARS-CoV-2 diagnostic testing workflow to support surge capacity. J. Clin. Virol. 2020, 128, 104469. [Google Scholar] [CrossRef] [PubMed]
  51. Zhang, J.; Fang, X.; Mao, Y.; Qi, H.; Wu, J.; Liu, X.; You, F.; Zhao, W.; Chen, Y.; Zheng, L. Real-time, selective, and low-cost detection of trace level SARS-CoV-2 spike-protein for cold-chain food quarantine. npj Sci. Food 2021, 5, 12. [Google Scholar] [CrossRef] [PubMed]
  52. Chang, H.-H.; Meyerhoefer, C.D. COVID-19 and the Demand for Online Food Shopping Services: Empirical Evidence from Taiwan. Am. J. Agric. Econ. 2021, 103, 448–465. [Google Scholar] [CrossRef]
  53. Bracale, R.; Vaccaro, C.M. Changes in food choice following restrictive measures due to COVID-19. Nutr. Metab. Cardiovasc. Dis. 2020, 30, 1423–1426. [Google Scholar] [CrossRef]
Figure 1. Study framework of COVID-19 transmission risk by CCF.
Figure 1. Study framework of COVID-19 transmission risk by CCF.
Foods 11 00908 g001
Figure 2. Spatial patterns of CCF consumption level (a) and urbanization level (b).
Figure 2. Spatial patterns of CCF consumption level (a) and urbanization level (b).
Foods 11 00908 g002
Figure 3. Probability distribution of strength-out degree (a) and BC (b).
Figure 3. Probability distribution of strength-out degree (a) and BC (b).
Foods 11 00908 g003
Figure 4. Community structure and flow patterns in the CCFTN.
Figure 4. Community structure and flow patterns in the CCFTN.
Foods 11 00908 g004
Figure 5. National of COVID-19 intensity transmission by CCF.
Figure 5. National of COVID-19 intensity transmission by CCF.
Foods 11 00908 g005
Figure 6. Spatial pattern of internal diffusion risk (a) and external input risk (b).
Figure 6. Spatial pattern of internal diffusion risk (a) and external input risk (b).
Foods 11 00908 g006
Figure 7. Spatial pattern of comprehensive risk (a) and dominant risk type (b).
Figure 7. Spatial pattern of comprehensive risk (a) and dominant risk type (b).
Foods 11 00908 g007aFoods 11 00908 g007b
Table 1. Some case of imported cold food supply chain-related coronavirus disease cold chain practitioners.
Table 1. Some case of imported cold food supply chain-related coronavirus disease cold chain practitioners.
Outbreak TimeOutbreak CityPlace Cases (Tested)Zero PatientSARS-CoV-2 Source
11 June 2020BeijingAgricultural products in 335 (>11,000,000) Xinfadi Wholesale marketEmployeeEnvironmental swab samples related to imported salmon
22 July 2020Dalian, LiaoningKaiyang company (Seafood 92 (>6,000,000) processing enterprises)
Dock
DockworkerOuter packaging of imported fish
24 September 2020Qingdao, ShandongQingdao port, Dock 2 (5781)Stevedores, DockworkerOuter packaging surface of imported frozen cod
11 October 2020Qingdao, ShandongDarang company of Qing-12 (10,920,000)
dao port
Stevedore DockworkerOuter packaging of imported frozen cod
25 October 2020Kashgar, XinjiangKashgar Airport 138 (4,746,500)StevedoreContainer from abroad
7 November 2020TianjinHailian Frozen Food Co., 12 (1,030,000)
Ltd Tianjin Hailian cold storage
Stevedore DockworkerOuter packaging of imported pig food packaging of frozen pork
9 November 2020ShanghaiPudong Airport 4 (>14,000)StevedoreAirborne container from North America
15 December 2020Dalian, LiaoningDalian Port Yidu cold chain 83 (6,379,000) Co., Ltd.StevedoreEnvironmental swab samples related to imported cold chain food
Sources: Health Times (www.jksb.com.cn, accessed on 10 September 2021); Chinese Center for Disease Control and Prevention (www.chinacdc.cn, accessed on 10 September 2021).
Table 2. Criteria for classification of different elements.
Table 2. Criteria for classification of different elements.
LevelJenks BreaksEqual Interval
IDRIEIRICTRICCF ConsumptionUR
Very Low0–8.150–1.051–40–3860–20
Low8.15–19.671.05–3.195–8386–60520–40
Medium19.67–35.783.19–8.149–12605–83440–60
High35.78–56.188.14–15.9913–16834–115460–80
Very high56.18–91.3415.99–42.8117–201154–187080–100
Note: The original values of IDRI and EIRI are very large, so the IDRI and EIRI are divided by 105 and 1410.
Table 3. Characteristics of the top 20 countries in the CCFTN.
Table 3. Characteristics of the top 20 countries in the CCFTN.
Export QuantityOut-StrengthBetweenness Centrality
ISO3Value
(1010 kcal)
Proportion (%)OutdegreeISO3Value
(1010 kcal)
Proportion (%)ISO3ValueOutdegree
IDN43.577.34173IDN7537.799.39USA2037.71187
ARG36.386.13134NLD5316.556.63NLD1368.40186
UKR34.285.78150UKR5142.206.41GBR1327.77161
NLD28.584.82186ARG4875.426.08MYS1217.67184
RUS26.924.54104USA4476.275.58CAN1082.32152
CAN25.624.32152MYS4367.645.44FRA1072.42179
USA23.944.03187DEU4184.125.21BEL1059.53183
MYS23.744.00184CAN3894.324.85CHN1005.16187
DEU23.643.98177ESP3034.503.78DEU863.48177
ESP17.242.91176RUS2799.303.49ESP815.76176
BRA14.662.47181BRA2653.963.31SAU813.87123
BEL12.762.15183BEL2334.452.91NZL752.39161
PNG12.232.0615CHN2127.722.65AUS719.76133
FRA11.751.98179FRA2104.122.62ITA697.17177
CHN11.381.92187NZL1384.961.73ZAF626.26156
BOL10.241.7335ITA1315.441.64ARE615.10141
BLR9.131.5465POL1304.791.63IND534.72171
PHL8.901.50120TUR1100.141.37UKR488.39150
NZL8.601.45161DNK1087.441.36CHE478.08119
POL8.361.41156PHL1068.041.33TUR477.70170
Top20391.9466.06/Top2062,109.1877.40///
Global593.33100.00/Global80,249.05100.00///
Note: ISO3 is the nation code, see the Table S4 for the full country name. “/” indicates no corresponding data.
Table 4. Characteristics of the CCFTN community structure.
Table 4. Characteristics of the CCFTN community structure.
NameIntra-Community TradeNodeEdgeCRIIDRIEIRICTRI
Quantity
(1010 kcal)
Proportion
(%)
NumberProportion
(%)
Number
1-NLD-DEU22,817.7328.435125.7613957.2938.1636.5711.88
2-IDN-UKR12,089.9915.076030.309723.2311.1219.046.23
3-USA-CAN6854.808.543015.153716.1325.9323.819.67
4-ARG-RUS5080.576.333517.683945.3123.3215.659.57
5-TUR-ARE1485.121.85178.591375.5325.1615.439.88
6-SRB- MKD119.150.1552.53208.6048.204.6311.40
Global80,249.05/198100.00119245.4224.6323.019.24
Note: The original values of IDRI and EIRI are very large; as such, the IDRI and EIRI are divided by 105 and 1410.
Table 5. Classification statistics of different COVID-19 transmission risks.
Table 5. Classification statistics of different COVID-19 transmission risks.
LevelInternal Diffusion RiskExternal Input RiskComprehensive Risk
NumberProportion (%)NumberProportion (%)NumberProportion (%)
Very low5326.7712965.15405.56
Low4924.753316.675122.73
Medium4020.202412.125125.76
High3718.6973.544525.76
Very high199.6052.531120.20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, C.; Yang, Y.; Feng, Z.; Xiao, C.; Liu, Y.; Song, X.; Lang, T. Cold Chain Food and COVID-19 Transmission Risk: From the Perspective of Consumption and Trade. Foods 2022, 11, 908. https://doi.org/10.3390/foods11070908

AMA Style

Zhang C, Yang Y, Feng Z, Xiao C, Liu Y, Song X, Lang T. Cold Chain Food and COVID-19 Transmission Risk: From the Perspective of Consumption and Trade. Foods. 2022; 11(7):908. https://doi.org/10.3390/foods11070908

Chicago/Turabian Style

Zhang, Chao, Yanzhao Yang, Zhiming Feng, Chiwei Xiao, Ying Liu, Xinzhe Song, and Tingting Lang. 2022. "Cold Chain Food and COVID-19 Transmission Risk: From the Perspective of Consumption and Trade" Foods 11, no. 7: 908. https://doi.org/10.3390/foods11070908

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