Rural-urban differentials in the determinants of under-five mortality in Bhutan

Tashi Dendup (School of Health and Society, Faculty of Arts, Social Sciences and Humanities, University of Wollongong, Wollongong, Australia)
Yun Zhao (School of Public Health, Faculty of Health Sciences, Curtin University, Perth, Australia)
I Gusti Ngurah Edi Putra (School of Health and Society, Faculty of Arts, Social Sciences and Humanities, University of Wollongong, Wollongong, Australia)

Journal of Health Research

ISSN: 2586-940X

Article publication date: 25 September 2020

Issue publication date: 20 April 2021



The differences in the distribution of factors associated with under-five mortality (UFM) can help explain the rural-urban inequities in UFM. The determinants contributing to UFM in rural and urban areas have not been previously explored in Bhutan. This study examined the factors associated with UFM in rural and urban Bhutan and the role of the factors in explaining UFM disparity.


The dataset of 6,398 single births (4,999 in rural and 1,399 in urban areas) from the 2012 Bhutan National Health Survey was analyzed. Logistic regression analysis accounting for the complex survey design was performed to investigate the determinants.


The UFM rate was 2.75 times higher in rural than in urban Bhutan. In rural communities, children of younger mothers, born in households without safe sanitation and electricity, and central and eastern regions had increased UFM odds. Whereas, children born to working mothers and educated fathers, and born in households with non-working household heads had lower UFM odds in urban areas. A higher number of births and smaller household size was associated with an increased UFM odds irrespective of rural-urban residence. Environmental factors were attributable for the largest portion of rural UFM disadvantage.


This study helps to understand the rural-urban differences in the factors influencing UFM in Bhutan. The findings suggest that policies aimed to improve environmental and socioeconomic conditions, women empowerment, and those aimed to enhance health utilization can help reduce the rural-urban child survival disparity and accelerate the achievement of the Sustainable Development Goal target.



Dendup, T., Zhao, Y. and Putra, I.G.N.E. (2021), "Rural-urban differentials in the determinants of under-five mortality in Bhutan", Journal of Health Research, Vol. 35 No. 3, pp. 226-239.



Emerald Publishing Limited

Copyright © 2020, Tashi Dendup, Yun Zhao and I Gusti Ngurah Edi Putra


Published in Journal of Health Research. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at


Under-five mortality (UFM) remains a pressing public health issue in developing countries. Globally, 5.6 million children died before reaching the age of 5 years in 2016, of which countries in the sub-Saharan and South Asian regions contributed to more than 80% of the deaths [1]. Given its close association with poverty, education, healthcare, nutrition, and the environment, child mortality is not only a good marker of children's health but also reflects a nation's wellbeing [1]. The under-five mortality rate (UFMR) has been considered as an important socioeconomic measure under the Sustainable Development Goals (SGDs) with an ambitious target to reduce it by ≤25 deaths per 1,000 live births by 2030 [2]. Despite the significant progress made by many countries in improving child mortality, child survival differentials by place of residence persist [3].

The risk of UFM is generally higher in rural than urban areas in developing nations [410]. Higher child deaths in rural areas can indicate an uneven distribution of adequate public resources and services. Some studies attribute low household income and wealth, mother's education, contraceptive use and vaccination to rural disadvantage in infant and child mortality [68, 11, 12]. Likewise, poor access to safe drinking water, electricity, health services and poor housing conditions contributed to high child mortality in rural settings [4, 9]. In contrast, better community factors besides socioeconomic and behavioral characteristics were shown to be the reasons for urban advantage in child survival [7]. A study also identified household wealth in addition to the child's gender and maternal obesity in urban areas, and region of residence and mother's marital status in rural areas as risk factors for infant mortality [5]. Rural-urban social and resource inequalities may lead to different behaviors among mothers and families that can potentially contribute to UFM disparities [3]. Differences in the distribution of factors influencing child mortality have been purported to explain the gap in mortality between urban and rural areas [4, 13, 14].

Bhutan is a small, landlocked Himalayan country in South Asia known for measuring its wellbeing and development using the Gross National Happiness Index. In Bhutan, health service is accessible free of cost and is provided by the State through its three-tiered health care system built upon primary health care principles. The UFMR declined from 134 in 1990 to 37 per 1,000 live births in 2012 [15], and further to 32 per 1,000 live births in 2016 [16]. The fertility rate has also declined from 5.6 in 1994 to 2.1 in 2012 [15], and further down to 1.9 in 2017 [17] reflecting the commendable socio-economic and health progress that Bhutan has achieved over the last two decades. However, the current UFMR is higher compared to some other developing countries in the region, such as Indonesia (26 per 1,000), Sri Lanka and Maldives (9 per 1,000 each), and Thailand (12 per 1,000) [16]. Bhutan is currently positioned 66th out of more than 192 countries in terms of global UFMR when ranked from highest to lowest [16]. Targeted interventions to improve child survival may not only help Bhutan achieve a UFMR level beyond that of the SDG target but may also positively impact other related socioeconomic indicators.

Around 65% of the population live in rural Bhutan and most of the population practice sustenance farming [18]. Rural Bhutan is characterized by scattered settlements, rugged terrain, inadequate public services including health and motorable roads, stronger traditional beliefs, and a poor socioeconomic and environmental conditions. These characteristics place rural residents, especially women and children, at a much higher risk of poor health. Socioeconomic transitions have also led to a rising trend of migration and urbanization, lifestyle changes, and health issues. Bhutan's projected urban population growth rate was among the highest in South Asia [19], which is attributable to the rapid rural-urban migration resulting in an increasing underserved migrant population [20]. The growing population can exert pressure on the inadequate services and the environment, which can affect health, especially that of women and children. Seasonal shortages and pressure on the already limited water supply in urban areas are being reported [21]. Although health care is free, stark rural-urban health disparities remain [15, 20, 22, 23]. Studies investigating specific factors influencing health in rural and urban areas are scarce.

Our previous study found mother's age, household size, parity, access to electricity and sanitation, and the region of residence were the predictors of UFM [24]. National surveys show higher rates of UFM in rural than urban Bhutan portraying a significant rural-urban health inequity [15, 25]. The factors contributing to this disparity have not been previously explored. Interventions directed to reduce child mortality in rural areas may yield a greater impact on the overall UFM, given the high UFM in rural settings. Such focused efforts to improve health and narrow disparities by understanding the health needs of the specific population is essential for a resource-constrained country like Bhutan. Therefore, this study examined the determinants of UFM in rural and urban Bhutan using the most current dataset of the Bhutan National Health Survey (BNHS) [15]. The potential factors contributing to the rural-urban differences in UFM were also explored.


Study data

This study used the datasets of the Bhutan National Health Survey (BNHS) that was conducted jointly by the Ministry of Health and the National Statistical Bureau of Bhutan in 2012 [15]. The survey aimed to generate health indicators and to observe the trend of health status. A two-stage sampling design that defined rural and urban areas in each of the 20 districts as sampling strata was used. The data were collected using five interviewer-administered questionnaire sets, namely the household, individual, women, immunization and domestic violence from November 2012 until February 2013. This survey successfully interviewed 13,256 (97%) households and 45,635 eligible individuals (87%). The response rate was ≥90% for both the domestic violence and questionnaires targeting women in the reproductive age range of 10-49 years. The study report of the 2012 BNHS has been previously published [15].

All the data files, such as household, women, birth history, immunization and domestic violence, were merged. In maintaining consistency with the national report, all live births in the last five years or more prior to the survey (November 2002 to October 2007) were extracted. This allowed for the youngest children to attain their first five years at the beginning of the BNHS in November 2012 [15]. For this study, only single births were selected to prevent potential confounding due to multiple pregnancies. After excluding eight observations without mortality information, the final sample used in the present study was 6,389.

Study variables

Under-five death was defined as a death that occurred before the children attained their fifth birthday. To enable logistic regression analysis, under-five death was coded as “1” and those who survived as “0”. The selection of explanatory variables was made by reviewing the current literature and considering the availability of variables in the dataset. Mosley and Chen's framework for analyzing UFM determinants in developing countries was used to group the explanatory variables into three categories [26]. Parent's education, wealth index that was created using the principal components analysis [15], household head and mother's working status, household size, and domestic violence were grouped as the socioeconomic determinants. The bio-demographic variables included the current age of mother, age when mother got first married and pregnant, marital status, birth order, birth interval, total births, the household head's sex and the child's sex. The World Health Organization's recommendation of a minimum of 33 months between two consecutive births was used to define birth interval [27]. The birth interval was also assessed as a continuous variable and using other thresholds (≤36 and >36 months, and ≤24 and >24 months). The region of residence, drinking water source, access to safe sanitation, the use of solid fuel for cooking, and the availability of electricity comprised the environmental determinants.

Statistical analysis

Cross-tabulations assessed the distribution of UFM by the explanatory variables. The unadjusted or crude odds ratio (COR) along with the 95% confidence intervals (CIs) were generated using simple logistic regression (SLR) that examined the association of each explanatory variable with the likelihood of UFM.

Determinants of UFM were then identified by using multiple logistic regression (MLR) analysis with a stepwise hierarchical method [28]. Varying significance levels with a backward elimination method were applied to minimize the risk of missing important determinants [29]. Initially, socioeconomic variables with a p-value of ≤0.2 in the SLR were included in an MLR model to assess their association with UFM, and variables associated with UFM at a 10% level (p < 0.1) were retained for consequent modeling. In the next model, bio-demographic variables with a p-value of ≤0.2 in the SLR were included and assessed their association with UFM in the presence of the retained socioeconomic variables. The variables with a p-value of <0.1 were kept again in the model. The final model combined the socioeconomic and bio-demographic variables retained in the preceding models with the environmental variables having a p-value ≤0.2 in SLR. The final model retained only those variables that were significant at 5% (p < 0.05). The adjusted odds ratio (AOR) and their 95% CIs were reported.

The above analyses were conducted separately for rural and urban areas. Some variables that are found to significantly predict child mortality in the literature, such as age and education of mother, sex of the child, wealth index, and use of solid fuel, were again included in the final model to reexamine their influence. Furthermore, using the full sample, separate models that adjusted for each group of variables were also fitted to assess their role in explaining the effect of place of residence (rural-urban areas) on UFM. For this, only those variables found significant in the bivariate analysis in each group were considered. The complex-sample analysis was performed in the data analysis to take into account the sampling weight and clustering effects emanating from the multistage sampling. All analyses were conducted using the STATA version 14 package [30].

Ethical approval

This study obtained ethical approval from the Bhutan Research Ethics Board of Health (Approval number: REBH/Approval/2016/031).


The total number of single births in the five-year period (November 2002 to October 2007) was 1,399 in urban and 4,990 in rural areas. Table 1 presents the distribution of UFM by determinants. The overall weighted UFMR was 37 per 1,000 live births, and the weighted UFMR in rural and urban areas was 44 per 1,000 and 16 per 1,000 live births, respectively. In urban areas, the UFMR was higher among boys, second child, among children born with a birth interval of ≥33 months, born to mothers aged 29–34 years, who married before 18 years of age, who were married, who experienced domestic violence and had >2 children. The prevalence was also higher among children born in male-headed households, born to uneducated parents and unemployed mothers, born in households with <5 members, and uneducated and employed household heads. The UFMR was also greater among children born in central Bhutan, and in households with access to safe drinking water, sanitation, electricity, and that used solid fuel.

The rates of UFM in rural areas was similar to urban areas for most of the reported characteristics. In contrast, the prevalence of UFM in rural areas was higher among third or higher-order children, children born to mothers who got pregnant before 20 years old and did not experience violence, and among children born with a birth interval of <33 months. The UFMR was greater among children born in households with educated and unemployed household heads and born in poor to middle-income households. The UFMR was also higher among those born in the eastern region, and born in households with safe drinking water, without safe sanitation and electricity, and that utilized solid fuel.

Table 2 provides the unadjusted odds ratio for rural and urban areas. Under-five children born to mothers with >2 children had a higher likelihood of UFM regardless of the residential area. Specifically, in urban areas, children born to educated mothers and fathers, working mothers, born in households headed by females, with educated and non-working household heads had significantly reduced odds of mortality. In rural areas, under-five children born in wealthier households had reduced odds while children living in eastern and central regions, born in households without safe sanitation and electricity, and that used solid fuels for cooking had a significantly greater likelihood of UFM.

The adjusted odds ratio estimated in the final model is presented in Figures 1 and 2 for urban and rural areas, respectively. Compared to bivariate analysis in Table 2, mother's education, and household head's sex, education, and working status in urban areas, and wealth index and solid fuel use in rural areas turned out to be not statistically significant in the multivariate analysis. However, the mother's age in rural areas and household size in urban areas emerged significant. Regardless of rural-urban residence, under-five children in households with >4 members had significantly reduced odds of UFM (by 91% in urban and 65% in rural areas), while those born to mothers who have >2 children had higher UFM odds. Independently, in urban areas, the likelihood of mortality among children of educated fathers was lower by 70% (AOR = 0.30, p = 0.043) than their counterparts. Under-five children of working mothers also had decreased odds of dying by 79% (AOR = 0.21, p = 0.016). Whereas, children in households with non-working household heads had lower odds of dying (AOR = 0.20, p = 0.032). The pseudo-R2 of the final urban model was 0.20.

Specifically, in rural areas (Figure 2), children born to mothers who were 29–34 (AOR = 0.77, p = 0.092) and >34 (AOR = 0.60, p = 0.002) years old had reduced odds of UFM. Compared to children born in the western region, the odds of UFM was 2.19 (p = 0.001) and 1.60 (p = 0.097) times higher among children born in the eastern and central regions, respectively. Children living in households without electricity (AOR = 1.72, p = 0.046) and unsafe sanitation (AOR = 1.50, p = 0.049) also had a higher likelihood of dying. The pseudo-R2 for the final rural model was 0.08.

The results showed that environmental factors were attributable to 46% of the effect of place of residence on UFM (Table 3). The effect of rural residence attenuated but remained statistically robust when adjusted for socioeconomic and bio-demographic factors. Although the association was not statistically strong in the models that adjusted for all groups of variables, the effect size was >1, suggesting the role of other factors not examined in this study. The associations between other birth interval measures and UFM were not significant (the results are not presented).


To our knowledge, this is the first study that examined the UFM determinants by rural-urban residence using a nationally representative dataset. Our study found that the weighted UFMR in rural Bhutan was 2.75 times higher than in urban Bhutan. Parity and household size were the common factors influencing UFM, indicating higher parity and smaller household size were associated with greater odds of mortality in both rural and urban areas. The findings suggest that environmental factors (electricity availability, the region of residence and safe sanitation) are important in influencing child health in rural areas. Whereas, some socioeconomic factors (employment status and education level) could play a stronger role in affecting children's health in urban Bhutan. The largest part of rural-urban disparities in UFM was attributable to the environmental factors, followed by socioeconomic factors. Our findings point to the need for targeted policies to reduce rural-urban health inequities in Bhutan that can also help accelerate the achievement of the SDG target.

The higher UFMR in rural areas in this study is consistent with findings in previous national surveys [15, 25] and other studies that showed higher mortality risk among rural children [46, 8]. This can be attributable to poor public service and health infrastructure, greater distance to health centers, poor socioeconomic and living conditions, inadequate sanitation, and safe drinking water in rural Bhutan [15, 25, 31, 32]. The results also show a greater risk of mortality among children born in households of rural communities with no safe sanitation and no access to electricity. The coverage of improved sanitation in 2012 was still low in rural (58%) compared to urban areas (93%) [15]. The high occurrence of diarrhea, malnutrition, and respiratory infection in rural Bhutan depicts the direct effect of poor sanitation and hygiene [15, 17, 25]. Studies have also demonstrated that access to electricity is associated with reduced odds of child mortality [33, 34] in rural areas [4]. Access to electricity can be crucial in improving hygiene, reducing indoor pollution from solid fuel use, enhancing health knowledge through access to radio and television [35], and may potentially increase income contributing towards child health.

Defying expectations, the impact of socioeconomic factors such as wealth index, parent's education, and working status was not evident in the adjusted analysis, especially in rural areas. Studies suggest that improved community-based health care services may buffer the adverse impact of socioeconomic factors [36, 37]. This can be particularly true in the Bhutanese context where access to health care has improved substantially, as reflected in high immunization (∼95%) and antenatal care coverage [15]. However, some environmental factors such as electricity, sanitation and solid fuel use were associated with UFM in rural areas suggesting the role of individual components used to construct the wealth index that may be useful for channeling interventions. This finding supports that environmental factors were responsible for the large part of the observed rural-urban UFM gap. Wealth index also emerged to be a significant UFM predictor in the unadjusted model.

Furthermore, the findings that children born to non-working mothers and uneducated fathers were more likely to experience UFM may indicate poor socioeconomic status to be an important factor contributing to health disparities in urban Bhutan. This is plausible given the rapid rural-urban migration leading to significant growth of the urban population in Bhutan, leaving a migrant population that may be socioeconomically disadvantaged. Poor urban areas were also found to experience higher UFM rates and stunting than in rural areas [11]. Addressing the health of urban residents, therefore, becomes equally important. More studies are needed to understand the health disparities in urban Bhutan. The results also show that urban under-five children of educated fathers had reduced odds of dying. While the benefit of the mother's education on their children's health is well documented [38], father's education can also be important in determining the child's health [39] through improving household economic conditions [40], care for children, and health utilization.

This study found that under-five children born in households with >4 members had decreased odds of dying independent of rural-urban residence. This is in accordance with findings from previous research [35, 41]. Larger family size may have more members to take care of children [42] and higher household income. In Bhutan, the practice of extended family living is still prevalent [42], and one in five people (20%) surveyed were found to live in an extended family [31]. Another finding in this study was that children in households with non-working household heads had reduced mortality risk. The household head is usually a male [15] and the oldest person in the house and may not be economically active; instead, he/she may be engaged in looking after the children at home while the parents of the child are away working. This situation is plausible, particularly in urban Bhutan, where the rising living costs may compel both the couples to work. Besides, older parents might also support the family financially. The actual mechanism of the effect of household size and economic activity of household heads requires further investigation.

Similar to the findings in other studies [43, 44], higher parity was correlated with a higher propensity of UFM both in rural and urban areas. The decline in the mother's health status and possible competition for attention and restricted resources among siblings may be some explanations for this. The multivariable analysis also showed that children from older mothers had reduced odds of dying in rural areas. Some previous work also found similar findings [8, 35, 41]. Better economic status and the ability to cope with pregnancy and related situations among older mothers can explain this attenuated risk. In agreement with the findings of a national survey [25], this study found that under-five children from rural areas of central and eastern regions had an increased probability of dying. Comparatively, western Bhutan is more developed with better public infrastructure and services. Malnutrition and stunting in children, inadequate access to safe sanitation, high prevalence of anemia in pregnant women, poor health literacy and poverty in rural communities of these regions are some plausible reasons for the higher UFM risk [17, 25, 31, 32].

Strengths and limitations

The use of complex samples analysis to prevent biased estimation arising from sampling strategy adopted in the 2012 BNHS and the stepwise regression informed by a widely used conceptual framework are some of the strengths of this study. However, there are some limitations. First, since this was a cross-sectional study, the directionality of the associations identified cannot be made. Owing to the huge missing observation for health-related factors and the non-availability of information on birth weight, these factors were not examined. The influence of health-related variables, including breastfeeding on UFM is indisputable [5, 9, 33, 43], and can be more important than socioeconomic and other distal factors with improvement in socioeconomic conditions. The effect of these factors needs to be investigated in future studies. The possible underestimation of UFMR from interviewing only those mothers who were alive and bias arising from recalling birth and mortality inherent to demographic and health surveys cannot be ruled out. Furthermore, to cover all possible under-five deaths of children being studied, all births from November 2002 to October 2007 were extracted, which allowed for five full calendar years of life (exposure). However, socio-economic and other characteristics were only collected in 2012. Given the time elapsed between birth extraction and socioeconomic variables collection periods, the associations identified in this study could be imprecise and need to be interpreted with an understanding of the changes in socioeconomic conditions in Bhutan over time. Future studies with matched measurement periods are warranted. Finally, the small number of under-five deaths for some categories of the variables examined in urban areas could have potentially influenced the results.


The data of BNHS 2012 suggests that UFMR was higher in rural than urban Bhutan, and the determinants influencing UFM varied by rural and urban areas. Older mother's age, living in the western region, and access to safe sanitation and electricity were significantly associated with lower odds of UFM among rural children. Whereas in urban areas, children of unemployed mothers and uneducated fathers had higher UFM odds, while those born in households with unemployed household heads had lower odds of dying. Lower parity and larger household size were associated with reduced UFM risk regardless of rural-urban residence. Similar studies with better designs, including those that investigate infant mortality and urban health disparities, are required to inform focused policies.

Focused policy interventions are needed to reduce the rural-urban disparities in UFM that can further help Bhutan hasten the attainment of the SDG target. Specifically, increasing health access and utilization, enhancing health knowledge, programs to improve access to electricity and sanitation and hygiene, and empowering women through education can potentially help improve child survival in rural areas. Targeting socioeconomically disadvantaged urban populations with health and socioeconomic improvement programs, and women empowerment through education can be cost-effective in urban areas. Additionally, equitable regional socioeconomic development policies are essential to reduce rural-urban health inequities. These may be undertaken in collaboration with relevant stakeholders such as the education sector, urban municipalities, primary health centers and the local governments and non-government organizations working at the grass-root level, including the community. These cross-sectoral interventions can not only help accelerate the reduction of UFMR beyond that of the SDG target but may also improve other socioeconomic and health indicators that in turn can positively impact development.


Adjusted association between under-five mortality and the determinants in urban Bhutan

Figure 1

Adjusted association between under-five mortality and the determinants in urban Bhutan

Adjusted association between under-five mortality and the determinants in rural Bhutan

Figure 2

Adjusted association between under-five mortality and the determinants in rural Bhutan

Characteristics and distribution of under-five mortality rate in urban and rural Bhutan

VariablesUrban (N = 1399, 27 died)Rural (N = 4990, 231 died)
Nn1 (%1)UFMRnn1 (%1)UFMR
(a) Bio-demographic determinants
Mother's age
≤28 years393417 (28.0)912021215 (25.7)46
29–34 years608653 (43.8)2317601621 (34.3)44
>34 years398420 (28.2)1120281884 (39.9)42
Mother's age when first married
<18 years445484 (32.5)2017551728 (32.5)52
≥18 years936982 (65.9)1428992713 (65.9)40
Not reported/missing1824 (1.6)0336279 (5.9)31
Moth ers age when first pregnant
<20 years563603 (40.5)1723182231 (47.3)47
≥20 years747790 (53.0)1624782286 (48.4)38
Not reported/missing8997 (6.5)7194202 (4.3)73
Mother's marital status
Married13241408 (94.6)1644184270 (90.5)45
Not married7581 (5.4)0572450 (9.5)31
Sex of the child
Girl685715 (48.0)1025072365 (50.1)38
Boy714774 (52.0)2124832354 (49.9)50
Total births
≤2 births664794 (53.3)614081317 (27.9)15
>2 births735695 (46.7)2735813402 (48.03)55
Not reported/missing 110
Sex of the household head
Female369403 (27.0)421911811 (38.4)43
Male10301086 (73.0)2027992909 (61.6)44
Birth order
First572658 (44.2)913961348 (28.6)39
Second413420 (28.2)2212712088 (24.7)38
≥Third414411 (27.6)1923231284 (46.7)50
Birth interval
<33 months10101,001 (67.2)1639283659 (77.5)45
≥33 months3650 (3.4)3999111 (2.4)19
Not reported/missing353438 (29.4)12963949 (20.1)40
(b) Socioeconomic determinants
Mother's education
No education583610 (40.9)3032233052 (64.6)45
With education814876 (58.8)617511660 (35.2)41
Not reported/missing24 (0.3)0168 (0.2)37
Education of the husband
No education363382 (25.7)3723872322 (49.2)47
With education1,0131077 (72.3)922662122 (45.0)43
Not reported/missing2330 (2.0)0337276 (5.8)22
Education of the household head
No education434452 (30.4)3132513015 (63.9)42
With education9611036 (69.5)917261698 (36.0)47
Not reported/missing41 (0.1)0137 (0.1)41
Household size
≤4583691 (46.4)2011821189 (25.2)61
>4816798 (53.6)1138083531 (74.8)38
Mother's working status
Working438488 (32.8)514681422 (30.1)34
Not working947992 (66.6)2134703267 (69.2)48
Not reported/missing149 (0.6)05230 (0.6)37
Household head's working status
Working11961276 (85.7)1731153168 (67.1)43
Not working200211 (14.2)418581537 (32.6)46
Not reported/missing32 (0.2)01714 (0.3)0
Wealth index
Poor/second/middle192194 (13.0)1439343712 (78.6)50
Fourth/richest12071295 (87.0)1610561008 (21.4)20
Domestic violence
Yes142124 (8.3)23729646 (13.7)36
No12571365 (91.7)1542594073 (86.3)45
Not reported/missing 21 (0.03)0
(c) Environmental determinants
Region of residence
Western6351037 (69.6)1312781563 (33.1)23
Central488226 (15.2)2518971364 (28.9)45
Eastern276227 (15.2)1818151792 (38.0)61
Place of residence
Urban13991489 (24.0)16
Rural 49904720 (76.0)44
Safe drinking water
Yes13941482 (99.5)1648404561 (96.6)44
No57 (0.5)0144151 (3.2)27
Not reported/missing 69 (0.2)220
Safe sanitation facilities
Yes12791386 (93.1)1628582681 (56.8)32
No116100 (6.7)1421262035 (43.1)59
Not reported/missing43 (0.2)064 (0.1)0
Solid fuel use
Yes1620 (1.3)4417841807 (38.3)56
No13831469 (98.7)1532052913 (61.7)36
Not reported/missing 11 (0.01)0
Electricity availability
Yes13941485 (99.7)1642084034 (85.5)38
No44 (0.3)0778677 (14.3)78
Not reported/missing11 (0.01)049 (0.2)0

Note(s): UFMR: weighted estimates of under-five mortality rate as deaths per 1000 live births; n: sample size; n1: weighted sample size; %1: weighted percentage

Unadjusted association between under-five mortality and determinants in urban and rural Bhutan

COR95% CIp-valueCOR95% CIp-value
(a) Bio-demographic determinants
Mother's age (Ref: ≤28 years) 0.246 0.844
29–34 years2.630.53–13.160.2340.970.67–1.400.864
>34 years1.220.30–4.930.7740.910.63–1.290.564
Mother's age when first married (Ref: <18 years) 0.469 0.066
≥18 years0.700.26–1.870.4690.760.57–1.020.066
Mothers age when first pregnant (Ref: <20 years) 0.917 0.148
≥20 years0.950.39–2.360.9170.810.60–1.080.148
Mother's marital status (Ref: Married) 0.103
Not marriedNA0.690.44–1.090.103
Sex of the child (Ref: Girl) 0.132 0.109
Total births (Ref: ≤2 births) <0.001 0.001
>2 births4.932.30–10.56<0.0013.711.87–7.350.001
Sex of the household head (Ref: Female) 0.017 0.755
Birth order (Ref: First) 0.355 0.255
Birth interval (Ref: <33 months) 0.415 0.262
≥33 months2.540.26–24.710.4150.400.08–2.120.262
(b) Socioeconomic determinants
Mother's education (Ref: No education) 0.035 0.623
With education0.190.04–0.880.0350.900.57–1.420.623
Education of the husband (Ref: No education) 0.027 0.594
With education0.230.06–0.840.0270.920.67–1.270.592
Education of the household head (Ref: No education) 0.048 0.559
With education0.270.08–0.990.0481.120.76–1.650.559
Household size (Ref: ≤4) 0.273 0.022
Mother's working status (Ref: Not working) 0.040 0.071
Household head's working status (Ref: Working) 0.015 0.732
Not working0.240.08–0.750.0151.090.64–1.860.732
Wealth index (Ref: Poorest/Second/Middle) 0.815 <0.001
Domestic violence (Ref: No) 0.594 0.364
(c) Environmental determinants
Region of residence (Ref: Western) 0.312 0.002
Safe drinking water (Ref: Yes) 0.555
Safe sanitation facilities (Ref: Yes) 0.856 0.002
Solid fuel use (Ref: No) 0.235 0.008
Electricity availability (Ref: Yes) 0.011

Note(s): COR: crude/unadjusted odds ratio; CI: confidence interval; Ref: reference group; NA: not available

Change in effect of place of residence (rural-urban) after adjusting for other groups of independent variables

Place of residence (ref: Urban)OR for ruralProportion change in OR
OR95% CIp-value
Model 1 (unadjusted model)2.891.74–4.81<0.001
Model 22.211.32–3.690.003−24%
Model 31.711.01–2.890.044−41%
Model 41.560.90–2.710.111−46%
Model 51.170.68–2.030.547−60%

Note(s): OR: Odds ratio; Model 2: Adjusted for significant bio-demographic variables from bivariate analysis (mother's age when first married, total births, birth order); Model 3: Adjusted for significant socioeconomic variables from bivariate analysis (husband's education, mother's working status, wealth index); Model 4: Adjusted for significant environmental variables from bivariate analysis (region, safe sanitation facilities, solid fuel use, electricity availability); Model 5: Adjusted for only for significant bio-demographic, socioeconomic, and environmental determinants


1United Nations Children's Fund [UNICEF]. Levels and trends in child mortality report 2017. New York, NY: UNICEF; 2017.

2United Nations. Sustainable Development Goals: 17 goals to transform our world. New York: United Nations; 2015.

3Tabutin D, Masquelier B. Mortality inequalities and trends in low-and middle-income countries, 1990-2015. Population. 2017; 72(2): 227-307. doi: 10.3917/popu.1702.0227.

4Van de Poel E, O'donnell O, Van Doorslaer E. What explains the rural-urban gap in infant mortality: household or community characteristics?. Demography. 2009; 46(4): 827-50.

5Adewuyi EO, Zhao Y, Lamichhane R. Risk factors for infant mortality in rural and urban Nigeria: evidence from the national household survey. Scand J Public Health. 2017 Jul; 45(5): 543-54. doi: 10.1177/1403494817696599.

6Saikia N, Singh A, Jasilionis D, Ram F. Explaining the rural-urban gap in infant mortality in India. Demogr Res. 2013 Sep; 29(18): 473-505. doi: 10.4054/DemRes.2013.29.18.

7Sastry N. What explains rural-urban differentials in child mortality in Brazil?. Soc Sci Med. 1997; 44(7): 989-1002. doi: 10.1016/s0277-9536(96)00224-9.

8Ettarh RR, Kimani J. Determinants of under-five mortality in rural and urban Kenya. Rural Rem. Health. 2012; 12: 1812.

9Yi B, Wu L, Liu H, Fang W, Hu Y, Wang Y. Rural-urban differences of neonatal mortality in a poorly developed province of China. BMC Publ. Health. 2011 Jun; 11(1): 477. doi: 10.1186/1471-2458-11-477.

10Günther I, Harttgen K. Deadly cities? Spatial inequalities in mortality in sub-Saharan Africa. Popul Dev Rev. 2012; 38(3): 469-86. doi: 10.1111/j.1728-4457.2012.00512.x.

11Van de Poel E, O'Donnell O, Van Doorslaer E. Are urban children really healthier? Evidence from 47 developing countries. Soc Sci Med. 2007 Nov; 65(10): 1986-2003. doi: 10.1016/j.socscimed.2007.06.032.

12Cleland J, Bicego G, Fegan G. Socioeconomic inequalities in childhood mortality: the 1970s to the 1980s. Health Transit Rev. 1992 Apr; 2(1): 1-18.

13Bocquier P, Madise NJ, Zulu EM. Is there an urban advantage in child survival in sub-saharan Africa? Evidence from 18 countries in the 1990s. Demography. 2011 May; 48(2): 531-58. doi: 10.1007/s13524-011-0019-2.

14Lalou R, LeGrand TK. Child mortality in the urban and rural Sahel. Popul. 1997; 9: 147-68.

15Ministry of Health. National health survey report 2012. Thimphu: Ministry of health, Royal Government of Bhutan; 2012.

16United Nations Children's Fund [UNICEF]. The state of the world's children 2017: children in a digital world. New York, NY: UNICEF; 2017.

17Ministry of Health. Annual health bulletin 2018. Thimphu: Ministry of Health, Royal Government of Bhutan; 2018.

18National Statistics Bureau. Statistical year book of Bhutan 2018. Thimphu: National Statistics Bureau, Royal Government of Bhutan; 2018.

19The World Bank. Leveraging urbanization in Bhutan. [updated 2015 September; cited 2019]. Available from:

20Asian Development Bank. Proposed grants and technical assistance grant Kingdom of Bhutan: health sector development program. (Project Number 51141-002). Asian Development Bank; 2018.

21Jamtsho T. Urbanization and water, sanitation and hygiene in Bhutan. Regional Health Forum. 2010; 14(1): 1-7.

22Sharma J, Zangpo K, Grundy J. Measuring universal health coverage: a three-dimensional composite approach from Bhutan. WHO South East Asia J Public Health. 2014 Jul-Dec; 3(3): 226-37. doi: 10.4103/2224-3151.206745.

23Jamtsho S, Corner R, Dewan A. Spatio-temporal analysis of spatial accessibility to primary health care in Bhutan. ISPRS Int J Geoinf. 2015 Sep; 4(3): 1584-604. doi: 10.3390/ijgi4031584.

24Dendup T, Zhao Y, Dema D. Factors associated with under-five mortality in Bhutan: an analysis of the Bhutan national health survey 2012. BMC Publ. Health. 2018 Dec; 18(1): 1375. doi: 10.1186/s12889-018-6308-6.

25National Statistics Bureau. Bhutan multiple indicator survey 2010. Thimphu: National Statistics Bureau; 2010.

26Mosley WH, Chen LC. An analytical framework for the study of child survival in developing countries. Popul Dev Rev. 1984; 10: 25-45.

27World Health Organization [WHO]. Report of a WHO technical consultation on birth spacing: Geneva, Switzerland 13-15 June 2005. Geneva: WHO; 2007.

28Victora CG, Huttly SR, Fuchs SC, Olinto MT. The role of conceptual frameworks in epidemiological analysis: a hierarchical approach. Int J Epidemiol. 1997 Feb; 26(1): 224-7. doi: 10.1093/ije/26.1.224.

29Bursac Z, Gauss CH, Williams DK, Hosmer DW. Purposeful selection of variables in logistic regression. Source Code Biol Med. 2008 Dec; 3(1): 17. doi: 10.1186/1751-0473-3-17.

30StataCorp LLC. Stata: data analysis and statistical software [Computer software]. Special Edition Release. StataCorp LLC; 2015.

31National Statistics Bureau. Bhutan living standards survey 2012. Thimphu: National Statistics Bureau and Asian Development Bank; 2013.

32National Statistics Bureau. Bhutan poverty assessment 2014. Thimphu: National Statistics Bureau; 2014.

33Charmarbagwala R, Ranger M, Waddington H, White H. The determinants of child health and nutrition: a meta-analysis. Washington, DC: World Bank; 2004.

34Karmaker SC, Lahiry S, Roy DC, Singha B. Determinants of infant and child mortality in Bangladesh: time trends and comparisons across South Asia. Bangladesh J. Med. Sci. 2014; 13(4): 431-7.

35Pandey MK. Maternal health and child mortality in rural India. (ASARC Working Paper 2009/12). Canberra: Australia South Asia Research Centre (ASARC); 2009.

36Welaga P, Moyer CA, Aborigo R, Adongo P, Williams J, Hodgson A, et al. Why are babies dying in the first month after birth? A 7-year study of neonatal mortality in northern Ghana. PLoS ONE. 2013; 8(3): e58924. doi: 10.1371/journal.pone.0058924.

37Phillips JF, Bawah AA, Binka FN. Accelerating reproductive and child health programme impact with community-based services: the Navrongo experiment in Ghana. Bull World Health Organ. 2006 Dec; 84(12): 949-55. doi: 10.2471/blt.06.030064.

38Cleland JG, Van Ginneken JK. Maternal education and child survival in developing countries: the search for pathways of influence. Soc Sci Med. 1988; 27(12): 1357-68.

39Breierova L, Duflo E. The impact of education on fertility and child mortality: Do fathers really matter less than mothers?. Cambridge, MA: National bureau of economic research; 2004.

40Lindeboom M, Llena-Nozal A, van der Klaauw B. Parental education and child health: evidence from a schooling reform. J Health Econ. 2009 Jan; 28(1): 109-31. doi: 10.1016/j.jhealeco.2008.08.003.

41Singh R, Tripathi V. Under-five mortality among mothers employed in agriculture: findings from a nationally representative sample. PeerJ. 2015 Jan; 2: e710. doi: 10.7717/peerj.710.

42Barth F, Wikan U. Situation of children in Bhutan: an anthropological perspective. Thimphu: Centre for Bhutan Studies; 2011.

43Adhikari R, Podhisita C. Household headship and child death: evidence from Nepal. BMC Int Health Hum Rights. 2010 Jun; 10(1): 13. doi: 10.1186/1472-698X-10-13.

44Lee HY, Van Do D, Choi S, Trinh OT, To KG. Trends and determinants of infant and under-five childhood mortality in Vietnam, 1986-2011. Glob. Health Action. 2016; 9: 29312. doi: 10.3402/gha.v9.29312.


The authors would like to thank the Ministry of Health of Bhutan for giving permission and access to use the dataset of the 2012 BNHS. The authors thank all the participants of the 2012 BNHS.Funding: No funding was received for the conduct of this study. The study used a secondary dataset.Conflict of interest: The authors have no conflicts of interest to declare for this study.

Corresponding author

Tashi Dendup can be contacted at: