Skip to main content

Correlations between geomagnetic field and global occurrence of cardiovascular diseases: evidence from 204 territories in different latitude

Abstract

Background

The correlation between stable geomagnetic fields and unstable geomagnetic activities with mortality, incidence, and prevalence of cardiovascular diseases (CVDs) remains ambiguous.

Method

To investigate the correlations between geomagnetic field (GMF) intensity and geomagnetic disturbance (GMD) and CVDs events in global, long-period scale, global and 204 countries and territories were included on the base of 2019 Global Burden of Disease study (GBD 2019). Data of GMF intensity, GMD frequency, CVDs events, weather and health economic indicators from 1996 to 2019 of included locations were collected. Linear regression and panel data modelling were conducted to identify the correlations between GMF intensity and CVDs events, multi-factor panel data analysis was also generated to adjust the effect of confounding factors.

Results

For the average data during 1996–2019, linear regression model revealed consistent positive correlations between total GMF (tGMF) intensity and mortality of total CVDs [coef = 0.009, (0.006,0.011 95%CI)], whereas negative correlations were found between horizonal GMF (hGMF) intensity and total CVD mortality [coef = -0.010 (-0.013, -0.007 95%CI)]. When considering the time trend, panel data analysis still demonstrated positive correlation between tGMF and total CVDs mortality [coef = 0.009, (0.008,0.009 95%CI)]. Concurrently, the hGMF negatively correlated with total CVDs mortality [coef = -0.008, (-0.009, -0.007 95%CI)]. When the panel models were adjusted for confounding factors, no reverse of correlation tendency was found between tGMF, hGMF and CVDs events. In high-income territories, positive correlation was found between geomagnetic storm (GMS) frequency and mortality of total CVDs [coef = 14.007,(2.785, 25.229 95%CI)], however, this positive trend faded away gradually with the latitude decreasing from polar to equator.

Conclusions

Stable and long-term horizontal component of GMF may be beneficial to cardiac health. Unstable and short-term GMF called GMD could be a hazard to cardiac health. Our results suggest the importance of regular GMF in maintaining cardio-health state and the adverse impacts of GMD on cardiac health.

Peer Review reports

Introduction

The Earth's magnetic field, also referred to as the geomagnetic field (GMF), is composed of two components: the Internal Field and the External Field [1]. Internal Field or Main Field is generated by the geomagnetic dynamo process in the Earth’s core, facilitateing the transportation of solar wind particles from the diurnal to nocturnal side of our planet and enabling further spatial redistribution. The Earth’s internal magnetic field is a vector field, and the total GMF intensity (tGMF, F) can be divided into two components: horizontal intensity (hGMF, H) and vertical intensity (vGMF, Z) (Fig. 1). The horizontal components can be described by orthogonal components X (directed towards geographic North), Y (directed towards geographic East). These five components can be interconnected through the following equation:

Fig. 1
figure 1

Components of geomagnetic field vector. *Geomagnetic field is described by five parameters, including horizontal intensity (H), vertical intensity (Z), total intensity (F) and the north (X) and east (Y) components of the horizontal intensity. X is the component along the Geographic North direction; Y is the component along the Geographic East direction; Z is the vertical component pointing vertically downwards; H is the magnitude of the horizontal component; F is total geomagnetic field intensity

$$F=\sqrt{{H}^{2}+{Z}^{2}}$$

Where H is given by

$$H=\sqrt{{X}^{2}+{Y}^{2}}$$

The average total intensity above Earth’s surface ranges from 20,000nT to 65,000nT, depending on latitude and altitude. Considering the vector of the GMF, the horizontal field (H) displays opposite trends compared to the total intensity, being greatest in the equatorial region and weakening towards the Earth's poles [2] (Fig. 2). The external field is highly dynamic and arises from interactions between the solar wind and Earth's internal magnetic field. Some of these interactions can lead to the generation of short-term geomagnetic disturbance (GMD), also referred to as geomagnetic storm (GMS) [3, 4]. These magnetic field disturbances are observable across the Earth's surface and can induce currents, referred to as geomagnetically induced currents (GICs) [5]. Prior research has established that short-term GMDs could have adverse effects on human health, including an increased suicide rate, higher hospitalization rates for acute coronary syndrome, and an elevated risk of stroke occurrences [6,7,8,9]. Observations have indicated a strong correlation between solar events and variations in the horizontal geomagnetic component. This correlation suggests that the horizontal geomagnetic field GMF may play a stabilizing role in relation to the charges generated by geomagnetic disturbances GMDs [10,11,12,13]. Hence, exploring the connection between environmental physical activities and diverse medical events, as well as understanding how these activities are implicated in specific diseases, can prove valuable for certain aspects of preventive measures. Cardiovascular diseases (CVDs), including conditions like ischemic heart disease (IHD) and stroke, constitute a primary cause of global mortality and a significant contributor to disability [14]. With an electrical conduction system and cardiac muscle cells capable of transmitting electrical stimuli, the cardiovascular system is considered susceptible to the influence of external electromagnetic fields [15,16,17,18]. A growing body of research has explored the link between solar and geomagnetic activity (GA) and CVD mortality and morbidity. Regarding short-term exposure to the fluctuating GMF, as mentioned earlier, GMDs have been associated with increased occurrences of total deaths, CVDs, myocardial infarction (MI), cardiac arrhythmia, and stroke events [6, 7, 19,20,21,22,23]. Importantly, certain ground-based investigations have revealed that a hypo-magnetic field (HMF), simulating the low-magnetic environment of outer space, can disrupt circadian rhythms and contribute to the development of mental and physiological disorders [24,25,26]. Hence, it is plausible to infer that elevated levels of consistent GMF)activity might have beneficial impacts on cardiac health by safeguarding organisms against irregular and transient solar activities.

Fig. 2
figure 2

Map of main field total intensity and horizontal intensity from World Magnetic Map (WMM) model for 2020.0. From https://www.ncei.noaa.gov/products/world-magnetic-model

CVDs exhibited a geographic pattern, characterized by a greater overall incidence of conditions like MI and hypertension in regions situated at higher latitudes [27,28,29]. While these geographic patterns could be partly attributed to factors such as temperature, ultraviolet exposure, and genetic variations, it remains uncertain whether the intensity of the GMF plays a role as a determinant of CVDs. Furthermore, as previously noted, the horizontal intensity of the GMF is highest at the equatorial region and gradually diminishes towards the Earth's poles.. Consequently, in this original investigation, we aimed to explore the link between stable GA, particularly the horizontal component of the GMF, and CVD risks across 204 countries with varying latitudes. Additionally, we examined the impact of unstable geomagnetic activity, specifically GMDs from solar events, on the occurrence of CVD events.

Method

Study design and index definition

As the GMF undergoes periodic fluctuations aligned with the solar cycle (approximately every 12 years), we selected an entire solar magnetic period (1996–2019, encompassing solar cycles 23 and 24) as the timeframe for our investigation [30]. Data of GMF and GMD, Cardiovascular events, meteorological factors, and economic indices from 204 countries or territories globally were gathered from authoritative databases. Both the average values for the entire 24-year period and annual data for each indicator were compiled, enabling separate cross-sectional and panel data analyses.

The countries or territories included in our study consist of 193 sovereign states and 2 quasi-sovereign states, which have been admitted by the United Nations and are members of the World Health Organization (WHO), along with 9 significant economies as defined by the World Bank (https://www.worldbank.org/en/home). Intensity, direction and fluctuation of GMF were all considered in our study, in order to thoroughly investigate the association between GMF and CVD. Apart from total intensity, vertical and horizonal intensity, which were defined as positive in the upward and southward directions respectively, were also taken into account and measured in nanotesla (nT) [31]. Due to the orthogonality relation among the intensities of tGMF, vGMF, and hGMF, our study exclusively analyzed tGMF and hGMF. The frequency of annual strong magnetic disturbances (defined as level G3-G5 by the National Oceanic and Atmospheric Administration, NOAA) was collected for the purpose of fluctuation analysis. This information can be accessed from the NOAA's official website: https://www.swpc.noaa.gov/noaa-scales-explanation [32]. The entire range of CVD, IHD, and stroke were selected as the subjects based on the International Statistical Classification of Diseases and Related Health Problems 10th Revision (ICD-10 2019).

Indicators related to weather and economy were introduced as covariates to mitigate bias. In terms of weather indicators, the maximum and minimum temperatures were measured in degrees Fahrenheit, and the occurrence of specific weather events was quantified using the fractions associated with six distinct conditions (FRSHTT: Fog/Rain/Snow/Hail/Thunder/Typhoon). All countries and territories were categorized based on their gross national income (GNI) per capita into higher income, upper middle income, lower middle income, and low income groups (The World Bank 2020) [33].

Data source

Data of GMF were obtained based on the 13th generation of the International Geomagnetic Reference Field (IGRF, http://www.geomag.bgs.ac.uk/research/modelling/IGRF.html). For each country or territory, the capital or political center was standardized as the location where the point GMF intensity was recorded, serving as a representation of the general GMF in that region. The specific latitude and longitude information for each location was also obtained from the IGRF model. Data for each individual day spanning from 1996 to 2019 was recorded sequentially and in batches using Node-JavaScript-Runtime, which operated on the server side and facilitated high concurrency. Subsequently, the total and annual GMF intensities were computed for the 204 regions encompassing the 24-year period. Further, the annual frequency of strong GMD was obtained from Space Weather Data Portal (https://lasp.colorado.edu/space-weather-portal/data), employing a filter criteria of G index ≥ 3. Data regarding the FRSHTT ratio, maximum and minimum temperatures were extracted from the National Oceanic and Atmospheric Administration (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/). The daily weather particulars for each region were aggregated into both total and annual data using the same methodology as applied to GMF intensity. The global burden of diseases (GBD) results tool (https://ghdx.healthdata.org/gbd-results-tool) was utilized to access data related to CVDs. This tool serves as a user-friendly platform founded on GBD 2019, offering annual data encompassing a diverse array of measures for all GBD causes, risks, impairments and injuries, across nearly every region globally. The core dataset utilized in this study encompasses 204 countries/territories observed from 1996 to 2019. The selected causes of interest include stroke, IHD, and the total CVDs. Measures of cause include mortality, incidence, and prevalence. For age standardization, the ratio per million people was chosen as the age and metric reference. For each sample, its income level was corresponded according to the GNI. All data underwent independent searching and processing by two researchers, followed by a cross-check to ensure accuracy.

Statistical analysis

Three main approaches were adopted to ascertain and estimate the effects of GMF on the outcomes. Initially, we present the correlation outcomes yielded by the ordinary least squares (OLS) method, as displayed in Table 1. While it is improbable that our independent variables are influenced by dependent variables, potential unobserved factors that are simultaneously associated with both GMF intensity and CVD events may still introduce disturbances to the results. To address this issue of omitted variables, we utilized both parametric and non-parametric methods to reevaluate our theoretical assumption. The parametric method involved a country fixed-effect model. The merit of this approach lies in its ability to exclude consistent unobserved factors, such as endowments and environmental factors, at the country level through the incorporation of a set of dummy variables. Additionally, mean value imputation, cluster robust standard error and time dummies were utilized to consummate the fixed-effect model. Alternative measurements of GMF (global annual GMF intensity, average GMF intensity of 204 locations) and CVD (total CVDs, IHD and stroke’s mortality, incidence and prevalence) were applied to improve robustness of the estimator in all three approaches mentioned above.

Table 1 The correlation results between total intensity (A), horizonal intensity (B) and difference measures of CVD, IHD and stroke in average, single factor panel and multiple-factor panel data analysis

To investigate the impact of GMD, linear regression was applied to explore the relationship between global GMS frequency and CVD mortality. Furthermore, our panel analysis exclusively encompassed high-income countries and territories. This choice was made to mitigate potential reporting biases stemming from backward development. Panel data was structured based on location and year. Additionally, we conducted stratified analysis according to latitude (0–15° near-equatorial region, 15–55° mid-latitude region, 55–90° near-polar region) was also conducted.

All the data analysis was processed using STATA (StataSE 16 64-bit), and statistical significance was defined as P < 0.05.

Result

For the cross-sectional study, annual data on GMF intensity and CVD events from 204 included territories, as well as global GMS data, were fully recorded without any gaps occurring during processing. In the panel data analysis, each indicator had 4896 (24*204) data points. Data on CVD events and GMF were complete, while 958 data points were missing for minimum and maximum air temperatures, and 1267 data points were missing for special weather events. All these gaps were filled by using the average values. Additionally, the GNI level was unavailable for 6 territories due to lack of recognized sovereignty.

Correlation between GMF intensity and CVDs events

In the cross-sectional data analysis using a linear regression model, a significant relationship was observed between GMF intensity and CVD events in terms of mean values. For total CVDs, evident positive correlations were found between tGMF intensity and total CVDs mortality [coef = 0.009, (0.006,0.011 95%CI)], incidence [coef = 0.022, (0.017,0.026 95%CI)] and prevalence [coef = 0.160 (0.122,0.198 95%CI)]. On the contrary, a contrasting pattern emerged with hGMF intensity, showing notable negative correlations with CVD mortality [coef = -0.010 (-0.013, -0.007 95%CI)], incidence [coef = -0.024, (-0.029,-0.018 95%CI)] and prevalence [coef = -0.198, (-0.243,-0.152 95%CI)]. For both IHD and stroke, consistent and conspicuous positive correlations were discovered between tGMF intensity and the rates of mortality, incidence, and prevalence. Conversely, negative correlations were identified between hGMF intensity and the rates of mortality, incidence, and prevalence. The relationships between the mortality rates of CVD, IHD, and stroke, and the intensities of tGMF and hGMF, were effectively illustrated using scatter plots and fitted lines in Fig. 3

Fig. 3
figure 3

Correlations between mortality rate (per 100000) and geomagnetic field intensity (nT) with total (A) and horizonal (B) component in CVD, IHD and S.*CVD = cardiovascular diseases, IHD = ischemic heart diseases, S = stroke

The panel data was configured by year and GNI level. In the single-factor panel regression, a positive correlation was observed between tGMF and the mortality [coef = 0.009, (0.008,0.009 95%CI)], incidence [coef = 0.019, (0.018,0.020 95%CI)] and prevalence [coef = 0.122 (0.112,0.131 95%CI)] of total CVDs. On the contrary, hGMF exhibited a negative correlation with total CVDs mortality [coef = -0.008, (-0.009, -0.007 95%CI)], incidence [coef = -0.018, (-0.020, -0.017 95%CI)] and prevalence [coef = -0.142, (-0.153, -0.132 95%CI)]. The correlation trends observed between GMF intensity (tGMF and hGMF) and IHD and stroke events were consistent with the correlation trends seen between GMF intensity and total CVDs. The concurrent association was also evident in multi-factor panel analysis. After accounting for confounding factors, tGMF continued to be identified as a risk factor for total CVDs, IHD and stroke, across mortality, incidence, and prevalence measures. Conversely, hGMF was indicated as a protective factor in all total CVDs events, as well as in the incidence [coef = -0.002, (-0.003, -0.001 95%CI)] and prevalence [coef = -0.010, (-0.016, -0.005 95%CI)] of IHD. Nonetheless, no statistically significant correlation was observed between hGMF and IHD mortality [coef = -0.0004, (-0.0008, 0.0001 95%CI)], stroke mortality [coef = -0.0001, (-0.0003, 0.0002 95%CI)] and prevalence [coef = 0.001, (-0.001,0.002 95%CI)]. No transformation in significance and tendency occurred after the adjustment of time dummies, the comprehensive results of correlations between GMF intensity and CVDs events were illustrated in Table 1.

Correlation between GMD frequency and CVDs events

Between the years 1996 and 2019, there was no noteworthy linear association detected between global GMD frequency and global CVDs mortality [coef = -0.226 (-0.478,0.170 (95%CI)] However, distinct positive and negative correlations were observed between GMD frequency and stroke mortality [coef = 0.339, (0.185,0.493 95%CI)] and IHD mortality [coef = -0.338, (-0.596, -0.080 95%CI)] respectively.

In our panel analysis, a total of 63 high-income countries or territories were encompassed within our study. Among these entities, 9 were allocated to the near-equatorial region group, 46 were classified under the mid-latitude group, and 8 were segregated into the near-polar region group.. Across all the included territories, GMD exhibited a significant association with elevated mortality rates for CVDs [coef = 14.007,(2.785, 25.229 95%CI)], IHD [coef = 13.688, (6.040, 21.335 95%CI)] and stroke [coef = 8.891, (4.928, 12.855 95%CI)]. For stratification analysis, significant correlations were observed in near-polar territories between GMD frequency and mortality rates of IHD [coef = 39.703, (7.315, 72.092 95%CI)] and stroke [coef = 19.461, (2.314, 36.607 95%CI)]. Comparable findings were obtained for mid-latitude regions, where the presence of GMS was associated with increased mortality rates for total CVDs ([coef = 15.452, (3.131, 27.775 95%CI)], IHD [coef = 14.556, (6.998, 22.114 95%CI)] and stroke [coef = 9.212, (4.853, 13.572 95%CI)]. However, with decreasing latitude and the concurrent decline in GMD frequency alongside the gradual increase in horizontal GMF, the trend of positive correlation between GMD and CVDs events diminished in near-equatorial regions. No significant relationship was observed between GMD frequency and mortality for any type of CVD. This suggests that GMD and horizontal GMF might respectively function as a risk factor and a protective factor for CVD events from an alternative perspective (Table 2).

Table 2 Correlations between geomagnetic storm frequency in developed territories with different latitude and difference measures of CVD, IHD and stroke

Discussion

In this international long-term analysis, we studied how tGMF, hGMF, and GMD) frequency affect CVDs. By including countries at different latitudes, we found evidence supporting hGMF's protective role against CVDs. As latitude decreases, GMD risk also reduces. Notably, our study is the first to explore the impact of various Geomagnetic Field components, including vector facets (horizontal and vertical) and contribution facets (internal and external), on CVD events.

Overall, our study demonstrated that both stable and unstable GMF are correlated to CVDs events. As for total intensity of stable GMF, our cross-sectional data showed clear positive correlation between tGMF and CVDs events, and this association remained similar and statistically significant in panel data analysis, which took GNI level and temporal components into consideration. To our knowledge, this study is the first to evaluate the impact of the Earth’s main field intensity on CVD. We used International Geomagnetic Reference Field (IGFR) Calculator to get large-scale, time-varying portion of Earth’s internal magnetic field from 1996 to 2019. Although our results about total GMF are similar to previous epidemiological studies which showed the positive relationship between total GA and CVDs, the meaning of geomagnetic indexes in those studies is totally different from us [19, 20, 34, 35]. In those studies, they used Kp or Ap index, which represents the disturbance of Earth’s magnetic field whereas our study was based on general internal magnetic field in different areas [19,20,21, 34, 36,37,38,39]. GMF protects Earth from solar wind and cosmic rays, stabilizing the ionosphere where Earth’s atmosphere meets the space [9]. However, in our study, total GMF positively related to CVDs events whereas there is an inverse correlation between horizontal GMF (hGMF) and CVDs. It can be due to the deflection of charged particles influenced by different component of GMF [40]. The horizontal GMF produces vertical electric field and exerts a vertical force on atmospheric charged particles, which impedes further separation of the charges. On the contrary, the vertical GMF produces a horizontal electric field and causes charges separated. Our findings about horizontal GMF are consistent with previous study showing the negative relationship between horizontal GMF and multiple sclerosis prevalence [41]. In the present study, we also evaluated the impact of irregular natural GMD. Defined as G-scale ≥ 3 (based on Kp index), GMD showed positive correlations with CVDs in all latitude regions, which are consistent with previous epidemiological studies. Those studies have demonstrated that intense geomagnetic activities (GMA) or GMD increases all-cause, CVD and MI deaths in specific regions, such as Mexico, Russia, Israel and U.S. [20, 34,35,36, 42]. It is noteworthy that, we included 63 developed countries at multiple latitudes and found those positive correlations diminished as latitude drops and no significant association was observed between GMD frequency and any type of CVDs in subequatorial regions. As shown in Fig. 2, horizontal GMF is greatest around equators and weakest towards poles while total and vertical GMF display opposite arrangement. It is consistent with our main findings that horizontal GMF negatively correlated with CVDs events and supports the idea that horizontal GMF may decrease the risk of adverse CVDs events in intense GMD or GMS.

Understanding the mechanisms by which stable GMF and unstable GMD influence the cardiovascular system is necessary for explaining the results in the present study. The GMF protects Earth’s atmosphere against charged particles from the solar wind, wthereby reducing ionization caused by cosmic rays [43]. Moreover, the notable correlation between horizontal or overall GMF and CVDs, while not robust enough to imply a causal relationship, does necessitate the examination of a potential third variable that could independently influence both factors. GMS is a major disturbance of Earth’s magnetosphere caused by energy exchange between the solar wind and Earth's surrounding atmosphere, resulting in changes in currents, plasmas, and fields within the magnetosphere (NOAA 2022). The perturbation of environmental electric and magnetic fields has the potential to influence the autonomic nervous system (ANS) and circadian rhythm. Several limited-scale cohort studies have showcased a connection between short-term geomagnetic disturbances (GMD) and reduced heart rate variability (HRV). This reduction in HRV serves as an indicator of ANS dysregulation, standing as an established marker for heightened cardiovascular risk [44,45,46]. Recently, a large cohort of elderly men over a 17-year follow-up demonstrated a significant adverse impact of GA and intense GMD on HRV and the association remained after adjusting for air pollution [19]. Reduced HRV could be the one of the biological mechanisms between GMD and CVDs in our study, however, it remains unclear as to how the body detects GMA variation and GMD. The photo/magnetic-reception system, including cryptochrome protein in retina and ferromagnetic receptors in magnetite-contained cell, is proposed as a mediator enabling organisms to sense variations in the environmental GMF related to GMD [19, 24]. Overstimulating this system disrupts the balance between sympathetic and parasympathetic activities in ANS, which further affects HRV. Furthermore, this system mediates light-dependent stimulation of GMD and influences the melatonin secretion in the pineal gland [47, 48]. Functioning as both an antioxidant and a circadian rhythm regulator, melatonin assumes a vital role in CVDs, encompassing conditions such as MI, ischemic strokes, and arrhythmias [49]. Hence, the mechanisms through which GMD heighten the cardiovascular risk of events entail the dysregulation of the ANS and the reduction of melatonin levels (Fig. 4).

Fig. 4
figure 4

The mechanistic basis of Internal Geomagnetic field and External Geomagnetic disturbance (GMD). *Geomagnetic field protects Earth from ionizing radiation coming from Solar Wind which causes disturbance of electric and magnetic field. Short-term geomagnetic disturbance modulates autonomic nervous system (ANS) and subsequently decreases heart rate variability (HRV). GMD overstimulates photo/magnetic-reception system in intrinsically photosensitive retinal ganglion cells (ipRGCs) and impacts circadian rhythm and melatonin secretion in pineal gland. GMD can alter the toxicity of air pollution which indirectly enhance the association between air pollutants and CVDs events

Another potential mechanism for interaction between GMD and CVDs is alterations in the toxicity of air pollution. During periods of GMD, charged particles caused by solar wind can create air ionization and produce electric currents that alter the chemistry of ambient aerosols and droplets [50]. While the precise impact of GMD on altering the electrochemical characteristics and toxicity of ambient particles requires further elucidation, numerous investigations have noted that heightened GMA and pronounced GMD could potentially amplify the connection between air pollutants (such as PM2.5, NO2, black carbon (BC), and ultrafine particles (PN)) and cardiovascular events [7, 23, 50, 51]. Under protection of GMF, these reactions take place on a much smaller scale, which is supportive of our results about the beneficial effects of hGMF on CVDs [9].

There are some limitations in our study. Firstly, the utilization of point data for weather and GMF to represent a country or territory's overall data may introduce biases, particularly in countries with large areas or latitude spans. Secondly, despite the high accuracy verification of the IGRF-13 model (Alken, Thebault et al. 2021), deviations between the model and actual GMF strength do exist. Moreover, the G-index, which is utilized as the sole metric for measuring geomagnetic storms in our study, may introduce potential biases in subsequent analyses due to its inherent limitations and lack of precision. Besides, similar to other studies of this type, the null hypothesis is not independent of the data used for testing. Additionally, it may be difficult to account for all confounding factors and missing data can introduce measurement bias in multi-factor panel analysis. Lastly, although we attested the significant correlations between GMF intensity and CVDs events, the data were not sufficient to demonstrate dose–effect relationship, this is also partly due to inherent limitations in ecological research, such as ‘lack of data on the joint distribution of exposure and disease’, ‘inability to control for potential confounding factors’, and ‘approximate exposure levels in the data’. Future prospective studies addressing such limitations would be needed to substantiate the causal relationship between GMF and CVDs events.

Conclusion

In conclusion, our study revealed a positive correlation between elevated levels of total GMF intensity and CVDs, alongside a negative correlation between the horizontal component of GMF and CVDs. Furthermore, we identified a positive correlation between unstable GMF activity—GMD—and CVD events. This hypothesis-generating investigation underscores the significance of consistent GMF exposure for maintaining cardiovascular well-being, while highlighting the detrimental effects of GMD on cardiac health. Further prospective studies are imperative to establish a causal relationship between GMF/GMD and CVD events.

Availability of data and materials

Data of GMF were obtained based on the 13th generation of the International Geomagnetic Reference Field (IGRF, http://www.geomag.bgs.ac.uk/research/modelling/IGRF.html). The annual frequency of strong GMD was obtained from Space Weather Data Portal (https://lasp.colorado.edu/space-weather-portal/data). Data regarding the FRSHTT ratio, maximum and minimum temperatures were extracted from the National Oceanic and Atmospheric Administration (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/). The GBD results tool (https://ghdx.healthdata.org/gbd-results-tool) was utilized to access data related to CVDs. All the included territories’ GNI levels were obtained from the World Bank 2020 (https://www.worldbank.org/en/home).

Abbreviations

ANS:

Autonomic nervous system

BC:

Black carbon

CVDs:

Cardiovascular diseases

FRSHTT:

Fog/rain/snow/hail/thunder/typhoon

GA:

Geomagnetic activity

GBD:

Global burden of disease

GICs:

Geomagnetically induced currents

GMD:

Geomagnetic disturbance

GMF:

Geomagnetic field

tGMF:

Total geomagnetic field

hGMF:

Horizonal geomagnetic field

vGMF:

Vertical geomagnetic field

GMS:

Geomagnetic storm

GNI:

Gross national income

HMF:

Hypo-magnetic field

HRV:

Heart rate variability

ICD-10:

International statistical classification of diseases and related health problems 10th-revision

IGFR:

International geomagnetic reference field

IHD:

Ischemic heart disease

MI:

Myocardial infarction

NOAA:

National oceanic and atmospheric administration

OLS:

Ordinary least squares

PN:

Ultrafine particles

WHO:

World health organization

References

  1. Glassmeier K-H, Vogt J. Magnetic Polarity Transitions and Biospheric Effects. Space Sci Rev. 2010;155(1):387–410.

  2. The Earth's main field. In: Introduction to Geomagnetic Fields. 2 edn. Edited by Campbell WH. Cambridge: Cambridge University Press; 2003: 1–66.

  3. Encyclopedia of Geomagnetism and Paleomagnetism: Springer Netherlands; 2007.

  4. Roberts PH, King EM. On the genesis of the Earth’s magnetism. Rep Prog Phys. 2013;76(9):096801.

    Article  PubMed  Google Scholar 

  5. Fadhil A. A comparison study between the internal and external geomagnetic Field at the location (3.5oN, 96oE). Diyala J Human Res. 2009;32:1.

    Google Scholar 

  6. Feigin VL, Parmar PG, Barker-Collo S, Bennett DA, Anderson CS, Thrift AG, Stegmayr B, Rothwell PM, Giroud M, Bejot Y, et al. Geomagnetic storms can trigger stroke: evidence from 6 large population-based studies in Europe and Australasia. Stroke. 2014;45(6):1639–45.

    Article  PubMed  Google Scholar 

  7. Vencloviene J, Grazuleviciene R, Babarskiene R, Dedele A, Grazulevicius T. Short-term nitrogen dioxide exposure and geomagnetic activity interaction: contribution to emergency hospitalization for acute coronary syndrome. Int J Environ Health Res. 2011;21(3):149–60.

    Article  CAS  PubMed  Google Scholar 

  8. Nishimura T, Tada H, Nakatani E, Matsuda K, Teramukai S, Fukushima M. Stronger geomagnetic fields may be a risk factor of male suicides. Psychiatry Clin Neurosci. 2014;68(6):404–9.

    Article  PubMed  Google Scholar 

  9. Erdmann W, Kmita H, Kosicki JZ, Kaczmarek Ł. How the geomagnetic field influences life on earth - an integrated approach to geomagnetobiology. Orig Life Evol Biosph. 2021;51(3):231–57.

    Article  PubMed  Google Scholar 

  10. Syed Zafar S, Sabri N, Jusoh M, Dagang AN, Adzni M, Umar R. Behaviour of Geomagnetic storm, horizontal geomagnetic field and solar wind parameters during solar flare and CMEs event. J Phys. 2021;1768:012003.

    Google Scholar 

  11. Rastogi RG. Magnetic storm effects in components of the geomagnetic field at low and middle latitudes. J Atmospher Solar-Terrestrial Phys. 2005;67(7):665–75.

    Article  Google Scholar 

  12. Roslan U: The Behavior Of Solar Wind Parameters And Geomagnetic Activity Indices For Geomagnetically Induced Current Observations; 2018.

  13. Agbo G, Azi A, Nnabo P, Onuorah L, Abbey T: Association Of Residual Horizontal Geomagnetic Field With Interplanetary Magnetic Field Parameters. 2020:15.

  14. Roth GA, Mensah GA, Johnson CO, Addolorato G, Ammirati E, Baddour LM, Barengo NC, Beaton AZ, Benjamin EJ, Benziger CP, et al. Global burden of cardiovascular diseases and risk factors, 1990–2019: update from the GBD 2019 Study. J Am Coll Cardiol. 2020;76(25):2982–3021.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Wang S, Zhou X, Huang B, Wang Z, Zhou L, Wang M, Yu L, Jiang H. Noninvasive low-frequency electromagnetic stimulation of the left stellate ganglion reduces myocardial infarction-induced ventricular arrhythmia. Sci Rep. 2016;6:30783.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wu F, Wang C, Xu Y, Ma J. Model of electrical activity in cardiac tissue under electromagnetic induction. Sci Rep. 2016;6(1):28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Jeong JH, Kim JS, Lee BC, Min YS, Kim DS, Ryu JS, Soh KS, Seo KM, Sohn UD. Influence of exposure to electromagnetic field on the cardiovascular system. Auton Autacoid Pharmacol. 2005;25(1):17–23.

    Article  CAS  PubMed  Google Scholar 

  18. Soltani D, Samimi S, Vasheghani-Farahani A, Shariatpanahi SP, Abdolmaleki P, Madjid Ansari A. Electromagnetic field therapy in cardiovascular diseases: a review of patents, clinically effective devices, and mechanism of therapeutic effects. Trends Cardiovasc Med. 2021;33(2):72–8.

    Article  PubMed  Google Scholar 

  19. Zilli Vieira CL, Chen K, Garshick E, Liu M, Vokonas P, Ljungman P, Schwartz J, Koutrakis P. Geomagnetic disturbances reduce heart rate variability in the normative aging study. Sci Total Environ. 2022;839:156235.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Dimitrova S, Babayev E, Georgieva K, Obridko V, Mustafa F. Possible effects of solar and geomagnetic activity on sudden cardiac death in middle latitudes. Sun and Geosphere. 2009;4:84–8.

    Google Scholar 

  21. Samsonov SN, Manykina VI, Kleimenova NG, Parshina SS, Strekalovskaya AA, Petrova PG. The HELIO-geophysical storminess health effects in the cardio-vascular system of a human in the middle and high latitudes. Wiad Lek. 2016;69(3 pt 2):537–41.

    PubMed  Google Scholar 

  22. Cornelissen Guillaume G, Gubin D, Beaty LA, Otsuka K. Some near- and far-environmental effects on human health and disease with a focus on the cardiovascular system. Int J Env Res Public Health. 2020;17(9):3083.

    Article  Google Scholar 

  23. Wang VA, Zilli Vieira CL, Garshick E, Schwartz JD, Garshick MS, Vokonas P, Koutrakis P. Solar activity is associated with diastolic and systolic blood pressure in elderly adults. J Am Heart Assoc. 2021;10(21):e021006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Xue X, Ali YF, Luo W, Liu C, Zhou G, Liu NA. Biological effects of space hypomagnetic environment on Circadian Rhythm. Front Physiol. 2021;12:643943.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Stoupel E, Tamoshiunas A, Radishauskas R, Bernotiene G, Abramson E, Israelevich P. Acute myocardial infarction (AMI) (n-11026) on days of zero geomagnetic activity (GMA) and the following week: differences at months of maximal and minimal solar activity (SA) in solar cycles 23 and 24. J Basic Clin Physiol Pharmacol. 2012;23(1):5–9.

    Article  PubMed  Google Scholar 

  26. Zhang HT, Zhang ZJ, Mo WC, Hu PD, Ding HM, Liu Y, Hua Q, He RQ. Shielding of the geomagnetic field reduces hydrogen peroxide production in human neuroblastoma cell and inhibits the activity of CuZn superoxide dismutase. Protein Cell. 2017;8(7):527–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Bhatnagar A. Environmental determinants of cardiovascular disease. Circul Res. 2017;121(2):162–80.

    Article  CAS  Google Scholar 

  28. Cabrera SE, Mindell JS, Toledo M, Alvo M, Ferro CJ. Associations of blood pressure with geographical latitude, solar radiation, and ambient temperature: results from the chilean health survey, 2009–2010. Am J Epidemiol. 2016;183(11):1071–3.

    Article  PubMed  Google Scholar 

  29. Mackay DF, Clemens TL, Hastie CE, Cherrie MPC, Dibben C, Pell JP. UVA and seasonal patterning of 56 370 myocardial infarctions across Scotland, 2000–2011. J Am Heart Assoc. 2019;8(23):e012551.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Kane R. Some implications using the group sunspot number reconstruction. Sol Phys. 2002;205(2):383–401.

    Article  Google Scholar 

  31. Alken P, Thébault E, Beggan CD, Amit H, Aubert J, Baerenzung J, Bondar T, Brown W, Califf S, Chambodut A. International geomagnetic reference field: the thirteenth generation. Earth Planets Space. 2021;73(1):1–25.

    Google Scholar 

  32. Gonzalez W, Joselyn J-A, Kamide Y, Kroehl HW, Rostoker G, Tsurutani B, Vasyliunas V. What is a geomagnetic storm? J Geophys Res. 1994;99(A4):5771–92.

    Article  Google Scholar 

  33. Ram R. PPP GDP per capita for countries of the world: a comparison of the new ICP results with World Bank data. Soc Indic Res. 2016;127(3):1057–66.

    Article  Google Scholar 

  34. Vencloviene J, Babarskiene R, Milvidaite I, Kubilius R, Stasionyte J. The effect of solar-geomagnetic activity during and after admission on survival in patients with acute coronary syndromes. Int J Biometeorol. 2014;58(6):1295–303.

    Article  PubMed  Google Scholar 

  35. Stoupel E. Cardiac arrhythmia and geomagnetic activity. Indian Pacing Electrophysiol J. 2006;6(1):49–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Zilli Vieira CL, Alvares D, Blomberg A, Schwartz J, Coull B, Huang S, Koutrakis P. Geomagnetic disturbances driven by solar activity enhance total and cardiovascular mortality risk in 263 U.S. cities. Environ Health. 2019;18(1):83.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Villoresi G, Ptitsyna NG, Tiasto MI, Iucci N. Myocardial infarct and geomagnetic disturbances: analysis of data on morbidity and mortality. Biofizika. 1998;43(4):623–31.

    CAS  PubMed  Google Scholar 

  38. Podolská K. The impact of ionospheric and geomagnetic changes on mortality from diseases of the circulatory system. J Stroke Cerebrovasc Dis. 2018;27(2):404–17.

    Article  PubMed  Google Scholar 

  39. Vencloviene J, Babarskiene R, Slapikas R, Sakalyte G. The association between phenomena on the sun, geomagnetic activity, meteorological variables, and cardiovascular characteristic of patients with myocardial infarction. Int J Biometeorol. 2013;57(5):797–804.

    Article  PubMed  Google Scholar 

  40. . In: Britannica, T Information Architects of Encyclopaedia "geomagnetic field" Encyclopedia Britannica, July 17, 2022 https://www.britannica.com/facts/geomagnetic-field.

  41. Wade B. Correlations of the earth′s magnetic field intensity with global prevalence of multiple sclerosis. American J Life Sci. 2013;1:31.

    Google Scholar 

  42. Mendoza B, Diaz-Sandoval R. Effects of solar activity on myocardial infarction deaths in low geomagnetic latitud regions. Nat Hazards. 2004;32:25–36.

    Article  Google Scholar 

  43. Aplin KL, Harrison RG, Bennett AJ. Effect of the troposphere on surface neutron counter measurements. Adv Space Res. 2005;35(8):1484–91.

    Article  CAS  Google Scholar 

  44. Costa MD, Redline S, Davis RB, Heckbert SR, Soliman EZ, Goldberger AL. Heart rate fragmentation as a novel biomarker of adverse cardiovascular events: the multi-ethnic study of atherosclerosis. Front Physiol. 2018;9:1117.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Berger M, Pichot V, Solelhac G, Marques-Vidal P, Haba-Rubio J, Vollenweider P, Waeber G, Preisig M, Barthélémy JC, Roche F, et al. Association between nocturnal heart rate variability and incident cardiovascular disease events: The HypnoLaus population-based study. Heart Rhythm. 2022;19(4):632–9.

    Article  PubMed  Google Scholar 

  46. Alabdulgader A, McCraty R, Atkinson M, Dobyns Y, Vainoras A, Ragulskis M, Stolc V. Long-term study of heart rate variability responses to changes in the solar and geomagnetic environment. Sci Rep. 2018;8(1):2663.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Cornélissen G, Halberg F, Sothern RB, Hillman DC, Siegelová J. Blood pressure, heart rate and melatonin cycles synchronization with the season, earth magnetism and solar flares. Scr Med (Brno). 2010;83(1):16–32.

    PubMed  Google Scholar 

  48. Burch JB, Reif JS, Yost MG. Geomagnetic disturbances are associated with reduced nocturnal excretion of a melatonin metabolite in humans. Neurosci Lett. 1999;266(3):209–12.

    Article  CAS  PubMed  Google Scholar 

  49. Dominguez-Rodriguez A, Abreu-Gonzalez P, Sanchez-Sanchez JJ, Kaski JC, Reiter RJ. Melatonin and circadian biology in human cardiovascular disease. J Pineal Res. 2010;49(1):14–22.

    CAS  PubMed  Google Scholar 

  50. Zilli Vieira CL, Koutrakis P. The impact of solar activity on ambient ultrafine particle concentrations: an analysis based on 19-year measurements in Boston, USA. Environ Res. 2021;201:111532.

    Article  CAS  PubMed  Google Scholar 

  51. Zilli Vieira CL, Link MS, Garshick E, Peralta AA, Luttmann-Gibson H, Laden F, Liu M, Gold DR, Koutrakis P. Solar and geomagnetic activity enhance the effects of air pollutants on atrial fibrillation. EP Europace. 2022;24(5):713–20.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Natural Science Foundation of China (81970325, and 82170375); Key Research and Development Project of Science & Technology Department of Sichuan Province (2022ZDZX0020); Chinese Medical Association Cardiovascular Branch (CSC) Clinical Research Special Fund Project (CSCF2020B04); “1· 3· 5 project for disciplines of excellence– Clinical Research Fund, West China Hospital, Sichuan University”.

Disclosures

There are no financial conflicts of interest to disclose.

Funding

This work was supported by the National Natural Science Foundation of China (81970325, and 82170375); Key Research and Development Project of Science & Technology Department of Sichuan Province (2022ZDZX0020); Chinese Medical Association Cardiovascular Branch (CSC) Clinical Research Special Fund Project (CSCF2020B04); “1· 3· 5 project for disciplines of excellence– Clinical Research Fund, West China Hospital, Sichuan University”.

Author information

Authors and Affiliations

Authors

Contributions

Professor Mao Chen set the original thought and topic about this study. Zheng Chai and Yan Wang wrote the main manuscript and prepared the figures and tables. Zheng Chai and Yi-Ming Li conducted the statistical analysis, Zhen-Gang Zhao and Yi-Ming Li provided suggestions for modification of the manuscript. All authors participated in the design, modification and optimization of this study.

Corresponding author

Correspondence to Mao Chen.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chai, Z., Wang, Y., Li, YM. et al. Correlations between geomagnetic field and global occurrence of cardiovascular diseases: evidence from 204 territories in different latitude. BMC Public Health 23, 1771 (2023). https://doi.org/10.1186/s12889-023-16698-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-023-16698-1

Keywords