Original Research

Sociodemographic disadvantage and low healthcare access as contributing factors to high maternal death in Papua, Indonesia: multilevel analysis of census population data

AUTHORS

name here
Danik Iga Prasiska
1,2 MPH ORCID logo

name here
Suk-Yong Jang
3,4 PhD ORCID logo

name here
Min Jin Ha
5,6 PhD, Assistance Professor * ORCID logo

CORRESPONDENCE

*Dr Min Jin Ha

AFFILIATIONS

1 Department of Public Health, Graduate School, Yonsei University, Seoul, Republic of Korea

2 Department of Public Health, Regional Health Office, Nganjuk, Indonesia

3 Department of Health Policy and Management, Graduate School of Public Health, Yonsei University, Seoul, Republic of Korea

4 Department of Biohealth Industry, Policy Analysis Division, Graduate School of Transdisciplinary Health Science , Yonsei University, Seoul, Republic of Korea

5 Department of Health Informatics and Biostatistics, Graduate School of Public Health, Yonsei University, Seoul, Republic of Korea

6 Biohealth Data Science, Graduate School of Transdisciplinary Health Science, Yonsei University, Seoul, Republic of Korea

PUBLISHED

4 November 2025 Volume 25 Issue 4

HISTORY

RECEIVED: 12 July 2024

REVISED: 6 September 2025

ACCEPTED: 24 September 2025

CITATION

Prasiska D, Jang S, Ha M.  Sociodemographic disadvantage and low healthcare access as contributing factors to high maternal death in Papua, Indonesia: multilevel analysis of census population data. Rural and Remote Health 2025; 25: 9349. https://doi.org/10.22605/RRH9349

AUTHOR CONTRIBUTIONSgo to url

This work is licensed under a Creative Commons Attribution 4.0 International Licence


Abstract

Introduction: Maternal mortality remains a major public health challenge in Indonesia, with Papua province experiencing the highest provincial maternal mortality ratio. Evidence on how both individual- and community-level factors contribute to maternal deaths in this context remains limited.
Methods: This study analyzed data from the 2020 Indonesian Population Census, including 18,886 women aged 15–49 years in Papua, of whom 169 experienced maternal death. Maternal deaths were identified using the census mortality module. Descriptive statistics and χ2 tests examined bivariate associations. Multilevel logistic regression models were applied to estimate the effects of individual- and community-level determinants, accounting for district-level clustering. Sensitivity analysis using modified Poisson regression with robust variance was conducted to assess robustness of estimates
Results: Maternal mortality was significantly associated with advanced maternal age (>35 years, adjusted odds ratio (aOR) 2.92, 95% confidence interval (CI): 2.45–3.47), nulliparity (aOR 5.94, 95%CI: 4.74–7.44), unmarried status (aOR 2.79, 95%CI: 2.19–3.56), low education (aOR 2.03, 95%CI: 1.38–3.00), and rural residence (aOR 1.95, 95%CI: 1.56–2.45). None of the community-level variables were statistically significant. Substantial between-regency variation persisted in the final multilevel model (intraclass correlation coefficient = 0.61; median odds ratio = 8.95), indicating strong contextual disparities.
Conclusion: Maternal mortality in Papua is driven primarily by sociodemographic disadvantage and rural residence, while community-level health system factors alone did not explain between-regency variation. Policies should prioritize improving access and quality of care for first-time and unmarried mothers, enhancing youth-friendly health services, and addressing rural barriers through referral strengthening, maternity waiting homes, and integration of traditional birth attendants. Long-term investments in education and health workforce retention in remote areas remain critical to reducing maternal deaths in Papua.

Keywords

health inequity, health system access, maternal mortality, Papua Indonesia, primary health care, rural health services, sociodemographic factors.

Introduction

Maternal mortality is a critical global health challenge and a central indicator of health system performance. In 2023 more than 260,000 women died from pregnancy and childbirth-related causes, equivalent to nearly one death every two minutes1. Sustainable Development Goal (SDG) 3.1 aims to reduce the global maternal mortality ratio (MMR) to fewer than 70 per 100,000 live births by 20302. Although the global MMR has declined by about 40% since 2000, progress has slowed in recent years, and over 90% of maternal deaths continue to occur in low- and middle-income countries (LMICs)1. Biomedical complications such as hemorrhage, eclampsia, and sepsis are the direct causes of maternal death, but broader determinants including socioeconomic inequities, cultural practices, and weaknesses in health systems strongly influence survival3.

Indonesia reflects these global challenges. The national MMR decreased from 305 in 2015 to 189 per 100,000 live births in 2020, yet it remains among the highest in South-East Asia and far from the SDG target4,5. National averages conceal extreme disparities: the 2020 Long-Form Population Census reported provincial MMRs ranging from 48 in Jakarta to 565 in Papua5. Such inequities underscore the need to examine maternal mortality in subnational contexts.

Access to maternal health services is deeply unequal across regions. Women in Java and Bali are more than three times as likely to complete four or more antenatal visits compared with women in Papua (odds ratio (OR): 3.6, 95% confidence interval (CI): 2.7–4.7)6. Skilled birth attendance and facility-based deliveries are also low in Papua7,8. Geographic barriers further exacerbate inequities: nationally the average distance to a health facility is 5 km, but in Papua and other eastern provinces it exceeds 30 km. Nearly 28% of Papuans, predominantly in rural areas, require more than an hour to reach a health center compared with fewer than 11% in many western provinces9. These rural–urban disparities delay access to emergency obstetric care and contribute to higher maternal mortality10.

Papua represents one of the most disadvantaged provinces in Indonesia, combining high maternal mortality with remote geography, limited health infrastructure, and entrenched cultural practices that hinder care-seeking11,12. Many districts are mountainous and sparsely populated, complicating the retention of skilled health workers13 and delivery of essential maternal health services. For these reasons, Papua provides a critical setting to investigate how structural and individual factors interact to shape maternal survival.

Previous studies on maternal mortality in Indonesia have mostly emphasized biomedical complications1,3 or relied on national-level averages, without accounting for sharp structural differences across regions. Few have applied a multilevel framework to capture both individual- and community-level determinants simultaneously. This study addresses that gap by analyzing maternal deaths in Papua using the 2020 Indonesia Population Census. Particular attention is given to rural–urban residence, as it reflects systemic inequities. By integrating sociodemographic characteristics with community-level health system indicators, this study provides new evidence to inform equitable maternal health strategies in rural and remote settings.

Conceptual framework

This study adopts the maternal health and mortality determinants framework proposed by Souza et al3, which organizes risk factors across individual, community, and structural levels. The framework emphasizes that maternal deaths arise from an interaction of biological, social, and health system determinants, and that the dominant drivers shift as countries transition from high to lower mortality contexts. Applying this perspective to Papua, where maternal mortality remains high, highlights the importance of both individual sociodemographic factors and systemic barriers such as healthcare access and workforce distribution.

Methods

The study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines. Data were drawn from the 2020 Indonesian Population Census, conducted by the Indonesian Statistical Agency (Badan Pusat Statistik, BPS). Access to anonymized data is available through the BPS research service portal. Maternal deaths were identified using the census mortality module, which asked households whether, since 1 January 2017, any woman who was pregnant or had ever been pregnant had died and, if so, whether the death occurred during pregnancy or within 42 days postpartum14,15. This approach, comparable to simplified verbal autopsy methods, enables identification of maternal deaths in contexts where medical certification is limited. Figure 1 illustrates the sample selection process. From 471,794 individuals in the 2020 Papua Province Census, we excluded males, women outside the reproductive age group of 15–49 years, and women aged 15–49 years without a pregnancy record since 2017, as maternal mortality can only be assessed among those with a pregnancy during the reference period. We also removed duplicate records to avoid double-counting. After these exclusions, the final analytic sample consisted of 18,886 women aged 15–49 years with a pregnancy record.

The outcome variable was maternal death, defined according to WHO guidelines as the death of a woman while pregnant or within 42 days of termination of pregnancy, excluding accidental or incidental causes15. The variable was coded as binary (maternal death and survival). Predictor variables were chosen based on prior literature3 and census data availability. Individual-level variables included age, parity, marital status, educational attainment, and place of residence.

Age and parity capture reproductive risk, while marital status and education proxy for social and economic support. Parity was categorized into three groups: nulliparity, defined as a first pregnancy; low multiparity, defined as one to three previous births; and grand multiparity, defined as four or more previous births. Educational attainment was classified into three categories based on the highest level of schooling completed. Low education included respondents who had never attended school, did not complete primary school, or completed only primary education. Middle education comprised those who had completed junior high school senior high school or a diploma. High education referred to those with a bachelor’s degree or higher. Place of residence was classified as rural or urban using pre-coded criteria, which consider population density, infrastructure, and administrative status. Rural residence in Papua is strongly associated with longer travel distances to health facilities, limited availability of skilled health workers, and reliance on traditional practices, making it a central variable of interest.

Community-level variables included the density of health workers per 1000 population, the number of hospitals, and the number of community health centers per regency. In Indonesia, a regency (kabupaten) represents the second-level administrative division beneath a province, governing both rural and semi-urban areas. It functions at a scale broadly equivalent to a county-level jurisdiction in the US or UK, or a local government area in Australia, and is administratively parallel to an urban city (kota). Community-level data were obtained from the Indonesian Statistical Agency and are publicly accessible through its official website.

table image Figure 1: Sample selection process based on data from 2020 Indonesian Population Census.

Statistical analysis

Descriptive statistics summarized the distribution of individual- and community-level variables. Bivariate associations between maternal death and covariates were examined using χ2 tests. Considering the hierarchical nature of the data, with individuals nested within regencies, multilevel logistic regression analyses were performed using the generalized linear mixed model procedure (PROC GLIMMIX)16. The modeling was conducted sequentially: beginning with an empty (null) model estimating baseline regency-level variance in maternal mortality; followed by a model adjusting for individual-level characteristics; and culminating in a fully adjusted model incorporating both individual- and community-level determinants.

Model fit and performance were evaluated by comparing −2 log likelihood values across sequential models. To quantify the proportion of variance attributable to contextual factors, the intraclass correlation coefficient (ICC) was calculated at the regency level. Additionally, the median odds ratio (MOR) was computed to provide an interpretable measure of heterogeneity in maternal mortality risk across regencies, representing the median increase in odds when moving between regencies with differing risk profiles.

For sensitivity analysis, modified Poisson regression with robust variance estimation was conducted through PROC GLIMMIX, aligning with the main approach by modeling cluster-specific random effects. Consistency of estimates across models supported the robustness of the results. All analyses were conducted using Statistical Analysis v9.4 (SAS Institute; https://www.sas.com), with a two-sided p-value <0.05 considered statistically significant.

Ethics approval

The study was based on anonymized secondary data from the 2020 Indonesia Population Census. No personally identifiable information was used, and institutional review board exemption was obtained from the Public Health Faculty, Jember State University, Indonesia, No. 541/KEPK/FKM-UNEJ/VIII/2024. The study followed the STROBE guidelines, and a completed checklist is provided in Supplementary table 1 to support transparency and reproducibility.

Results

Table 1 presents the distribution of sociodemographic and healthcare-related factors according to maternal survival status. Maternal deaths were significantly more frequent among adolescents (<19 years, 2.1%) and women older than 35 years (0.8%) compared with those aged 19–35 years (0.5%; p<0.001). Nulliparous women experienced the highest proportion of deaths (2.3%) relative to women with low multiparity (0.4%) or grand multiparity (0.2%; p<0.001). Similarly, single or divorced women had a substantially higher risk of death (2.3%) compared with married women (0.4%, p<0.001).

Education showed a clear gradient, with deaths most common among women with low education (0.9%), compared with middle (0.5%) and high education (0.3%; p<0.001). Rural residence was also associated with higher maternal mortality (0.8%) than urban residence (0.4%, p<0.001).

At the community level, maternal deaths were more common in districts with a health worker density less than 2.5 per 1000 population (1.0%) compared with those at or greater than this threshold (0.3%, p<0.001). On average, regencies reported 2.06 hospitals (standard deviation (SD) 1.75) and 18.71 community health centers (SD 6.25).

Table 2 shows the variance components and model fit statistics from the multilevel logistic regression analysis. In the null model, the variance at the regency level was 5.50 (2.52), indicating substantial heterogeneity across districts. The corresponding ICC was 0.62, suggesting that 62% of the total variation in maternal deaths was attributable to differences between regencies. The MOR was 9.37, meaning that, on average, the odds of maternal death could increase almost 10-fold if a woman moved from a lower risk to a higher risk regency.

When individual-level covariates were added (model 1), the variance remained high (5.57), with an ICC of 0.62 and MOR of 9.51, indicating persistent contextual disparities even after adjusting for personal characteristics. In the fully adjusted model (model 2), which included both individual- and community-level factors, the variance decreased slightly to 5.27, with the ICC at 0.61 and MOR at 8.95. This suggests that community-level health system factors only modestly reduced unexplained regional variation. The model fit improved considerably from the null model (−2 log likelihood reduced from 9021.14 to 7909.43), confirming that the inclusion of covariates substantially enhanced explanatory power.

Figure 2 presents the community-level random intercepts with 95% confidence intervals, illustrating the heterogeneity in maternal mortality risk across regencies. Several regencies had significantly higher or lower baseline risks even after adjusting for individual- and community-level covariates, consistent with the high ICC and MOR values. This visualization reinforces the finding that large unexplained disparities persist between regencies, highlighting the importance of contextual factors beyond individual sociodemographic characteristics and facility availability.

Estimates for individual-level variables are reported from model 1, while contextual factors are presented from the full multilevel model 2. This presentation allows comparison between individual- and community-level effects on maternal mortality. As shown in Table 3, several individual-level characteristics were significantly associated with maternal deaths. Compared with women aged 19–35 years, those older than 35 years had almost threefold higher odds of death (aOR: 2.92, 95%CI: 2.45–3.47), while adolescents (<19 years) did not differ significantly. Nulliparous women had markedly elevated odds of death compared to women with low multiparity (aOR: 5.94, 95%CI: 4.74–7.44), whereas grand multiparity was not associated with increased risk (aOR: 0.90, 95%CI: 0.68–1.19). Single or divorced women faced almost three times higher odds of maternal death compared with married women (aOR: 2.79, 95%CI: 2.19–3.56). Educational attainment was inversely associated with maternal mortality: women with low education (aOR: 2.03, 95%CI: 1.38–3.00) and middle education (aOR: 1.58, 95%CI: 1.08–2.31) had significantly greater risk compared with those with high education. Residence also played an important role: women in rural areas were almost twice as likely to die as their urban counterparts (aOR: 1.95, 95%CI: 1.56–2.45). At the community level, however, none of the contextual variables showed statistically significant associations with maternal death. Sensitivity analyses using modified Poisson regression with robust variance produced similar estimates, confirming the robustness of the main findings (see supplementary table 3).

Table 1: Sociodemographic and healthcare-related factors of maternal deaths in Java province, Indonesia (2020 Indonesian Population Census)

Factor Characteristic Variable Total Survived (n=18,717) Died (n=169) p-value
(weighted χ2)
n n Weighted % n Weighted %
Sociodemographic Age (years) <19 1924 1875 97.9 49 2.1 <0.001
19–35 12,146 12,070 99.5 76 0.5
>35 4816 4772 99.2 44 0.8
Parity Nulliparity 3964 3853 97.7 111 2.3 <0.001
Low multiparity 10,619 10,577 99.6 42 0.4
Grand multiparity 4303 4287 99.8 16 0.2
Marital status Single/divorced 3407 3311 97.7 96 2.3 <0.001
Married 15,479 15,406 99.6 73 0.4
Education Low 9159 9054 99.1 105 0.9 <0.001
Middle 8618 8558 99.5 60 0.5
High 1109 1105 99.7 4 0.3
Residence Urban 5723 5685 99.6 38 0.4 <0.001
Rural 13,163 13,032 99.2 131 0.8
Healthcare-related Health worker density <2.5/1000 people 7712 7596 99.0 112 1.0 <0.001
≥2.5/1000 people 11,193 11,121 99.7 57 0.3
Number of hospitals: 2.06±1.75  
Number of community health centers: 18.71±6.25  

Table 2: Estimates of multilevel logistic regression model of maternal deaths in Java province, Indonesia (2020 Indonesian Population Census)

Parameter Null model Model 1 Model 2
Variance (standard error) 5.50 (2.52) 5.57 (2.24) 5.27 (2.15)
Intra-class correlation 0.62 0.62 0.61
Median odds ratio 9.37 9.51 8.95
–2 log likelihood 9021.14 7916.46 7909.43

Table 3: Assessment results of sociodemographic and healthcare-related factors associated with maternal deaths in Java province, Indonesia (2020 Indonesian Population Census)

Level Characteristic Variable Null model Model 1 Model 2
aOR (95%CI) p-value aOR (95%CI) p-value aOR (95%CI) p-value
Individual Age (years) <19 0.91 (0.74–1.12) 0.3603 0.91 (0.74–1.12) 0.3605
19–35 1.00   1.00  
>35 2.92 (2.45–3.47) <0.0001 2.92 (2.45–3.47) <0.0001
Parity Nulliparity 5.94 (4.74–7.44) <0.0001 5.94 (4.74–7.44) <0.0001
Low multiparity 1.00   1.00  
Grand multiparity 0.90 (0.68–1.19) 0.4625 0.90 (0.68–1.19) 0.4595
Marital status Single/divorced 2.79 (2.19–3.56) <0.0001 2.79 (2.19–3.55) <0.0001
Married 1.00      
Education Low 2.03 (1.38–3.00) 0.0004 2.03 (1.08–3.00) 0.0004
Middle 1.59 (1.09–2.32) 0.0172 1.58 (1.08–2.31) 0.0176
High 1.00      
Residence Rural 1.60 (1.30–2.00) <0.0001 1.95 (1.56–2.45) <0.0001
Urban 1.00      
Community Health worker density <2.5/1000   1.20 (0.15–9.34) 0.8612
    1.00  
Number of hospitals   1.66 (0.83–3.34) 0.153
Number of community health centers   1.14 (0.97–1.33) 0.1161

aOR, adjusted odds ratio. CI, confidence interval.

table image Figure 2: Caterpillar plot of community-level random intercepts from multilevel logistic regression.

Discussion

This study applied a multi-level analysis to examine both individual- and community-level factors associated with maternal death events, aiming to inform evidence-based policy strategies to reduce maternal mortality in Papua province. While the data were collected in 2020, they remain relevant for analysis of maternal mortality because the census provides the most comprehensive source of maternal death events in Papua. We also acknowledge that subsequent events, particularly the COVID-19 pandemic and recent health policy reforms, may have influenced maternal health outcomes. This limitation should be considered when interpreting the findings. Nevertheless, the 2020 Indonesian Population Census remains the only province-wide dataset capturing maternal deaths in Papua, and thus provides uniquely valuable evidence despite the time gap.

At the individual level, sociodemographic characteristics of women were significantly associated with maternal deaths. Pregnancy in women aged more than 35 years was associated with a higher risk of maternal death, whereas pregnancy in women aged less than 19 years was not statistically significant compared with those aged 19–35 years. Nevertheless, the proportion of maternal deaths among adolescents was highest (Table 1). Evidence from the US consistently shows elevated risks beginning at age 35 years17,18, likely due to chronic comorbidities and physiological changes that reduce adaptability to pregnancy, compounded by social determinants such as limited healthcare access and sociodemographic disadvantage. In contrast, the risk associated with adolescent pregnancy is more variable and context-dependent. Some studies demonstrate that very young maternal age, particularly in resource-limited settings, increases maternal mortality through biological immaturity and restricted access to care19, while others find no significant or only modestly elevated risk compared with women aged 20–35 years. This suggests a J- or U-shaped relationship between maternal age and mortality, with the lowest risk observed among women aged 20–29 years20.

Globally, WHO identifies complications in pregnancy and childbirth as the leading cause of death among girls aged 15–19 years1,21. In Indonesia, child marriage is often influenced by sociodemographic and cultural factors and remains a critical contributor to adolescent pregnancy22-24. Although the legal minimum age for marriage has recently been raised, the dispensation system still allows early marriage in practice23. Stronger policy enforcement, together with comprehensive reproductive health education and community engagement, are essential to reduce the risks associated with early marriage and adolescent pregnancy.

Nulliparous women were found to be at higher risk of maternal death compared with women who had previous childbirth experience, a finding consistent with systematic reviews from LMICs25. In our analytic sample, nulliparity accounted for approximately 21% of pregnancies, underscoring the importance of first pregnancies as a vulnerable group in Papua. Nationally, the 2017 Indonesian Demographic and Health Survey reported that about 7% of girls aged 15–19 years had already begun childbearing, and the prevalence is higher in rural and eastern provinces, including Papua26. These patterns indicate that both adolescent and first pregnancies remain important contributors to maternal risk in the region. Evidence from China indicates that advanced maternal age elevates the risk of gestational hypertension and pre-eclampsia, particularly among nulliparous women27.

Marital status also emerged as an important determinant. Women without partners experienced a substantially higher risk of maternal death compared with married women. Indonesian evidence shows that unmarried women and adolescents with unintended pregnancies are less likely to complete the recommended four antenatal care visits, reflecting both stigma and barriers within the health system. Although programs such as the youth care health program Pelayanan Kesehatan Peduli Remaja (Adolescent-Friendly Health Services) are designed to provide confidential and supportive reproductive health services, coverage and utilization remain limited in Papua28. This is consistent with findings from the US, where married women demonstrated better maternal outcomes, partly due to greater social and financial support29. In Malaysia, unmarried women were more likely to underutilize prenatal services, engage in unhealthy behaviors such as smoking and alcohol consumption, and report a higher prevalence of depressive symptoms30. Such factors increase vulnerability during pregnancy and childbirth. Collectively, these findings underscore the importance of supportive programs that enable both nulliparous and unmarried women to access timely and high-quality prenatal care without stigma.

Education and place of residence were also significant contributors to maternal death. Women with lower levels of education and those living in rural areas faced a higher risk, consistent with findings from previous studies 3,25,31,32. Higher educational attainment is known to enhance reproductive decision-making, promote positive health-seeking behaviors, strengthen women’s autonomy, and facilitate the timely use of maternal health services33. In contrast, rural residence often reflects structural barriers, including limited availability of healthcare facilities, slower adoption of medical technologies, and greater travel distances to reach obstetric care, all of which contribute to delays in receiving appropriate treatment34,35.

At the community level, more than half of the variability in maternal deaths was attributable to contextual factors, underscoring the importance of applying a multilevel analytical approach. This finding should be interpreted with caution. Health workforce availability emerged as an important factor. Health worker density below WHO’s recommended threshold of 2.5 per 1000 population was not statistically associated with higher maternal mortality in this study. This finding is consistent with a previous published study in Indonesia36. Previous research from Sub-Saharan Africa and other settings has demonstrated that greater health workforce density correlates with improved maternal health outcomes37,38. Retention of healthcare workers in remote and rural areas of Indonesia remains a national challenge39, particularly in Papua, where difficult geographical access and limited career opportunities hinder long-term retention. Strengthening incentives, professional development, and community-based deployment strategies will be essential to improve health worker availability in these underserved regions.

Although the number of hospitals and community health centers in a regency was included as community-level variables, neither showed a statistically significant association with maternal mortality. This finding suggests that the mere availability of facilities may not be sufficient to reduce maternal deaths. Similar results have been reported in other low-resource settings, where hospital proximity or distance alone was not necessarily associated with maternal mortality40. Several explanations are possible. First, facility presence does not guarantee functionality, quality of care, or readiness to manage obstetric emergencies41. Many rural health centers in Indonesia lack trained personnel, essential supplies, and reliable referral systems42,43, which limits their ability to prevent maternal deaths even when facilities are physically available. Second, geographic barriers remain substantial, and women living in remote areas may still face long travel times to reach hospitals, particularly in emergency situations44,45. Third, under-reporting and misclassification of deaths outside health facilities may distort the apparent association46,47. As reporting systems improve and more deaths are captured within facilities, stronger associations may emerge.

Our results should also be interpreted within the specific context of Papua. Previous work has highlighted that maternal health services in Papua are constrained by systemic challenges, including shortages of skilled personnel, weak referral networks, and the difficulty of delivering services across mountainous and remote areas48. These constraints may explain why the number of hospitals and community health centers was not significantly associated with maternal deaths in our analysis: the presence of facilities does not guarantee functionality, quality of care, or timely access in emergencies. In addition, qualitative research from Jayawijaya district, Papua, found that recognition of maternal and newborn complications was often delayed, and that cultural norms, reliance on traditional birth attendants, and transport barriers limited care-seeking even when services were available49. Together, these findings suggest that reducing maternal mortality in Papua requires strategies that go beyond increasing facility counts to address both structural and sociocultural barriers. Practical approaches include improving the integration of traditional birth attendants into referral pathways, strengthening emergency transport systems, and expanding maternity waiting homes in remote districts as an interim strategy to improve timely access to life-saving care.

Strengths and limitations

This study is the first to examine the combined effects of individual- and community-level healthcare-related factors on maternal deaths in Indonesia. A major strength lies in the use of census data, which provides nationwide coverage and captures a large number of maternal death events across diverse populations, thereby enhancing the representativeness and generalizability of the findings.

Several limitations should also be acknowledged. First, under-reporting of maternal deaths by census respondents is possible, which may have led to an underestimation of the true burden. Second, the census relied on retrospective reporting of deaths that occurred within the previous 3 years, raising the possibility of recall bias, especially for sensitive events such as maternal deaths. Third, the data did not capture pregnancies that did not progress to a live birth, thereby excluding an important subgroup of high-risk women. For example, women with a history of miscarriage, abortion, or stillbirth often face elevated risks in subsequent pregnancies, yet these events were not recorded in the census.

In addition, the census lacked detailed biomedical and clinical variables, which limited the ability to identify the direct medical causes of maternal deaths including antenatal or postnatal care utilization. Important contextual factors such as distance to the nearest health facility, travel time, transportation barriers, number of private clinics, presence of a traditional healer, and referral pathways were also unavailable. A traditional healer is a non-medical practitioner who relies on herbal, spiritual, and belief-based practices, and is distinct from a traditional birth attendant, who assists with childbirth, although both are non-medical providers. These factors are well-recognized determinants of maternal survival, particularly in geographically remote areas such as Papua, and their absence may have led to underestimation of health system barriers. Consequently, our analysis may not fully reflect the role of health system performance in shaping maternal outcomes.

Future research should address these gaps by linking census data with geospatial information, health facility assessments, and biomedical data. Such integrated approaches would provide a more comprehensive understanding of how sociodemographic, cultural, and health system factors interact to influence maternal mortality in Indonesia.

Conclusion

This study demonstrates that maternal mortality in Papua province is shaped by a combination of individual and community factors. Advanced maternal age, nulliparity, unmarried status, low education, and rural residence were all strongly associated with higher risk of maternal death. Although community-level health system indicators such as health worker density and the number of facilities were not significantly associated with maternal deaths, the large contextual variation observed between districts highlights the influence of systemic and geographic barriers.

The findings underscore that reducing maternal mortality in Papua requires interventions that go beyond increasing the number of health facilities or workers. Practical strategies should prioritize improving access and quality of care for vulnerable groups, particularly first-time and unmarried mothers, through adolescent- and youth-friendly health services. For rural and remote districts, strengthening emergency referral systems, investing in maternity waiting homes near referral hospitals, and integrating traditional birth attendants into formal referral pathways are feasible measures that can provide immediate impact. Long-term efforts to retain health workers in remote areas and address structural inequities in education and health access remain essential.

By situating individual risks within the broader rural and remote context of Papua, this study highlights the importance of targeted, context-sensitive strategies to accelerate progress toward reducing maternal mortality and achieving national and global goals.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Data availability

The Indonesian Statistical Agency (Badan Pusat Statistik, BPS) provided consent for their data to be used for research and analysis purposes and publish the findings. BPS-provided anonymous census data is available for the public by request at https://pst.bps.go.id.

Conflicts of interest

The authors declare that they have no relevant financial or non-financial interest to disclose.

References

1 World Health Organization. Maternal mortality. Available: https://www.who.int/news-room/fact-sheets/detail/maternal-mortalityweb link (Accessed September 2025).
2 United Nations Department of Economic and Social Affairs. Goals: 3 Ensure healthy lives and promote well-being for all at all ages. Available: https://sdgs.un.org/goals/goal3#targets_and_indicatorsweb link (Accessed September 2025).
3 Souza JP, Day LT, Rezende-Gomes AC, Zhang J, Mori R, Baguiya A, et al. A global analysis of the determinants of maternal health and transitions in maternal mortality. The Lancet Global Health 2024; 12(2): e306-e316. https://doi.org/10.1016/S2214-109X(23)00468-0DOI, https://www.ncbi.nlm.nih.gov/pubmed/38070536\PMid:38070536
4 Statistics Indonesia (BPS). Profile of the Indonesian Population: Results of the 2015 Intercensal Population Survey. [In Indonesian]. Jakarta: Statistics Indonesia (BPS), 2016.
5 Statistics Indonesia (BPS). Mortality in Indonesia: Results of the 2020 Population Census Long Form. [In Indonesian]. Jakarta: Statistics Indonesia (BPS), 2024.
6 Laksono AD, Rukmini R, Wulandari RD. Regional disparities in antenatal care utilization in Indonesia. PLoS ONE 2020; 15(2): e0224006. https://doi.org/10.1371/journal.pone.0224006DOI, https://www.ncbi.nlm.nih.gov/pubmed/32053621\PMid:32053621
7 Laksono AD, Wulandari RD. Regional disparities of facility-based childbirth in Indonesia. Trends in Sciences 2021; 18(21): 387. https://doi.org/10.48048/tis.2021.387DOI
8 UNICEF Indonesia. Improving maternal and newborn health services in eastern Indonesia. Findings from an external review. Jakarta: UNICEF, 2017.
9 Mahendradhata Y, Trisnantoro L, Listyadewi S, Soewondo P, Marthias T, Harimurti P, et al. The Republic of Indonesia health system review, vol. 7. New Delhi: World Health Organization, 2017.
10 Rajkotia Y, Gergen J, Djurovic I, Koseki S, Coe M, Jobarteh K, et al. Re-envisioning maternal and newborn health in Indonesia: How the private sector and civil society can ignite change. Washington DC: USAID, 2016.
11 Wahyudi R, Sondakh F, Permatasari, VNAR. Social protection inclusion in Indonesia's remote areas: Identifying and addressing gaps and challenges. 2024. Available: https://documents1.worldbank.org/curated/en/099062624071020814/pdf/P177341-b77c6fc9-f2db-4e5c-b24f-c41e1ce945ca.pdfweb link (Accessed 15 October 2025).
12 Wahyuni RN, Ikhsan M, Damayanti A, Khoirunurrofik K. Inter-district road infrastructure and spatial inequality in rural Indonesia. Economies 2022; 10(9): 229. https://doi.org/10.3390/economies10090229DOI
13 Muharram FR, Sulistya HA, Swannjo JB, Firmansyah FF, Rizal MM, Izza A, et al. Adequacy and distribution of the health workforce in Indonesia. WHO South-East Asia Journal of Public Health 2024; 13(2): 45-55. https://doi.org/10.4103/WHO-SEAJPH.WHO-SEAJPH_28_24DOI, https://www.ncbi.nlm.nih.gov/pubmed/39995001\PMid:39995001
14 World Health Organization. WHO guidance for measuring maternal mortality from a census. Geneva: World Health Organization, 2013.
15 World Health Organization. International classification of diseases. 11th revision. 2018. Available: https://www.who.int/standards/classifications/classification-of-diseasesweb link (Accessed 15 October 2025).
16 SAS Institute Inc. The GLIMMIX Procedure. SAS/STAT® 15.2 User’s Guide. 2022. Available: https://documentation.sas.com/doc/en/statug/15.2/statug_glimmix_toc.htmweb link (Accessed 21 October 2025).
17 MacDorman MF, Thoma M, Declercq E, Howell EA. Causes contributing to the excess maternal mortality risk for women 35 and over, United States, 2016-2017. PLoS ONE 2021; 16(6): e0253920. https://doi.org/10.1371/journal.pone.0253920DOI, https://www.ncbi.nlm.nih.gov/pubmed/34185810\PMid:34185810
18 Huang RS, Spence AR, Abenhaim HA. Age-related disparities in national maternal mortality trends: A population-based study. PLoS ONE 2025; 20(1): e0316578. https://doi.org/10.1371/journal.pone.0316578DOI, https://www.ncbi.nlm.nih.gov/pubmed/39752379\PMid:39752379
19 Mayor S. Pregnancy and childbirth are leading causes of death in teenage girls in developing countries. BMJ 2004; 328(7449): 1152. https://doi.org/10.1136/bmj.328.7449.1152-aDOI, https://www.ncbi.nlm.nih.gov/pubmed/15142897\PMid:15142897
20 Nove A, Matthews Z, Neal S, Camacho AV. Maternal mortality in adolescents compared with women of other ages: evidence from 144 countries. The Lancet Global Health 2014; 2(3): e155-e164. https://doi.org/10.1016/S2214-109X(13)70179-7DOI, https://www.ncbi.nlm.nih.gov/pubmed/25102848\PMid:25102848
21 World Health Organization. Maternal mortality: evidence brief. Geneva: World Health Organization, 2019.
22 Ratnaningsih M, Wibowo HR, Goodwin NJ, Rezki AAKS, Ridwan R, Hadyani RN, et al. Child Marriage Acceptability Index (CMAI) as an essential indicator: an investigation in South and Central Sulawesi, Indonesia. Global Health Research and Policy 2022; 7(1): 32. https://doi.org/10.1186/s41256-022-00252-4DOI, https://www.ncbi.nlm.nih.gov/pubmed/36163286\PMid:36163286
23 Grijns M, Horii H. Child marriage in a village in West Java (Indonesia): Compromises between legal obligations and religious concerns. Asian Journal of Law and Society 2018; 5(2): 453-466. https://doi.org/10.1017/als.2018.9DOI
24 Judiasih SD, Rubiati B, Yuanitasari D, Salim EF, Safira L. Efforts to eradicate child marriage practices in Indonesia: Towards Sustainable Development Goals. Journal of International Women's Studies 2020; 21(6): 135-149.
25 Symonds NE, Vidler M, Wiens MO, Omar S, English LL, Ukah UV, et al. Risk factors for postpartum maternal mortality and hospital readmission in low- and middle-income countries: a systematic review. BMC Pregnancy and Childbirth 2023; 23(1): 303. https://doi.org/10.1186/s12884-023-05459-yDOI, https://www.ncbi.nlm.nih.gov/pubmed/37120529\PMid:37120529
26 Devi YP, Ekoriano M, Sari DP, Muthmainnah M. Factors associated with adolescent birth in Indonesia: a national survey. Rural and Remote Health 2022; 22(2): 7219. https://doi.org/10.22605/RRH7219DOI, https://www.ncbi.nlm.nih.gov/pubmed/35598832\PMid:35598832
27 Luo J, Fan C, Luo M, Fang J, Zhou S, Zhang F. Pregnancy complications among nulliparous and multiparous women with advanced maternal age: a community-based prospective cohort study in China. BMC Pregnancy and Childbirth 2020; 20(1): 581. https://doi.org/10.1186/s12884-020-03284-1DOI, https://www.ncbi.nlm.nih.gov/pubmed/33008331\PMid:33008331
28 Andriana A, Kamil I, Budiarty E. Adolescent-friendly health services in Jakarta and Makassar: A gender and cultural analysis. 2019. Available: https://www.alignplatform.org/resources/adolescent-friendly-health-services-jakarta-and-makassar-gender-and-cultural-analysisweb link (Accessed 15 October 2025).
29 Fafard St-Germain AA, Kirby RS, Urquia ML. Reproductive health among married and unmarried mothers aged less than 18, 18-19, and 20-24 years in the United States, 2014-2019: A population-based cross-sectional study. PLoS Medicine 2022; 19(3): e1003929. https://doi.org/10.1371/journal.pmed.1003929DOI, https://www.ncbi.nlm.nih.gov/pubmed/35271581\PMid:35271581
30 Mohd Zain N, Low WY, Othman S. Impact of maternal marital status on birth outcomes among young Malaysian women: a prospective cohort study. Asia Pacific Journal of Public Health 2015; 27(3): 335-347. https://doi.org/10.1177/1010539514537678DOI, https://www.ncbi.nlm.nih.gov/pubmed/25005933\PMid:25005933
31 Tajvar M, Hajizadeh A, Zalvand R. A systematic review of individual and ecological determinants of maternal mortality in the world based on the income level of countries. BMC Public Health 2022; 22(1): 2354. https://doi.org/10.1186/s12889-022-14686-5DOI, https://www.ncbi.nlm.nih.gov/pubmed/36522731\PMid:36522731
32 Khalili M, Mashrouteh M, Haghdoost A, Torkian S, Chegeni M. Systematic review and meta-analysis of maternal mortality ratio and related factors in the Islamic Republic of Iran. Eastern Mediterranean Health Journal 2023; 29(5): 380-401. https://doi.org/10.26719/emhj.23.063DOI, https://www.ncbi.nlm.nih.gov/pubmed/37306175\PMid:37306175
33 Yakubu Y, Mohamed Nor N, Abidin EZ. A systematic review of micro correlates of maternal mortality. Reviews on Environmental Health 2018; 33(2): 147-161. https://doi.org/10.1515/reveh-2017-0050DOI, https://www.ncbi.nlm.nih.gov/pubmed/29729149\PMid:29729149
34 Harrington KA, Cameron NA, Culler K, Grobman WA, Khan SS. Rural-urban disparities in adverse maternal outcomes in the United States, 2016-2019. American Journal of Public Health 2023; 113(2): 224-227. https://doi.org/10.2105/AJPH.2022.307134DOI, https://www.ncbi.nlm.nih.gov/pubmed/36652639\PMid:36652639
35 Kozhimannil KB, Interrante JD, Henning-Smith C, Admon LK. Rural-urban differences in severe maternal morbidity and mortality in the US, 2007-15. Health Affairs (Millwood) 2019; 38(12): 2077-2085. https://doi.org/10.1377/hlthaff.2019.00805DOI, https://www.ncbi.nlm.nih.gov/pubmed/31794322\PMid:31794322
36 Utomo B, Romadlona NA, Naviandi U, BaharuddinNur RJ, Makalew R, Liyanto E, et al. Census block based loglinear regression analysis of health and social determinants of maternal mortality in Indonesia 2010-2021. Scientific Reports 2025; 15(1): 9397. https://doi.org/10.1038/s41598-025-91942-9DOI, https://www.ncbi.nlm.nih.gov/pubmed/40102491\PMid:40102491
37 Rosser JI, Aluri KZ, Kempinsky A, Richardson S, Bendavid E. The effect of healthcare worker density on maternal health service utilization in Sub-Saharan Africa. American Journal of Tropical Medicine and Hygiene 2022; 106(3): 939-944. https://doi.org/10.4269/ajtmh.21-0727DOI, https://www.ncbi.nlm.nih.gov/pubmed/35026729\PMid:35026729
38 Liu J, Eggleston K. The association between health workforce and health outcomes: a cross-country econometric study. Social Indicators Research 2022; 163(2): 609-632. https://doi.org/10.1007/s11205-022-02910-zDOI, https://www.ncbi.nlm.nih.gov/pubmed/35310535\PMid:35310535
39 Kurniati A, Efendi F, Ismawiningsih I, Mulyani N, Zakaria Z, Ambarwati R, et al. Retention of doctors and dentists to serve in remote areas in Indonesia: A discrete choice experiment. Journal of Multidisciplinary Healthcare 2024; 17: 2215-2225. https://doi.org/10.2147/JMDH.S459158DOI, https://www.ncbi.nlm.nih.gov/pubmed/38741922\PMid:38741922
40 Hanson C, Cox J, Mbaruku G, Manzi F, Gabrysch S, Schellenberg D, et al. Maternal mortality and distance to facility-based obstetric care in rural southern Tanzania: a secondary analysis of cross-sectional census data in 226 000 households. The Lancet Global Health 2015; 3(7): e387-e395. https://doi.org/10.1016/S2214-109X(15)00048-0DOI, https://www.ncbi.nlm.nih.gov/pubmed/26004775\PMid:26004775
41 Mweemba C, Mapulanga M, Jacobs C, Katowa-Mukwato P, Maimbolwa M. Access barriers to maternal healthcare services in selected hard-to-reach areas of Zambia: a mixed methods design. Pan African Medical Journal 2021; 40: 4. https://doi.org/10.11604/pamj.2021.40.4.28423DOI, https://www.ncbi.nlm.nih.gov/pubmed/34650654\PMid:34650654
42 Rudiyanti N, Utomo B. Challenges of health workers in primary health facilities in implementing obstetric emergency referrals to save women from death in Indonesia: A qualitative study. Belitung Nursing Journal 2024; 10(6): 644-653. https://doi.org/10.33546/bnj.3525DOI, https://www.ncbi.nlm.nih.gov/pubmed/39601029\PMid:39601029
43 Mahmood MA, Hendarto H, Laksana MAC, Damayanti HE, Suhargono MH, Pranadyan R, et al. Health system and quality of care factors contributing to maternal deaths in East Java, Indonesia. PLoS One 2021; 16(2): e0247911. https://doi.org/10.1371/journal.pone.0247911DOI, https://www.ncbi.nlm.nih.gov/pubmed/33635928\PMid:33635928
44 Rizkianti A, Saptarini I, Rachmalina R. Perceived barriers in accessing health care and the risk of pregnancy complications in Indonesia. International Journal of Women’s Health 2021; 13: 761-772. https://doi.org/10.2147/IJWH.S310850DOI, https://www.ncbi.nlm.nih.gov/pubmed/34429661\PMid:34429661
45 Rambu Ngana F, Eka Karyawati AAIN. Scenario modelling as planning evidence to improve access to emergency obstetric care in eastern Indonesia. PLoS ONE 2021; 16(6): e0251869. https://doi.org/10.1371/journal.pone.0251869DOI, https://www.ncbi.nlm.nih.gov/pubmed/34106942\PMid:34106942
46 Ahmed SMA, Cresswell JA, Say L. Incompleteness and misclassification of maternal death recording: a systematic review and meta-analysis. BMC Pregnancy and Childbirth 2023; 23(1): 794. https://doi.org/10.1186/s12884-023-06077-4DOI, https://www.ncbi.nlm.nih.gov/pubmed/37968585\PMid:37968585
47 Qomariyah SN, Sethi R, Izati YN, Rianty T, Latief K, Zazri A, et al. No one data source captures all: A nested case-control study of the completeness of maternal death reporting in Banten Province, Indonesia. PLoS ONE 2020; 15(5): e0232080. https://doi.org/10.1371/journal.pone.0232080DOI, https://www.ncbi.nlm.nih.gov/pubmed/32379774\PMid:32379774
48 Marthias T, Trisnantoro L. Health system for maternal health - a case study from Papua, Indonesia. BMC Public Health 2012; 12(Suppl 2): A24. https://doi.org/10.1186/1471-2458-12-S2-A24DOI
49 Rosales A, Sulistyo S, Miko O, Hairani LK, Ilyana M, Thomas J, et al. Recognition of and care-seeking for maternal and newborn complications in Jayawijaya district, Papua province, Indonesia: a qualitative study. Journal of Health, Population and Nutrition 2017; 36(1): 44. https://doi.org/10.1186/s41043-017-0122-0DOI, https://www.ncbi.nlm.nih.gov/pubmed/29297380\PMid:29297380