Global COVID-19 pandemic trends and their relationship with meteorological variables, air pollutants and socioeconomic aspects.

Meteorological variables, air pollutants, and socioeconomic factors are associated with COVID-19 transmission. However, it is unclear what impact their interactions have on COVID-19 transmission, whether their impact on COVID-19 transmission is linear or non-linear, and where the inflexion points are. This study examined 1) the spatial and temporal trends in COVID-19 monthly infection rate of new confirmed cases per 100,000 people (Rn) in 188 countries/regions worldwide from March to November 2020; 2) the linear correlation between meteorological variables (temperature (T), rainfall (R), wind speed (WS), relative humidity (RH), air pressure (AP)), air pollutants (nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), ozone (O3)) and socioeconomic aspects (population density (PD), gross domestic product per capita (GDP), domestic general government health expenditure per capita (GHE)) and Rn, and 3) the interaction and non-linear effects of the different variables on Rn, based on GeoDetector and Boosted regression tree. The results showed that the global Rn had was spatially clustered, and the average Rn increased From March to November 2020. Global Rn was negatively correlated with meteorological variables (T, R, WS, AP) and positively correlated with air pollutants (NO2, SO2, O3) and socioeconomic aspects (GDP, GHE). The interaction of SO2 and O3, SO2 and RH, and O3 and T strongly affected Rn. The variables effect on COVID-19 transmission was non-linear, with one or more inflexion points. The findings of this work can provide a basis for developing a global response to COVID-19 for global sustainable development.

Introduction

COVID-19 pandemic has become an important challenge to global sustainable development (Naidoo and Fisher, 2020; Yin et al., 2021). On January 27, 2021, the global number of confirmed COVID-19 cases had surpassed 100 million (https://coronavirus.jhu.edu/map.html), representing more than 1/70th of the global population. COVID-19 pandemic created an adverse effect on human life, economy, environment, energy, and transportation (Nundy et al., 2021; Yao et al., 2020), which made sustainable development goals are being derailed (SDGs) (Naidoo and Fisher, 2020). Faced with the continuous transmission of COVID-19, the key to reducing virus spread and promoting SDGs back on track is establishing the drug and non-drug prevention and control strategies, which are suitably adapted to each country's reality (Zhao et al., 2020; Li et al., 2020). Establishing a non-drug prevention and control strategy relies on identifying variables that may affect the spatial and temporal COVID-19 transmission, especially identifying the effect of natural conditions and socioeconomic aspects on COVID-19 transmission (Han et al., 2021).

Meteorological variables, air pollutants, and socioeconomic aspects are critical links affecting COVID-19 transmission (Bolaño-Ortiz et al., 2020; Bontempi et al., 2020; Copiello and Grillenzoni, 2020; Fu et al., 2021; Van Doremalen et al., 2020). Temperature (T), rainfall (R), relative humidity (RH), wind speed (WS), and air pressure (AP) affect the virus activity in the air and objects surface (Chin et al., 2020), and affect human activities (Van Doremalen et al., 2020). Air pollution increases human vulnerability (e.g., respiratory problems) to the virus effects and affects their transmissibility (Andrée, 2020; Brandt et al., 2020; Conticini et al., 2020). Socioeconomic aspects have also been identified as linked to the COVID-19 pandemic (Copiello and Grillenzoni, 2020). In addition, the interaction between meteorological variables, air pollutants and socioeconomic aspects further complicates the transmission mechanism of COVID-19 (Coccia, 2020; Lin et al., 2020a; Srivastava, 2021). For example, among meteorological factors, wind speed may affect air humidity, further restricting SARS-CoV-2 virus activity (Chin et al., 2020; Lin et al., 2020a), affecting the stagnation and diffusion of particulate matters mixed with the SARS-CoV-2 virus (Rendana, 2020). The temperature may interact with air pollutants to affect human susceptibility to the virus (Sillman and Samson, 1995; Domingo and Rovira, 2020). Air pollutants affect not only human susceptibility to the SARS-CoV-2 virus (Brandt et al., 2020; Conticini et al., 2020; Frontera et al., 2020a; Marquès and Domingo, 2022), but also affects regional meteorological conditions and human activity intensity (Fenger, 2009). Moreover, socioeconomic factors are closely related to regional environmental conditions and medical levels (Copiello and Grillenzoni, 2020). The complex interaction of variables makes slowing COVID-19 transmission be a challenge.

Therefore, to establish a non-drug prevention and control strategy, it is necessary first to identify the linkages between the COVID-19 transmission and meteorological variables, air pollutants and socioeconomic aspects (Qu et al., 2020; Zhao et al., 2020). Determine whether the linkage between them is linear or non-linear, and whether there is an inflexion point at which the impact of variables on COVID-19 transmission shifts rapidly (Chien and Chen, 2020; Xie and Zhu, 2020). More importantly, the impact of interactions between variables on COVID-19 transmission needs to be analyzed to clarify the transmission path of COVID-19 and develop countermeasures (Lin et al., 2020a; Srivastava, 2021). Scholars need to pay more attention to these two aspects. This study aims to understand the evolution of global COVID-19 transmission from March to November 2020 and to analyze the interaction and non-linear effects of variables that can affect it - meteorological variables (T, R, WS, RH, AP), air pollutants (nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), ozone (O3)) and socioeconomic aspects (population density (PD), gross domestic product per capita (GDP), domestic general government health expenditure per capita (GHE)).

Materials and methods

Data acquisition and processing

In order to eliminate the effect of vaccination campaigns on identifying COVID-19 transmission mechanisms, which began in many countries in December 2020 (https://covid19.who.int/), this study focused on the period from March (COVID-19 be declared a pandemic) to November 2021. The data of COVID-19 cases, basic data, meteorological variables, air pollutants, and socioeconomic used in this study are described in Table 1 . The number of newly confirmed cases in a given period only reflects the absolute number of newly confirmed cases. It does not strip the effect of the total population on newly confirmed cases (Jahangiri et al., 2020). Therefore, this study used the rate of new confirmed cases per 100,000 people (R ) to reflect the COVID-19 transmission status. R was obtained by dividing the number of new confirmed cases per month by the country's total population and multiplying by 100,000 (Copiello et al. 2020). Meteorological and air pollutants data were monthly averaged (March to November 2020) at the country level using the zonal statistics tool in ArcGIS 10.5. In order to reflect the socioeconomic situation of each country, the 5-year average (2015–2019) was used for GDP per capita. The same process was applied concerning domestic general government health expenditure per capita (2013–2017).

Table 1

Data of COVID-19 cases, basic data, meteorological variables, air pollutants, and socioeconomic used in the study.

Dimension	Variable	Units	Period	Data type	Resolution	Source
COVID-19 cases	Newly confirmed cases	Number	March 2020–November 2020	–	Country	World Health Organization (https://covid19.who.int/)
Basic data	Administrative boundaries	–	2020	Vector	Country	Resources and Environmental Science and Data Center of Chinese Academy of Sciences (http://www.resdc.cn/)
Meteorological variables	Temperature	°C	March 2020–November 2020	Raster	0.1°	ERA5-Land monthly average data (https://cds.climate.copernicus.eu/)
	Rainfall	mm
	Wind speed	m/s
	Relative humidity	%
	Air pressure	Pa
Air pollutants	Nitrogen dioxide	mol/m2	March 2020–November 2020	Raster	7 km × 3.5 km	Sentinel-5P image in the Google Earth Engine platform
	Sulfur dioxide
	Carbon monoxide
	Ozone
Socioeconomic aspects	Population	Number	2019	–	Country	World Bank public database (https://data.worldbank.org.cn/)
	GDP per capita	US$/capita	2015–2019	–	Country
	Domestic general government health expenditure per capita	US$/capita	2013–2017	–	Country
	Population density	Persons/km2	2020	Raster	Year, 1 km	Gridded Population of the World (GPW) v4.0 (https://sedac.ciesin.columbia.edu/data/collection/gpw -v4)

Analysis methods

Spatial autocorrelation analysis

Spatial autocorrelation analysis can reveal the R pattern and COVID-19 transmission (Han et al., 2021), which can be applied globally and locally. Global spatial autocorrelation assessed the R spatial pattern in all the studied areas. Global Moran's I is one of the most commonly used indices for global spatial autocorrelation analysis (Cliff and Ord, 1981). It is calculated according to the following formula:where W represents the weight of the spatial relationship between i and j, (X -‾X) represents the deviation of the attribute value of element i from its mean, and n is the total number of elements. The value of Global Moran's I ranges from [−1,1]. If the index is close to 0, it indicates that the distribution of R has a random pattern. If the index is close to −1, the R distribution is dispersed. If the index is close to 1, the pattern of R is clustered.

Global spatial autocorrelation analysis can reveal the degree of spatial dependence in the distribution of R but cannot identify the local differences. Therefore, we applied Local Moran's I (Anselin, 1995), which measure the Local Indicators of Spatial Association (LISA) to explore the correlation of R between the different countries. Local Moran's I was calculated according to the formula. (2):where W represents the weight of the spatial relationship between i and j, (X - ) represents the deviation of the attribute value of element i from its mean, and n is the total number of elements. Based on Local Moran's I calculations, the LISA aggregation can classify R distribution into two types of positive correlations: High-High and Low-Low; two types of negative correlations: High-Low and Low-High; and non-significant correlations. Positive correlations indicate that R is in line with the surrounding region's trend observed. High-High type indicates high in the middle high periphery, and the Low-Low type indicates low in the middle low periphery. Negative correlations indicate that R is opposite to the trend of the surrounding region. High-Low type indicates low in the middle high periphery, and low-high type indicates high in the middle low periphery. Non-significant correlations represent that R has a random pattern. Significant correlations were considered at a P < 0.05. Both Global and Local Moran's I was computed using ArcGIS 10.5.

Statistical analysis

The Spearman correlation was applied to identify the linear correlation between variables and R . P < 0.01 and P < 0.05 indicated a significant correlation at 99% and 95% confidence levels. The multicollinearity test was applied to eliminate the effect of redundant variables on R . The variables were eliminated if the variance inflation factor (VIF) was greater than 10. Both Spearman correlation analysis and multicollinearity test were implemented in SPSS 24.0.

GeoDetector

GeoDetector (http://www.geodetector.org/) is a statistical method used to detect the driving forces affecting the spatial differentiation of geographic phenomena (Wang et al., 2010). It reflects the explanatory ability of independent variables to dependent variables through the sum of variances ratio under different classifications of independent variables to the sum of variances of dependent variables. Factor detector and interaction detector is the sub-detectors of GeoDetector. Factor detector was applied to detect the explanatory ability of independent variables to the spatial differentiation of the dependent variables. Their magnitude can be measured by q value (Wang et al., 2010). The calculation formula of q value is as follows: where h represents the category of the independent variable (h = 1, 2 …, L); N and N are the number of units in category h and the entire region, respectively; σ2 and σ2 represent the variance of the dependent variable in the category h and the entire region, respectively (Eq. (3)). SSW represents the sum of the category variances of the independent variable, and SST represents the total variance of the dependent variable in the entire region (Eq. (4)). q indicates the ability of independent variables to explain dependent variables, q∈[0,1], the larger the q value, the greater the explanatory ability of the independent variable. In addition, the GeoDetector can test the significance of the q value.

The Interaction detector has a key advantage that can recognize the interaction between two independent variables compared to conventional statistical methods. The q value reflects the ability of bivariate interactions to explain the dependent variable. The direction and mode of interaction between the two variables can be assessed by comparing the single variable q value and the double variables q value. Since the GeoDetector requires the input data of the independent variables to be categorical data, we used the natural breakpoint method to classify the selected independent variables into five categories for modelling (Peng et al., 2019).

Boosted regression tree

Boosted regression tree (BRT) (Elith et al., 2008) is a regression algorithm built on a traditional classification regression tree algorithm (CART). BRT builds multiple regression trees based on self-learning and multiple random selections. The process of multiple fitting can gradually reduce the error of model fitting and then steadily improve the simulation accuracy of regression trees. During operation, BRT will be iterated and repeated several times. Some data will be randomly selected to analyze the effect degree of independent variables on dependent variables. The remaining data will be used for cross-validation. Finally, the average value of the generated multiple regression trees will be calculated. In the study, BRT was assessed using R version 3.6.1 (R Core Team, Vienna, Austria) using the “gbm” package and BRT equation package (Elith et al., 2008). The distribution function was set as “Gaussian” in the parameter setting, and the learning rate was 0.05. 80% of the data was extracted for training and 20% for validation. A Pearson correlation between observed values and fitted ones of the training set and validation set. Also, a Root mean square error (RMSE) obtained from the cross-validation procedure was calculated to test the fitting effect of the model. The closer the Pearson correlation coefficient (P < 0.01) is to 1, and the closer the RMSE is to 0, the better is the model fitting.

The results of the BRT model can reflect the relative contribution of the variables to R through response curve plots. The horizontal coordinate of the response curve plots shows the variables change, and the vertical coordinate is the marginal effect on R with the variables change. When the marginal effect value is less than 0, effect of increasing variable on R is weakened. When the marginal effect value is equal to 0, the effect of the variable increasing on R is unchanged, and when the marginal effect value is greater than 0, the effect of increasing variable on R is enhanced.

Results

Spatial and temporal transmission of COVID-19

From March to November 2020, the global average R was an increase (Table 2 ). The average R was 28.62 in March and 394.59 in November. Between March and November 2020, the countries with the high level of R were San Marino, Qatar, Bahrain, Aruba, Andorra and Luxembourg. Andorra has the highest R , which was 3498.74. The coefficient of variation of R ranged between 155.38% and 276.62%. Before July 2020, there were large differences in R among different countries, with the coefficient of variation of R higher than 200%. From July to November, the difference was reduced.

Table 2

Descriptive statistic of global R from March to November 2020.

		March	April	May	June	July	August	September	October	November
Minimum	Rn	0.01	0.02	0.01	0.03	0.01	0.03	0.01	0.01	0.01
	Country	Papua New Guinea	Yemen	China	Viet Nam	Laos	Laos	Laos	Laos	Yemen
	Continent	Oceania	Asia	Asia	Asia	Asia	Asia	Asia	Asia	Asia
Maximum	Rn	664.50	983.46	1477.82	1406.89	884.49	1626.31	1903.79	3498.74	2716.42
	Country	San Marino	San Marino	Qatar	Qatar	Bahrain	Aruba	Aruba	Andorra	Luxembourg
	Continent	Europe	Europe	Asia	Asia	Asia	South America	South America	Europe	Europe
Mean		28.62	51.63	56.21	74.71	96.82	123.84	145.17	290.71	394.59
Standard deviation		79.17	105.90	139.55	170.30	164.65	202.27	241.93	485.14	613.12
Coefficient of variation		276.62%	205.11%	248.27%	227.95%	170.06%	163.33%	166.65%	166.88%	155.38%

As shown in Fig. 1 , in March 2020, the high level of R was mainly distributed in European countries and the United States, while in April, R was also observed to increase in Brazil, Peru, and Chile. In May, R increased in Central Asian countries and some African countries. There was an R increase in South America, North America, Europe, and some African countries from June to November.

Fig. 1

Global spatial distribution of R from March to November 2020.

Spatial pattern of COVID-19 transmission

In all the studied months (Table 3 ), R was significantly clustered. The Global Moran's I of R ranged from 0.12 to 0.71 (P < 0.01). However, there was a decrease of R spatial aggregation between March and May and increased from June to August. High Global Moran's I were identified in October (0.50) and November (0.71), while the lowest was registered in May (0.12).

Table 3

Global Moran'I of the spatial distribution of R from March to November 2020.

	March	April	May	June	July	August	September	October	November
Moran'I	0.34	0.29	0.12	0.21	0.25	0.30	0.23	0.50	0.71
Z score	6.98	6.04	2.85	4.34	4.75	5.89	4.46	9.72	13.27
P value	<0.001	<0.001	0.004	<0.001	<0.001	<0.001	<0.001	<0.001	<0.001

Note: Temperature (T), rainfall (R), wind speed (WS), relative humidity (RH), air pressure (AP), nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), ozone (O3), population density (PD), gross domestic product per capita (GDP), domestic general government health expenditure per capita (GHE).

The R spatial pattern varied with the virus transmission (Fig. 2 ). Local Moran's I results showed that the R spatial distribution was mainly High-High and Low-Low clustering types. The High-High cluster types were distributed in Europe from March to April and October to November. From May to September, the High-High cluster pattern was and mostly observed in South America. Low-Low clusters types were identified between March and November in Africa and Asia. Among other clustering types with less distribution, High-Low clustering types were only observed in South Africa in May and Russia and Libya in August. Low-High clustering types were sporadically distributed in Europe, Africa and South America.

Fig. 2

Local aggregation of R from March to November 2020.

Correlation between variables and COVID-19 transmission

Spearman correlation (Table 4 ) showed that T (−0.27), R (−0.09), WS (−0.07), AP (−0.07) had a significant negative correlation with R (P < 0.01). Contrary, NO2 (0.34), SO2 (0.24), O3 (0.29), GDP (0.45) and GHE (0.43) were significantly positive correlated with R (P < 0.01). RH, CO, and PD were not significantly correlated with R .

Table 4

Spearman correlation coefficients results.

	Rn	T	R	WS	RH	AP	NO2	SO2	CO	O3	PD	GDP	GHE
Rn	1.00
T	−0.27**	1.00
R	−0.09**	0.08**	1.00
WS	−0.07**	0.12**	−0.34**	1.00
RH	−0.03	−0.01	0.78**	−0.21**	1.00
AP	−0.07**	0.33**	0.08**	0.58**	0.31**	1.00
NO2	0.34**	−0.21**	−0.36**	−0.13**	−0.36**	−0.11**	1.00
SO2	0.24**	−0.43**	−0.13**	0.05	−0.03	0.04	0.35**	1.00
CO	0.03	0.04	−0.18**	−0.06*	−0.12**	0.08**	0.36**	0.19**	1.00
O3	0.29**	−0.48**	−0.22**	0.15**	−0.24**	−0.03	0.45**	0.34**	0.16**	1.00
PD	0.02	0.15**	0.09**	0.06*	0.15**	0.31**	0.32**	0.06*	0.05*	0.03	1.00
GDP	0.45**	−0.42**	−0.07**	0.27**	0.03	0.20**	0.26**	0.40**	−0.10**	0.48**	0.07**	1.00
GHE	0.43**	−0.36**	−0.03	0.28**	0.09**	0.27**	0.23**	0.37**	−0.08**	0.44**	0.11**	0.96**	1.00

Note: *P < 0.05, **P < 0.01, Temperature (T), rainfall (R), wind speed (WS), relative humidity (RH), air pressure (AP), nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), ozone (O3), population density (PD), gross domestic product per capita (GDP), domestic general government health expenditure per capita (GHE).

Among meteorological variables, T, R, RH and AP were negatively correlated with air pollutants, while T and R were negatively correlated with socioeconomic aspects. RH and AP were positively correlated with socioeconomic aspects. NO2, SO2 and O3 were positively correlated with socioeconomic aspects but negatively correlated with meteorological variables. CO was negatively correlated with meteorological and socioeconomic variables. The multicollinearity test (Table S1) showed that the VIF were between 1.14 and 7.61 (VIF<10). Therefore, all the variables were considered in the subsequent analysis.

Single variable effects on COVID-19 transmission

Factor detector results (Table 5 ) showed that all variables had a significant effect on R . Nevertheless, SO2, GDP, NO2, and GHE had a relatively high effect. They could explain 15.70%, 9.20%, 8.90% and 8.20% of the variance, respectively.

Table 5

q value of single variable effect on R.

	T	R	WS	RH	AP	NO2	SO2	CO	O3	PD	GDP	GHE
q value	0.09**	0.01**	0.01**	0.02**	0.00	0.09**	0.16**	0.02**	0.04**	0.01*	0.08**	0.09**
P value	<0.001	0.002	0.006	<0.001	0.523	<0.001	<0.001	<0.001	<0.001	0.014	<0.001	<0.001

Note: *P < 0.05,**P < 0.01, Temperature (T), rainfall (R), wind speed (WS), relative humidity (RH), air pressure (AP), nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), ozone (O3), population density (PD), gross domestic product per capita (GDP), domestic general government health expenditure per capita (GHE).

Interaction effects on COVID-19 transmission

The interaction detector can reflect the interaction of two variables on R (Table 6 ). The q value of interaction was higher in O3∩SO2 (0.33), SO2∩RH (0.30), O3∩T (0.28), NO2∩RH (0.27) and GHE∩SO2 (0.26).

Table 6

q value statistics of the effect of variables interaction effect on R.

	T	R	WS	RH	AP	NO2	SO2	CO	O3	PD	GDP	GHE
T	0.09
R	0.11	0.01
WS	0.14	0.04	0.01
RH	0.23	0.16	0.04	0.02
AP	0.12	0.03	0.02	0.04	0.00
NO2	0.21	0.10	0.12	0.27	0.12	0.09
SO2	0.22	0.19	0.24	0.30	0.20	0.23	0.16
CO	0.15	0.04	0.05	0.05	0.06	0.14	0.26	0.02
O3	0.28	0.06	0.06	0.07	0.05	0.17	0.33	0.06	0.04
PD	0.13	0.03	0.03	0.04	0.02	0.12	0.17	0.04	0.06	0.01
GDP	0.17	0.11	0.14	0.15	0.15	0.16	0.23	0.10	0.19	0.11	0.08
GHE	0.18	0.11	0.13	0.17	0.14	0.17	0.26	0.11	0.22	0.11	0.10	0.09

Note: Temperature (T), rainfall (R), wind speed (WS), relative humidity (RH), air pressure (AP), nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), ozone (O3), population density (PD), gross domestic product per capita (GDP), domestic general government health expenditure per capita (GHE).

Non-linear effects on COVID-19 transmission

Accuracy verification results of BRT (Table S2) showed that the Pearson correlation coefficient between the observed and fitted values of the training set and verification set was 0.86 and 0.71 (P < 0.01), respectively. The RMSE was 0.03. The correlation coefficients were close to 1, and RMSE was close to 0. The model had a good fitting effect.

Results of the BRT (Fig. 3 ) showed that the average explanatory power of air pollutants, meteorological variables and socioeconomic aspects were 10.23%, 9.12%, and 4.57%, respectively. The contribution of 12 selected variables to the change in R , were in the decreasing order: SO2 (21.80%), RH (17.20%), WS (15.80%), O3 (8.40%), NO2 (8.10%), T (6.90%), GHE (6.40%), PD (4.10%), R (3.40%), GDP (3.20%), CO (2.60%), AP (2.30%). The variables effect on R was non-linear. The SO2, RH, NO2, GHE, PD, R and GDP effect on R changed from negative to positive with the increasing values. On the contrary, the effect of T, WS, O3, CO and AP on the R changed from positive to negative with the increasing values of variables.

Fig. 3

Marginal effect of variables change on R.

When the atmospheric SO2 concentration was less than 3.5 × 10−4 mol/m2, SO2 had a negative effect on R . However, when the concentration was higher than 3.5 × 10−4 mol/m2, the effect was positive. At concentrations higher than 6.5 × 10−4 mol/m2, the effect was stable. When the atmospheric O3 concentration was less than 0.133 mol/m2, O3 had a positive effect on R , while concentration higher than 0.133 mol/m2, the effect of O3 turned negatively. When atmospheric NO2 concentration was less than 0.4 × 10−5 mol/m2, NO2 negatively affected R , but the negative effect turned to positive when concentration higher than 0.4 × 10−5 mol/m2. When atmospheric CO concentration was less than 0.034 mol/m2, CO positively affected R . However, when the concentration was higher than 0.034 mol/m2, CO had a slightly negative effect on COVID-19 transmission.

Among meteorological variables, when the RH was lower than 80%, RH negatively affected R , and the negative effect was relatively smooth. While atmospheric RH was higher than 80%, the effect turned positive until it stabilizes higher than 90% after the RH. When WS was slower than 3.5 m/s, WS had a positive effect on R , but when the WS was faster than 3.5 m/s, WS had a negative effect on R and turned stable. There were multiple inflexion points for the effect of T on R . When T below 6 °C, between 9 °C and 12 °C, or higher than 33 °C, T had a positive effect on R . However, when T between 6 °C and 9 °C, or between 12 °C and 33 °C, T shifted to a negative effect on R . When R was less than 100 mm, R has a negative effect on R , which gradually decreased with the increase of R. When R was higher than 100 mm, the marginal effect of R on R was close to 0. Overall, R was slightly affected by changes in R. When AP was less than 8300 Pa, R was positively affected by AP, and the positive effect reached the maximum when AP was about 8200 Pa. In contrast, with an AP higher than 8300 Pa, the effect of AP on R turned negatively.

When GHE was less than 3500 US$/capita, the effect on R was negative, while GHE was more than 3500 US$/capita, the effect on R was positive and stabilized when GHE was more than 4700 US$/capita. When PD was less than 700 persons/km2, the negative effect of PD on R gradually weakens, and after PD was higher than 700 persons/km2, PD been a smooth, positive effect. When GDP lower than 40,000 US$/capita, GDP had a negative effect on R , while GDP higher than 40,000 US$/capita, GDP had a slightly positive effect.

Discussion

Spatial and temporal characteristics of the global COVID-19 pandemic

Understanding the spatial distribution and monthly characteristics of R is one way to understand the transmission mechanism of COVID-19. Regarding spatial distribution characteristics, the countries with the highest R between March and November 2020 were mainly San Marino, Qatar, Bahrain, Aruba, Andorra and Luxembourg (Table 2). These countries with limited resources and are more dependent on exchanges with other countries, making it challenging to adopt effective border control measures (Bontempi, 2020a). As a result, they are more vulnerable to COVID-19 transmission, resulting in a high number of new confirmed cases. Countries with low R were mainly located in East Asia and Africa. In East Asia countries such as China, South Korea and Japan, the COVID-19 pandemic was quickly controlled due to the strict measures imposed (Shaw et al., 2020). In addition, in the East Asian culture that emphasizes collectivism, the social isolation policy implemented by the government is also easier to be implemented (An and Tang, 2020). Africa is the second most populous continent globally, but only a few countries showed high R . Limited testing, a much younger population, climatic differences may affect the trajectory of the COVID-19 pandemic in Africa (Kandel et al., 2020; Maeda and Nkengasong, 2021). The spatial distribution of the global R varies monthly (Fig. 1). R observed a decline in the summer months in European countries, while they increased in the winter months in South America and South Africa in the southern hemisphere. Seasonal variations in R may be related to seasonal variations in human susceptibility to infectious diseases (Dowell, 2001).

Variables affecting global COVID-19 transmission

Effect of meteorological variables on COVID-19 transmission

Meteorological factors play a crucial role in the transmission and control of COVID-19 (Fu et al., 2021; Guo et al., 2021; Islam et al., 2021). The effects of T and RH on COVID-19 transmission are of great interest among meteorological variables (Auler et al., 2020; Ma et al., 2020). Our study showed that the effects of T and RH on COVID-19 transmission were non-linear. Below 12 °C, T increased COVID-19 transmission, while T between 12 °C and 33 °C reduced virus transmission. The inflexion point of RH effect is 80%, RH greater than 80% promotes the transmission of COVID-19 and vice versa. Some studies observed that T and RH were negatively correlated with COVID-19 transmission (Guo et al., 2021; Sarkodie and Owusu, 2020) and suggested that low T and dry environments increased COVID-19 transmission. The reason is that low T and humidity create conditions for the virus survival (Chin et al., 2020). Reduced T also inhibits the phagocytic function of human alveolar macrophages (Ma et al., 2020), making the human body more vulnerable to the COVID-19 virus. Several studies (Auler et al., 2020; Pani et al., 2020; Prata et al., 2020) hold a different view. They suggested that T or RH was positively correlated with the transmission of COVID-19. Two different viewpoints are due to the non-linear effects of T and RH on COVID-19 transmission. Lin et al. (2020a) suggest that the effect of T and humidity on COVID-19 is affected by the interaction between T and RH. High humidity promotes the transmission of COVID-19 when the T is reduced and decrease the transmission when the T is high.

WS and AP also affect COVID-19 transmission. Our results showed that when the WS was slower than 3.5 m/s, the transmission increased, while when the WS was faster than 3.5 m/s, there was a decrease in virus transmission. Wind may affect the time that virus is suspended and its dispersal distance (Bolaño-Ortiz et al., 2020). Droplets with viruses are more likely to be suspended in the air at reduced WS or settle on the surfaces that humans can contact. Contrary, high WS may blow away droplets and reduce their concentration in suspended particles (Rendana, 2020). High AP may shorten droplet suspension and reduce COVID-19 transmission (Lin et al., 2020b). When AP is lower than 8300 Pa, we found that AP may favour the atmosphere's droplets' suspension, increasing the likelihood of humans being infected. The reverse situation is observed if after AP is higher than 8300 Pa.

Our results showed that when monthly R was lower than 100 mm, there was a reduction in COVID-19 transmission. R may lead to a decrease in T, which in turn reduces the body blood supply, slowing the body to supply immune cells to the nasal mucosa (Shahzad et al., 2020; Wu et al., 2020). At the same time, R may also reduce the ability of solar radiation to kill the COVID-19 virus affecting the regional humidity, aggravating the virus transmission (Biasin et al., 2021).

Effect of air pollutants on COVID-19 transmission

Air pollution is a consequence of human activities, weather and topographical conditions. Air pollution is thought to be a key factor in coronavirus transmission and infection in two ways (Zhang et al., 2020): 1) COVID-19 is a respiratory infectious disease spread through the air. After COVID-19 is adsorbed onto atmospheric particles, it can spread and contribute to the increase of infected people; and 2) Air pollution increases the population vulnerability to virus effects since it weakens human immunity. Inhabitants that live near polluted areas have an immune system weak and are more vulnerable to respiratory or cardiovascular diseases. Therefore, the population exposed to high air pollution have a high probability of contracting and be affected by the virus (Brandt et al., 2020; Conticini et al., 2020).

High concentrations of SO2, O3, NO2, and CO are common in urban areas (Fenger, 2009). The correlation result of this study found a significant positive correlation between O3, NO2 and CO and R . O3 causes oxidative stress, and exposure to high concentrations reduces lung-lining fluid's antioxidant levels (Domingo and Rovira, 2020). PM2.5 decline caused by COVID-19 lockdown (He et al., 2020) may lead to the reactive uptake of HO2 radicals by aerosol particles, thus stimulating O3 production (Li et al., 2019). New produced O3 may contribute to COVID-19 transmission. Exposure to NO2 may induce airway epithelial cells to synthesize pro-inflammatory cytokines, leading to respiratory tract inflammation and reducing viruses' immunity (Persinger et al., 2002). High ambient NO2 may also be responsible for the extensive lung injury and acute effects of COVID-19 (Frontera et al., 2020a). The effects of CO on the human body stem from its ability to bind haemoglobin. High CO concentrations reduce the reactive oxygen species production and impair myoglobin function (Burnett et al., 1998), resulting in asphyxia-related death or hypoxic tissue damage (Zhao et al., 2019) and high human susceptibility COVID-19 virus.

However, unlike the research of Bashir et al. (2020) and Zhu et al. (2020), we found that SO2 is positively correlated with R . The SO2 pollution effect on COVID-19 transmission may have been overlooked. Long-term exposure to SO2 leads to cardiovascular and respiratory disease (Cole et al., 2020). The response of COVID-19 transmission to SO2 is critical since it has a high concentration in the atmosphere (Fenger, 2009).

Effect of socioeconomic aspects on COVID-19 transmission

Demographic and socioeconomic factors may also affect COVID-19 transmission (Bolaño-Ortiz et al., 2020; Bontempi et al., 2020). COVID-19 pandemic has challenged the disaster response capabilities of governments (Shaw et al., 2020). It is widely accepted that good socioeconomic and health conditions increase infectious diseases prevention (Fu et al., 2021). Coccia (2021) found that high health expenditures can effectively reduce the fatality rate of the COVID-19 pandemic. Nevertheless, this study found that GDP and GHE can reduce the infection rate of COVID-19 when they are at a low level, but a high level not (Fig. 3). Indicated that better medical facilities could avoid worsening symptoms in people after infection but can not contain the virus from spread effectively. Individual and social behavioural sciences, which can be very heterogeneous in different countries, may have affected this phenomenon (Bavel et al., 2020). Mendenhall (2020) argued that, compared with simple socioeconomic conditions, other factors such as hypertension, diabetes, respiratory diseases, systemic racism, distrust of science and leadership, and decentralized health care system in some countries could be behind the high COVID-19 transmission.

Previous works found a positive correlation between population density and the number of COVID-19 infections in Beijing, China (Han et al., 2021). High population density will increase the possibility of person-to-person contact, thereby aggravating the transmission of COVID-19. However, here we found no significant correlation between R and population density at a global scale, although the interaction of population density with T and SO2 significantly affected COVID-19. The effects of population density on the COVID-19 transmission is not consistent over different geographical areas, which are likely affected by local environmental, cultural, and policy factors (Bavel et al., 2020).

Non-linear relationship between environmental variables and COVID-19 transmission

Linear correlations analysis may be limited to identify correctly the variables that affect CIVID-19 transmission. For example, a negative correlation between T and humidity and COVID-19 has been observed in previous works (Guo et al., 2021; Sarkodie and Owusu, 2020), while others reported the opposite trend (Doğan et al., 2020; Islam et al., 2020). In this context, the relations between COVID-19 and meteorological variables, air pollutants and socioeconomic aspects in different regions may be affected by the non-linear effect.

Unlike linear correlation analysis, GeoDetector has no linear assumptions (Wang et al., 2010). Our result of the factor detector was different from the correlation analysis. For instance, SO2 had the strongest explanatory power but had a low correlation coefficient (Table 4 and Fig. 3). This show that SO2 had a non-linear effect on COVID-19 transmissibility. Previous works based on linear correlations may have overlooked the SO2 effect. In addition, although PD and RH not significantly correlated R , they significantly correlated most meteorological and socioeconomic variables. Therefore, PD and RH may indirectly affect the transmission of COVID-19 by affecting the likelihood of human exposure to the virus and its activity (Han et al., 2021; Lin et al., 2020a).

The effect of interactions between variables on COVID-19 transmission cannot be ignored. Effects of bivariate interactions on COVID-19 transmission increased compared to single variables, with SO2 and O3, SO2 and RH, O3 and T being the most critical influences (Table 6). As atmospheric oxidants, O3 can efficiently promote SO2 oxidation in aerosol media and result in a high level of particulate sulfate formation particle contaminant in the air (Wang et al., 2016). Particle contaminants can promote the spread of the SARS-CoV-2 virus and reducing human immunity (Bontempi, 2020b, Marquès and Domingo, 2022). The interaction between SO2 and RH also significantly affected COVID-19 transmission. RH affected heterogeneous reactions of SO2 with soot particles. The reaction progress generated by the ketone and sulfate species can lead to inflammation and cellular damage (He et al., 2017) and thus increased human susceptibility to COVID-19 (Cole et al., 2020; Marquès and Domingo, 2022). The interaction of O3 and T also affect COVID-19 transmission. O3 is produced in the troposphere by photochemical oxidation of CO, methane, and non-methane volatile organic compounds by the hydroxyl radical in the presence of reactive nitrogen oxides (Park, 2019). High T promotes the production of O3 (Sillman and Samson, 1995), which in turn forms O3 pollution and affects human health (Domingo and Rovira, 2020). In addition, Rendana (2020) indicated that low WS promote stagnation of particles, thus facilitating COVID-19 transmission. This study also confirmed that the interaction between WS and air pollutants strongly affects COVID-19 transmission (Table 6).

The effects of meteorological variables, air pollution, and socioeconomic aspects on COVID-19 transmission have an inflexion point. Some scholars tried to identify the T threshold effect on COVID-19 transmission. Chien and Chen (2020) found that the effect of T on COVID-19 transmission has a peak at 20–25 °C. Xie and Zhu (2020) observed that when the T is below 3 °C, the number of confirmed cases increased by 4.86% for every 1 °C decrease. Hoang and Tran (2021) identified a 9% increase in confirmed cases of COVID-19 when the T was below 8 °C. BRT results showed that all variables, not only T, have a non-linear effect on R , with positive and negative effects, and their effects have one or more inflexion points (Fig. 3). Identifying the inflexion points of variables effect on COVID-19 created the possibility of a more targeted response to COVID-19.

Potential pathways of variables effect on global COVID-19 transmission

Meteorological factors, air pollutants and socioeconomic factors affect COVID-19 transmission in three aspects: virus activity, virus transmission and human susceptibility.

COVID-19 virus activity is sensitive to changes in T and humidity. A laboratory study published by Chin et al. (2020) has reported that the virus is very stable at 4 °C but is sensitive to heat. The virus survival time decreases with increasing T. High T destroys the lipid layer of the virus, reducing its activity, stability and infection potential, and therefore, the transmission. Chan et al. (2011) reported that SARS-CoV could survive for more than 5 days on smooth surfaces at T between 22 and 25 °C and 40%–50% RH. However, when the T or RH increased, the viability of SARS-CoV rapidly decreased. Therefore, SARS-CoV-2 may be unstable in a high T and humidity environment and affect the human infection rate. The main routes of COVID-19 transmission are 1) direct or indirect contact with infected subjects, (2) large droplets dispersed by coughing or sneezing, and (3) inhalation of small airborne particles remaining in the air (Domingo et al., 2020; Qu et al., 2020).

The likelihood of becoming infected after exposure to the virus may be affected by T and air pollutants. Human immunity decreases at reduced T (Ma et al., 2020). In addition, exposure to respiratory inflammation and cardiovascular disease caused by air pollutants may increase the probability of COVID-19 infection and the virus effects on human health (Brandt et al., 2020; Conticini et al., 2020).

Gaps in research, limitations, and uncertainties

Although our study provides some evidence for the relationship between the non-linear and interactive effects of meteorological variables, air pollutants, and socioeconomic effects on global COVID-19 transmission, there are important limitations (e.g.,.data availability, resolution, and quality). This represents an impediment to conducting conclusive studies on factors affecting COVID-19 transmission (Sarkodie and Owusu, 2021; Zhou et al., 2020). Currently, publicly available COVID-19 case data are only region-level counts. The COVID-19 case data we used is the national scale and cannot accurately reflect the spatial distribution characteristics of COVID-19 cases, and some information might have been lost. This occurred as well in meteorological, air pollutants and socioeconomic data. For socioeconomic aspects, we only selected PD, GDP and GHE, and there was a lack of consideration of national policies regarding the lockdown measures, cultural background, and personal background of the case (presence of historical illness, activity status, and smoking and alcohol use) (Hu et al., 2021). Second, to better understand COVID-19 dispersion patterns, an interdisciplinary, multidimensional approach should be encouraged to draw more definitive conclusions (Bontempi et al., 2020). Our study analyzed the effects of meteorological, air pollutants, and socioeconomic variables on COVID-19 using only geographic and statistical methods. Behavioural and social sciences that consider social context, scientific communication, and ethical norms are also important to assess COVID-19 transmission (Bavel et al., 2020). Third, as the COVID-19 vaccine becomes widespread worldwide, virus transmission becomes more complex (Wouters et al., 2021). We did not consider vaccination data, which may pose some restrictions to our analysis. Nevertheless, although our study has some limitations, it contributes to understanding the effects of environmental and socioeconomic variables on COVID-19 transmission and can inform decision-making on COVID-19 prevention and control at a global scale.

Conclusion

During the study period, global R increased, and COVID-19 was spreading around the world. High R is mainly distributed in North America, Europe, and South America. Global R is spatially correlated, and the correlation increased with the COVID-19 pandemic intensification. Global R was negatively correlated with meteorological variables (T, R, WS, AP) and positively correlated with air pollutants (NO2, SO2, O3) and socioeconomic aspects (GDP, GHE). However, the linear correlation method cannot fully reveal the relationship between variables and COVID-19 transmission. Interactions between variables play a key role in affecting COVID-19 transmission. The interaction of SO2 and O3, SO2 and RH, and O3 and T strongly affected R , more than each variable individually. Consider the non-linear effect, the average explanatory power of meteorological variables, air pollutants, and socioeconomic aspects for COVID-19 transmission in BRT were 9.12%, 10.23%, and 4.57%, respectively. Despite limitations in data availability, behavioural and social science analysis, and vaccination data prevented us from drawing more definitive conclusions, the variables interaction and non-linear effect on COVID-19 transmission we studied to some extent shed light on the complex mechanisms of COVID-19 transmission and provide some basis for developing a global policy response. Based on the potential transmission routes of COVID-19, social distancing and environmental disinfection are currently effective means to slow COVID-19 transmission, and controlling environmental pollution can improve people's overall ability to resist respiratory infections.

Author contributions section

Yi Han: Methodology, Data Formal analysis, Visualization, Writing – original draft. Wenwu Zhao: Conceptualization, Formal analysis, Writing – review & editing. Supervision. Paulo Pereira: Writing – review & editing.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.