Skip to main content

Spatial patterns and the associated factors for breast cancer hospitalization in the rural population of Fujian Province, China



Despite the known increasing incidence of breast cancer in China, evidence on the spatial pattern of hospitalization for breast cancer is scarce. This study aimed to describe the disparity of breast cancer hospitalization in the rural population of Southeast China and to explore the impacts of socioeconomic factors and heavy metal pollution in soil.


This study was conducted using the New Rural Cooperative Medical Scheme (NRCMS) claims data covering 20.9 million rural residents from 73 counties in Southeast China during 2015–2016. The associations between breast cancer hospitalization and socioeconomic factors and soil heavy metal pollutants were evaluated with quasi-Poisson regression models and geographically weighted Poisson regressions (GWPR).


The annual hospitalization rate for breast cancer was 101.40/100,000 in the studied area and the rate varied across different counties. Overall, hospitalization for breast cancer was associated with road density (β = 0.43, P = 0.02), urbanization (β = 0.02, P = 0.002) and soil cadmium (Cd) pollution (β = 0.01, P = 0.02). In the GWPR model, a stronger spatial association of Cd, road density and breast cancer hospitalization was found in the northeast regions of the study area while breast cancer hospitalization was mainly related to urbanization in the western regions.


Soil Cd pollution, road density, and urbanization were associated with breast cancer hospitalization in different regions. Findings in this study might provide valuable information for healthcare policies and intervention strategies for breast cancer.

Peer Review reports


Female breast cancer is a leading cause of cancer in the world and in China [1]. Nowadays, the increasing health burden of breast cancer has become a critical public health issue in China, with cases accounting for 9.1% of all newly diagnosed cancers and 3.9% of all cancer deaths in China [2]. More importantly, the distribution of breast cancer in China varies greatly [3], contributing to a difference in hospitalization for breast cancer [4, 5]. While such geographic disparities might be influenced by the differences in the distribution of environmental exposures and socioeconomic factors [6,7,8], previous evidence was, however, usually limited by the assumption of constant estimates over the entire geographic area and unable to capture the geographical variations in the relationships [9, 10].

Socioeconomic factors could influence women’s lifestyle [11], including dietary pattern, physical activity and reproductive behavior [12], which are known factors for breast cancer risk [1, 6, 13, 14] and hospitalization [15]. In addition, recent evidence suggests that hospitalizations for breast cancer might also be influenced by travel distances [16]and social security [17]. However, no study has systemically assessed the potential role of indicators for socioeconomic development and breast cancer hospitalization.

Moreover, special attention should be paid to the impact of environmental exposures, especially soil heavy metal pollution, considering the biological [18, 19] carcinogenesis effects from heavy metals. As heavy metals in soil are difficult to break down [20], soil heavy metals could serve as an indicator of long-term exposure to heavy metals [21]. While previous studies mainly focused on the associations between dietary heavy metals and breast cancer [22,23,24,25,26], there is less information on the association between chronic exposure to heavy metals in soil and breast cancer hospitalization, especially in the rural population who mainly face to the soil heavy metals.

Assessment of the spatial pattern and associated factors for breast cancer hospitalization is important for the planning of cancer-related health resources and identifying priority regions for possible prophylactic interventions. In this study, we focus on the spatial autocorrelation of breast cancer hospitalization in the rural population of Southeast China and explore the association between socioeconomic factors, soil heavy metals and breast cancer hospitalization by using general spatial autocorrelation. We further used the geographically weighted Poisson regression (GWPR) model to explore spatial variations of the associations between socioeconomic factors, soil heavy metals and breast cancer.

Materials and methods

Data sources

This study used data from the New Rural Cooperative Medical Scheme (NRCMS) in Fujian Province, which is located in the southeast of China (ranging from 23°33′ ~ 28°20′ N to 115°50′~ 120°40′ E). It covers an area of about 124,000 km 2, with a total population of 37.2 million in 2015, of which 20.9 million were rural. By the end of 2014, 98.9% of the rural population ( participated in the NRCMS. Therefore, we extracted medical insurance claims data of hospitalized female breast cancer patients within the NRCMS from January 1, 2015, to December 31, 2016. Due to the small number of rural residents in Taijiang County, Gulou County, Licheng County, Fengze County, Cangshan County, and Xiamen City, people in these counties were covered by the Urban Basic Medical Insurance and were excluded from the analysis, leaving 73 counties covering about 98% of the rural population in the study area.

Breast cancer hospitalization was identified using the International Classification of Diseases (ICD) 10 code C50. The hospitalization data consisted of age, sex, county of registered residence, year of hospitalization, and cause of hospitalization. Cases were selected based on the county of registered residence rather than hospital attended, which ensures a better representation of the spatial distribution of the exposed population to the environmental and socioeconomic factors. The rural female population data was obtained from the Fujian Statistical Yearbook ( All data were aggregated at the county level. The rate of breast cancer hospitalization was then calculated as the number of female breast cancer hospitalizations divided by female rural population in 73 counties.

Socioeconomic factors

Socioeconomic factors included population density, urbanization level, road density, per capita GDP, medical institutions, and health workers at the county level. The population density was defined as the number of people per square kilometer in each county. Urbanization level was defined as the proportion of people who live in urban areas among the entire population in each county. Road density was defined as the length of road per square kilometer in each county. Per capita GDP was defined as the ratio of the gross domestic product realized to the resident population in each county. Medical institutions were defined as the number of medical institutions per 1000 population in each county. Health workers were defined as the number of health workers per 1000 population in each county. We used the average level of these covariates during 2015–2016 to estimate their association with breast cancer hospitalization (supplementary Data 1 and Supplementary Table 1).

Soil heavy metals

In this study, the data on soil heavy metals was obtained from the Central Station of Environmental Monitoring of Fujian Province. The details of the sample collection and testing can be found in Chen’s study [27]. Briefly, all samples were from the forest land, grassland, orchard, tea garden, paddy field, dry land and tidal flat land in the map of land-use types in Fujian Province, and the grid was divided based on the geographic information system platform. The tea, orchards and arable land were laid out with a grid density of 8 km×8 km. For arable land in some key agricultural areas, 4 km×4 km of densely laid points were used. For tidal flats, one point was laid every 16 km along the coastline. Specifically, a total of 1114 surface soil samples were collected from the study area. The spatial distributions of heavy metals were estimated in terms of the ordinary kriging technique based on point data. ArcMap software was used for image registration and calculating pollutant regions of soil heavy metals at the county level. Supplementary Table 2 presents the descriptive statistics of soil heavy metals.

Statistical analysis

Spatial autocorrelation analysis

To explore the spatial distribution pattern of hospitalizations for breast cancer, the global Moran’s I index and the local Moran index were used in this study [28]. The global Moran’s I index evaluated the spatial autocorrelation degree of the entire study area, while the local Moran statistics result, namely, local indicators of spatial association (LISA), was used to explore the existence of anomalies in local areas [29]. The corresponding formulas [30] are as follows:

$$I = \frac{{\Sigma \begin{array}{*{20}{c}}k\\{i = 1}\end{array}\Sigma \begin{array}{*{20}{c}}k\\{j = 1}\end{array}({X_i} - \bar x)({X_j} - \bar x)}}{{{S^2}\Sigma \begin{array}{*{20}{c}}k\\{i = 1}\end{array}\Sigma \begin{array}{*{20}{c}}k\\{j = 1}\end{array}{W_{ij}}}}$$
$${I^*} = \Sigma \begin{array}{*{20}{c}}m\\{p \ne q}\end{array}{W_{pq}}{Z_p}{Z_q}$$

In the formulation, I represents the global Moran’s I index; k represents the number of grids. Xi and Xj stand for the values of grid i and grid j, and Wij is the spatial weight matrix between grid i and grid j (Wij = 1 if grid i is adjacent with grid j, and Wij = 0 otherwise). S2 is the variance of the whole grid. I* is the local Moran’s I index; Wpq is the standardized spatial weight matrix; and Zp and Zq represent standardized values of grid p and q.

Associations with socioeconomic factors and soil metal pollution

Quasi Poisson regression model was used to evaluate the overall association between the socioeconomic variables (Urbanization, Road density, Population density, Health workers, Medical institutions and Per capita GDP) and breast cancer hospitalization in women. In the univariate model, we tested the association with each of these variables, while the multivariable model mutually adjusted for the statistically significant variables identified in the univariate model.

Then, we selected statistically significant socioeconomic variables in the multivariable model, and incorporated each soli heavy metal pollutant into the model, separately [31]. This implementation allowed a global assessment of the existence of relationships, as well as their directions, between the different covariates and the rate of breast cancer hospitalization in women before analyzing the effect of the association by county.

In the next step, geographically weighted Poisson regression was used to identify the factors associated with breast cancer hospitalization. The geographical weight regression model allows the association between two factors to vary spatially, which has been used in several previous environmental studies [32,33,34,35]. The model of GWPR is shown in Eq. (2),

$$\begin{array}{l}{y_i} \sim Poisson\left[ {{N_i}exp\left( {{u_i}} \right)} \right]\\{u_i}(\{ {x_{k,i}}\} ) = {\Sigma _k}{\beta _k}({u_i},{v_i}){x_{k,i}}\end{array}$$

where yi, ui, xk, i and Ni are, respectively, dependent variable (the number of women who were hospitalized for breast cancer), the linear predictor, kth independent variable including the constant term, and the offset variable corresponding to population size at risk (defined as the average rural female population in 2015 and 2016) at the location i. It should be noted that the estimated risk of hospitalization for breast cancer at location i is given by the term of exp(ui). (ui, vi) is the x-y coordinate of the ith location; and coefficients βk(ui, vi) are assumed to be smoothly varying conditional on the location.

The estimates of the local coefficients at the location i are obtained by fitting a usual Poisson regression model to the data subset around the regression point i with a geographical weighting function. The standard errors of estimated local coefficients can be derived using local regression theory. A descriptive measure of goodness-of-fit for Poisson regression is the percent of deviance explained:

$${\rm{ pde}}{\upsilon _i}{\rm{ = 1 - de}}{\upsilon _i}{\rm{/nullde}}{\upsilon _i}$$

The data were screened and cleaned using exploratory data analysis using SPSS26.0 and Microsoft Excel. Global Moran’s I, LISA values, were measured using ArcMap software (version 10.6). We used GWR software (version 4.0.80) to implement a GWPR to explore spatial non-stationarity in our data. All vector maps were created using ArcMap software (version 10.6). The study was approved by the ethical committee at Fujian Medical University (committee’s reference number 2019-27).


Spatial distribution of breast cancer hospitalizations and socioeconomic factors

Table 1 illustrates the demographic characteristics of hospitalized breast cancer patients. The mean rate of breast cancer hospitalization in the rural population of Fujian Province was 101.40/105. The spatial distribution map (Fig. 1, left) shows that the values of breast cancer hospitalization rates in each county in the study area have obvious spatial aggregation.

Table 1 Demographic characteristics of breast cancer hospitalizations
Fig. 1
figure 1

Breast cancer hospitalization rate and LISA cluster map of 73 counties in southeast China

The Moran Index was calculated by using edges and corners as adjacent elements, and the global Moran index was 0.40 (p < 0.01), which suggested significant clustering. We further employed local Moran statistics to observe the local aggregation of breast cancer hospitalization rates. According to the result of local Moran statistics, two ‘high-high’ cluster regions, two ‘low-low’ cluster regions, as well as one ‘low-high’ cluster region were detected as statistically significant clusters (Fig. 1, right). From the LISA cluster map, the significant ‘high-high’ clusters A and B were located in southeastern Fujian. Additionally, another ‘high-high’ cluster C was located in central Fujian. The large clusters a and b with low risk existed in southern and mid-eastern Fujian, respectively. Moreover, there was ‘low-high’ cluster L between ‘high-high’ clusters A and B.

Table 2 displays the results of the Quasi-Poisson Regression model with socioeconomic factors. In univariate analyses, urbanization, population density and road density were significantly associated with breast cancer hospitalization rate. In the multivariable model, urbanization rate (β = 0.02, P = 0.002) and road density (β = 0.43, P = 0.02) were significantly associated with breast cancer hospitalizations. However, medical institutions, health workers and per capita GDP were not significantly associated with breast cancer hospitalization rate at the 0.05 level.

Table 2 Summary of parameters in the Quasi-Poisson regression model in socioeconomic variables

Spatial relationship between soil heavy metals and breast cancer hospitalization

In this study, ArcMap software was used to extract the polluted areas of arsenic (As), cadmium (Cd), chromium (Cr), copper (Cu), mercury (Hg), nickel (Ni), lead (Pb), and zinc(Zn) in these 73 counties. Cd was the most serious pollutant covering 54.27% of the total area, followed by Hg, As, and Pb. To eliminate the influence of the area factor of each county, the polluted area of each county was converted into the percentage of soil polluted area of each county. Distribution of percentage of polluted areas in 73 counties of Fujian province was visualized by color in Fig. 2.

Fig. 2
figure 2

Percentage of the contaminated area of soil heavy metals in Fujian Province

After adjustment for urbanization and road density, Cd was the only heavy metal significantly associated with the hospitalization rate of breast cancer (β = 0.005, P = 0.016; Table 3). Using the GWPR analysis, the spatial variation of the estimated effects from urbanization, road density and Cd over the 73 counties is shown in Fig. 3. The association between breast cancer hospitalization and urbanization was the highest in the north and the west of Fujian Province, while the association with Cd was the highest in the northeast part. The strongest association between the hospitalization rate of breast cancer and road density was found in counties in Northeast and Southwest Fujian.

Table 3 Summary of parameters in the Quasi-Poisson regression model in soil heavy metals
Fig. 3
figure 3

Parameters of predicting variables by county in the GWPR

The red color on the map indicates a higher value of a local parameter estimate


Key results

In this study, we investigated the association of breast cancer hospitalization rates with socioeconomic factors and the spatial distribution of soil heavy metals at the county level in the rural population of Fujian Province, China. Urbanization and road density were positively associated with the hospitalization rate of breast cancer. We further used geographically weighted Poisson regression analysis to explore the spatial relationship between Cd, road density, urbanization and the hospitalization rate of breast cancer. We found a stronger spatial association of Cd, road density and breast cancer in the northeast of Fujian, while the spatial relationship between urbanization and breast cancer was in western Fujian.

Due to the different spatial distributions of environmental and socioeconomic factors, the distribution of disease presents dramatic spatial heterogeneity in geographic regions. Our study found a distinct pattern of breast cancer hospitalization in Fujian Province, characterized by significantly higher clusters in the southeast of Fujian versus lower clusters in southern regions. The spatial autocorrelation analysis further showed three high clusters (A, B, C), two low clusters (a,b), and one high-low cluster(L) for the hospitalization of breast cancer. This spatial distribution may be partly because of the uneven distribution of medical resources for the rural population in Fujian province, with most provincial and high-quality hospitals located in relatively developed coastal cities. However, it’s worth noting that the developing area (C) had a high hospitalization rate, suggesting stronger demand for medical resources and unequal resource distribution in this area. Recently, due to the rapid economic development in China [36], even people in developing areas in Fujian province have increased their incomes and improved health awareness of cancer. In addition, women have more frequent exposure to environmental pollutants (such as heavy metals, noise and dioxins) [37], resulting in an increased incidence of breast cancer and demand for related health services. Therefore, sufficient attention should be paid to the areas in need and immediate action should be taken to improve the fairness of medical resources allocation in Fujian province.

In this study, we found a positive association between urbanization and breast cancer hospitalization in the majority of the counties, which was consistent with the previously reported association between urbanization and breast cancer risk [38,39,40,41,42,43]. In recent years, more resources have shifted in China towards the development of small towns, especially in southeastern coastal areas of China, implying that these areas are experiencing a rapid development period [44]. Women in these areas may have a higher living standard, increased health awareness, and were more likely to afford breast cancer treatments. In addition, the development of economy may also lead to more life pressures, higher body mass index(BMI) [45], and changes in women’s lifestyle and reproductive behaviors [12, 46], which are known risk factors for breast cancer and may be responsible for the observed associations.

Together with the increase of urbanization, road density is also increasing in China recently. Positive association between road density and breast cancer hospitalization was found in many counties, possibly due to several reasons. First, higher road density means that there may be more traffic pollution. Studies on traffic-related pollutants have suggested that NO2 might be associated with a higher risk of breast cancer [47, 48]. PAHs, mainly originating from traffic, were also associated with breast cancer [49]. What’s more, one recent study from a Danish nationwide cohort found that road traffic noise was related to an increased risk of breast cancer [50]. Second, a higher road density also means more convenient transportation to medical institutions [51], which would consequently increase the number of hospitalizations. The strongest association between the low hospitalization rate of breast cancer and low road density in Southwest Fujian further supported this possibility, considering the quite limited medical resources in western Fujian. However, the strong association in Northeast Fujian might be the result of environmental pollution.

We observed substantial geographic heterogeneity in the association between Cd and breast cancer hospitalization, with the strongest association observed among women in the northeast region (region D). Cd has been classified as a known carcinogen [52]. The increased hospitalization for breast cancer in regions with higher Cd pollution is biologically plausible as Cd is involved in the carcinogenesis of breast cancer by activating ER [53, 54], inducing cell proliferation, differentiation, and apoptosis [55, 56], and modulating gene expression. Epidemiological studies have also indicated that urinary Cd levels were higher in breast cancer patients than in normal controls [57, 58], while breast cancer tissue had four times higher Cd levels than surrounding healthy tissue [31].

Considering the geographical disparity of risk factors for breast cancer hospitalization, different measures should be taken to improve women’s health. In western areas, which are in the rapid development stage, urbanization is the main risk factor for breast cancer and women have more demand for high quality medical resources. Health education to increase awareness of socioeconomic and lifestyle risk factors might reduce their future risk of breast cancer. In addition, regular screening services could be provided to detect breast cancer earlier and avoid hospitalizations caused by severe diseases. What’s more, for those coastal east areas, road density and Cd are the main risk factors. It is therefore vital to strengthen environmental governance, such as limiting the discharge of industrial heavy metals, especially Cd, and promoting the use of electric vehicles.

The major strength of our study is the use of geographically weighted regression models, which are able to detect the spatially varied magnitude of the associations and provide insight into priority regions for health resource planning and allocation. Moreover, we had access to a variety of data on socioeconomic factors and soil heavy metals, allowing us to explore the potential influencing factors.

Despite these, our findings should be interpreted in the context of several limitations. First, the design of ecological studies limits the ability to prove causality. The observed associations may be influenced by many unobserved confounding factors, such as other cultural factors, policy factors, and other natural factors. Second, the research unit of this paper is the county level divided by administrative regions. However, the distribution of environmental factors may be related to the natural geographical environment, and the corresponding effect might be attenuated. Subsequent studies are, therefore, needed to confirm the relationship between Cd and breast cancer using regions divided by the natural environment.


Soil Cd pollution, road density, and urbanization were associated with breast cancer hospitalization in southeast China. The findings in this paper will help to understand the regional impact of Cd, road density and urbanization on the hospitalization rate of breast cancer, and provide valuable reference for the prevention strategy for breast cancer in southeast China and similar regions.

Data availability

The hospitalization data that support the findings of this study are applied from the New Rural Cooperative Medical Scheme (NRCMS) database, and is therefore not publicly available. Part of the population dataset analyzed in this study can be found at, and provided in supplementary data 1. Data is, however, available from the corresponding author upon reasonable request and with permission of the Fujian Medical Security Bureau.


  1. Sung H, Ferlay J, Siegel RL, et al. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and Mortality Worldwide for 36 cancers in 185 Countries[J]. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Y Zhang, Z Lyu Fsong, et al. Trends of incidence and mortality of breast cancer worldwide and in China[J]. J Multidisciplinary Cancer Management(Electronic Version). 2021;7(02):14–20.

    Google Scholar 

  3. Fei X, Lou Z, Christakos G, et al. A Geographic Analysis about the Spatiotemporal pattern of breast Cancer in Hangzhou from 2008 to 2012[J]. PLoS ONE. 2016;11(1):e147866.

    Article  Google Scholar 

  4. Q Xu Xzeng. Analysis of 30075 cancer cases covered by New Rural Cooperative Medical Scheme in Jinjiang city,Fujiang province during 2009–2013[J]. J Practical Oncol. 2015;30(05):477–82.

    Google Scholar 

  5. Y Xue, L Li Lzhen. Analysis of clinical characteristic and incidence tend on inpatients with breast cancer from 1990~1999[J]. Chinese Journal of Primary Medicine Pharmacy, 2002(05):13–15.

  6. Gentil J, Dabakuyo T S, Ouedraogo S, et al. For patients with breast cancer, geographic and social disparities are independent determinants of access to specialized surgeons. A eleven-year population-based multilevel analysis[J]. BMC Cancer. 2012;12:351.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Merletti F, Galassi C, Spadea T. The socioeconomic determinants of cancer[J]. Environ Health. 2011;10(Suppl 1):7.

    Article  Google Scholar 

  8. Robert SA, Strombom I, Trentham-Dietz A, et al. Socioeconomic risk factors for breast cancer: distinguishing individual- and community-level effects[J]. Epidemiology. 2004;15(4):442–50.

    Article  PubMed  Google Scholar 

  9. Tu J, Xia Z G. Examining spatially varying relationships between land use and water quality using geographically weighted regression I: model design and evaluation[J]. Sci Total Environ. 2008;407(1):358–78.

    Article  CAS  PubMed  Google Scholar 

  10. Wang N, Mengersen K, Tong SL, et al. County-level variation in the long-term association between PM2.5 and lung cancer mortality in China[J]. SCIENCE OF THE TOTAL ENVIRONMENT; 2020. p. 738.

  11. Fei X, Wu J, Kong Z, et al. Urban-rural disparity of breast cancer and socioeconomic risk factors in China[J]. PLoS ONE. 2015;10(2):e117572.

    Article  Google Scholar 

  12. Linos E, Spanos D, Rosner B A, et al. Effects of reproductive and demographic changes on breast cancer incidence in China: a modeling analysis[J]. J Natl Cancer Inst. 2008;100(19):1352–60.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Adams J, White M. Are there socioeconomic gradients in stage and grade of breast cancer at diagnosis? Cross sectional analysis of UK cancer registry data[J]. BMJ. 2004;329(7458):142.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Harper S, Lynch J, Meersman S C, et al. Trends in area-socioeconomic and race-ethnic disparities in breast cancer incidence, stage at diagnosis, screening, mortality, and survival among women ages 50 years and over (1987–2005)[J]. Cancer Epidemiol Biomarkers Prev. 2009;18(1):121–31.

    Article  PubMed  Google Scholar 

  15. Silva D. Hospitalization and ambulatory costs related to breast cancer due to physical inactivity in the Brazilian state capitals[J].PLoS One, 2022,17(1):e261019.

  16. de Souza B C, Dos S F F SL, et al. Regional disparities in the flow of access to breast cancer hospitalizations in Brazil in 2004 and 2014[J]. BMC Womens Health. 2020;20(1):137.

    Article  Google Scholar 

  17. Ventura-Alfaro C E, Torres-Mejia G, Avila-Burgos L S. Hospitalization and mortality in Mexico due to breast cancer since its inclusion in the catastrophic expenditures scheme[J]. Salud Publica Mex. 2016;58(2):187–96.

    Article  PubMed  Google Scholar 

  18. Some metals and metallic compounds[J]. IARC Monogr Eval Carcinog Risk Chem Hum. 1980;23:1–415.

    Google Scholar 

  19. Isidori M, Cangiano M, Palermo F A, et al. E-screen and vitellogenin assay for the detection of the estrogenic activity of alkylphenols and trace elements[J]. Volume 152. Comparative Biochemistry And Physiology C-Toxicology & Pharmacology; 2010. pp. 51–6. 1.

  20. Burges A, Epelde L. Impact of repeated single-metal and multi-metal pollution events on soil quality[J]. Chemosphere. 2015;120:8–15.

    Article  CAS  PubMed  Google Scholar 

  21. Chen K, Liao Q L, Ma Z W, et al. Association of soil arsenic and nickel exposure with cancer mortality rates, a town-scale ecological study in Suzhou, China[J]. Environ Sci Pollut Res Int. 2015;22(7):5395–404.

    Article  CAS  PubMed  Google Scholar 

  22. Adams S V, Newcomb P A, White E. Dietary cadmium and risk of invasive postmenopausal breast cancer in the VITAL cohort[J]. Cancer Causes Control. 2012;23(6):845–54.

    Article  Google Scholar 

  23. Adams S V, Quraishi S M, Shafer M M, et al. Dietary cadmium exposure and risk of breast, endometrial, and ovarian cancer in the Women’s Health Initiative[J]. Environ Health Perspect. 2014;122(6):594–600.

    Article  PubMed  Google Scholar 

  24. Garcia-Esquinas E, Pollan M. Arsenic exposure and cancer mortality in a US-based prospective cohort: the strong heart study[J]. Cancer Epidemiol Biomarkers Prev. 2013;22(11):1944–53.

    Article  PubMed  Google Scholar 

  25. Itoh H, Iwasaki M, Sawada N, et al. Dietary cadmium intake and breast cancer risk in japanese women: a case-control study[J]. Int J Hyg Environ Health. 2014;217(1):70–7.

    Article  CAS  PubMed  Google Scholar 

  26. NÚÑEZ, O, FERNÁNDEZ-NAVARRO P, MARTÍN-MÉNDEZ I, et al. Arsenic and chromium topsoil levels and cancer mortality in Spain[J]. Environ Sci Pollut Res. 2016;23(17):17664–75.

    Article  Google Scholar 

  27. Z Chen. Geo-accumulation index and potential ecological risk on Soil Heavy Metals:an evaluation of case in Fujian[J]. J Subtropical Resour Environ. 2016;11(04):37–45.

    Google Scholar 

  28. Soleimani M. Spatial and temporal analysis of myocardial infarction incidence in Zanjan province, Iran[J]. BMC Public Health. 2021;21(1):1667.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Xie Z, Qin Y, Li Y, et al. Spatial and temporal differentiation of COVID-19 epidemic spread in mainland China and its influencing factors[J]. Sci Total Environ. 2020;744:140929.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. L Ding, Y Xiao-Mei W, Sheng-Lan, et al. Research on evolutionary process and driving forces of spatially economic disparity within Counties based on technology of ESDA—GIS[J]. Econ Geogr. 2013;33(05):31–6.

    Google Scholar 

  31. Andelkovic M, Djordjevic A B, Miljakovic E A, et al. Cadmium tissue level in women diagnosed with breast cancer - a case control study[J]. Environ Res. 2021;199:111300.

    Article  PubMed  Google Scholar 

  32. Mor Builv et al. Z, Chemtob D,. Use of Geographically Weighted Poisson Regression to examine the effect of distance on Tuberculosis incidence: A case study in Nam Dinh, Vietnam[J]. PLOS ONE, 2018,13(11).

  33. Liu Y, Ji Y J, Shi Z B et al. The Influence of the Built Environment on School Children’s Metro Ridership: An Exploration Using Geographically Weighted Poisson Regression Models[J]. SUSTAINABILITY, 2018,10(12).

  34. Ribeiro M C, Pereira MJ. Modelling local uncertainty in relations between birth weight and air quality within an urban area: combining geographically weighted regression with geostatistical simulation[J]. Volume 25. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH; 2018. pp. 25942–54. 26.

  35. Wang JY, Wang S J, LI SJ. Examining the spatially varying effects of factors on PM2.5 concentrations in chinese cities using geographically weighted regression modeling[J]. Environ Pollut. 2019;248:792–803.

    Article  CAS  PubMed  Google Scholar 

  36. Y Yang YXU. Research on spatial and temporal pattern of Urban Expansion in Fujian based on DMSP/OLS Nighttime Light Data[J]. Urbanism and Architecture. 2021;18(30):17–9.

    Google Scholar 

  37. Brody JG, Moysich K B, Humblet O, et al. Environmental pollutants and breast cancer: epidemiologic studies[J]. Cancer. 2007;109(12 Suppl):2667–711.

    Article  CAS  PubMed  Google Scholar 

  38. Baquet C R. Socioeconomic factors and breast carcinoma in multicultural women[J]. Cancer. 2000;88(5 Suppl):1256–64.

    Article  PubMed  Google Scholar 

  39. Dano H, Andersen O. Socioeconomic status and breast cancer in Denmark[J]. Int J Epidemiol. 2003;32(2):218–24.

    Article  PubMed  Google Scholar 

  40. Faggiano F, Partanen T, Kogevinas M et al. Socioeconomic differences in cancer incidence and mortality[J].IARC Sci Publ, 1997(138):65–176.

  41. Pollan M. High-risk occupations for breast cancer in the swedish female working population[J]. Am J Public Health. 1999;89(6):875–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Robert SA, Strombom I, Trentham-Dietz A, et al. Socioeconomic risk factors for breast cancer: distinguishing individual- and community-level effects[J]. Epidemiology. 2004;15(4):442–50.

    Article  PubMed  Google Scholar 

  43. Sarfati D, Blakely T, Shaw C, et al. Patterns of disparity: ethnic and socio-economic trends in breast cancer mortality in New Zealand[J]. Cancer Causes Control. 2006;17(5):671–8.

    Article  PubMed  Google Scholar 

  44. Y Yang YXU. Research on spatial and temporal pattern of Urban Expansion in Fujian based on DMSP/OLS Nighttime Light Data[J]. Urbanism and Architecture. 2021;18(30):17–9.

    Google Scholar 

  45. Liu W, Wang L J, Qi JL, et al. [Disease burden of breast cancer in women in China, 1990–2017][J]. Zhonghua Liu Xing Bing Xue Za Zhi. 2021;42(7):1225–30.

    CAS  PubMed  Google Scholar 

  46. Howell A, Anderson A S, Clarke R B, et al. Risk determination and prevention of breast cancer[J]. Breast Cancer Res. 2014;16(5):446.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Crouse D L, Goldberg M S, Ross N A, et al. Postmenopausal breast cancer is associated with exposure to traffic-related air pollution in Montreal, Canada: a case-control study[J]. Environ Health Perspect. 2010;118(11):1578–83.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Reding K W, Young M T, Szpiro A A, et al. Breast Cancer risk in relation to Ambient Air Pollution exposure at residences in the Sister Study Cohort[J]. Cancer Epidemiol Biomarkers Prev. 2015;24(12):1907–9.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Mordukhovich I, Beyea J, Herring A H, et al. Vehicular traffic-related polycyclic aromatic hydrocarbon exposure and breast Cancer incidence: the long island breast Cancer Study Project (LIBCSP)[J]. Environ Health Perspect. 2016;124(1):30–8.

    Article  CAS  PubMed  Google Scholar 

  50. Sorensen M, Poulsen A H, Kroman N, et al. Road and railway noise and risk for breast cancer: a nationwide study covering Denmark[J]. Environ Res. 2021;195:110739.

    Article  PubMed  Google Scholar 

  51. Fei X, Chen W, Zhang S, et al. The spatio-temporal distribution and risk factors of thyroid cancer during rapid urbanization-A case study in China[J]. Sci Total Environ. 2018;630:1436–45.

    Article  CAS  PubMed  Google Scholar 

  52. IARC. Cadmium and cadmium compounds[J]. Rep Carcinog. 2011;12:80–3.

    Google Scholar 

  53. Ali I, Damdimopoulou P, Stenius U, et al. Cadmium-induced effects on cellular signaling pathways in the liver of transgenic estrogen reporter mice[J]. Toxicol Sci. 2012;127(1):66–75.

    Article  CAS  PubMed  Google Scholar 

  54. Ali I, Hurmerinta T, Nurmi T, et al. From pure compounds to complex exposure: Effects of dietary cadmium and lignans on estrogen, epidermal growth factor receptor, and mitogen activated protein kinase signaling in vivo[J]. Toxicol Lett. 2016;253:27–35.

    Article  PubMed  Google Scholar 

  55. Kluxen F M, Diel P, Hofer N, et al. The metallohormone cadmium modulates AhR-associated gene expression in the small intestine of rats similar to ethinyl-estradiol[J]. Arch Toxicol. 2013;87(4):633–43.

    Article  PubMed  Google Scholar 

  56. Huff J, Lunn R M, Waalkes M P, et al. Cadmium-induced cancers in animals and in humans[J]. Int J Occup Environ Health. 2007;13(2):202–12.

    Article  CAS  PubMed  Google Scholar 

  57. Mcelroy J A, Shafer M M, Trentham-Dietz A, et al. Cadmium exposure and breast cancer risk[J]. J Natl Cancer Inst. 2006;98(12):869–73.

    Article  PubMed  Google Scholar 

  58. Nagata C, Nagao Y, Nakamura K, et al. Cadmium exposure and the risk of breast cancer in japanese women[J]. Breast Cancer Res Treat. 2013;138(1):235–9.

    Article  CAS  PubMed  Google Scholar 

Download references


Grateful acknowledgment is extended to Zengwen Chen for his offering the original data of soil heavy metals in Fujian province.


Open access funding provided by Karolinska Institute. HY is supported by the National Natural Science Foundation of China [grant no: 82204132], Natural Science Foundation of Fujian Province [grant no: 2021J01721], the Startup Fund for High-level Talents of Fujian Medical University [grant no: XRCZX2020007], and Startup Fund for Scientific Research, Fujian Medical University [grant no: 2019QH1002].

Author information

Authors and Affiliations



XH, ZH and HY were responsible for the study concept and design. MS, XW, XT, and ZR did the data and project management. MS did the data cleaning and analysis. MS, XW, and HY interpreted the data. MS and HY drafted the manuscript. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

Corresponding author

Correspondence to Haomin Yang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Fujian Medical University (committee’s reference number: 2019-27). All methods were carried out in accordance with the principles stated in the Declaration of Helsinki. As the database for NRCMS was designed for administration purposes rather than academic research, patients were given oral informed consent when involved in the NRCMS. The oral informed consent procedure was approved by the Ethics Committee of Fujian Medical University.

Consent for publication

Not applicable.

Competing interests

Not applicable.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Table 1 source of the used dataset

Supplementary Table 2 contaminated area of soil metal elements (km2)

Supplementary Table 3 Summary of local parameters in the GWPR

Supplementary Data 1 Socioeconomic data used in the analysis

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 The Creative Commons Public Domain Dedication waiver ( 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

Song, M., Huang, X., Wei, X. et al. Spatial patterns and the associated factors for breast cancer hospitalization in the rural population of Fujian Province, China. BMC Women's Health 23, 247 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Soil heavy metals
  • Socioeconomic factors
  • Breast cancer hospitalization
  • Spatial heterogeneity
  • GWPR