How do temperature, humidity, and air saturation state affect the COVID-19 transmission risk?

Environmental parameters have a significant impact on the spread of respiratory viral diseases (temperature (T), relative humidity (RH), and air saturation state). T and RH are strongly correlated with viral inactivation in the air, whereas supersaturated air can promote droplet deposition in the respiratory tract. This study introduces a new concept, the dynamic virus deposition ratio (α), that reflects the dynamic changes in viral inactivation and droplet deposition under varying ambient environments. A non-steady-state-modified Wells-Riley model is established to predict the infection risk of shared air space and highlight the high-risk environmental conditions. Findings reveal that a rise in T would significantly reduce the transmission of COVID-19 in the cold season, while the effect is not significant in the hot season. The infection risk under low-T and high-RH conditions, such as the frozen seafood market, is substantially underestimated, which should be taken seriously. The study encourages selected containment measures against high-risk environmental conditions and cross-discipline management in the public health crisis based on meteorology, government, and medical research.

Introduction

The novel coronavirus disease 2019 (COVID-19) has stimulated unprecedented public health concerns (Azman and Luquero 2020; Liu et al. 2020b; Morawska et al. 2020). Two years have passed since its early outbreak, COVID-19 prevention has destined to be a prolonged battle. It is impracticable for all human active areas to take long-term disinfection and isolation of activities. Therefore, our attention needs to be directed to high-risk areas (Coccia 2021b). Without considering other artificial factors, an urgent need is to identify the exposure scenarios that may promote COVID-19 transmission.

It has been reported that seasonal cyclicity is a ubiquitous feature of acute infectious diseases, which is also commonly observed in respiratory viral diseases, such as SARS and MERS (Killerby et al. 2018; Wang et al. 2020). Therefore, it is reasonable to speculate that environmental and meteorological factors affect the COVID-19 (Xie and Zhu 2020). The analysis of how certain meteorological conditions may have affected the initial spread of COVID-19 at the country, city, or regional level has become an important line of research recently (Li et al. 2022; Rahimi et al. 2021). To control the pandemic, many studies have examined meteorological/weather conditions that might influence the spread of coronaviruses by exploring the association between these factors and the COVID-19 cases over different worldwide locations (see Appendix A: Table A.1).

Some scholars only selected T and RH to analyze the correlation between meteorological parameters and the spread of COVID-19 (absolute humidity and rainfall are also considered as a form of humidity). The preliminary results showed T and RH significantly related to the transmission risk (Haque and Rahman 2020; Mendez-Arriaga 2020; Pani et al. 2020; Prata et al. 2020; Sobral et al. 2020; Tosepu et al. 2020; Wu et al. 2020).

Some scholars not only focused on T and RH but also comprehensively analyzed the other meteorological parameters, such as solar radiation, wind speed, and air pollutant (Diao et al. 2021; Ma et al. 2020; Menebo 2020; Rosario et al. 2020; Runkle et al. 2020; Şahin 2020; Xu et al. 2020a). The solar radiation showed a negligible on the COVID-19 transmission, even in the comparison of high- and low-altitude regions (Song et al. 2022). The cities with low wind speed and high levels of air pollution exceeding safe levels of particulate matter had higher numbers of COVID-19-related cases and deaths (Coccia 2020c, a, 2021a, c, d, b, 2022b).

In addition, many scholars did their epidemic research from many different angles and combined meteorological parameters (T and RH) with population mobility. The conclusion is obvious that the higher population density and migration scale tends to lead to prolong the spread (Ahmadi et al. 2020; Diao et al. 2021; Lin et al. 2020; Liu et al. 2020a).

In brief, T and RH are not only the main meteorological parameters widely analyzed by scholars but also significantly correlated with the COVID-19 pandemic, even with the incorporation of additional variables and controls (solar radiation, wind speed, and air pollutant, and population migrated).

However, the association between meteorological parameters and COVID-19 cases also would be affected by local artificial factors, including social, political, and economic (Lai et al. 2020). The results in quantifying environmental effects on transmission risk are limited and equivocal (Briz-Redón and Serrano-Aroca 2020).

Therefore, to reject the effect of other social factors, quantifiable mathematical models of transmission risk are needed (Bin et al. 2018; Borro et al. 2020; Qian et al. 2012; Sharma and Balasubramanian 2020). Wells-Riley model (W-R model) is one of the most popular models for quantitatively assessing the infection risk of airborne diseases (Keene 1955; Riley et al. 1978; Vuorinen et al. 2020; Wagner et al. 2009; Zhang et al. 2020). Understandably, the W-R model has been further improved and extended by many latest studies to include more realistic factors and considerations since the COVID-19 outbreak. For instance, Nordsiek et al. developed an optimization W-R model for computational risk assessment suitable for mono/poly-pathogen aerosols (Nordsiek et al. 2021). Andrews et al. used a modified W-R model to estimate the impact of ventilation on infection probability for inmates sharing a cell with an infector (Urrego et al. 2015). Zhou and Koutsopoulos proposed a modified W-R model for risk analysis in public transportation systems, which integrated with a simulation model of subway operations (Zhou and Koutsopoulos 2021). Zhai et al. and Zhang et al. developed a modified W-R model for integrating the spatial distribution of pathogen concentrations (Guo et al. 2021; Zhai and Li 2021). Fierce et al. tried to develop the quadrature-based model of droplet risk by introducing a dose-response framework, which could reflect the evolution of particles after they are expelled, their deposition in the airways, and the subsequent risk of initial infection (Fierce et al. 2021).

Although the various studies above used different methods for risk analysis of airborne diseases, they collectively treated the shared air space as a close-to-greenhouse environment. The T, RH, and P are assumed to be constant. However, the epidemic outbreak in different regions keeps transmitting under varying environmental parameters. It can be seen from Table 1 that those environmental parameters were believed to significantly promote or inhibit the spread of COVID-19. Therefore, a non-steady-state-modified Wells-Riley model of the dynamic change of environmental parameters is required.

Table 1

Control parameters for different exposure scenarios

Environment parameters	Case 1	Case 2	Case 3	Case 4	Case 5	Case 6
Temperaturea (T, °C)	40	60	20	24	20	5
Relative humidityb (RH, %)	90	50	30	60	80	95
Viral inactivation (λ, h−1)	0.8	0.99	0.4	0.9	0.2	0.15
Saturated vapor pressure (Pv, hPa)	66.37	99.64	7.01	17.9	18.7	8.3
Vapor pressure differencec (|Pv−Pvo|, hPa)	4.27	37.5	55.1	44.2	43.4	53.8

aHigh temperature: T > 30 °C; middle temperature: 10 °C< T < 30 °C; low temperature: T < 30 °C

bHigher humidity: RH > 80%; middle humidity: 40% < RH < 80%; lower humidity: RH < 30%

cAbsolute value of the saturated vapor pressure difference between ambient air and respiratory tract air

Respiratory boundary conditions: T = 37 °C, RH = 99%, and Pvo = 62.1 hPa

This study aims to answer the following questions: How are T, RH, and P linked to infection risk? Which exposure scenarios promote or inhibit COVID-19 transmission? The answers to these questions will help local healthcare policymakers grade human activity regions according to the COVID-19 transmission risk.

Materials and methods

Sample and data

A thorough understanding of infectious droplets is a primary factor in the study of COVID-19 airborne transmission risks, such as the diameter distribution, SARS-CoV-2 viral load, viral inactivation in the air, and deposition fraction in the respiratory tract. This study summarizes the number distribution of exhaled droplets from different respiratory activities and the viral load in the body fluids of infected persons (see Appendix A: Table A.2 and Table A.3, respectively).

Notably, the viral inactivation in the air and the droplet deposition in the respiratory tract are significantly influenced by ambient environmental factors (T, RH, and P) and should be discussed with emphasis (Ishmatov 2020; Longest et al. 2011; Xi et al. 2015; Xie and Zhu 2020).

Viral inactivation characteristics in the air

The significant changes in virus survival balance after being exhaled from the infector lead to differences in the disease infection risk in different environments (Pica and Bouvier 2012; van Doremalen et al. 2020; Wei and Li 2016). Therefore, several scholars have conducted studies on aerosolized viral inactivation under varying T and RH and aerosolized the viruses into a rotating drum, where the aerosols were held at the desired T and RH. Following this, the samples of air from the drum were collected using an impinger (Doremalen et al. 2013; Pyankov et al. 2018; Sattar et al. 1984). This study summarizes viral inactivation characteristic of respiratory disease under varying T and RH, as shown in Appendix A: Figure A.1 and Table A.4

When T was at a specific value (22–25 °C), the viruses were found to survive better at low-RH levels (<33%) and high-RH levels (>85%). The middle-RH level (50–75%) was found to be the least favorable for the survival of the viruses. And, within the range of certain RH, there was a significant negative correlation between the viral inactivation rate and T (Colas et al. 2014; Harper 1961; Lin and Marr 2020b; Noti et al. 2013). Therefore, the viral inactivation rate in the air (λ) is significantly affected by the ambient T and RH, which further inhibits or promotes the spread of an epidemic in local areas (Mao et al. 2020; Walsh et al. 2020).

Droplet deposition in the respiratory tract

Diameter growth of inhaled droplets in the airways

To analyze pulmonary drug delivery, plenty of studies on particle deposition in the airways have been conducted by medical scholars (Austin et al. 2010; Gralton et al. 2011; Knight 2010; Morrow 2010; Nicas et al. 2005; Yeh and Raabe 1976). During inhalation, when the respiratory tract air (T = 37 °C, RH = 99% (Winkler-Heil et al. 2017)) is mixed with the inhaled ambient air (relatively low RH and T), the small particle (<10 μm) can grow owing to condensation. Therefore, the concept of enhanced condensational growth (ECG effect) was proposed by medical scholars in which a drug aerosol is inhaled in combination with cold-saturated air; this would cause the small particle diameter to increase, thereby promoting lung drug deposition (Deng et al. 2020; Longest and Hindle 2010, 2011; Longest et al. 2011; Tian et al. 2011). The diameter growth factors of inhaled supersaturation droplets in the airways are summarized (Kim et al. 2013; Kreyling 1984; Longest et al. 2011; Sarangapani and Wexler 1996; Xi et al. 2013; Xi et al. 2015) (see Appendix A: Table A.5, and Figure A.2). Understandably, the ECG effect would also increase the deposition risk of infectious droplets.

Droplet deposition fraction in the respiratory tract

It is known that viruses can only cause infection when inhaled by susceptible populations and when successfully deposited in the respiratory tract (Akhbarizadeh et al. 2021; Pu et al. 2020). Furthermore, the deposition fraction of infectious droplets is dependent on the aerosol droplet diameter distribution (Asgari et al. 2019; Fröhlich-Nowoisky et al. 2016). Figure 1 summarizes the distribution of exhaled droplets from different respiratory activities and the corresponding deposition fractions in the respiratory tract (ICRP 1994).

Fig. 1

Number distribution of exhaled droplets from different respiratory activities and the corresponding deposition fraction in the respiratory tract

Large particles are mainly deposited through inertial impaction, whereas small particles are mainly deposited through molecular diffusion (Fröhlich-Nowoisky et al. 2016). It is noteworthy that supersaturated air can lead to significant growth in inhaled small particle diameter in the airways due to the ECG effect. For example, the total deposition of 0.3-μm particles in the airways may rise from 13% (when supersaturation is not considered) to 90% (under supersaturated conditions). It is evident that the infection risk is set to surge.

Measures of variables

Viral inactivation rate in the air (λ)

The viral inactivation rate in the air (λ) is significantly affected by the ambient T and RH, which further inhibits or promotes the spread of an epidemic in local areas. The value of λ was selected from existing studies on aerosolized viral inactivation under varying T and RH (Appendix A: Table A.4)

Dynamic virus deposition ratio (α)

Considering the dynamic changes in particle size and viral deposition in different ambient environments, our study proposes a new term, the dynamic virus deposition ratio (α) (Eq. (1)), which is defined as the ratio between the viral load deposited in the respiratory tract under varying environments (Qdeposition) and the total viral load inhaled by a susceptible population (Qtotal).

where Qdeposition is the viral load deposited in the respiratory tract under varying ambient environments (RNA copies), and Qtotal is the total viral load inhaled by a susceptible population (RNA copies).

where is the viral concentration of droplets with I diameter after considering the ECG effect on particle size (RNA copies/mL), is the number of droplets with i diameter after considering the ECG effect on particle size (part. m−3), is the volume of a single droplet (mL) as a function of the droplet diameter after considering the ECG effect on particle size (d), and is the deposition fraction of droplets with i diameter after considering the ECG effect on particle size (the values were selected according to the ICRP data (ICRP 1994) ).

The supersaturation diameter growth factors for inhaled air with different saturated vapor pressures (d/d0) can be calculated from the following equation (Xi et al. 2015):

where P is the saturated vapor pressure of the ambient air (hPa), and p is the saturated vapor pressure of air in the respiratory tract (T = 37 °C, RH = 99%, P = 62.1 hPa).

However, droplet condensation in the airways (ECG effect) changes their particle size (), which modifies the concentration of viral charges in the droplets (, copies/mL) due to a higher water fraction. Meanwhile, it should not change the absolute number of viral copies in the droplet. (For example, a 1-μm droplet has the same number of viral charges after condensation, whereas their concentration is reduced owing to a higher liquid volume.) Therefore, the absolute number of viral copies in a droplet with i diameter after considering the ECG effect can only be a dynamic cumulative value (Eq. (5)), not directly calculated by the original volume and concentration (Eq. (2)).

where is the total viral load of droplets with diameter i after considering the ECG effect on the particle size (RNA copies).

Model and data analysis procedure

Modified Wells-Riley model

Quanta emission

A quantum is defined as the dose of airborne droplet nuclei required to cause infection in 63% of susceptible persons. In the original W-R model, the quanta emission rate remains uncertain and highly variable (Gba et al. 2020; Harrichandra et al. 2020). To estimate the spatiotemporal variation in virion concentrations, an Italian scholar, Buonannoa, developed a modified W-R model based on a dose-response method to calculate the virion emission rate through a mass balance, which allows consideration of per-particle variation in viral load and the airborne transmission risk between scenarios (Buonanno et al. 2020). In particular, the viral load emitted was expressed in terms of the quanta emission rate (ER).

where ER is the quanta emission rate from different respiratory activities (quanta h−1), j indicates the different respiratory activities considered (breath = 1, voice = 2, cough = 3), C is the viral concentration in the throat of the asymptomatic infector (RNA copies/mL), C is a conversion factor defined as the ratio between one infectious quantum and the infectious dose expressed in viral RNA copies (Li et al. 2010; Watanabe et al. 2010; Yu et al. 2004) (0.01–0.1), p is the breathing rate per person (Marmett et al. 2020) (m3/s), N is the number of droplets with i diameter from j respiratory activities, and V(d0) is the volume of a single droplet (mL) as a function of the droplet diameter (d0); the droplets are assumed to be standard spheres.

Quanta concentration of the shared air space

After being expelled from the infector, the quanta concentration at time t (Gammaitoni and Nucci 1997), n(t), as follows:

where n(t) is the quanta concentration from different respiratory activities at time t, I is the number of asymptomatic infectors, AER is the air exchange rate via ventilation (natural ventilation (Alfano et al. 2012, Stabile et al. 2017), 0.2 h−1), K is the particle deposition on surfaces (released from 1.5 m at a speed of 1 × 10−4 m/s (Chatoutsidou and Lazaridis 2019), 0.24 h−1), λ is the viral inactivation rate of respiratory disease (the values were selected according to Appendix A: Table A.4), V is the volume of the shared air space(m−3), and t is the exposure time (s).

Infection risk under varying environments

Thereafter, the dynamic virus deposition ratio, α, is used to correct the breathing rate of the susceptible population, p, and the modified dimensionless inhalation rate, p′, is expected to reflect the dynamic risk due to the ECG effect more accurately.

Infection risk (R, %) as a function of exposure time (t) of susceptible people.

In the end, the basic reproductive number (R0) could be calculated if we multiply the peak of infection risk by the number of all exposed persons during the exposure time (Rothman et al. 2011).

Data analysis procedure

Example case

ECG effect is a concept of pulmonary drug delivery in which the drug aerosol is inhaled combination with cold-supersaturated air. The subsequent small particle condensation growth will promote lung drug deposition (Xi et al. 2015). Similarly, the ECG effect would increase the virus deposition risk of airborne diseases.

To illustrate the impact of ECG effect on respiratory particle deposition and Virion exposure, this paper presents the calculation process of dynamic virus deposition ratio, α, in typical cold-saturated conditions (example case: T = 5 °C, RH = 95%), as shown in Fig. 4.

Fig. 4

Infection risk trends as a function of time (an asymptomatic infector remained inside the space for 10 min). The difference (△R) between RLTS (T = 10 °C, RH = 100%) and RMTUS (T = 25 °C, RH = 60%) for three respiratory activities (cough, voice, breath) are compared, respectively

According to Eq. (4), the droplet diameter would change during inhalation. The smaller the particle size, the bigger the growth factor, and only a small particle (<10 μm) shows significantly grow (see Fig. 2(a)). It is assumed that these droplets are deposited after condensation, and the final droplet distribution ranged from 3 to 100 μm. This change is unfortunate to reach the high-risk deposition peak in the airways (see Fig. 2(b)). As a result, the number of deposited particles increased by 3 times (see Fig. 2(c)), and the dynamic virus deposition ratio, α, of this condition could go as high as 51.94% (see Fig. 2(d)). Thus, it is evident that the virion exposure and infection risk increases dramatically (example case: T = 5 °C, RH = 95%).

Fig. 2

Calculation process of dynamic virus deposition ratio, α, in typical cold-saturated condition (example case: T = 5 °C, RH = 95%)

Results

Quanta emission rate from different respiratory activities

Figure 3 shows the ER (quanta.h−1) trends as a function of time and viral load in the throat (C, Log10 RNA copies. mL−1) and the quanta-RNA copies correction factor (C, 0.01-0.1) from three respiratory activities (breath, voice, cough) and light exercise activity level (p, 1.38 m3/h). For the sole purpose of simplifying the discussion, zones representative of low (<1 quantum h−1) and high (>1 quantum h−1) quanta emission are separated by the red dotted line indicated (ER, 1 quantum h−1).

Fig. 3

ERj (quantum. h−1) trends as a function of the time and the viral concentration in the throat (C, Log10 RNA copies. mL−1) and quanta-RNA copies correction factor (C, 0.01-0.1) for the three respiratory activities (breath, voice, cough). Only the first 15 days of data are calculated after onset. Zones representative of low (<1 quantum h−1) and high (>1 quantum h−1) quanta emission are separated by the red dotted line indicated

The trends of the quanta emission rate are similar to that of the viral load of SARS-CoV-2 in a patient’s body fluids, peaking approximately 4 to 6 days after onset and then decreasing over time. In the case involving typical clinical symptoms of the COVID-19 patient (cough condition), the quanta emission rate covers the range of 10−4 to 105 quanta.h−1, up to 5.55 × 104 quanta.h−1. A high emission in the case of coughing was achieved at 11 days out of 15. Speech is one of the most important communication methods between people (voice condition). The quanta emission rate covering the range of 10−5 to 104 quanta.h−1 is approximately 1/10 of the coughing condition. A high emission in the case of speaking can be achieved only at five days out of 15. It is noteworthy that breathing is a continuous behavior of an asymptomatic person during his/her daily activities (breath condition). Although the quanta emission rate fails to achieve a high emission zone in the breath condition, more viruses might accumulate over time due to having a higher frequency than cough and voice conditions.

Infection risk for different environments

This section aims to show how the environmental parameters (T, RH, and P) are linked to the disease transmission risk in two different microenvironments, low-temperature-saturated condition (low-T: T = 10 °C, saturated: RH = 100%, R) and middle-temperature-unsaturated condition (middle-T: T = 25 °C, unsaturated: RH = 60%, R).The following input data are defined: C is 108 RNA copies/mL, AER is 0.2 h−1 (natural ventilation (Buonanno et al. 2020)), K is 0.24 h−1 (Chatoutsidou and Lazaridis 2019), λ is 0.4 h−1 and 0.8 h−1 respectively for R and R (according to Appendix A: Table A.4), Ci is 0.02 (Watanabe et al. 2010), V is 8 m*6 m*3.5 m = 168 m3, and p is 1.38 m3.h−1 (light exercise activity level (Adams 1993)). It is assumed that only an asymptomatic infector remained inside the environment for 10 min; the infection risk is simulated for up to 6 h, as shown in Fig. 4.

The trends highlight that the brief stopover of the asymptomatic infector in the microenvironments leads to a peak risk of infection at 10 min. The more intense the respiratory activity, the higher the quanta emission rate. Therefore, the highest disease infection risk was 50% for cough, 6.2% for voice, and 1.2% for breath. In particular, it is known from the literature mentioned above that the viruses are found to survive better in low-temperature air, and the subsequent particle growth of inhaled supersaturated air promotes lung deposition. Therefore, it is understandable that a higher peak infection risk is reached at R (T = 10 °C, RH = 100%) rather than R (T = 25 °C, RH = 60%). Cough, ∆R = 15±3.3%, p < 0.01; Voice, ∆R = 2±0.5%, p < 0.01; Breath, ∆R = 0.0022±0.0006%, p < 0.01.

Typical scenarios studies

The research results in the “Infection risk for different environments” section support the argument that infection risk varies with the environment. Therefore, six typical exposure scenarios from daily life were chosen to analyze the effects of environmental parameters on disease transmission risk. The control parameters for the different exposure scenarios are summarized in Table 1.

First, the exposure scenarios of the shower room (case 1, high-T and high-RH conditions (China 1998)), and the sauna (case 2, high-T and middle-RH conditions (China 1998)) were selected to analyze the infection risk under high temperature. Second, the most common exposure scenario during daily life is selected as a comparison condition (case 4, middle-T and middle-RH conditions (China 2017)). As the flu season approaches, the indoor exposure scenarios during the rainy condition in autumn (case 5, middle-T and high-RH conditions (Nan et al. 2009)) and during heating conditions in winter (case 3, middle-T and low-RH conditions (Yan 2015)) were selected to analyze the infection risk under middle temperature. Finally, based on the conclusions in the “Model and data analysis procedure” section, the infection risk in the scenario (cold and damp) must be taken seriously. Therefore, the extreme typical exposure scenario of the seafood cold store (case 6, low-T and high-RH conditions (China 2010)) is analyzed.

For all the typical exposure scenarios considered in the simulations, the asymptomatic infector remains inside the environment for 5 min, as shown in Fig. 5. The aim is to compare the infection risk during speaking (voice, j = 2) in the shared air space. For this reason, the following assumptions for all scenarios are defined: 4 persons are always present; 1 new visitor every 5 min enters; every visitor remains inside for 5 min; thus, 5 persons are simultaneously present. The other input data is consistent with the “Infection risk for different environments” section.

Fig. 5

Infection risk (Rmax) and number of the exposed persons (Nmax) for all scenarios

Although the quanta emission rate is consistent in each scenario, viral inactivation in the air (λ) and the deposition fraction in the respiratory tract (ω) are significantly influenced by the environment. Thus, remarkable differences in disease transmission between the scenarios are noted. The brief stopover of the asymptomatic infector in the microenvironments leads to a peak risk of infection at 5 min (Rmax, the peak of infection risk). After the asymptomatic infector leaves, the quanta concentration slowly drops to a safe level (Tterminal, the high infection risk duration of the microenvironment). During this period, the number of all the persons who visited the microenvironment is counted (Nmax, the number of all exposed persons). In the case of middle-T conditions, such as case 3, case 4, and case 5 (20 °C, 24 °C, 20 °C), as the RH increases gradually (30%, 60%, 80%), the Rmax decreases first and then increases later (3.6%, 2.1%, 4.8%); this may be because the middle RH (45–60%) is found to be the least favorable for viral survival. Furthermore, in the case of high-RH conditions such as case 1, case 5, and case 6 (90%, 80%, 95%), as the ambient T reduces gradually (40 °C, 20 °C, and 5 °C), the P difference in vivo and in vitro gradually increases (4.27 hPa, 43.4 hPa, 53.8 hPa) and Rmax also shows a rising trend (2.3%, 4.8%, 5.2%). The large value for case 6 is evidently due to the low viral inactivation in the air (SARS-CoV-2 was found to survive better in low-temperature conditions) and the high deposition fraction in the respiratory tract (the supersaturated air could promote lung deposition due to the ECG effect). This could be another explanation as to why COVID-19 transmission was enhanced in the frozen seafood market, in which the environment is usually cold and damp.

The graph shows the brief stopover of the asymptomatic infector (first 5 min) and the high infection risk duration of the microenvironment (Tterminal). Although the quanta emission rate is consistent in each scenario, viral inactivation in the air (λ) and the deposition fraction in the respiratory tract (ω) are significantly influenced by the environment. The highest Rmax and the longest Tterminal might be an explanation as to why COVID-19 transmission was enhanced in the frozen seafood market, in which the environment is usually cold and damp (case 6).

Discussion and limitations

Seasonality of COVID-19

The thermal properties of ambient air, as well as relative humidity, may affect the transmissibility and viability of the virus (Rahimi et al. 2021). With the following changes in the seasons toward more temperate weather, it is of pivotal importance to know the role of environmental conditions in the transmission of the virus to raise awareness on the prevention of disease spread, and especially in a specific period (e.g., autumn, winter, spring, or summer) of a given geographical area (Coccia 2021b).

Seasonal cyclicity is a ubiquitous feature of acute respiratory infectious diseases, such as SARS and MERS (Killerby et al. 2018; Wang et al. 2020). For COVID-19, Mario Coccia (2022) found that the arithmetic mean of confirmed cases, hospitalizations of people, and admissions to intensive care units (ICUs) is significantly equal (p <0.01), despite the different interventions on COVID-19 in 2020 (strong lockdowns measures) and 2021 (health policy of vaccinations). Hence, findings here suggest that the COVID-19 pandemic is driven by seasonality and environmental factors, and the intervention strategies should be set up in summer and fully implemented during autumn and winter period (Coccia 2022b). Li et al. (2022) employed a meta-analysis to address these questions using results from 2813 published articles and observed that a rise in T would promote virus transmission in cold seasons, while significant negative influences in the warm season. RH significantly inhibited the COVID-19 transmission on the national scale. They suggested that the control measures should be developed according to local meteorological properties for the individual city (Li et al. 2022).

In this paper, we explored the complex relationship between ambient environment (T and RH) and COVID-19 transmission risk by using a modified W-R model. R0 is the key index for evaluating the transmission risk of infectious disease (Yuan et al. 2020). Figure 6 shows the comparison between the R0 simulated by the original W-R model and the simulated by the modified W-R model. The original W-R model fails to consider the influence of the environment on viral inactivation and droplet deposition and shows the same result for R0 in each case, always 3.18. However, the values simulated by the modified W-R model changed with the environmental parameters. As mentioned earlier, the middle-RH (45–60%) and high-T (>30 °C) scenarios were found to be the least favorable to the survival of the viruses (Colas et al. 2014, Harper G. 1961, Lin and Marr 2020a, Noti et al. 2013), which might be the reason why the disease infection risk is reduced in cases 1, 2, and 4. The droplet deposition risk in the airways under low-T and high-RH conditions is ignored by the original W-R model (Ishmatov 2020; Longest et al. 2011; Xi et al. 2015; Xie and Zhu 2020), which might be the reason why the disease infection risk is substantially underestimated in cases 5 and 6. In particular, the infection risk in a cold-damp scenario, such as a seafood cold store, can reach 6 times higher than in other scenarios at most.

Fig. 6

Comparison between the R0 simulated by the original W-R model and the simulated by the modified W-R model

That is, there does exist a certain range of environmental conditions that can significantly promote or inhibit the spread of COVID-19. In particular, the argument was further supported by these news that employees who worked in cold chain transportation were more vulnerable to being diagnosed with COVID-19 (https://baijiahao.baidu.com/s?id=1728401608245389772&wfr=spider&for=pc), and the covering of imported frozen seafood such as salmon (http://finance.sina.com.cn/china/gncj/2020-06-13/doc-iirczymk6749091.shtml) and shrimp ( https://www.best73.com/news/27562.html) tested positive for SARS-CoV-2 too.

Limitations

However, several limitations of our study should be acknowledged. First, this research only focused on meteorological/weather conditions that might influence the spread of COVID-19. Further studies on the effects of artificial lockdown interventions on disease transmission, such as masks and social distance, are needed. Second, the limitations of space prevent us from covering other important meteorological factors in this study, such as air pollutant, solar radiation, and wind speed.

Xu et al. and Shen et al. found that the reduction of total aerosols in the air and the contribution of the increase in fine-mode particles during the lockdown in Wuhan (Shen et al. 2021; Xu et al. 2020b). A study from California observed that air pollutants have a significant correlation with the COVID-19 epidemic (Bashir et al. 2020). Another study also indicated that an increase in particulate matter concentration causes more COVID-19 cases and mortality (Coccia 2020a, 2021a, 2021c, 2021d; Srivastava 2021).

Despite the ultraviolet light was constantly being used for sterilization, Song et al. found that the solar radiation showed a negligibility on the COVID-19 transmission (Runkle et al. 2020), even in the comparison of high- and low-altitude regions of China (Song et al. 2022). However, a previous study reported that high solar radiation might be indicated as the main climatic factor for the tropical state that suppresses the spread of COVID-19 in Brazil (Rosario et al. 2020).

Like solar radiation, wind speed also has its two sides. Some scholars argue that the wind would accelerate the virus’ travel from one location to another. This perhaps explains why the wind speed is higher, the number of cases increases in their study (Şahin 2020; Sarkodie and Owusu 2020). However, another study detected that high wind speed can be indicated as the main factor that suppresses the spread of COVID-19 (Rosario et al. 2020). Considering the transmission dynamics of the SARS-CoV-2 can also be explained by stagnation of atmospheric pollutants and viral agents, low wind speeds may promote a longer permanence of viral particles (e.g., the SARS-CoV-2) in polluted air, thus favoring the spread of COVID-19 in society (Coccia 2020c, a, 2021a, c, d, b, 2022b). Due to the highly uncertain and complex effect of environmental factors on the pandemic, the influence of meteorological factors on the newly confirmed COVID-19 cases varied greatly among existing studies, and no consistent conclusion can be drawn, which only can be controlled in laboratory experiments or model simulation (Islam et al. 2021). In order to improve the accuracy of transmission risk prediction, a more comprehensive model to analyze environment-mediated COVID-19 transmission should be employed in further research work, especially ventilation and wind speed (closely related to the HVAC system of indoor spaces). These conclusions are only a theoretical attempt since the transmission of COVID-19 is influenced by multiplicate factors.

Conclusions and prospects

The original W-R model ignored the influence of the environment on viral inactivation and droplet deposition, which might lead to the inappropriate distribution of emergency aid. In this study, a new approach is proposed to fill the gaps in knowledge when quantitatively evaluating the influence of environmental parameters on the spread of respiratory tract infections. We analyzed the currently available data on the effects of the environment on disease infection risk and identified the most critical parameters (T, RH, and P), and developed an alternative modified W-R model to address the aforementioned problems by proposing the dynamic virus deposition ratio (α) to assess the infection risk under varying environments. Then, it is applied to six typical exposure scenarios in daily life, reflecting how the environmental parameters are linked to viral inactivation and particle deposition, thus quantitatively affecting the transmission risk. This modified approach will enable healthcare workers and disinfection equipment to be assigned to the high-risk environment for COVID-19 prevention, before the infection outbreak, rather than passively waiting.

The results revealed that the temperature-risk relationship was negatively linear when the mean temperature was below 24 °C and became flat above 24 °C, indicating that a rise in T would significantly reduce the transmission of COVID-19 in the cold season (case 3/4/5/6), while inconspicuously negative influences in the hot season (case 1/2/4). So, the public and governments could not expect the high T to eradicate the epidemic in the hot season. The environment with middle RH has the lowest disease transmission probability. In addition, the infection risk under low-T and high-RH conditions is substantially underestimated, which should be taken seriously. In other words, high-risk environments, such as the frozen seafood market, are ill-equipped to effectively control the epidemic due to inadequate information, attention, and action, which may lead to higher risk or something even more catastrophic.

COVID-19 pandemic has threatened the public health and national security of countries seriously, which asked each country to make critical decisions to take advantage of important opportunities in the presence of highly restricted time in turbulent situations of emergency (Coccia 2020b). These measures for constraining future pandemic threats can be divided into two main categories, pharmaceutical interventions (such as vaccinations (Coccia 2022a) and health expenditures (Coccia 2021c)) and non-pharmaceutical interventions (such as lockdown policies (Askitas et al. 2021, Flaxman et al. 2020) and meteorology (Xie and Zhu 2020)). International experience of the COVID-19 pandemic crisis has shown that well-planned public policies and scientific and economic coordination policies are effective interventions for reducing high-risk aspects before and during an outbreak (Coccia 2021e, 2022c).

Recent studies show that accelerated transmission dynamics of COVID-19 is due to mainly to the mechanism of air pollution-to-human transmission (airborne viral infectivity) rather than human-to-human transmission (Coccia 2020c). The spread of COVID-19 is due to systemic causes: disease commonness (viral characteristics, transmission route), environment individuality (meteorological conditions, social context), and individual difference (disease susceptibility, immune response) (Coccia 2020c). According to a 2021 research review of the COVID-19 environment, more than 36% of them were directly or indirectly related to the indoor and outdoor environment, 16% to meteorological factors, 14% to fomites, and 34% to other factors (Rahimi et al. 2021). The still-growing body of literature on COVID-19 suggests the importance of environmental conditions in transmitting the disease beyond geographical/interpersonal borders. Despite various complicated and co-dependent factors during the COVID-19 transmission process, focusing on high-risk conditions will assist in solving these complex decision-making problems.

In general, this study suggests that sustainable science for crisis management in the presence of COVID-19 pandemic threat is also based on the cross-discipline study of meteorology, government, and socioeconomic factors, and not only on medical research. The high-risk condition may generate most part of the infections and indicate that selected containment measures are more sustainable social policy than a blind full lockdown, which provides useful information for policymakers if the COVID-19 coexists with humans for a long time.

(DOCX 772 kb)