Bridging geographical disparities across 368 townships with healthcare system and socioeconomic factors in Taiwan

A universal health insurance program such as the National Health Insurance in Taiwan offers a wide coverage and increased access to healthcare services. Despite its ongoing efforts to enhance healthcare accessibility, differences in health for people living in urban and resource-deprived areas remain substantial. To investigate the longitudinal impact of the healthcare system and other potential structural drivers such as education and economic development on geographical disparities in health, we designed a panel study with longitudinal open secondary data, covering all 368 townships in Taiwan between 2013 and 2017. Our findings indicated higher mortality rates in the mountainous and rural areas near the east and south regions of the island in both years. Multivariate analyses showed an increase in the density of primary care physicians (PCP) was associated with lower all-cause mortality (β = − 0.72, p < 0.0001) and cardiovascular disease mortality (β = − 0.41, p < 0.0001). Effect of PCP is evident, but merely focusing on access to healthcare is still not enough. Additional measures are warranted to address the health disparities existing between urban and underprivileged areas.

For certain population groups, actual uptake of healthcare services can be often impeded by structural barriers despite universal availability of these resources within a health system.People living in different geographical areas may experience different levels of accessibility due to physical distances and barriers, which subsequently leads to geographical variation in health [1][2][3] .Numerous studies have highlighted the significant role of individuallevel factors in shaping geographic inequalities 4,5 .Other studies, however, have shown that the influence from the nature of places, or the place's structural elements, are much more important in our understanding of geographical variations in health.For example, a study suggested that regional healthcare markets, which are made up by predominantly groups of healthcare providers (supply), can generate significant geographic variation in healthcare utilization and spending after accounted for patient-level (demand) characteristics 6 .Juhn et al. (2021) 7 examined COVID-19 epidemiology and found that community factors such as income level and living arrangements are significant factors of COVID-19 burden disparities.Other studies have also found significant impact of structural determinants of health on health equity such as education, employment, and other social and public policies 8,9 .
Taiwan is an island country with a population of roughly 23 million people and a land area of just slightly greater than Belgium (~ 36,000 square kilometers).It is divided into two sides by a central mountain range that spans the length of the country: the western terrain, where 90% of the population lives, is mostly flat and the eastern topography with mostly mountains and natural undeveloped landscape.There are five major geographic regions: north, center, south, east, and offshore archipelago islands which include Penghu, Matsu, and Kinmen.These regions consist of 368 cities/townships and districts in total, with six major municipalities that are directly Measurements.Outcomes.Main outcome measures in this study were all-cause mortality (per 100,000), cardiovascular disease (CVD) mortality (per 100,000), and mortality attributable to unintentional injuries (UI) (per 100,000).All-cause mortality was computed from the total number of deaths from all causes (International Classification of Diseases, 10th Revision, Clinical Modification (ICD-10-CM): A00-Y98).CVD mortality was calculated using the number of deaths from heart-related diseases excluding hypertension (ICD-10-CM: I01-I02.0,I05-I09, I20-I25, I27, I30-I52).UI mortality accounted for the total number of deaths from accidents or unplanned events (ICD-10-CM: V01-X59, Y85-Y86).All outcome variables were measured at the townshiplevel and age-standardized with the population's age distribution.Our source of data was the Causes of Death by Township Statistics from the Ministry of Health and Welfare (MOHW), Taiwan.For UI deaths, misclassification bias is unlikely since the primary cause of death on the death registry is typically recorded based on the immediate cause of death.In the rare event of accidents or unplanned events, they should be appropriately recorded as the primary cause of death.
Healthcare system.To characterize the healthcare system at the township-level, we considered two main aspects under the NHI: workforce and institutional capacity.For workforce, primary care physician (PCP) density (per 100,000 population) was used to represent the availability of PC doctors in each township to offer generalist care.Density of specialists was not considered as statistics were categorized by medical specialty and a physician might be registered under multiple specialties, preventing an accurate estimate of specialists.
To sufficiently describe the institutional capacity of the hospitals, we used density of acute care (AC) beds (per 100,000 population), density of long-term care (LTC) beds (per 100,000 population), and density of intensivecare unit (ICU) beds (per 100,000 population).AC beds represent the number of hospital beds dedicated to short term care, not including those for mental conditions.LTC beds include the total number of hospital beds for chronic conditions, and ICU beds measure the number of hospital beds provided in intensive care wards, including neonatal ICUs.Although these three measures appear quite similar, they each capture distinct dimensions of a hospital's capacity for care delivery and are mutually exclusive.Socioeconomic structure.Aside from healthcare system attributes, we included township-level structural determinants of health like education, economic development, and demographics.For education, we summed up the proportions of populations having achieved each level of education and applied the following equation to estimate the education level of each township according to the number of years in education (starting from elementary school up to postgraduate studies; where 29 is the maximum number of years of education that can be achieved assuming two years of masters studies and five years of doctoral studies on average).Thus, a higher calculated score would indicate higher educational achievement for each township/city.For the level of economic development, we ranked the 368 townships according to their annual median household income and categorized them into tertiles: low, middle, and high.Cut-off points were made at 33% (New Taiwan (NT) $66,680; ~ US$2,245) and 66% (NT$133,359; ~ US$4,490) using 2013 as the reference year.If the township annual median household income was ≤ NT$66,680, its economic development was classified as 'low'; if it was between NT$66,681 and NT$133,358, its economic development was classified as 'middle'; and if it was ≥ NT$133,359, then its economic development was classified as 'high' .Income values were not adjusted since the inflation rate from 2013 to 2017 was acceptable at 2.92%.
Demographics.For the general demographics profile of each township, population density, sex ratio, and proportion of population aged 65 or over were measured for the two years of study.Population density was calculated by dividing population size by the geographical surface area (km 2 ).Sex ratio was the number of males per 100 females.Proportion of population aged over 65 years was the population size aged 65 years or over divided by total population in each township.These demographic attributes were considered because they might exert confounding influence on our outcomes of interest.

Data analysis.
For 2013 and 2017, we calculated descriptive statistics for population health outcomes, healthcare system, education, economic development, and demographics in the 368 townships.Geospatial distribution of health outcomes was mapped across townships using QGIS 3.26 which is an open source geographic information system for geographical mapping.Due to the non-normality nature of the study variables, bivariate analyses were performed using Kruskal-Willis test to compare distributions of the three mortality rates across different economic levels for each year, while Mann-Whitney U test was conducted to compare the health indicators between the two years.Spearman correlation metric was used to examine correlations between all study variables.Effect size (ES) was estimated using two different tests: Cohen's D for continuous variables and Cramer's V for categorical variables.
Fixed effects regression model estimated the impact of healthcare system and structural determinants on the change in population health over time, from 2013 to 2017.Fixed effects model was chosen in order to avoid bias due to unobserved township-level factors that vary over time or can be correlated with study variables.In the multivariate analyses, we weighted the models by population density, proportion of population aged over 65, sex ratio, and adjusted for all predictors.Multicollinearity diagnostics were performed to confirm no predictor variables were collinear.All statistical analyses were performed using SAS statistical software package, version 9.4.

Results
Table 1 shows the distribution of study variables for years 2013 and 2017.Of the three population health indicators, UI mortality showed the most obvious decline from 2013 to 2017 (ES: 0.0038), followed by all-cause mortality (ES: 0.0023), while CVD mortality marginally increased (ES: 0.0005), all of which were statistically insignificant (p > 0.05).Mean standardized all-cause mortality was 514.72 deaths per 100,000 in 2013 and 503.06 deaths per 100,000 in 2017; mean standardized CVD mortality was 56.73 deaths per 100,000 in 2013 and 57.20 deaths per 100,000 in 2017; and mean standardized UI mortality was 34.67 deaths per 100,000 in 2013 and 32.46 deaths per 100,000 in 2017.
For healthcare system, changes in all variables were also statistically insignificant.Mean PCP density increased from 54.13 physicians per 100,000 in 2013 to 55.24 physicians per 100,000 in 2017 (ES: 0.0005).Densities of AC, LTC, and ICU beds showed minimal changes.
Mean level of educational achievement improved significantly over the 4-year period, with an ES of 0.1766 (p < 0.01).Economic development level also shifted slightly towards the higher end, but without statistical significance (high economic development: 29.35% in 2013 and 33.97% in 2017, p = 0.9621).
Figure 1 presents the geographical distribution of the three health indicators separately for 2013 (left panel) and 2017 (right panel).Although all-cause and UI mortality rates have generally decreased over the 4-year period, higher mortality rates were still detected in the mountainous and rural areas near the east and south regions of the island such as Taoyuan District in Kaohsiung and Zhuoxi Township in Hualien.Conversely, lowest mortality rates were found in the metropolitan areas concentrated in the west and towards the country's capital in the north, irrespective of time.
When we stratify the townships according to their economic development, there was a significant difference (p < 0.0001) in all-cause mortality rates across different levels for 2013 and 2017 (Fig. 2).If we compare the indicator between 2013 and 2017 for each level, most substantial improvement was observed in townships with middle development, whose rate declined from 493.4 deaths per 100,000 in 2013 to 465.6 deaths per 100,000 in 2017, with an ES of 0.0119.For townships with high-and low-development levels, the decrease in all-cause mortality was to a lesser degree (ES: 0.0007).However, the declines were all statistically insignificant regardless of development level.
For CVD mortality, townships with high economic development exhibited significantly lower rates compared to both middle-and low-level townships (p < 0.0001 for 2013 and p < 0.001 for 2017) (Fig. 2).It is important to highlight that there was a rising (but statistically insignificant) trend in CVD mortality rate from 2013 to 2017, which was consistent across all townships irrespective of their level of economic development.
For UI mortality, a gradient emerged, with low-level economy townships showing the most significant decline and high-level economy townships displaying the least improvement: low-level (39.9 to 33.8 deaths per 100,000, ES: 0.0181, p = 0.0327), middle-level (30.2 to 28.5 deaths per 100,000, ES: 0.0041, p = 0.3227), and high-level (19.7 to 19.3 deaths per 100,000, ES: 0.0009, p = 0.6738).In both years, significant differences in UI mortality rates existed across townships of different economic development levels (p < 0.0001).In Tables 2 and 3, strength of correlation between continuous explanatory variables and the three health outcomes is presented separately for 2013 and 2017.For both years, densities of PCP, AC bed, LTC bed, and ICU bed presented significant negative association with all mortality rates.AC bed density and ICU bed density seemed to have the strongest negative correlation with all three mortality rates, followed by PCP density and then LTC bed density.When examining socioeconomic influence, township's educational achievement achieved even stronger negative correlation with all three mortality rates compared to healthcare system variables in both years.This suggests that improving the population's educational attainment may help prevent mortality risks from CVD, UI, or all causes.Townships with higher sex ratios were observed with higher CVD, UI, and all-cause mortality rates.
Table 4 shows the results from both crude and adjusted fixed effects models.When taking into account the change of township characteristics between 2013 and 2017, we see that an increase in PCP density was attributable to a decline in all-cause mortality (crude: β = − 0.72, p < 0.0001; adjusted: β = − 0.72, p < 0.0001) and CVD mortality (crude: β = − 0.41, p < 0.0001; adjusted: β = − 0.41, p < 0.0001) for the 368 townships.This indicates that, each unit increase in township's PCP density would cause its all-cause mortality and CVD mortality to fall by 0.72 and 0.41, respectively.Similarly, an increase in AC bed density would indicate a drop in UI mortality by 0.05 in the adjusted model (β = − 0.05, p = 0.01), while an increase in ICU bed density would have the opposite effect (β = 0.36, p = 0.005).The socioeconomic determinants, education or economic development, did not seem to exert a significant impact on any of the population health outcomes after adjusting for all covariates.

Discussion
Through our panel study of geographical disparities in Taiwan, we discovered that population health had improved overall from 2013 to 2017.Within the 4-year period, greatest improvements in UI and all-cause mortality were generally seen in areas with low and middle economic development.Nevertheless, significant   www.nature.com/scientificreports/however, still face rising health inequalities despite their sustained improvements in health.In Europe, for example, a review of health disparities in the continent suggested obvious patterns of within-country geographical differences in health, including the UK, Spain, France, and Germany 16 .Similarly, a repeated cross-sectional study from Canada revealed worsening inequalities in self-reported health among employed and unemployed citizens between 2000 and 2014 17 .Studies from Asia provided similar findings 18,19 .
There are several possible factors contributing to the geographical disparities in Taiwan.As our findings indicate, access to primary care (PC) may be one of them.In her prominent series of work on PC, Professor Barbara Starfield highlighted the positive impact of PC on health and the effect of the supply of PCPs on reducing incomeinequality on health [20][21][22] .This has also been the case for different geographical areas in Taiwan.By facilitating first-contact and people-centered care, in contrast to specialty care, PC places a greater focus on prevention and effectively increases access to health services for disadvantaged population groups 21 .A previous study of 102 countries found that expanded coverage of PC services is an affordable means that can be beneficial to population particularly in low and middle-income settings 23 .Evidence suggested that provision of PC in socially and economically deprived areas may decrease the likelihood of child mortality, potentially due to adequate referral mechanisms, continuity of care, and receiving PC in public facilities 24 .
Since PCP was found to be protective factor for all-cause mortality and CVD mortality in this study, but not UI mortality, it is possible to speculate that the former two health indicators included in part mortality amenable to healthcare.These are deaths from certain causes that are potentially preventable with timely and effective health care.This has often been used as a measure of health system performance in developed countries and includes deaths due to tuberculosis, diabetes mellitus, and ischemic heart disease 25 .Our evidence suggests the vital role of PCPs in sustaining the performance of Taiwan's healthcare system.
It is essential to promote an equitable distribution of PCP services among townships.Nonetheless, the impact of PCP density on UI mortality did not appear to be significant.Other factors, such as modes of transportation, street planning, and pedestrian-friendly infrastructure, may play a more important role in preventing unforeseen accidents and injuries [26][27][28] .Unfortunately, these data were mostly not available at the township level.We did observe, however, a social gradient in the decline of UI mortality, i.e., the largest drop in UI mortality in low economic development townships, followed by middle and high economic areas, respectively.While improved access to AC was linked to a reduction in UI mortality, we also noticed a higher incidence of UI mortality in areas with greater access to ICU facilities.This could be attributed to a higher number of critically ill patients receiving treatment in those areas.Furthermore, the reduction in UI and all-cause deaths may serve as contributing factors to the increased life expectancy during this period.
High economic development townships exhibited the lowest CVD mortality rate compared to middle and low development townships.There was, however, a slight increase in this health indicator across all townships from 2013 to 2017.Despite many high-income countries experiencing significant reductions in CVD deaths over the past decades, Taiwan, like some other high-income nations, is witnessing a slowdown in the rate of decline.This trend can be attributed to factors such as sedentary lifestyle and the increasing prevalence of obesity, diabetes, and hypertension 29,30 .From the 1980s to the 2000s, Taiwan experienced a significant reduction in CVD mortality.Nevertheless, dietary shifts among younger generations towards higher calorie intake and sedentary lifestyles led to health concerns like overweight and metabolic syndrome.This trend is further compounded by an aging population with a higher prevalence of comorbidities, all of which collectively contribute to the rising CVD mortality rates.
Our study has other limitations.We did not consider the impact of specialist clinicians due to several reasons; first, data availability of specialists practicing at the township level was scarce and; second, calculating physician density for each area was difficult due to the possibility of a physician being registered with two or more www.nature.com/scientificreports/specialties.Also to factor in the contribution by employment and access to work, we used data on the economic development level of each township as a proxy.The economic development level was determined by ranking annual median household income of all 368 townships.We believe that annual median household income can serve as a good representation of people's employment opportunities and condition, because it reflects how families are earning; higher incomes would not only indicate access to work, but it would also suggest having jobs with better wages.Results from our correlation analyses indicated that sex ratio was positively associated with all three mortality endpoints, which suggests that areas with more males had more deaths from CVD, UI, and all-causes.A potential gender disparity could exist, but because demographics was not our main variable of interest, future investigation is needed to examine this association more closely.Lastly, although a period of four years for a longitudinal analysis may not seem extensive, we were able to demonstrate the existing geographical variations in health and their significant determinants during this period.Another longitudinal analysis may be feasible should more updated and comprehensive data becomes available in the near future.
Since the government's release of the Health Inequalities in Taiwan report in 2016 15 , which was led by Sir Michael Marmot and his team from University College London, a certain level of awareness had been raised about the health disparities among different social groups and their possible underlying factors.Unfortunately, it appears that the issue of health inequality is still being viewed through the lens of 'health care' , rather than being approached from a broader scope, like 'health in all policies' (HiAP).Currently, there are incentives to recruit more junior and government-sponsored doctors to work in mountainous and rural areas.The Integrated Delivery System provides ongoing and routine specialist outpatient services in these areas 31 .Tertiary hospitals are also subsidized for allocating critical care and trauma personnel to these areas.Other related social and public policies include providing financial support for low-income households through social welfare schemes and programs aimed at preventing chronic diseases through regular public health screening 32 .While the positive impact of PCPs is evident, targeted interventions for disadvantaged areas are essential.We must acknowledge that health is shaped by individuals' lifestyles, which are, in turn, shaped by the social and economic environments in which we live in, such as the conditions in which we are born, grow, live, work, and age 10,33 .

Conclusions
Findings of this study highlight the potential of primary care access in promoting a more equitable distribution of health.Despite ongoing efforts to improve healthcare accessibility, our findings also underscore the persistence of health disparities between individuals residing in urban and resource-deprived areas over time.Further research and targeted interventions are warranted to bridge these gaps and ensure that all populations can access quality healthcare services and achieve improved health outcomes. https://doi.org/10.1038/s41598-023-42124-ywww.nature.com/scientificreports/

Figure 1 .
Figure 1.Geographical distribution of age-standardized mortality rates (all-cause, cardiovascular disease, unintentional injury) for 2013 and 2017 in Taiwan.Overall, higher mortality rates were found in the mountainous and rural areas near the east and south regions of the island in 2013 and 2017, such as Taoyuan District in Kaohsiung and Zhuoxi Township in Hualien.On the contrary, lowest mortality rates were detected in the metropolitan areas concentrated in the west and towards the country's capital in Taipei.Note: (A) Age standardized all-cause mortality per 100,000 people.(B) Age standardized CVD mortality per 100,000 people.(C) Age standardized unintentional injury mortality per 100,000 people.

Table 1 .
Distribution of study variables for 2013 and 2017.CVD cardiovascular disease; UI unintentional injury; PCP primary care physician; AC acute care; LTC long-term care; ICU intensive-care unit.p: p-value calculated using Mann-Whitney U test or Chi-square test; Effect size: Cohen's D for continuous variables and Cramer's V for categorical variables.
Vol.:(0123456789) Scientific Reports | (2023) 13:15007 | https://doi.org/10.1038/s41598-023-42124-ywww.nature.com/scientificreports/ Population health comparison by township economic development level and year.There was significant difference across townships in all three mortality rates for 2013 and 2017.Townships with low economic development generally had the highest mortality rates.For all-cause and UI mortality, improvements were observed in the 4-year period especially for townships with low and middle economic development.

Table 4 .
Fixed effects model for population health disparities.CVD cardiovascular disease; UI unintentional injury; PCP primary care physician; AC acute care; LTC long-term care; ICU intensive-care unit; SE standard error.*Model adjusted for sex ratio and proportion of population aged 65 + and weighted by population density.