Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Determinants and Causes of Neonatal Mortality in Jimma Zone, Southwest Ethiopia: A Multilevel Analysis of Prospective Follow Up Study

  • Gurmesa Tura Debelew ,

    gurmesatura@gmail.com

    Affiliation Department of Reproductive Health and Health Service Management, School of Public Health, Addis Ababa University, Addis Ababa, Ethiopia

  • Mesganaw Fantahun Afework,

    Affiliation Department of Reproductive Health and Health Service Management, School of Public Health, Addis Ababa University, Addis Ababa, Ethiopia

  • Alemayehu Worku Yalew

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Addis Ababa University, Addis Ababa, Ethiopia

Abstract

Background

Ethiopia is among the countries with the highest neonatal mortality with the rate of 37 deaths per 1000 live births. In spite of many efforts by the government and other partners, non-significant decline has been achieved in the last 15 years. Thus, identifying the determinants and causes are very crucial for policy and program improvement. However, studies are scarce in the country in general and in Jimma zone in particular.

Objective

To identify the determinants and causes of neonatal mortality in Jimma Zone, Southwest Ethiopia.

Methods

A prospective follow-up study was conducted among 3463 neonates from September 2012 to December 2013. The data were collected by interviewer-administered structured questionnaire and analyzed by SPSS V.20.0 and STATA 13. Verbal autopsies were conducted to identify causes of neonatal death. Mixed-effects multilevel logistic regression model was used to identify determinants of neonatal mortality.

Results

The status of neonatal mortality rate was 35.5 (95%CI: 28.3, 42.6) per 1000 live births. Though significant variation existed between clusters in relation to neonatal mortality, cluster-level variables were found to have non-significant effect on neonatal mortality. Individual-level variables such as birth order, frequency of antenatal care use, delivery place, gestation age at birth, premature rupture of membrane, complication during labor, twin births, size of neonate at birth and neonatal care practice were identified as determinants of neonatal mortality. Birth asphyxia (47.5%), neonatal infections (34.3%) and prematurity (11.1%) were the three leading causes of neonatal mortality accounting for 93%.

Conclusions

This study revealed high status of neonatal mortality in the study area. Higher-level variables had less importance in determining neonatal mortality. Individual level variables related to care during pregnancy, intra-partum complications and care, neonatal conditions and the immediate neonatal care practices were identified as determinant factors. Improving antenatal care, intra-partum care and immediate neonatal care are recommended.

Introduction

Globally, about 6.6 million children die before their 5th birthday each year. About 5 million of this occurs in the first year of life and nearly 3 million die within the first 28 days of birth. This indicates that about 44% of under-five deaths and 60% of infant deaths are accounted by the neonatal mortality. Moreover, the share of neonatal mortality from the under-five death rose from 37% in 1990 to 44% in 2012 [1]. This clearly points out that it is difficult to achieve the desired millennium development goal (MDG4) target for the two-thirds reduction of child mortality by 2015 without particular focus on neonatal mortality.

More than 98% of these deaths occur in developing countries. Sub-Saharan Africa has the highest risk of death in the first month of life and among the regions showing the least progress in reducing the neonatal mortality rate. Most of these deaths are caused by infectious diseases, pregnancy-related complications, delivery-related complications, including intra-partum asphyxia, birth trauma and premature birth which can easily be prevented [1], [2].

In Ethiopia, neonatal mortality rate (NMR) has long been very high. In spite of many efforts by the government and other stakeholders, non-significant and very sluggish decline has been achieved in the last 15 years. The neonatal mortality rate for the years 1991–1995, 1996–2000, 2001–2005 and 2006–2011 were 46, 42, 39 and 37 per 1000 live births, respectively [3], [4], [5]. Moreover, about 63% of infant deaths in the country occur during the first month of life. Thus, accelerated reduction in neonatal mortality is increasingly critical for progress towards the MDG4.

To do this, identifying the determinants and causes of neonatal mortality at the local context is very crucial and timely issue. However, many of the neonatal deaths happen at home and unrecorded that make obtaining sampling-frame and sources of data very difficult. As a result, studies on this topic including verbal autopsies (VAs) are limited in Ethiopia. Moreover, the very few available studies are only facility-based and cross-sectional in design, which are not the preferred designs to establish causal relationships necessitating community-based prospective study supported by verbal autopsy (VA).

Therefore, this longitudinal community based study aimed to fill these gaps by determining the status of neonatal mortality and identifying the determinants and causes at different levels by applying multilevel analysis. The findings of the study will also be used as inputs for policy makers and program implementers at national as well as regional levels to design evidence-based intervention strategies to tackle the problems of neonatal mortality.

Methods and Materials

Study design and setting

This study was a community-based prospective follow up conducted in Jimma Zone from September 2012-December 2013. Jimma Zone is one of the 17 Zones of the Oromia Regional State of Ethiopia having a total of 17 rural districts and two town administrations. According to the 2007 national population and housing census, the Zone has a total population of 2.6 million, of which 88.7% are rural residents [6], [7].

Sample size and sampling technique

The minimum required sample size for this study was determined by using Epi-Info V.3.5.1 by considering two sample comparisons of proportions based on the following assumptions. The outcome variable was neonatal mortality. Among all the determinants of neonatal mortality considered, educational status of mothers was found to give the largest sample size. Based on this, the prevalence of neonatal mortality among mothers having educational status of secondary or above was estimated to be 4.0% (P1 = 0.040) and among those who didn't attend secondary education was to be 8.1% (P2 = 0.081) [8]; 95% level of confidence and 80% power were considered. A ratio of 1∶3 was used (r = 3). As multistage-clustered sampling method was used, a design effect of 2 was considered. Finally, 10% was added for non-responses and miss-to-follow up and the final sample size became 3604.

Multistage-clustered sampling technique was used to identify a cohort of pregnant women to be enrolled in the follow up for the study. At first stage, the Zone was stratified as rural districts (17 in number) and town administrations (2 in number, Jimma and Agaro). Then, by considering time and logistics, 5 districts (30%) were selected by simple random sampling from the 17 districts. At second stage, all the selected 5 districts were clustered by ‘Kebeles’ (A ‘kebele’ is the smallest administrative unit having 5000 population in average) and stratified in to urban and rural ‘Kebeles’.

Then, by simple random sampling method, 9 rural ‘Kebeles’ and 2 urban ‘Kebeles’ were selected from each selected district. This number of clusters (‘kebeles’) was determined based on expected number of pregnant women per ‘Kebele’. Jimma town administration and Agaro town administration have 13 and 5 ‘Kebeles’, respectively and all were included purposefully. With this, a total of 73 Clusters (‘Kebeles’) were included in the study. Then, for all selected ‘kebeles’, pregnant women were enumerated by using house-to-house visit and all obtained were enrolled in the study (Figure S1).

Measurements

The dependent variable for this study was neonatal mortality and the independent variables were divided into two levels. Level 2 (higher-level variables) included community or cluster level variables such as place of residence, access to health centers and access to hospitals. Level 1 (lower-level variables) included individual and household characteristics such as: socio-demography, wealth quintiles, maternal obstetric factors, maternal health care use, conditions of labor, characteristics of the neonates and neonatal care practices. The detail descriptions and measurements are given below (Table 1).

thumbnail
Table 1. Description of variables and measurement for the study, Jimma Zone, Southwest Ethiopia, September 2012-December 2013.

https://doi.org/10.1371/journal.pone.0107184.t001

Instruments

The data were collected by using pre-tested interviewer administered structured questionnaires which were adapted from different literatures. The indicators for the wealth index were adapted from Ethiopian Demographic and Health Survey (EDHS) [5]. Indicators to measure birth preparedness and complication readiness (BP & CR) were adapted from the safe motherhood questionnaires developed by maternal and neonatal health program of Johns Hopkins Program for International Education in Gynecology and Obstetrics (JHPIEGO) [9]. Indicators for neonatal care practices were adapted from the World Health Organization (WHO) minimum neonatal care packages [10]. Data on causes of neonatal death were collected by using structured verbal autopsy questionnaire adapted from the standard VA questionnaire developed and validated by WHO, Johns Hopkins University (JHU) and London School of Hygiene and Tropical Medicine [11]. All the questionnaires were prepared in English, then translated to local languages ‘Afan Oromoo’ and Amharic and used to collect the data after back translating to English by different experts to check its consistency.

Data collection process

As this was prospective follow up study, data were collected in three phases. First, home-to-home visit was made to enumerate pregnant women from the selected 73 clusters. Then, all the identified pregnant women were enrolled in the study as a cohort. At a baseline, data on basic socio-demography, economy and birth preparedness and complication readiness were collected. Then, just at the end of neonatal period, maternal service use (antenatal care (ANC), delivery place and attendant and postnatal care), conditions of labor, neonatal characteristics and neonatal care practices were collected. For died neonates, VAs were conducted within 15–30 days of death.

Females, who had completed 10th grade or above were recruited, trained and collected the data. The VAs were conducted by two experienced females. The data collection process was supervised strictly by trained supervisors and principal investigators. To control the quality of data, in addition to training, pretest, supervision and use of local languages, the inter-item consistency of the indicators to measure the composite score of wealth index, BP & CR and neonatal care practices were checked by using Chronbach-alpha at 0.7 cut-off points.

Data management and analysis

The collected data were coded and entered into Epidata V.3.1 to minimize logical errors and design skipping patterns. Then, the data were exported to SPSS for windows version 20.0 for cleaning, editing and analysis. Descriptive analysis was done by computing proportions and summary statistics. Socioeconomic quintiles were determined by using Principal Component Analysis (PCA). Birth preparedness and complication readiness was computed by composite indicator of five items. Similarly, neonatal care practice was determined by composite variable of 12 items by using PCA. As Jimma and Agaro town administrations were both purposefully included, the status of neonatal mortality was estimated by calculating weighted percentage based on the complex sample survey procedure to avoid underestimation.

Bivariate analysis was done by using cross-tabulation to see associations between the dependent and independent variables. Then, all variables having P-value <0.25 were considered as candidates for the final model. As multistage-clustered sampling technique was used because of the different levels of factors, mixed-effects multilevel logistic regression model was used by using STATA 13. This model was preferred in order to avoid the clustering effects as well as ecological fallacy. ‘Kebeles’ were considered as clusters and ‘kebele’ level variables were taken as higher-level (level 2). Neonates were nested within their family and households. As a result, neonatal individual variables, delivery conditions and household characteristics were taken as lower level (level 1) variables (Table 1).

Goodness of fit of the multilevel model was tested by the log likelihood ratio (LR) test. To evaluate the extent of the cluster variation in influencing neonatal mortality, intercept-only model was fitted as Logit (pij)  = γ00+u0j and Interclass Correlation Coefficient (ICC = ρ) was determined by dividing variance between groups (δ2µ0) by variance between groups (δ2µ0) plus variance within groups (δ2e). However, within group variance (δ2e) can't be directly obtained for dichotomous outcome variables; instead, it was estimated by dividing δ2µ0 by δ2µ0 plus ∏2/3. Then, to identify the determinant factors, the full model was fitted as: Logit(pij)  = γ00 + γ01 Zj + γ10 Xij + u0j + u1j Xij. Where: Logit(pij)  =  dependent variable at unit i in cluster j, Xij =  individual explanatory variable in cluster j, Zj =  group level explanatory variable, γ00 =  fixed intercept, γ01 and γ10 =  fixed slopes and u0j and u1j =  random effects at level 2.

Multicollinearity between the independent variables was assessed by using variance inflation factors (VIF>10 considered as existence of multicollinearity) before interpreting the final output. However, only skill of delivery attendant (VIF = 10.9) had multicollinearity with place of delivery (VIF = 9.1, reduced to 1.8 when delivery attendant was dropped). As a result, they were included in the model alternatively by dropping the other. For the rest of the variables, the VIF was <3. In addition, cross-level two-way interactions were checked, particularly between place of residence, access to health facilities and maternal health care use (ANC, place of delivery and delivery attendants). Individual-level two-way interactions between prematurity (gestation age at birth), twin births and size of neonate at birth were checked. However, no significant interaction was detected (P>0.05 for each). The VAs were interpreted by two independent pediatricians and third pediatrician interpreted in case of disagreements.

Ethical consideration

Ethical approval was obtained from the Institutional Review Board (IRB) of College of Health Sciences of Addis Ababa University as well as IRB of Oromia Regional State Health Bureau. Following this, formal letters and permissions were secured from all respective local administrators. Written informed consent was obtained from each respondent before actual data collection. Issues of confidentiality were maintained by removing any identifiers from the questionnaire. To protect vulnerable group, data collectors were trained to maintain confidentiality and provide necessary health information based on the need of the participants and arrange referral to health facilities for sick neonates.

Results

Response rate

It was planned to include a sample of 3604 neonates. However, after excluding incomplete questionnaires, abortion cases, missed-to-follow up, maternal deaths and stillbirths, a total of 3463 live-births were included in the analysis for this study making a response rate to be 96.1%.

The detail process and flow of the study is indicated below (Figure 1).

thumbnail
Figure 1. Flow-diagram of the overall study process, Jimma Zone, Southwest Ethiopia, Sept 2012-Dec 2013.

This figure shows that a total of 3696 pregnant women were obtained during enumeration. After excluding 14 women because of different reasons, 3682 women were interviewed at a baseline. Again after excluding 70 incomplete and inconsistent questionnaires, 3612 were included in the analysis and enrolled in the follow up. After the follow up, total of 3472 deliveries happened; among which 38 were twins and 47 were stillbirths. From the total of 3463 live births, 110 died before 28 days.

https://doi.org/10.1371/journal.pone.0107184.g001

Socio-demographic characteristics

Of the total 3463 live-births included in the analysis, 2602(75.1%) were from urban residence. Majority, 2209(63.8%), of the mothers of the neonates were in the age group of 20–29 years with a mean and standard deviation of 26.6±5.0. Oromoo was the dominant ethnic group, 3033(87.6%) and Muslim was the leading religion, 3019(87.2%). More than half of the mothers, 1871(54.0%), didn't attend any formal education. The great majority, 3280(94.7%), of the mothers were housewives and farmer was the leading occupation of their husbands, 2459(71.0%). Nearly half, 1779(51.4%), of the neonates were males with male-to-female ratio of 1.06∶1.00. The rate of twin births was 76(2.2%) (Table 2).

thumbnail
Table 2. Socio-Demographic characteristics of Respondents, Jimma Zone, Southwest Ethiopia, September 2012-December 2013 (n = 3463).

https://doi.org/10.1371/journal.pone.0107184.t002

Status of neonatal mortality

From a total of 3463 live-births, 110 died before 28 days of birth making weighted neonatal mortality rate of 35.5 (95%CI: 28.3, 42.6) per 1000 live-births (urban  = 20.0 (95%CI: 9.0, 31.0) and rural 36.2 (95%CI: 28.7, 43.7). From these, 76(69.1%) died within the first week of life making weighted early neonatal mortality rate to be 23.7 (95%CI: 18.5, 30.3). In this study, 47 stillbirths happened making a stillbirth rate and perinatal mortality rate to be 16.5 (95%CI: 12.2, 22.4) and 39.8 (95%CI: 32.9, 48.1) per 1000 births, respectively (Table 3).

thumbnail
Table 3. The status of neonatal mortality in Jimma Zone, Southwest Ethiopia, September 2012-December 2013 (n = 3463).

https://doi.org/10.1371/journal.pone.0107184.t003

Determinants of neonatal mortality

To evaluate the applicability of the mixed-effects multilevel logistic regression model, the ICC (ρ) was calculated in the empty model and it was found to be 0.100 indicating that 10.0% of the variation is contributed by between-cluster variation. The test of the preference of log likelihood Vs logistic regression was also strongly significant (P = 0.0003). Then, the full model was run by including all the cluster level and individual level variables and the ICC (ρ) was increased to 0.174 suggesting model improvement. This again indicated that 17.4% of the variation is attributed to cluster level variables suggesting the preference of multilevel analysis. The preference of log likelihood Vs logistic regression was again improved and strongly significant (P = 0.0001) (Table S1).

After adjusting in the final two-level mixed-effects logistic regression model, cluster level variables (place of residence, access to BEmOC and CEmOC) had non- significant associations with neonatal mortality. Similarly, among the lower level variables maternal socio demography and economy such as age, education, occupation and wealth quintiles had non-significant associations. Whereas, maternal obstetric and health care during pregnancy and delivery, delivery conditions, neonatal conditions and neonatal care practice had significant associations.

First birth order (OR = 5.45; 95%CI: 1.81, 16.40) and birth order of 5th or above neonates (OR = 2.61; 95%CI: 1.43, 4.74) were more likely to die during neonatal period as compared to 2nd–4th order. Those whose mothers had 1–3 ANC visits (OR = 0.51; 95%CI: 0.28, 0.93) and 4 or more visits (OR = 0.35; 95%CI: 0.18, 0.68) were less likely to die during neonatal period as compared to those who had no ANC visit at all. Neonates born at health centers were 43% (OR = 0.43; 95%CI: 0.17, 0.99) less likely to die during neonatal period as compared to those who were born at home. However, hospital delivery (OR = 0.73; 95%CI: 0.31, 1.70) and skilled attendants (OR = 0.57; 95%CI: 0.28, 1.16) had non-significant association as compared to home delivery and non skilled attendants, respectively. Prematurity (GA at birth <37weeks) was found to increase the likelihood of neonatal death as compared to term births (OR = 2.09; 95%CI: 1.03, 4.22). Premature and prolonged rupture of membrane before the onset of labor had increased the likelihood of neonatal death. Rupture of membrane 1–12 hours (OR = 2.71; 95%CI: 1.13, 6.53) and >12 hours (OR = 7.74; 95%CI: 2.27, 26.38) before the onset of labor had significantly higher risk of neonatal death as compared to rupture of membrane after the onset of labor.

The occurrence of obstetric complications during labor (OR = 6.77; 95%CI: 3.82, 12.00) and twin births (OR = 8.21; 95%CI: 3.46, 19.47) were among the strong predictors of neonatal mortality. Similarly, small size (OR = 1.95; 95%CI: 1.11, 3.42) and big size (OR = 10.73; 95%CI: 5.65, 20.37) at birth were found to increase the likelihood of neonatal death as compared to average size neonates. Not having good comprehensive neonatal care practice was the other strong predictor of neonatal mortality (OR = 10.36; 95%CI: 5.13, 20.94) (Table 4).

thumbnail
Table 4. Multilevel analysis of factors associated with neonatal mortality, Jimma Zone, Southwest Ethiopia, Sept 2012-Dec 2013.

https://doi.org/10.1371/journal.pone.0107184.t004

Causes of neonatal mortality

Out of the total 110 neonatal deaths occurred, VAs were conducted for 99 cases. Initially, the two interpreting physicians agreed on the most probable cause of death for 86 out of 99 cases, which was 86%. For 6 (6.1%) of the cases, first and second most probable causes were exchanged between the physicians making the overall agreement rate to be 93%. As a result, third physician interpreted the 13 cases and agreed with either of the physicians in 11 cases and the rest 2 cases were not agreed upon and classified as unspecified cause of death. With this, birth asphyxia (47.5%), neonatal infections (34.3%) and prematurity (11.1%) were the three leading causes of neonatal mortality (Figure 2).

thumbnail
Figure 2. Causes of neonatal death, Jimma Zone, Southwest Ethiopia, Sept 2012–Dec 2013.

This figure shows that birth asphyxia (47.5%), neonatal infections (34.3%) and prematurity (11.1%) were the three leading cause of neonatal death accounting for nearly 93%.

https://doi.org/10.1371/journal.pone.0107184.g002

Discussions

The status of neonatal mortality rate in this study is 35.5 (95%CI: 28.3, 42.6) per 1000 live births. This is similar with the finding of the EDHS 2011 in which NMR was 37.0 (95%CI: 33.7, 40.3) [5]. However this is higher than the findings of other countries with high neonatal mortality like India (19.5) [12] and Indonesia (23.8) [13]. This clearly indicated that the situation of neonatal mortality is still high and non-progressing suggesting targeted interventions by all partners at different levels.

In this study, 63 (57.3%) and 76 (69.1%) out of the 110 died neonates occurred on the date of birth and in the 1st week of life, respectively making early NMR to be 23.7 per 1000 live-births. This is in line with other prior studies in which more than three-quarter of neonatal deaths occurred in the first week of life [13][17]. The reason could be majority of neonatal mortality in developing countries are related to conditions of labor, intra-partum and the immediate newborn care practices. Moreover, the major causes of death are birth asphyxia, early neonatal infections and prematurity. This again clearly pointed out that neonatal survival interventions have to target the intra-partum as well as immediate and early neonatal periods.

In this study, higher-level factors (place of residence, access to health centre and hospital) were found to be less important in predicting neonatal mortality. Similarly, lower level factors related to basic socio-demography and economy were less important. Instead, lower level factors related to maternal and neonatal complications and care before, during and after delivery were the most important predictors of neonatal mortality. This finding is consistent with some other prior studies in which higher level variables had less importance in determining neonatal mortality [14], [16], [18], [19]. However, this finding is in contrast to the findings of study done in rural India in which both the higher and lower level factors were equally important [12]. This pointed out that the higher-level factors that are less amenable to short term interventions, have to be considered as distant factors and special focus needs to be given to the immediate proximal factors of neonatal mortality which are more feasible to be intervened.

First birth order and birth order of five or above were found to increase the likelihood of neonatal mortality by more than five and two times respectively. This may be due to high risk of occurrence of complications during delivery among nuliparous and grandmultiparous mothers. In this study, having ANC visit and giving birth at health center were found to decrease neonatal mortality significantly. This is consistent with previous study conducted in Ethiopia and other countries [13], [20]. This may be because; during ANC visits, necessary health condition of mothers can be screened and treated earlier. Moreover, health facility delivery is very necessary in detecting complications earlier and providing clean and safe delivery.

In the contrary, delivery at hospital and skilled attendants at birth had no significant association with neonatal mortality. Similarly, in the systematic review and meta-analysis we conducted before this study, in 9 of the 19 studies included, health facility delivery had no significant effect on neonatal mortality [21]. Similar findings have been reported in some other studies in the country and abroad [15], [22]. This may be because, majority of hospital deliveries in developing countries, including Ethiopia, are as a result of self referral after all attempts failed. In this study for example, 43.5% of hospital deliveries were because of occurrence of problems during labor. This highlighted the importance of addressing the first and the second delays. This means, giving birth at hospital attended by skilled attendant may not be a guarantee to avert neonatal mortality unless the woman comes to the right health facility at the right time. The other possible explanation could be the low coverage of hospital delivery and skilled attendant at birth that lead to low power to detect significant effects.

Intra-partum conditions, neonatal conditions and immediate neonatal care practices were the most important determinants of neonatal mortality in this study. Prematurity (<37weeks of GA at birth) and the occurrence of maternal complications during labour were found to increase the likelihood of neonatal mortality by two times and nearly seven times, respectively. Similarly, rupture of membrane within 1–12 hours and >12 hours before the onset of labor had increased neonatal mortality by nearly three times and eight times, respectively. Twin births was also increased the risk of neonatal mortality by eight times as compared to singletons. Small size and big size at birth were increased the likelihood of neonatal mortality by about two and eleven times, respectively. Neonates having poor neonatal care were about ten times more likely to die during neonatal period as compared to those who received good comprehensive neonatal care.

These findings are consistent with other previous studies conducted in the country and abroad in which the intra-partum and neonatal conditions were found to be the important predictors of neonatal mortality [14], [15], [17], [18], [19], [22]. This may be explained by the fact that premature and twin births are more likely to be under-weight and more prone to complications and infections. Similarly, long staying PROM and intra-partum complications increase the risk of infections and birth asphyxia. Thus, provision of comprehensive neonatal care including clean cord care, thermal care and appropriate feeding have the potential to avert some of these risks significantly.

In line with this, the verbal autopsy identified birth asphyxia, neonatal infections and prematurity as the major causes of neonatal mortality accounting for 93%. This is also in line with other prior studies where the three causes accounted for more than four-fifth of neonatal mortality [14], [23][25]. While infection is the leading cause of neonatal mortality in many studies in developing countries [23], birth asphyxia was the leading cause of death (47.5%) in this study. This finding may be explained by the high proportion of home deliveries in the absence of skilled attendants. Knowledge and skill deficiencies in prevention, diagnosis and management by unskilled attendants might have contributed to perinatal complications associated with neonatal asphyxia in these communities. The association between elevated neonatal mortality and home deliveries by unskilled attendants highlights the importance of delivery care education to local service providers.

As a programmatic implication, this study pointed out that targeting individual level factors like intra-partum conditions and the immediate neonatal period can significantly improve neonatal survival. Moreover, targeting the three major causes of death (birth asphyxia, neonatal infections and prematurity) can avert more than nine in ten of neonatal deaths.

This study may have its own limitation in that some of the medical terms were difficult to translate exactly to local languages, which might have affected the respondents' understanding. To reduce this limitation, local language experts translated the instruments. In addition, local data collectors who were fluent in local languages and familiar with local terms collected the data. The verbal autopsy was based on mothers' report of signs and symptoms concerning the underlying causes of neonatal death. This may not be as specific as the clinical diagnosis in identifying the exact cause of death. Therefore, future researches need to consider additional clinical and laboratory-based identification of causes of neonatal death by including facility based data for better interventions. In addition, the verbal autopsies were conducted within 15–30 days of neonatal death which might lead to recall bias of the exact conditions of death of the neonates. To address this, standardized and validated questionnaire was adapted in order that the mothers remember the conditions and report more valid data.

Conclusions

This study revealed that the status of neonatal mortality in the study area is still high. The great majority of neonatal deaths occurred in the first week of life. Higher-level variables were found to be less important in determining neonatal mortality. Lower (individual) level variables were found to be important determinants of neonatal mortality. Birth order, ANC visits, place of delivery, prematurity, PROM, complication during labor, twin births, size of neonate at birth and neonatal care practices were identified as determinants of neonatal mortality. Birth asphyxia, neonatal infections and prematurity were identified as the leading causes of neonatal mortality accounting for more than 90%. Ensuring the continuum of care from pregnancy through delivery to the immediate postnatal period should be in place so as to address neonatal mortality in the study area. Specifically, increasing adequate ANC visits and health facility delivery by addressing the first and second delays are very crucial. Provision of comprehensive neonatal care such as cord care, thermal care and early initiation of breast-feeding are recommended.

Supporting Information

Figure S1.

Schematic presentation of sampling method. This figure shows the multistage clustered sampling methods based the proportional allocation to the size of the population.

https://doi.org/10.1371/journal.pone.0107184.s001

(TIF)

Table S1.

Parameter coefficients and test of goodness-of-fit of the mixed effect multilevel model, in Jimma Zone, Southwest Ethiopia, September 2012-December 2013. This table shows the parameter estimates of the multilevel logistic regression, including the fixed effects, random effect at level 2, the Infraclass correlation Coefficient, LR test and level of significance both in the empty-model and full model.

https://doi.org/10.1371/journal.pone.0107184.s002

(DOCX)

Acknowledgments

The authors would like to thank all the study participants, data collectors, supervisors and all other individuals involved in any process of this study.

Author Contributions

Conceived and designed the experiments: GTD MFA AWY. Performed the experiments: GTD MFA AWY. Analyzed the data: GTD MFA AWY. Contributed reagents/materials/analysis tools: GTD MFA AWY. Wrote the paper: GTD MFA AWY. Read the final manuscript and approved for submission: GTD MFA AWY.

References

  1. 1. UNICEF WHO, The World Bank, United Nations Population Division (2013) Levels and Trends in Child Mortality: Estimates Developed by the UN Inter-agency Group for Child Mortality Estimation (UNIGME). New York, USA: UNICEF, WHO, the World Bank, United Nations Population Division.
  2. 2. UNICEF (2009) The State of the World's Children: Maternal and Newborn Health. New York, USA: UNICEF.
  3. 3. Central Statistics Agency (CSA) [Ethiopia], ORC Macro (2001) Ethiopian Demographic and Health Survey 2000. Addis Ababa, Ethiopia, and Calverton, Maryland, USA: Central Statistical Authority and ORC Macro.
  4. 4. Central Statistics Agency (CSA) [Ethiopia], ORC Macro (2006) Ethiopian Demographic and Health Survey 2005. Addis Ababa, Ethiopia, and Calverton, Maryland, USA: Central Statistical Authority and ORC Macro.
  5. 5. Central Statistics Agency (CSA) [Ethiopia], ORC Macro (2011) Ethiopian Demographic and Health Survey 2011. Addis Ababa, Ethiopia, and Calverton, Maryland, USA: Central Statistical Authority and ORC Macro.
  6. 6. Jimma Zonal Health Desk (2011) Annual report of the Zonal Health Desk of Jimma Zone. Jimma, Ethiopia: Jimma Zonal Health Desk.
  7. 7. Central Statistical Agency (CSA) [Ethiopia] (2007) National Population Census Report of Ethiopia. Addis Ababa, Ethiopia: CSA [Ethiopia].
  8. 8. Kebede B, Gebeyehu A, Sharma H, Yifru S (2012) Prevalence and associated factors of neonatal mortality in North Gondar Zone, Northwest Ethiopia. Ethio J Health Dev 26: 66–71.
  9. 9. Maternal and Neonatal Health Program JHPIEGO (2004) Monitoring Birth Preparedness and complication readiness tools and indicators for maternal and newborn health. Maternal and Neonatal Health Program JHPIEGO.
  10. 10. World Health Organization (WHO) (2009) Managing programs to improve child health. Geneva, Switzerland: WHO. WHO website. Available: http://whqlibdoc.who.int/publications/2009/9789241598729_eng_Module1.pdf?ua=1 Accessed 2012 Feb 10.
  11. 11. Anker M, Black RE, Coldham C, Kalter HD, Quigley MA, et al. (1999) A standard verbal autopsy method for investigating causes of death in infants and children. Geneva, Switzerland: WHO/CDS/CSR/ISR/99.4. WHO website. Available: http://whqlibdoc.who.int/hq/1999/WHO_CDS_CSR_ISR_99.4.pdf Accessed 2012 March 20.
  12. 12. Upadhyay R, Dwivedi P, Rai S, Misra P, Kalaivani M, et al. (2011) Determinants of Neonatal Mortality in Rural Haryana: A Retrospective Population Based Study. Indian Pediatrics 49: 291–294.
  13. 13. Titaley CR, Dibley MJ, Agho K, Roberts CL, Hall J (2008) Determinants of neonatal mortality in Indonesia. BMC Public Health 8: 232–247.
  14. 14. Tariq P, Kundi Z (1999) Determinants of Neonatal Mortality. JPMA 49: 56–60.
  15. 15. Mekonnen Y, Tensou B, Telake DS, Degefie T, Bekele A (2013) Neonatal mortality in Ethiopia: trends and determinants. BMC Public Health 13: 483–487.
  16. 16. Paude D, Thapa A, Shedain PR, Paudel B (2013) Trends and Determinants of Neonatal Mortality in Nepal: Further Analysis of the Nepal Demographic and Health Surveys, 2001–2011. Calverton, Maryland, USA: Nepal Ministry of Health and Population, New ERA, and ICF International.
  17. 17. Ben Hamida Nouaili E, Chaouachi S, Ben Said A, Marrakchi Z (2010) Determinants of neonatal mortality in a Tunisian population. Tunis Med 88: 42–45.
  18. 18. Quamrul H, Islam R, Hossain K (2010) Effects of Demographic Characteristics on Neonatal, Post neonatal, Infant and Child Mortality. Current Research Journal of Biological Sciences 2: 132–138.
  19. 19. Eltinay S, Rahman S, Ansari WE, Nimeri N, Latiph E, et al. (2013) Fetal and perinatal determinants of neonatal mortality in the state of Qatar: A pilot phase univariate and multivariate analysis by PEARL study. J Postgrad Med Inst 27: 235–241.
  20. 20. Arokiasamy P, Gautam A (2008) Neonatal mortality in the empowered action group states of India: trends and determinants. J Biosoc Sci 40: 183–201.
  21. 21. Tura G, Fantahun M, Worku A (2013) The effect of health facility delivery on neonatal mortality: systematic review and meta-analysis. BMC Pregnancy and Childbirth 13: 18–27.
  22. 22. Mahmood A (2002) Determinants of Neonatal and Post-neonatal Mortality in Pakistan. The Pakistan Development Review 41: 723–744.
  23. 23. Lawn J, Kerber K, Enweronu-Laryea C, Cousens S (2010) 3.6 million neonatal deaths: What is progressing and what is not? Semin Perinatol 34: 371–386.
  24. 24. Manji K (2009) Situation analysis of newborn health in Tanzania: Current situation, existing plans and strategic next steps for newborn health. Dar es Salaam: Ministry of Health and Social Welfare, Save the Children.
  25. 25. Federal MOH [Ethiopia], SESOG, L10K (2009) Addressing Community Maternal and Neonatal Health in Ethiopia. Report from National Scoping Exercise and National Workshop to Increase Demand, Accesses and Use of Community Maternal and Neonatal Health Services. Addis Ababa, Ethiopia: Federal MOH [Ethiopia], SESOG, L10K.