Article Text

Regional differences in the utilisation of antenatal care and skilled birth attendant services during the COVID-19 pandemic in Nigeria: an interrupted time series analysis
  1. Rotimi Felix Afolabi1,
  2. Mobolaji M Salawu1,
  3. Eniola A Bamgboye1,
  4. Segun Bello1,
  5. Ayo Stephen Adebowale1,2,
  6. Magbagbeola David Dairo1,
  7. Steven N Kabwama3,
  8. Irene Wanyana4,
  9. Susan Kizito5,
  10. Rawlance Ndejjo5,
  11. Rhoda K Wanyenze5,
  12. Olufunmilayo Ibitola Fawole1
  1. 1 Department of Epidemiology and Medical Statistics, Faculty of Public Health, College of Medicine, University of Ibadan, Ibadan, Nigeria
  2. 2 Population and Health Research Entity, Faculty of Humanity, North-West University, Mafikeng, South Africa
  3. 3 Department of Community Health and Behavioural Sciences, School of Public Health, College of Health Sciences, Makerere University, Kampala, Uganda
  4. 4 Department of Epidemiology and Biostatistics, School of Public Health, College of Health Sciences, Makerere University, Kampala, Uganda
  5. 5 Department of Disease Control and Environmental Health, School of Public Health, College of Health Sciences, Makerere University, Kampala, Uganda
  1. Correspondence to Dr Rotimi Felix Afolabi; rotimifelix{at}yahoo.com

Abstract

Introduction The COVID-19 pandemic has had a substantial negative impact on the utilisation of essential health services (EHS) globally, especially in resource-limited settings such as Nigeria. High maternal deaths associated with low access to and utilisation of EHS such as antenatal care (ANC) and skilled birth attendants (SBAs) remain a concern during the COVID-19 era. The study assessed the COVID-19 pandemic effects on ANC and SBA utilisation across regions in Nigeria.

Methods Monthly data on ANC and SBA between January 2017 and July 2021 were obtained from the Federal Ministry of Health database. An interrupted time-series analysis, implemented using the Prophet model, was conducted to compare the regional variation of outcomes during the COVID-19 pandemic. Average percentage changes (PC) between the observed and predicted outcomes including their 95% CI were reported.

Results From March 2020 to July 2021, the number of ANC visits was significantly lower than expected by a 16%–43% change in five of the six regions in Nigeria. The highest significant reduction was in North-West (PC=−43.4; 95% CI: −52.6 to –34.1) and the least in South-West (PC=−15.5; 95% CI: −24.8 to –6.1), with no significant change in the South-East. The number of deliveries by SBA was significantly lower than expected by a 18%–43% change in all the regions (p<0.01). North-East (PC=−43.3; 95% CI: –51.7 to –34.9) and South-West (PC=−18.3; 95% CI: −25.2 to –11.5), respectively, had the highest and the least decline in SBA utilisation. Overall, ANC and SBA patterns of change were relatively similar across the north-south divide though the change effect was considerably pronounced in the north.

Conclusion There was a substantial reduction in ANC and SBA utilisation due to the COVID-19 pandemic in Nigeria, especially in the northern regions. Targeted and contextually relevant interventions should be implemented to alleviate the impact of emergency response on access to EHS and promote access to care during the pandemic.

  • COVID-19
  • Public Health

Data availability statement

Data are available on reasonable request. Data analysed in this study were obtained from DHIS2, domiciled with the FMoH Nigeria. Readers interested in replicating the analysis may contact the FMoH Nigeria for access to the data.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • The COVID-19 non-pharmaceutical interventions such as movement restrictions have been shown to disrupt routine essential health services (EHSs).

  • Policy-makers use research evidence to make informed decisions on maintaining EHS during public health emergencies. Such evidence has included the impact of pandemics on the utilisation of EHS such as antenatal care (ANC).

WHAT THIS STUDY ADDS

  • The COVID-19 pandemic substantially lowered ANC attendance and utilisation of skilled birth attendant (SBA) services in the six regions of Nigeria, although with some regional differences.

  • The northern regions witnessed a significant higher decline in both ANC and SBA utilisation relative to the southern regions.

  • The South-South region experienced the largest reduction in ANC utilisation, while the South-East region had the largest reduction in SBA at delivery among the southern regions.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • Implementing restrictions and precautions during public health emergencies in Nigeria requires careful consideration of the attendant impact on routine EHSs.

  • The evidence emphasises the need to strengthen region-specific policies in light of sociocultural contexts that promote sustained access to quality maternal and new-born care during emergencies such as COVID-19 in Nigeria and other similar settings.

Introduction

Public health emergencies including the COVID-19 pandemic often disrupt the delivery of essential health services (EHSs) including antenatal care (ANC) visits and facility deliveries by skilled birth attendants (SBA).1 2 Evidence suggests that some interventions in response to the COVID-19 pandemic especially the non-pharmaceutical interventions have considerably disrupted the provision of EHS, affecting the lives of millions of women and families.3 According to Riley et al,4 about a 10% reduction in pregnancy-related and neonatal healthcare coverage could increase maternal and neonatal mortalities, respectively, by 28 000 and 168 000 deaths worldwide. The disruptions have been attributed to distorted health workforce and supply chains, and decreased care-seeking behaviour for non-COVID-19 healthcare-related visits, among other factors.5–7 Besides, patients’ fear of acquiring an infection while in public places during the pandemic, especially in healthcare facilities reduces the demand for EHS and keeps people away from receiving healthcare.8 Worrisome, about a decade of progress in reproductive, maternal and child health could be stalled if not reversed.6 9

The COVID-19 pandemic has impacted on the safety of childbirth, access to treatment for complications in health facilities, and delivery of already-fragile maternal and child health services.10 The disruptions affected ANC, a major determinant of maternal mortality rate, and important components of maternal care11–13 including access to institutional-based deliveries.2 4 6 10 14 15 Proper utilisation of ANC considerably enhances safe motherhood with improved maternal outcomes. Skilled assistance at birth has been promoted over the years14 16 and has positively impacted on maternal health outcomes. However, ANC service utilisation may not necessarily suggest child delivery at an institutional-based facility.3 This implies that not only improved use of ANC services but also the presence of skilled birth assistants (SBA) at delivery should be targeted in the ongoing COVID-19 pandemic response, particularly in resource-limited settings such as Nigeria.

In 2018, a national survey in Nigeria reported that 67% of women who gave birth in the 5 years preceding the survey received ANC from a SBA at least once for their last birth.17 The onset of COVID-19 in Nigeria diverted human and financial resources from routine healthcare services to pandemic-response activities, particularly in the early part of 2020.15 18 Service utilisation was especially affected by the COVID-19 non-pharmaceutical interventions such as physical distance and movement restrictions (eg, curfews, partial and total lockdowns) and the fear of being infected with the virus at these health centres.7 Using data obtained from the health information system, studies have reported mixed results for the effect of COVID-19 pandemic-related disruptions on EHS.19–21 For instance, Arsenault et al 20 reported a range of 5%–33% significant decline in maternal health services including ANC visits and facility-based delivery in nearly half of the 10 countries studied. Similarly, a study among 12 sub-Saharan Africa countries observed about a 2%–6% modest reduction in the utilisation of maternal and child health services.21 These suggest varying impact of COVID-19 pandemic on EHS across countries and regions.

In Nigeria, researchers have reported a 26% reduction in health services access, a 16% drop in the number of ANC visits, and a 6% decrease in facility deliveries in the first half of 2020 relative to 2019 due to the COVID-19 pandemic.15 22 These disruptions could leave 700 000 women without access to facility-based delivery and consequently may result in childhood and maternal mortality rate increasing by 20% and 10%, respectively, over 1 year.15 Therefore, reducing preventable maternal and child deaths requires sustained access to quality maternal and newborn care. These previous studies were based on mathematical simulations or restricted to local or qualitative studies.15 22 This study analysed real-world national data which would provide more robust findings. Specifically, the study assessed the effect of the COVID-19 pandemic on ANC and SBA utilisation across the six regions in Nigeria using an interrupted time series (ITS) analysis. The findings of the study may highlight the importance of context-specific planning for robust maternity services in any public health emergency response.

Methods

Study design and setting

This cross-sectional study was conducted to assess the effect of the COVID-19 pandemic on some EHS indicators in Nigeria. Nigeria is Africa’s most populous country and has a large, socioculturally diverse population of about 216 million.23 Administratively, Nigeria has 36 states including a Federal Capital Territory zoned into six geopolitical groups: namely, North-Central, North-East, North-West, South-East, South-South and South-West as shown in figure 1. Northern regions (North-Central, North-East and North-West) predominantly house Hausa/Fulani ethnic group and practise Islam religion; southern regions (South-East, South-South and South-West) are dominated by Christians of Igbo ethnic group in the South-East and Yoruba in the South-West.

Figure 1

Map of Nigeria showing boundaries of six geopolitical zones.

Healthcare is being delivered at three levels which are tertiary, secondary and primary healthcare, respectively, managed by the federal, state and local governments. Nigeria has a low health workers-to-population ratio of 1.6 health workers including doctors, nurses and midwives per 1000 persons, below the WHO recommendation of 4.5 health workers for every 1000 persons.22 By 5 June 2022, Nigeria had recorded 256,250 COVID-19 cases and 3143 deaths with a 1.2% case fatality rate.24

Data source and extraction

The secondary data used for the study were obtained from the health management information system—an electronic database of the Federal Ministry of Health (FMoH), Nigeria. The data were retrospectively extracted by FMoH staff who had access to the electronic medical records (EMR) system which operates on the District Health Information Software 2 (DHIS2). The nationally representative dataset provides information on several DHIS2 indicators including EHS. Specifically, the state-level data contained monthly information regarding the number of ANC visits and the number of deliveries by SBA in Nigeria from January 2017 to July 2021.

Study variables

In this study, the exposure was the COVID-19 pandemic. Two of the key components of EHS were considered as outcome measures. These are the monthly number of ANC visits and deliveries by SBA across the six regions in Nigeria. Although the first COVID-19 case was confirmed in Nigeria on 27 February 2020, 1 March 2020 was used to define the outcome measures before and after the emergence of the disease—assuming a negligible change would have happened.

Data management and analysis

Data preparation

The obtained state-level data consisted of a few missing data (1.9%) on ANC utilisation in the months of September–December 2020, across the regions but South-South. Of the 39 missing data, December (41.0%) has the highest, followed by November (38.5%) and the least (10.3%) both in September and October. To make the data conform with ITS analysis and address the identified challenges including under-reporting and over-reporting, the missing data were imputed with the mean of non-missing values. Further, a Hampel (x, window_size=5, n=3, imputation=True) function in Python was applied to detect outlier (defined as any observed value within a rolling window of size five greater than three median absolute deviations) and replace such with its rolling window’s median. Only a few outliers in the state-level data on ANC utilisation (3.8%) and on deliveries by SBA (2.1%) was flagged and replaced with the rolling median (see online supplemental appendix 1). The data were then aggregated into six regions according to the administrative geopolitical zones in Nigeria. This was plotted using time series plots to visualise regional differences in trend and any unusual features of the data.

Supplemental material

Thereafter, the ITS analysis using Prophet model was employed to evaluate the potential impact of the COVID-19 pandemic related-disruption on the two EHS indicators with the assumption that the level and trend of each outcome would remain the same if there was no interruption (COVID-19 pandemic).

Interrupted time series

An ITS analysis is often used in public health to evaluate the impact of an interruption like an intervention or an exposure.25 26 In this study, the emergence of COVID-19 disease was regarded as an interruption. The ITS regression model was used to compare the regional variation of the predicted outcomes to the observed outcomes during the COVID-19 era. Modelling of the data in the pre-COVID-19 enables assessment of the underlying trend and yields a counterfactual for what could have happened in the absence of the COVID-19 pandemic when extrapolated into the post-COVID-19 onset. Provided the underlying trend has been accounted for, differences between counterfactual (predicted) and observed data at post-COVID-19 onset could be evaluated. The ITS analysis implemented using Facebook Prophet was employed. The data analysis was implemented using Python V.3.10 programming language and the level of significance was set at 5%.

The Prophet model approach

The Facebook’s Core Data Science team developed Prophet as a forecasting model.27 Prophet is robust to outliers and shifts in the trend and it can handle various seasons of historical time series data. Relative to other approaches, studies have recommended Prophet due to its open source algorithm, automated nature, accuracy, efficient and speedy time series analysis and forecasting.28–30 The forecasting model is an additive model where nonlinear trends are modelled with yearly, monthly, weekly or daily seasonality including holiday effects.

The Prophet model is described in equation (1)

Embedded Image (1)

where:

Embedded Image is the outcome measure.

Embedded Image is the trend function (models non-periodic changes).

Embedded Image represents seasonality (models periodic changes).

Embedded Image represents effects of holidays and events (models often non-periodic shocks; this component was not considered in this study as the data do not have an effect of the holiday term in trend forecasts).

Embedded Image is an error term (models idiosyncratic changes).

Of note, Prophet employs a Fourier series to estimate the seasonality effects, and the seasonality models are specified as the periodic functions of time t.

Model hyperparameter

Hyperparameter tuning is the process of obtaining the best hyperparameters to improve model performance. The process is crucial for the optimal data training process; thus, the best set of hyperparameters is obtained before data training. After the data were preprocessed, a preliminary forecast was done to tune hyperparameters using a grid search cross-validation approach. The approach used a 4×4×2 grid of three parameters (changepoint prior scale (0.001, 0.01, 0.05, 0.1); seasonality prior scale (0.01, 0.1, 1.0, 10.0); seasonality mode (additive, multiplicative) to choose the best set of hyperparameters. The changepoint prior scale defines the trend flexibility, the seasonality prior scale defines the seasonality flexibility, while the seasonality mode identifies seasonality as either additive or multiplicative. With prediction date ranging from February 2019 to October 2019 and parallelisation over different periods, the hyperparameters were evaluated on root mean square error averaged over a 3-month horizon. The best combination of tuned parameters with the minimum RSME for each of the outcomes by region is presented in table 1.

Table 1

The tuned best-value parameters

Modelling and forecasting

In practice, a model is often trained on observed data and used to predict the counterfactual.29 31 Therefore, the data were divided into two parts: training and testing sets. The training set was used to optimise and train a model such that the structured underlying pattern of the ordered data is obtained. The testing set was used to fit and evaluate a prediction model. Specifically, data up to February 2020 were used for training and validating the model using the best-tuned hyperparameters; thereafter, the model was applied to predict the outcome for the remaining months in the data (March 2020–July 2021). Predicted and observed data were plotted together to visualise the fitting accuracy of the model and quantified to assess the level of change due to COVID-19. The change was assessed using percentage change (PC), as expressed in equation (2)

Embedded Image (2)

where: Embedded Image = actual value at time t; Embedded Image = predicted value at time t; n=number of data points

Additionally, the test of significance for the differences between the observed and predicted outcome’s mean counts was conducted using the Wilcoxon signed-rank test. The 95%CIs for PC including the corresponding p values were reported.

Model prediction accuracy

The prediction accuracy of a model using a cross-validation approach was conducted. Such that predictions were made for each of the outcomes during COVID-19 pandemic on a 3-month horizon. The prediction performance was evaluated by comparing the trends in the actual and predicted outcomes using the mean absolute percentage error (MAPE). This was calculated on a rolling window of the predictions after sorting by horizon. The MAPE is a commonly used performance metric to evaluate modelling capability and predictive ability. It is a measure of percentage errors in prediction, expressed as in equation (3).

Embedded Image (3)

The MAPE is an easy-to-interpret scale-independent measure that compares prediction performance between various scaled datasets.28 31 32 A forecast’s MAPE value ≤10% is interpreted as highly accurate, 11%–20% good, 21%–50% reasonable and >50% inaccurate forecasting.33

Patient and public involvement statement

Patients or the public were not involved in the design, implementation, reporting or dissemination plans of this study

Results

The trend of ANC and SBA utilisation measures by region between January 2017 and July 2021 in Nigeria

The trends of the observed monthly counts of ANC visits and deliveries by SBA are presented in figure 2. The highest number of ANC visits was observed in North-West, followed by North-East, while the least count was observed in South-East though the trend was similar to that of South-South. Except for 2017 when South-West had relatively lower counts compared with North-Central, the two regions had relatively the same trend of ANC utilisation. From January 2020, a decrease in the count was observed in most regions, prominently in the North-West and North-East regions. The observed decline in counts was, however, reversed starting from October 2020.

Figure 2

Regional trend for monthly observed counts of ANC visits and deliveries by SBA. ANC, antenatal care; SBA, skilled birth attendant.

An increasing trend in the number of deliveries by SBA was observed prominently in the northern regions. This, however, was reversed at the onset of the pandemic with a decrease in the monthly counts for SBA utilisation at delivery. Similar to the trend in ANC utilisation, the North-West had the highest number of deliveries by SBA, followed by the North-East and South-West regions.

Visualisation of the fitting accuracy of the model across the regions in Nigeria

Predicted and observed counts of ANC visits comparison

The fitted accuracy of the model for ANC utilisation using both the observed and cumulative counts are shown in figure 3 and Online supplemental figure 1, respectively. The models, for each of the regions, perform well as most observed counts are within the 95% CI for the predicted values before the emergence of the COVID-19 pandemic. If there were no COVID-19 outbreak, the monthly reported number of ANC visits was expected to increase in all the regions, except for South-East and South-West regions that had downward trend movements. Additionally, the observed monthly ANC visits nosedived during COVID-19 in most regions, but South-East region (figure 3). Corroborating the fitting accuracy of the model, the predicted and observed cumulative trend values were approximately the same pre-COVID-19 emergence period in all the regions (see online supplemental figure 1). Reduced ANC utilisation was observed post-COVID-19 emergence in all regions, except the South-East region, as the expected monthly cumulative number of ANC visits was respectively higher than the observed values.

Supplemental material

Figure 3

ITS plots showing the observed and predicted monthly values of ANC visits according to regions in Nigeria. Note: The blue portion indicates the 95% CI for prediction, while the grey area indicates the COVID-19 era (March 2020–July 2021). ANC, antenatal care; ITS, interrupted time series.

Predicted and observed counts of deliveries by SBA comparison

The fitted accuracy of the models for SBA utilisation at delivery, comparing both the observed and cumulative counts with predicted values in each of the regions are respectively shown in figure 4 and online supplemental figure 2. Before the emergence of the COVID-19 pandemic, the models fit well as most observed counts for deliveries by SBA are within the 95% CI for the prediction (figure 4). If there was no COVID-19 outbreak, the monthly reported number of deliveries by SBA was expected to increase in all regions except South-West region. The observed monthly deliveries by SBA decreased sharply at the onset and during the COVID-19 pandemic in the northern regions relatively compared with the southern regions (figure 4). Similarly, predicted and observed cumulative trend values were approximately the same pre-COVID-19 emergence in all the regions (see online supplemental figure 2). The counterfactual monthly cumulative number of deliveries by SBA was respectively higher than the observed values in all regions.

Supplemental material

Figure 4

ITS plots showing the observed and predicted monthly values of deliveries by SBA according to regions in Nigeria. Note: The blue portion indicates the 95% CI for prediction, while the grey area indicates the COVID-19 era (March 2020–July 2021). ITS, interrupted time series; SBA, skilled birth attendant.

Assessment of the level of change in outcome measures due to COVID-19

The PC in average counts of ANC visits and deliveries by SBA, within March 2020 and July 2021, respectively, are shown in figure 5. The chart reveals negative PC between the observed and the predicted mean counts for both outcomes in most regions. Table 2 also corroborates these findings.

Figure 5

Percentage change in predicted and observed mean counts for ANC visits and deliveries by SBA according to regions in Nigeria (March 2020–July 2021). ANC, antenatal care; SBA, skilled birth attendant.

Table 2

Differences between the observed and predicted outcome’s mean counts by regions in Nigeria

ANC utilisation

The PCs in the number of ANC visits between observed and predicted values were significant in all regions but in South-East. The significant negative PC (p<0.001) in ANC utilisation observed in South-South (PC=−32.2; 95% CI: –37.6 to –26.8) was similar to that of the northern regions (North-East—PC=−37.4, 95% CI: –45.9 to –29.0; North-Central—PC=−30.4; 95% CI: –36.5 to –24.4; North-West—PC=−43.4; 95% CI: –52.6 to –34.1) as their 95% CI overlapped in figure 5. The change effect was, however, different from that of other southern regions.

SBA utilisation

In all the regions, a significant negative PC (p<0.001) in the mean difference of deliveries by SBA was observed with the highest decline in North-East (PC=−43.3; 95% CI: –51.7 to –34.9) and the least in South-West (PC=−18.3; 95% CI –25.2 to –11.5). Figure 5 reveals that the pattern of changes in observed and predicted mean counts for SBA utilisation was similar among the northern regions and the southern regions as their 95% CI overlapped, respectively.

Evaluation of model performance

The model performance metric—MAPE is presented in table 3. A high-accurate prediction of ANC utilisation was observed in all regions (MAPE<10%) except for the North-West (MAPE=15.5%) and South-South (MAPE=11.2%) regions with a good prediction for about 3 months. The prediction accuracy improved from a month to 3-month horizon forecast for ANC utilisation nearly in all the regions. The North-Central (MAPE=1.1%), South-East (MAPE=2.6%) and South-West (MAPE=1.4%), respectively, had a low error rate indicating a very good prediction (table 3).

Table 3

Model performance metric based on MAPE

A high-accurate prediction (MAPE<10%) of SBA utilisation at delivery for the quarterly predictions was observed in the South-West, South-South, North-Central and North-West regions, while a good prediction was observed in North-East (MAPE<13.3%) and South-East (MAPE<11.9%) regions. The error rate improved from being good to a high-accurate prediction for the quarterly prediction in North-West, while the reverse scenario was observed in North-East and South-East regions. A low error rate, indicating a very good prediction, was observed in both South-South (MAPE=0.8%) and North-West (MAPE=3.4%) in all the horizons (table 3).

Discussion

This study assessed the influence of the COVID-19 pandemic on the utilisation of ANC and SBA services across the six regions in Nigeria between March 2020 and July 2021. The significant negative PC in the utilisation of the EHS indicators in this study suggests that the COVID-19 pandemic substantially lowered ANC attendance and utilisation of SBA services in the six regions of Nigeria. Regional differences in the effects of COVID-19 on the utilisation of ANC and SBA services were observed with a reduction of service utilisation in the northern regions of Nigeria compared closely with each other. Due to a long-standing sociocultural belief system in most parts of northern Nigeria, the demand for health services had always been lower and more elastic than the demand from southern Nigeria. A typical example is reflected in the vaccination coverage which had been perennially lowered in the North with a DPT3 coverage as low as less than 20% in some northern states.34 It is not surprising that any disruptions in routine health services would produce a higher and longer-lasting impact in most parts of Nigeria’s northern regions.

Of note, the northern region with a lower COVID-19 pandemic prevalence35 compared with the southern region has the highest volume of EHS indicators utilisation as shown in this study. The highest volume recorded in the northern regions may be associated with the higher population in the regions relative to the southern regions. The northern regions, however, experienced a relatively substantial decline in both ANC and SBA utilisation. This suggests that the extent of EHS disruptions was not directly linked to the COVID-19 disease burden. North-south divide in Nigeria’s population could indeed be linked with the observed regional differences in EHS utilisation due to the COVID-19 pandemic and its policy responses. This emphasises the need to strengthen region-specific policies that promote sustained access to quality maternal and newborn care during emergencies such as the COVID-19 pandemic in Nigeria.

The health-system shocks such as the COVID-19 pandemic and its response have been connected to a reduction in EHS demand and delivery. In this study, the numbers of ANC visits and deliveries by SBA were substantially lower between March 2020 and July 2021 than the expected figures in the absence of the COVID-19 pandemic in five of the six regions of Nigeria. Nigeria, such as many countries across the globe, has experienced COVID-19-related morbidity and mortality among healthcare workers, leading to EHS delivery staff shortages.22 This shortage could cause delay in EHS access and poor usage of approved health facilities, among others. Before the COVID-19 pandemic, Nigeria’s health system was weak in terms of EHS delivery. It is possible that the pandemic, and perhaps the lockdown and other non-pharmaceutical interventions implemented to mitigate the spread of the disease further undermined the access to and utilisation of the EHS in Nigeria.36 Other studies7 11 37 38 have alluded to several other reasons such as drug stock-out, harassment by law enforcement workers, and transportation difficulties as attributable to a significant reduction in the utilisation of EHS during the COVID-19 pandemic.

In this study, ANC services utilisation declined by 16%–43%, while SBA at deliveries reduced by 18%–43%. This finding was in agreement with a study among 12 countries in sub-Saharan Africa including Nigeria21 that reported a reduction in the utilisation of EHS. Even though the impact of the COVID-19 pandemic on EHS access was comparable in the northern regions, South-South was the only southern region with a substantial impact on its ANC visits. This may be due to the lower levels of ANC utilisation in these regions compared with South-East and South-West as previously reported in studies conducted in Nigeria.39 40 As a matter of fact, the northern region is noted for poor healthcare-seeking behaviour which got worsened during the pandemic. The observed differences may further be compounded during the pandemic due to the higher incidence of poverty in the northern regions where healthcare systems are more likely to be inadequate or fragile. The link between access to care and poverty remains. This may critically impact on EHS affordability, availability and accessibility thereof. However, in the later part of the period when movement restriction including lockdowns were lifted, failure to seek care predominated.

The impact of COVID-19 related disruption on the use of SBA at delivery was substantially higher in the northern regions relative to the southern regions. The simultaneous reduction in ANC and SBA utilisation in the northern regions may suggest a decrease in demand for EHS. The associated barriers to decreased care-seeking behaviours such as distorted health workforce and supply chains, fragile health facilities, poor access and fear of infection have been reported in the literature globally.5–8 The varying degree of regional differences in the reduction of ANC utilisation may be compounded by cultural and socioeconomic differences such as educational attainment, religion and ethnicity, especially along the north-south divide of Nigeria, as suggested in previous studies.39 41 42 Additionally, instability and conflict such as Boko Haram insurgence (predominantly in North-East), Banditry (North-West), herdsmen and farmer clash (North-Central), and militancy (South-South) may further explain the reduction in EHS utilisation.43 44

Study limitations and strengths

A few limitations were observed in this study. First, the data from EMR could have provided an incomplete picture of the actual values of each outcome measure due to reporting bias and missing values. Second, the obtained data constrained the exploration of additional variables while assessing the impact of the pandemic on the EHS indicators. Of note, missing data issues are still a concern but gradually reducing. Even though health facility reports are known to have multiple data quality issues, evidence of improvements abounds.21 Third, the representativeness of the DHIS data is controversial, especially with limited clear-cut subnational denominators. Thus, it would be unclear whether the declines in service utilisation could be entirely ascribed to failure to seek care or whether there was a factual less need to seek care. Furthermore, the lack of subnational denominators restricted more granular analysis to clarify whether any increase or decrease in the number of pregnancies (and by extension, ANC service utilisation) by late 2020, was attributable to COVID-19. Nonetheless, the study has been strengthened using a large nationally representative dataset. Besides, the strength of the study includes the use of the Facebook Prophet model which is robust to outliers and shifts in the trend, and of high accuracy for prediction.22–24 These enhanced reliable inferences. Also, the quantification of the regional variations in the level of change due to the COVID-19 pandemic in the utilisation of the EHS indicators remains a strength of this study.

Conclusions

Regional variation and substantial decline in the utilisation of ANC and SBA were observed in Nigeria during the COVID-19 pandemic in March 2020–July 2021. While the northern regions witnessed a significant decline in both ANC and SBA utilisation due to the COVID-19 pandemic, South-South experienced the most reduction in ANC utilisation and South-East the most reduction in SBA at delivery among the southern regions. These findings highlight the need for context-specific interventions including tracking of disruptions of EHS and appropriate mitigation during disease outbreaks and pandemics. The forecasting results have the potential to assist governments and stakeholders in policy decisions to maintain progress towards achieving the Sustainable Development Goal health targets by 2030.

Data availability statement

Data are available on reasonable request. Data analysed in this study were obtained from DHIS2, domiciled with the FMoH Nigeria. Readers interested in replicating the analysis may contact the FMoH Nigeria for access to the data.

Ethics statements

Patient consent for publication

Ethics approval

This study used a secondary dataset obtained from FMoH that uses aggregated data and does not contain personal identifiers. Administrative authorisation to access and use the data was granted by FMoH Nigeria.

Acknowledgments

We appreciate the administrative staff (FMoH Nigeria) who helped facilitate the release of the data for this study.

References

Supplementary materials

Footnotes

  • Handling editor Seema Biswas

  • Twitter @salawubolaji, @wanyanairene, @rawleng

  • Contributors RFA conceptualised the study; MMS, EAB, SB and OF sourced data for the study; RFA led the descriptive and modeling analyses with substantive support from ASA, MMS, SB and EAB. RFA prepared the initial draft. All authors reviewed critically the manuscript, made substantial contribution to produce the final manuscript and approved the final manuscript. RFA acts as the guarantor of this study.

  • Funding Financial support for this study was obtained from the Bill & Melinda Gates Foundation (INV-019313).

  • Disclaimer The views, opinions, and content of this publication are those of the authors and do not necessarily reflect the views, opinions, or policies of the Bill and Melinda Gates Foundation.

  • Map disclaimer The inclusion of any map (including the depiction of any boundaries therein), or of any geographic or locational reference, does not imply the expression of any opinion whatsoever on the part of BMJ concerning the legal status of any country, territory, jurisdiction or area or of its authorities. Any such expression remains solely that of the relevant source and is not endorsed by BMJ. Maps are provided without any warranty of any kind, either express or implied.

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.