Geographic barriers to care persist at the community healthcare level: Evidence from rural Madagascar

Health


Abstract

Geographic distance is a critical barrier to healthcare access, particularly for rural communities with poor transportation infrastructure who rely on non-motorized transportation. There is broad consensus on the importance of community health workers (CHWs) to reduce the effects of geographic isolation on healthcare access. Due to a lack of fine-scale spatial data and individual patient records, little is known about the precise effects of CHWs on removing geographic barriers at this level of the healthcare system. Relying on a high-quality, crowd-sourced dataset that includes all paths and buildings in the area, we explored the impact of geographic distance from CHWs on the use of CHW services for children under 5 years in the rural district of Ifanadiana, southeastern Madagascar from 2018–2021. We then used this analysis to determine key features of an optimal geographic design of the CHW system, specifically optimizing a single CHW location or installing additional CHW sites. We found that consultation rates by CHWs decreased with increasing distance patients travel to the CHW by approximately 28.1% per km. The optimization exercise revealed that the majority of CHW sites (50/80) were already in an optimal location or shared an optimal location with a primary health clinic. Relocating the remaining CHW sites based on a geographic optimum was predicted to increase consultation rates by only 7.4%. On the other hand, adding a second CHW site was predicted to increase consultation rates by 31.5%, with a larger effect in more geographically dispersed catchments. Geographic distance remains a barrier at the level of the CHW, but optimizing CHW site location based on geography alone will not result in large gains in consultation rates. Rather, alternative strategies, such as the creation of additional CHW sites or the implementation of proactive care, should be considered.

Introduction

In 2012, the UN General Assembly declared universal health coverage key to sustainable human and economic development. However, the majority of the world’s population still lacks access to essential heath services, including primary care [1]. There are multiple barriers to primary care access, including financial costs, low provider to population ratios, and geographic distance [2, 3]. Geographic distance to primary health clinics (PHCs) is especially problematic in contexts where the transportation network is sparse and the primary mode of transportation is non-motorized [4]. In these primarily rural regions, the use of primary care decreases exponentially for populations living at increasing distances from PHCs, known as the “distance-decay” effect [59]. Low access in turn negatively affects population health metrics such as maternal mortality [10] or hospitalization rates for severe malaria [11].

Community health programs, often implemented via community health workers (CHWs), have emerged as a potential solution to increase healthcare access in geographically isolated communities [12]. Usually residents of these communities, CHWs are located closer to their patients than PHCs and provide basic maternal and child care, as well as a range of services such as HIV and TB treatment. CHW programs became a focus of international development funding organizations following the Alma Alta declaration in 1978 and were implemented in countries around the world. CHWs have since played an integral role in national health systems and several universal health coverage strategies [1315]. Although located in closer proximity to underserved communities than other levels of the health system (e.g. PHC, district hospital), CHWs can be responsible for large geographic catchments or isolated, rural populations, limiting their ability to adequately serve everyone. As such, the World Health Organization recommends that community health programs adapt CHW catchments according to the local geography, including the proximity of households and population density [16]. Yet, it is unknown whether a meaningful ‘distance-decay’ relationship—similar to that observed for health facilities—can persist at the community health level, where distances to seek care are significantly shorter. Whether this relationship exists has significant implications for how CHW systems are designed geographically, and if and how community health programs should be optimized to the local geography.

Measuring the impact of geographic barriers on community health poses particular challenges. Studies on barriers to healthcare typically use data from national surveys, such as Demographic Health Surveys or Multiple Indicators Cluster Surveys (e.g. [17, 18]). Household coordinates in these surveys are randomly displaced to protect participants’ anonymity [19], preventing the assessment of geographic barriers at small spatial scales [20]. Dedicated surveys have measured households’ reported time to different health providers and the impact on healthcare access [21, 22], but these rarely assess the impact on CHW access (but see Ahmed et al. [23] for a review of the few, existing studies) and could be subject to biases in people’s perception of time. An alternative is the use of geographic information from routine health system data. While PHCs typically keep records of a patient’s residence to the lowest administrative level [9], this information is not systematically recorded in CHW registries. Studying the impact of geographic barriers at the community level therefore requires new methods and data if we are to optimize community health programs for the local geography of CHW catchments.

Madagascar is illustrative of the challenges to and opportunities for achieving universal health coverage by improving the design of community health systems. Following a history of French colonial occupation and political instability post-independence, Madagascar has witnessed little economic growth over the past half-century, severely limiting investment in public health, with negative consequences for healthcare access [24, 25]. There are fewer than 2 medical doctors and fewer than 3 nurses or midwives per 10,000 people in the country [26], and 43% of the population lives more than an hour walking time from a PHC, a common threshold for defining adequate access [4]. In contrast, the WHO recommended index for achieving universal health coverage is 44.5 doctors, nurses, or midwives per 10,000 people [27]. Nationally, health professionals cite common barriers to healthcare access seen in other settings such as the cost of medicines, geographic distance, mistrust of the health system, and unsanitary conditions at health facilities [28]. In addition, studies of healthcare access in rural districts of Madagascar have identified geographic distance to PHCs as a key barrier to primary care, even when the health system is strengthened and financial barriers have been reduced [29, 30]. Consequently, healthcare seeking rates for children under 5 years old are low, ranging from 40–60% nationally [31], and there is an opportunity for a strengthened community health program to help meet these needs. The national community health policy requires that two volunteer CHWs be locally selected and responsible for each fokontany, the smallest administrative area in Madagascar (mean population of approx. 1200) [32]. However, fokontany range in size from 1.5 km2 to 3747 km2. Given this range of catchment sizes, populations in rural Madagascar may continue to face geographic barriers when seeking care at the CHW level, which could limit the Malagasy health system’s ability to provide adequate coverage.

Here, we examined the impact of geographic distance to CHWs on the use of CHW services for children under 5 years in the rural district of Ifanadiana, southeastern Madagascar. Ifanadiana district has been established as a model district for universal health coverage by the government of Madagascar, implemented in partnership with the healthcare non-governmental organization (NGO) Pivot (see the ‘Study area and intervention’ section), with a particular focus on strengthening community healthcare. We build on previous work where we obtained accurate estimates of distance to PHCs and CHWs for every household in the district using a high-resolution spatial dataset that includes all households and transportation networks used by residents in the district [33]. From this dataset, we calculated an aggregate measure of geographic dispersion from the CHW site for each fokontany in our study area, and we assessed its relationship with all-cause child consultations with CHWs, controlling for relevant factors. We then explored two potential solutions to geographic barriers for isolated communities in Ifanadiana: optimizing a single CHW site location with regards to distance and installing an additional second site per catchment.

Methods

Study area and intervention

Madagascar is an island nation of approximately 587,041 km2 with a population of 28.1 million people located off the eastern coast of the African continent. Sixty-one percent of the population lives in rural areas [34], a mosaic of agriculture lands, forest, and savanna. The road network is sparse with only 10% characterized as passable year-round [35], and non-motorized transport by foot or bike is the norm. The majority of the population are agricultural producers, growing primarily rice, and 81% of the population lives below the international poverty line [34]. National government expenditure on healthcare is 19.75 USD, less than half the expenditure of neighboring countries [34], and many indicators of child health are below Sustainable Development Goal 2030 targets (Table 1).

Following national policy, each commune (subdivision of a district with ∼15 000 people) in Madagascar has at least one primary health center (PHC) that provides primary care for people of all ages. Each fokontany (subdivision of a commune with ∼1200 people) has two community appointed CHWs that provide integrated community case management (iCCM) [32]. Recent national policy recommends one CHW per 50 households in rural settings, but the policy has not yet been implemented [36]. Community health care, including diagnoses and distribution of limited medicines, takes place at the community health site of each fokontany, a static location, although more recent policies (2023 onward) recommend home visits. iCCM focuses primarily on malaria, diarrheal disease, and respiratory infections [36]. Those with severe symptoms or other diseases are referred to PHCs. CHWs are unpaid volunteers, making only marginal profits from the sale of certain medicines, and receive monthly supervision at PHCs.

Ifanadiana is a rural district located in the Vatovavy region in southeastern Madagascar. It has a population of approximately 200,000 people, distributed in 13 communes and 195 fokontany. Fokontany in the district range in size from 1.94–138.95 km2, with an average size of 20.36 km2 (S1 Table). The district is topographically complex, with areas of steep terrain, and landcover is a mosaic of rice fields and savanna. The transportation network consists of one paved road (<1% of the total road network in the district) that travels east-west and a secondary, unpaved road that follow a north-south axis [33]. The rest of the transportation network is tertiary roads or footpaths. Most of the population practices agricultural and nearly three-quarters live in extreme poverty [37].

In 2014, the Madagascar Ministry of Public Health (MMoPH) and the NGO Pivot began a partnership to strengthen the community health system as part of a broader intervention at all levels of care to achieve universal health coverage in Ifanadiana District (S2 Table). Public health facilities provide over 80% of healthcare in the district, primarily at PHCs, but utilization rates decrease by 50% for every increase in travel time of 1 hour to the PHC [9]. The strengthened community health program aimed to address this and revolved around an enhanced standard of care model which included providing CHWs with additional training; frequent on-site supervision by a team of nurses; the construction of a physical building for consultations in partnership with local communities; support to the supply chain, with medicines free of charge; and a modest stipend [32]. As of February 2022, this intervention had been implemented in seven communes, equivalent to about half the district, with plans to expand to the entire district by 2023.

Data collection

Distance to primary health centers and community healthcare workers.

For each fokontany, we estimated the average distance of all households in the fokontany to two types of health facilities, primary health centers (PHCs) and community health sites where CHWs are based. For this, we used a high-resolution spatial dataset available on OpenStreetMap, which was collected via participatory mapping (described in detail in Ihantamalala et al. [33]). This dataset included all paths, roads, buildings, and residential areas (defined as groups of four or more buildings) within Ifanadiana District. Overall, the dataset included over 100,000 buildings, 20,000 km of footpaths and more than 5000 residential areas. For each residential area, we estimated the network distance to the nearest health facility of each type using the Open Source Routing Machine algorithm (http://project-osrm.org/), along our mapped transport network of roads and paths. We then aggregated these values to the fokontany level by taking a weighted mean for all residential areas within a fokontany (residential areas were weighted by the number of buildings they comprised), resulting in an average household distance to the nearest PHC and community sites. We refer here to the average household distance to the fokontany’s community health site as geographic dispersion.

Survey data.

Household wealth is strongly associated with higher access to healthcare in Ifanadiana [38]. Therefore, we calculated the average wealth score for each fokontany using data from the IHOPE cohort, which includes longitudinal surveys at 80 clusters across Ifanadiana district, conducted every two years from 2014–2018 [37]. The sampling scheme was a two-stage cluster sampling design, where 40 clusters were sampled at random within each of two strata, within the Pivot initial catchment and the rest of the district. Twenty households per cluster were then randomly selected to be surveyed. The IHOPE cohort was implemented by the Madagascar National Institute of Statistics and was based primarily on the internationally validated Demographic and Health Surveys (DHS). Further details on the longitudinal survey can be found in Miller et al. [37] and Ezran et al. [39]. A wealth score was calculated for each of the 1600 households through a principal components analysis of household assets following standard DHS methods [40]. Because household wealth scores in Ifanadiana are very stable over time (S1 Fig), we summarized this value as the mean for each cluster across survey years. We then used inverse-distance weighting to interpolate these values from 80 clusters across the study area [41], using leave-one-out cross-validation to optimize the power used to weight the distance, via the spatstat package [42]. To aggregate this value to the scale of the fokontany, we extracted the wealth score at the coordinates of the four primary residential areas of each fokontany and assigned the mean of these values to the fokontany.

Community health system data.

To study the impact of geographic barriers on CHW utilization, we obtained the monthly number of CHW consultations for children under 5 years old from the CHW monthly reports to the District Office. All CHWs in these fokontany receive additional training and supervision from Pivot, including training on data reporting methods. In addition, mobile teams provide regular supervision and trained staff evaluate subsets of data for accuracy. These CHW reports, therefore, are a unique source of high quality data at the community level. We included data from 80 fokontany in seven communes in Ifanadiana, from May 2018—May 2021, except for one commune (Ambohimanga du Sud), where CHW support began in June 2020. Consultations here refer to all-cause consultations for children under 5 years old who are residents of that fokontany. Per national policy, CHW perform iCCM, and therefore only treat children under 5 years old who present with malaria, diarrheal disease, or respiratory infection. To obtain per capita utilization rates (per child under 5), we used population data collected by Pivot in 2021 as part of the proactive CHW program. We interpolated population estimates to other years assuming a 2.8% growth rate, following MMoPH methods.

In addition, we collected information to control for characteristics of the CHWs that may influence healthcare seeking behavior. This included the average age of the CHWs (median: 43.5, range: 20–63) and the gender of the CHWs in a community site (all men, all women, mixed). This data is routinely collected by the Pivot Community Health team. Because CHW composition was stable over the time period of this study (May 2018—May 2021), we assumed that these CHW characteristics were constant for the duration of our study. We also accounted for the introduction of a pilot initiative of proactive care in one of the communes (Ranomafana) in October 2019. Pivot teams also conduct monthly on-site supervisions of CHWs, observing multiple consultations and recording the percentage of cases of malaria, pneumonia, and diarrhea correctly managed (diagnosis, treatment, recommendations, etc.). From these measurements, we calculated the proportion of consults correctly treated across all three diseases as a metric of quality of care. This variable was used in a supplemental analysis to examine the impact of whether the CHW quality of care modified the relationship between our geographic proxies and healthcare access.

Potential solutions to geographic barriers

Following the statistical analysis estimating the relationship between geographic dispersion and consultation rates, we considered two solutions to geographic barriers to CHW access: optimizing the location of a single community health site and installing a second community health site per fokontany in an optimized location. We assessed whether relocating the community health site or optimizing the location of a second site to minimize geographic dispersion could substantially improve CHW consultation rates in the study area. To optimize the location of one CHW site, we identified the location within an existing residential area in that fokontany that resulted in the lowest geographic dispersion value. The optimal location was identified by calculating the geographic dispersion value for the fokontany from each residential area, following the method above, and selecting the residential area with the lowest value. To optimize the location of a second site, we calculated the geographic dispersion value for each fokontany with one site in the current location and the second site within an existing residential area. We assumed each household would use the community site closest to their residence within their fokontany, and calculated the average household distance to the closest community site for each potential new second location. The optimal second location was the residential area that resulted in the lowest average geographic dispersion value for the fokontany. We then used the fit statistical model to predict the change in consultation rates given new optimized site locations and reduced estimates of geographic dispersion, and compared predicted consultation rates at optimized locations to predicted rates at present locations.

We created an online dashboard to share the results of these optimization exercises for the full district of Ifanadiana (https://research.pivot-dashboard.org). This interactive cartographic tool allows public health stakeholders to identify current, optimal single, and optimal second CHW site locations for all 195 fokontany in Ifanadiana district. In addition, they can explore the predicted change in CHW consultation rates for both potential solutions for the 80 fokontany included in the statistical analysis. To ensure adoption of this tool by CHW supervisors at Pivot, we conducted a training on tool use following its establishment. The tool is available in English and French, the latter of which is spoken by CHW supervisors.

Results

Geographic barriers to community healthcare

CHWs conducted 73,857 consults from May 2018—May 2021, with an average annual consultation rate of 1.73 (1.51 sd) consultations per capita. These consultations are all-cause consultations for children under-5 seen by a CHW. Data was missing for 62 out of 2436 month-fokontany combinations, with one fokontany missing ten months of data, but the majority of fokontany (56) missing none. Geographic dispersion estimates ranged from 0.21 to 4.12 km, with a mean of 2.06 km (Fig 1). Geographic dispersion was a function of both the spatial pattern of residential zones relative to the community health site and the overall absolute size of the fokontany (Fig 1B and 1C). The fokontany with the lowest geographic dispersion had 93.8% of the population living within 1 km of the community health site, with a maximum distance of 1.65 km. In contrast, the fokontany with the highest geographic dispersion had 14.8% of the population living within 1 km of the community health site, with a maximum distance of 8.34 km.

Fig 1. Spatial patterns of geographic access to health care at community health sites in Ifanadiana.

A) Geographic dispersion of the fokontany population from the CHW (in km), based on the transportation network and location of all residential areas. Panels B) and C) show the routes from residential areas to CHW sites in a fokontany with high geographic dispersion (B) and low geographic dispersion (C). Blue points represent the community health site for that fokontany and shaded squares represent residential areas, with larger squares representing areas with more buildings. CHW: community health worker, PHC: primary health center. Data Sources: MMoPH (administrative boundaries), Pivot (health center locations), OpenStreetMap (transportation routes, buildings). All maps have been created from spatial datasets using the R programming language.


https://doi.org/10.1371/journal.pgph.0001028.g001

Our statistical model suggested that the geographic dispersion of a fokontany is associated with CHW consultation rates. Specifically, consultation rates were lower in fokontany where residential areas are more geographically dispersed from CHW sites (Table 2 and Fig 2A). An increase in geographic dispersion of 1 km corresponded to a decrease in monthly consultation rates of 28.1% (95% CI: 19.0% – 36.2%, Table 2). For example, while keeping all other variables unchanged and at their mean value, an increase in geographic dispersion from 1 km to 4 km corresponded to a predicted decrease in annual consultation rates from 2.18 to 0.81 per capita. We also found a trend for higher CHW consultation rates in fokontany that are further from PHCs (Table 2 and Fig 2B), but this trend was not significant (IRR = 1.019, 95% CI = 0.988–1.052). Wealthier fokontany had lower consultation rates than fokontany with lower socio-economic levels (IRR = 0.686, 95% CI = 0.533–0.882, Table 2). Consultation rates were higher in fokontany that had mixed gender CHWs (IRR = 1.514, 95% CI = 1.154–1.984) or all men teams (IRR = 1.423, 95% CI = 1.06–1.919), compared to teams of all women (Table 2). Our supplemental analysis found that these results were robust to the inclusion of the quality of care provided by CHWs in the model (S2 Fig).

thumbnail

Fig 2. CHW consultation rates for children under 5 years decrease with increasing CHW geographic dispersion and increase with increasing distance to the nearest PHC.

Lines represent the estimated marginal mean effects of each variable with 95% confidence intervals from the fitted model, holding all other variables constant at their mean value. Points represent the mean annual CHW consultation rate for each fokontany from May 2018—May 2021. CHW: community health worker, PHC: primary health center.


https://doi.org/10.1371/journal.pgph.0001028.g002

thumbnail

Table 2. Results of statistical model estimating effects of geographic barriers on CHW consultation rates.

Incidence rate ratios (IRR) of variables included in the model. Following Wasserstein et al. [48], we have chosen to provide 95% confidence intervals and not to dichotomize these continuous measures into thresholds of significance (e.g. p-values). CHW: community health worker, PHC: primary health center.


https://doi.org/10.1371/journal.pgph.0001028.t002

Impact of single optimized community health site placement on CHW consultations

Given the association between CHW geographic dispersion in a fokontany and CHW consultation rates, we explored whether optimized CHW site placement could be a potential solution to reduce geographic barriers to community health access. Forty-six fokontany (of 80) had an existing community site that was within 500m euclidean distance of the optimal site placement (Fig 3), resulting in an average reduction in geographic dispersion of 0.062 km. For this exercise, we considered these forty-six instances as fokontany where the CHW site was already optimally located. Of the 34 fokontany with an optimal site placement further than 500m from the existing site, four had a PHC located within their borders, and the optimal location was in the same residential area as the PHC. For the remaining 30 fokontany, the optimal CHW site placement was in a different geographic location than the current CHW site or PHC (Fig 3). The quartile of fokontany with the highest geographic dispersion (above 3.3 km) experienced an average reduction in geographic dispersion of 19.1% and absolute change in geographic dispersion of 0.71 km when calculating geographic dispersion using the optimal site placement compared to the current CHW site location. In contrast, the quartile of fokontany with the lowest geographic dispersion (below 0.90 km) experienced an average reduction in geographic dispersion of 15.3% and absolute reduction in geographic dispersion of 0.31 km when calculating geographic dispersion using the optimal site placement compared to the current CHW site location. In some instances, particularly for large fokontany, the distance between optimal and current site locations was several kilometers, but the change in geographic dispersion was small (Fig 3). This was because the size of the fokontany and road network necessarily placed a limit on the lowest geographic dispersion possible when only placing one CHW site.

thumbnail

Fig 3. The majority of existing CHW sites are located near the optimal location.

Absolute change in CHW geographic dispersion based on optimal CHW site location compared to existing CHW site location. The fill color represents the reduction in geographic dispersion and points represent existing site locations (gray) and optimal site locations (black). CHW: community health worker, PHC: primary health center. Data Sources: MMoPH (administrative boundaries), Pivot (health center locations), OpenStreetMap (transportation routes, buildings). All maps have been created from spatial datasets using the R programming language.


https://doi.org/10.1371/journal.pgph.0001028.g003

Using the fitted GLMM, the predicted change in annual consultation rates given optimal CHW site location was a modest increase of 0.13 consultations per year per child, or a 7.4% increase in annual consultation rates, across all 80 fokontany. Focusing on only the 30 fokontany with an optimal site placement more than 500m from the existing site, we predicted a 15.8% increase in consultation rates with optimized CHW site placement. For the quartile of fokontany that had the highest geographic dispersion, optimal CHW site placement resulted in a 15.2% increase in consultation rates compared with current CHW site placement (Fig 4). For the quartile of fokontany with the lowest geographic dispersion, optimal CHW site placement resulted in a 12.6% increase in consultation rates compared with current CHW site placement (Fig 4).

thumbnail

Fig 4. Both proposed solutions have the strongest predicted improvement in consultation rates in the most geographically dispersed fokontany.

The predicted change in annual CHW consultation rates following the optimization of a single CHW site (top row) or the addition of a second CHW site in an optimized location (bottom row). Each line represents the predicted change for one fokontany and fokontany are grouped into quartiles determined by their geographic dispersion estimates (Q1: least dispersed, Q4: most dispersed).


https://doi.org/10.1371/journal.pgph.0001028.g004

Impact of a second community health site on CHW consultations

The second solution we explored was the addition of a second community site to each fokontany in a location optimized to reduce the geographic dispersion of that fokontany. The optimal location of the second community site in four fokontany was the same location as a PHC (Fig 5), and so we did not include them in further analyses. The remaining 76 fokontany had a second site on average 2.29 km (0.98 km sd) from the existing site. On average, 42.2% of the population of each fokontany lived closer to the second site than the current site, representing the potential population served by a second community health site. The quartile of sites with highest geographic dispersion experienced an average reduction in geographic dispersion of 38.2%, or 1.27 km, following the addition of a second site (Fig 5). The quartile of sites with the lowest geographic dispersion experienced and average reduction in geographic dispersion of 48.8%, or 0.44 km, following the addition of a second site (Fig 5).

thumbnail

Fig 5. The addition of a second community health site reduces geographic dispersion by over 500m for the majority of fokontany.

Absolute change in CHW geographic dispersion following addition of a second site in an optimal location compared to current geographic dispersion. The fill color represents the reduction in geographic dispersion and points represent existing site locations (gray) and second site locations (black). CHW: community health worker, PHC: primary health center. Data Sources: MMoPH (administrative boundaries), Pivot (health center locations), OpenStreetMap (transportation routes, buildings). All maps have been created from spatial datasets using the R programming language.


https://doi.org/10.1371/journal.pgph.0001028.g005

After applying the fitted GLMM, we found that the addition of a second site was predicted to increase the average annual consultation rate from 1.81 consults per capita to 2.38 consultations per capita, an increase of over 30%, for all 76 fokontany. Considering the quartile of fokontany that had the highest geographic dispersion, the creation of a second site resulted in a predicted 49.1% mean increase in consultation rates (Fig 4). For the quartile of fokontany with the lowest geographic dispersion, the creation of a second site resulted in a predicted 27.3% mean increase in consultation rates (Fig 4).

Discussion

Geographic distance to health facilities is one of the main barriers to healthcare access, particularly in rural communities that often live several hours travel time from the nearest clinic. Community health workers (CHWs) have worked for decades within these isolated communities to provide basic health services, especially to children under five years old, and to increase geographic access to primary care. However, even at this finer spatial scale, distance to CHWs may continue to negatively impact healthcare access. A renewed push for community health has translated into updated international guidelines to optimize community health programs in order to maximize their impact, including recommendations to adapt programs to the local geography of CHW catchments. Using a comprehensive transport network dataset, we examined whether geographic barriers persist at the CHW level for child consultation rates in a rural district of Madagascar with a strengthened community health program. We found that consultation rates were lower in fokontany where households lived further from the CHW site (i.e. higher geographic dispersion), even after controlling for factors such as wealth and healthcare provider characteristics. Optimizing the location of these CHW sites to minimize geographic dispersion alone had minimal impact on predicted consultation rates. However, the addition of a second community site was predicted to increase consultation rates by nearly 50% in the most geographically dispersed fokontany. While geographic barriers clearly play a role in limiting healthcare access in this setting, our results suggest that optimizing the location of a single community site will not reduce geographic barriers. Rather, alternative solutions, such as installing additional sites in more dispersed fokontany or the implementation of proactive care programs, may be needed to significantly reduce geographic barriers to community health access.

We found that geographic barriers to healthcare access persist in Ifanadiana, even at the level of community healthcare. Other studies in Ifanadiana have found that, despite the installation of an intervention that has virtually removed user fees and improved the quality of care at all levels of the health system, geographic barriers persist and negatively impact hospital referrals [49], use of primary health centers [9], child vaccination coverage [50], and malaria case ascertainment [51]. In addition, results from the IHOPE cohort in Ifanadiana showed that, while financial inequalities in the intervention catchment had been significantly reduced between 2014 and 2018, geographic inequalities persisted for most coverage indicators evaluated [29]. This trend is not unique to Madagascar, and negative impacts of geographic accessibility on primary and secondary healthcare access are well studied (see Guagliardo [52] for a review). The findings of this study go one step further and show that geographic isolation remains a barrier for CHW access, where distances to seek care are much smaller. This contributes to a nascent, but growing literature studying the geography of community healthcare [53, 54], with the general consensus that geographic barriers also exist at this much finer spatial scale. Often, a threshold of 5km distance or 1 hour travel time is used to define spatial coverage of healthcare access (e.g. [4]), but we find that, for community care, a meaningful effect of distance on healthcare utilization remains even within that definition of “coverage”.

Our results have important implications for the design of national community health programs in Madagascar and elsewhere. The community health policy of Madagascar in recent years required two CHWs per fokontany regardless of their size or geographic dispersion of the population. In late 2022, national community health guidelines were updated to recommend one CHW for every 300 households in a fokontany [36], although this policy has not yet been implemented. However, the guidelines do not include recommendations for the distribution of CHWs or health sites with respect to spatial patterns in population densities or transport networks. Our results suggest that in addition to the recommended ratios, the location of CHWs should be adapted to the local geography of CHW catchments, and additional interventions may be necessary to reach isolated communities. Our results also support the need for proactive community case management programs (pro-CCM), which have been implemented elsewhere [55], and were recently recommended in the MMoPH’s community health policy. Under pro-CCM, CHWs visit each household within a catchment at regular intervals to provide basic services and conduct educational campaigns in addition to other elements, such as cost-free care, CHW professionalization and case follow-up regimes [56, 57]. Given the strong association between high geographical dispersion and low consultation rates, the implementation of pro-CCM programs in Madagascar has great potential for improving health outcomes.

Geographic optimization has been proposed as one potential solution to reduce barriers to healthcare [33, 53], but our optimization exercise predicts that optimizing the location of single community health sites would only marginally improve consultation rates. There are several potential explanations for this. First, the majority of fokontany already had CHW sites in geographically optimal locations. Pivot staff reports that CHW site locations in Ifanadiana are usually chosen by the fokontany mayor, and are often in one of the larger villages in the fokontany, if not the largest. Because of this, a significant proportion of the population by default lives within the same residential area as the CHW site, resulting in lower geographic dispersion. This implies that communities are able to identify optimal locations without the use of a complex geographic algorithm, and should be included in decision making regarding community health programs. Second, in Ifanadiana, communities in larger fokontany are necessarily further from community health sites (Fig 1). Given the larger area of these fokontany, there is a limit to how much optimizing placement can reduce geographic dispersion. Installing a second community health site in a geographically optimized location, however, resulted in much higher predicted increases in consultation rates than optimizing one site alone. This effect was greatest in the most dispersed fokontany, and such a solution could be targeted to large, dispersed fokontany where one site cannot adequately serve the whole population [58]. Another alternative solution is the creation of pro-CCMs, where CHWs visit all households in their catchment on a regular basis [32, 55]. Proactive community care programs have been associated with reductions in child mortality rates and increases in access to primary healthcare [55]. A pilot proactive community care program is currently being implemented in one commune in Ifanadiana district, with positive initial impacts on consultation rates and quality of care [32]. While the limited geographic scale of the program precludes formal statistical analysis, preliminary results show that this program is qualitatively increasing consultation rates in the more dispersed fokontany (S3 Fig).

In addition to geographic variables, we included several confounding variables describing characteristics of the population and CHWs in our model that warrant discussion. We found that CHW consultation rates were lowest for wealthier populations, counter to the general findings that wealthy populations have higher access to care [59]. Indeed, evidence for this positive relationship exists in Ifanadiana at the level of PHCs [29]. However, because we focused only on consultation rates with CHWs, we captured healthcare access for lower income populations that may not use PHCs because of financial barriers, and therefore rely more on CHWs for primary care. This suggests that CHWs are meeting their goal of serving vulnerable populations missed by other levels of the health system. In addition, we found that the gender make-up of CHWs influenced consultation rates, with consultation rates highest for all-men or mixed-gender teams. Historically, the majority of CHWs in Ifanadiana have been men, although there has been a focus on increasing gender equality among CHWs following international guidance from the WHO [60]. Following this guidance, there has been some concern that established gender norms in a field previously dominated by men could limit the effectiveness of women CHWs [6163]. Our finding that mixed gender teams have the highest consultation rates suggests that including CHWs of multiple genders at each site is a potential way to increase gender equality without negatively impacting healthcare seeking behaviors, as has been suggested elsewhere [61]. Additional study, particularly using qualitative methods, is needed to further explore patient perceptions of CHW gender in Ifanadiana.

Here, we advance new methods to study the impact of geographic barriers at the community health level, a topic that is currently under-researched due to intrinsic challenges in the spatial resolution of data at these finer scales. Our use of a comprehensive geographic dataset with all the transportation routes and buildings in the district, publicly available on OpenStreetMap, allowed us to examine spatial patterns of healthcare access within CHW catchments with remarkable accuracy. Ihantamalala et al. [33] found that these methods more accurately represented residents’ travel time than cost-distance maps often used to estimate potential access [4]. Recent studies of healthcare access [53, 64] have made use of fine-scale geographic data from OpenStreetMap or mobile phone data to obtain more accurate estimations of populations’ geographical isolation. Several of these data sources are open-source and available online. For example, OpenStreetMap is an open-source dataset of transportation networks that can be viewed via a website (https://www.openstreetmap.org/) or downloaded via the Overpass API (https://overpass-turbo.eu/). In addition, Maina et al. [65] collated the location of over 98,000 public health facilities from 50 sub-Saharan African countries, which is available through the WHO’s Global Malaria Programme (https://www.who.int/malaria/areas/surveillance/public-sector-health-facilities-ss-africa/en/). At the primary healthcare level, Weiss et al. [4] estimated travel times to these health facilities. Given the increasing availability of these data, future studies should leverage this fine-scale data to examine spatial patterns in healthcare inequalities, identifying underserved communities that may be otherwise missed in analyses that aggregate spatial data to coarser scales [66].

Our study had several limitations. CHW consultation data were reported at the level of the fokontany and individual visits are not georeferenced to each patient’s residence. Therefore, our conclusions are limited to inference at the scale of the fokontany, not the individual. Collecting patient-level data would allow for the consideration of more detailed spatial patterns of healthcare utilization within the fokontany, such as identifying a distance threshold for catchment populations where CHW utilization is lower than average. Additionally, we infer the wealth of 80 fokontany from surveys in a sample of 1600 households that was not necessarily representative of those fokontany. Nevertheless, the distribution of wealth in the eighty spatial clusters of the IHOPE cohort had a very consistent and strong spatial structure (e.g. higher wealth in areas near the paved road and around the larger towns) that justifies our spatial interpolation methods. Finally, we only examined one barrier to healthcare: geographic distance. There are many other factors that influence care-seeking behavior and access to care, including patient perceptions of care and treatment [67], costs of care [59], and CHW gender [68]. We partially addressed these other factors in our model by including variables for CHW age and gender and in the supplemental analysis that includes a measure of quality of care. However, the objective of this analysis was an in-depth consideration of geographic barriers, and future studies should consider these other factors in equal detail.

In conclusion, our results support WHO’s recent guidelines for community health programs, which suggest that the local geography of CHW catchments be taken into account in program design and implementation. We show that geographic barriers to care can persist at the community health level and we advance new methods to better understand these barriers in local settings. More research on geographic access to healthcare at multiple levels of the health system, particularly at finer spatial-scales, is needed to progress towards the collective goal of reaching “everyone, everywhere” included in international UHC policies.

References

  1. 1.
    Sacks E, Schleiff M, Were M, Chowdhury AM, Perry HB. Communities, universal health coverage and primary health care. Bull World Health Organ. 2020 Nov 1;98(11):773–80. pmid:33177774
  2. 2.
    O’Donnell O. Access to health care in developing countries: breaking down demand side barriers. Cad Saúde Pública. 2007 Dec;23:2820–34. pmid:18157324
  3. 3.
    Peters DH, Garg A, Bloom G, Walker DG, Brieger WR, Hafizur Rahman M. Poverty and Access to Health Care in Developing Countries. Ann N Y Acad Sci. 2008 Jul 25;1136(1):161–71. pmid:17954679
  4. 4.
    Weiss DJ, Nelson A, Vargas-Ruiz CA, Gligorić K, Bavadekar S, Gabrilovich E, et al. Global maps of travel time to healthcare facilities. Nat Med. 2020 Dec;26(12):1835–8. pmid:32989313
  5. 5.
    Khan AA. An integrated approach to measuring potential spatial access to health care services. Socioecon Plann Sci. 1992 Oct 1;26(4):275–87. pmid:10123094
  6. 6.
    Noor AM, Zurovac D, Hay SI, Ochola SA, Snow RW. Defining equity in physical access to clinical services using geographical information systems as part of malaria planning and monitoring in Kenya. Trop Med Int Health TM IH. 2003 Oct;8(10):917–26. pmid:14516303
  7. 7.
    Feikin DR, Nguyen LM, Adazu K, Ombok M, Audi A, Slutsker L, et al. The impact of distance of residence from a peripheral health facility on pediatric health utilisation in rural western Kenya. Trop Med Int Health TM IH. 2009 Jan;14(1):54–61. pmid:19021892
  8. 8.
    McLaren ZM, Ardington C, Leibbrandt M. Distance decay and persistent health care disparities in South Africa. BMC Health Serv Res. 2014 Nov 4;14(1):541. pmid:25367330
  9. 9.
    Garchitorena A, Ihantamalala FA, Révillion C, Cordier LF, Randriamihaja M, Razafinjato B, et al. Geographic barriers to achieving universal health coverage: evidence from rural Madagascar. Health Policy Plan. 2021 Nov 11;36(10):1659–70. pmid:34331066
  10. 10.
    Gabrysch S, Campbell OM. Still too far to walk: Literature review of the determinants of delivery service use. BMC Pregnancy Childbirth. 2009 Aug 11;9(1):34. pmid:19671156
  11. 11.
    O’Meara WP, Noor A, Gatakaa H, Tsofa B, McKenzie FE, Marsh K. The impact of primary health care on malaria morbidity—defining access by disease burden. Trop Med Int Health TM IH. 2009 Jan;14(1):29–35. pmid:19121148
  12. 12.
    Berman PA, Gwatkin DR, Burger SE. Community-based health workers: Head start or false start towards health for all? Soc Sci Med. 1987 Jan 1;25(5):443–59. pmid:3118477
  13. 13.
    World Health Organization, PEPFAR, UNAIDS. Task shifting: rational redistribution of tasks among health workforce teams: global recommendations and guidelines. World Health Organization; 2007.
  14. 14.
    Rifkin SB. Alma Ata after 40 years: Primary Health Care and Health for All—from consensus to complexity. BMJ Glob Health. 2018 Dec 1;3(Suppl 3):e001188. pmid:30622747
  15. 15.
    Henry B Perry, editor. Health for the people: National community health worker programs from Afghanistan to Zimbabwe. USAID MCHIP; 2020.
  16. 16.
    World Health Organization. WHO guideline on health policy and system support to optimize community health worker programs. Geneva; 2018.
  17. 17.
    Feldacker C, Emch M, Ennett S. The who and where of HIV in rural Malawi: Exploring the effects of person and place on individual HIV status. Health Place. 2010 Sep;16(5):996–1006. pmid:20598623
  18. 18.
    Lohela TJ, Campbell OMR, Gabrysch S. Distance to care, facility delivery and early neonatal mortality in Malawi and Zambia. PloS One. 2012;7(12):e52110. pmid:23300599
  19. 19.
    Rutstein S, Rojas G. Guide to DHS Statistics. Calverton, Maryland; 2006. (Demographic and Health Surveys Methodology).
  20. 20.
    Warren JL, Perez-Heydrich C, Burgert CR, Emch ME. Influence of Demographic and Health Survey Point Displacements on Distance-Based Analyses. Spat Demogr. 2016 Jul;4(2):155–73. pmid:27453935
  21. 21.
    Chayovan N, Hermalin AI, Knodel J. Measuring Accessibility to Family Planning Services in Rural Thailand. Stud Fam Plann. 1984;15(5):201–11. pmid:6495360
  22. 22.
    Thin Zaw PP, Liabsuetrakul T, Htay TT, McNeil E. Equity of access to reproductive health services among youths in resource-limited suburban communities of Mandalay City, Myanmar. BMC Health Serv Res. 2012 Dec 15;12:458. pmid:23241510
  23. 23.
    Ahmed S, Chase LE, Wagnild J, Akhter N, Sturridge S, Clarke A, et al. Community health workers and health equity in low- and middle-income countries: systematic review and recommendations for policy and practice. Int J Equity Health. 2022 Apr 11;21(1):49. pmid:35410258
  24. 24.
    Razafindrakoto M, Roubaud F, Wachsberger JM. Puzzle and Paradox: A Political Economy of Madagascar. 1st ed. Cambridge University Press; 2020.
  25. 25.
    World Health Organization. Global Health Expenditure Database [Internet]. 2022. Available from: https://apps.who.int/nha/database
  26. 26.
    World Health Organization. Global Health Workforce Statistics. Geneva, Switzerland; 2022.
  27. 27.
    World Health Organization. Health workforce requirements for universal health coverage and the sustainable developments goals. Geneva, Switzerland; 2016. (Human Resources for Health Obvserver Series). Report No.: 17.
  28. 28.
    Razafison R. Primary health care: back to basics in Madagascar. Bull World Health Organ. 2008 Jun 1;2008(6):421–2. pmid:18568266
  29. 29.
    Garchitorena A, Miller AC, Cordier LF, Randriamanambintsoa M, Razanadrakato HTR, Randriamihaja M, et al. District-level health system strengthening for universal health coverage: evidence from a longitudinal cohort study in rural Madagascar, 2014–2018. BMJ Glob Health. 2020 Dec 1;5(12):e003647. pmid:33272943
  30. 30.
    Sayre D, Steinhardt LC, Irinantenaina J, Dentinger C, Rasoanaivo TF, Kapesa L, et al. Baseline malaria prevalence and care-seeking behaviours in rural Madagascar prior to a trial to expand malaria community case management to all ages. Malar J. 2021 Oct 26;20(1):422. pmid:34702255
  31. 31.
    UNICEF. Madagascar Multiple Indicator Cluster Survey version 6 [Internet]. 2018. Available from: https://mics.unicef.org/
  32. 32.
    Razafinjato B, Rakotonirina L, Andriantahina JB, Cordier LF, Andriamihaja R, Rasoarivao A, et al. Evaluation of a novel approach to community health care delivery in Ifanadiana District, Madagascar. bioRxiv. 2020 Dec 14;2020.12.11.20232611.
  33. 33.
    Ihantamalala FA, Herbreteau V, Révillion C, Randriamihaja M, Commins J, Andréambeloson T, et al. Improving geographical accessibility modeling for operational use by local health actors. Int J Health Geogr. 2020 Jul 6;19(1):27. pmid:32631348
  34. 34.
    The World Bank. World Development Indicators [Internet]. 2022. Available from: https://data.worldbank.org/
  35. 35.
    World Food Program. Logistics Capacity Assessment [Internet]. 2022 [cited 2022 Oct 25]. Available from: https://dlca.logcluster.org/
  36. 36.
    Ministère de la Santé Publique. Guide de mise en oeuvre harmonisée du programme de santé communautaire à Madagascar. Antananarivo; 2022.
  37. 37.
    Miller AC, Ramananjato RH, Garchitorena A, Rabeza VR, Gikic D, Cripps A, et al. Baseline population health conditions ahead of a health system strengthening program in rural Madagascar. Glob Health Action. 2017 Jan 1;10(1):1329961. pmid:28621206
  38. 38.
    Garchitorena A, Miller AC, Cordier LF, Rabeza VR, Randriamanambintsoa M, Razanadrakato HTR, et al. Early changes in intervention coverage and mortality rates following the implementation of an integrated health system intervention in Madagascar. BMJ Glob Health. 2018 Jun 1;3(3):e000762. pmid:29915670
  39. 39.
    Ezran C, Bonds MH, Miller AC, Cordier LF, Haruna J, Mwanawabenea D, et al. Assessing trends in the content of maternal and child care following a health system strengthening initiative in rural Madagascar: A longitudinal cohort study. PLOS Med. 2019 Aug 20;16(8):e1002869. pmid:31430286
  40. 40.
    Rutstein S, Johnson K. The DHS Wealth Index. Calverton, Maryland: ORC Macro; 2004 Jan. (DHS Comparative Reports). Report No.: 6.
  41. 41.
    Shepard D. A two-dimensional interpolation function for irregularly-spaced data. In: Proceedings of the 1968 23rd ACM national conference. New York, NY, USA: Association for Computing Machinery; 1968. p. 517–24. (ACM ‘68).
  42. 42.
    Baddeley A, Rubak E, Turner R. Spatial Point Patterns: Methodology and Applications with R. London: Chapman and Hall/CRC Press; 2015.
  43. 43.
    Zuur A, Ieno EN, Walker N, Saveliev AA, Smith GM. Mixed Effects Models and Extensions in Ecology with R. Springer Science & Business Media; 2009. 579 p.
  44. 44.
    Brooks M E, Kristensen K, van Benthem K J, Magnusson A, Berg C W, Nielsen A, et al. glmmTMB Balances Speed and Flexibility Among Packages for Zero-inflated Generalized Linear Mixed Modeling. R J. 2017;9(2):378.
  45. 45.
    Hartig F. DHARMa: Residual Diagnostics for Hierarchical (Multi-Level / Mixed) Regression Models [Internet]. 2019. Available from: http://florianhartig.github.io/DHARMa/
  46. 46.
    R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021. Available from: https://www.R-project.org/
  47. 47.
    Morton B, Vercueil A, Masekela R, Heinz E, Reimer L, Saleh S, et al. Consensus statement on measures to promote equitable authorship in the publication of research from international partnerships. Anaesthesia. 2022;77(3):264–76. pmid:34647323
  48. 48.
    Wasserstein RL, Schirm AL, Lazar NA. Moving to a World Beyond “p < 0.05.” Am Stat. 2019 Mar 29;73(sup1):1–19.
  49. 49.
    Ihantamalala FA, Bonds MH, Randriamihaja M, Rakotonirina L, Herbreteau V, Révillion C, et al. Geographic barriers to establishing a successful hospital referral system in rural Madagascar. BMJ Glob Health. 2021 Dec 1;6(12):e007145. pmid:34880062
  50. 50.
    Rajaonarifara E, Bonds MH, Miller AC, Ihantamalala FA, Cordier L, Razafinjato B, et al. Impact of health system strengthening on delivery strategies to improve child immunisation coverage and inequalities in rural Madagascar. BMJ Glob Health. 2022 Jan 1;7(1):e006824. pmid:35012969
  51. 51.
    Hyde E, Bonds MH, Ihantamalala FA, Miller AC, Cordier LF, Razafinjato B, et al. Estimating the local spatio‐temporal distribution of malaria from routine health information systems in areas of low health care access and reporting. Int J Health Geogr. 2021 Feb 12;20(1):8. pmid:33579294
  52. 52.
    Guagliardo MF. Spatial accessibility of primary care: concepts, methods and challenges. Int J Health Geogr. 2004 Feb 26;3(1):3. pmid:14987337
  53. 53.
    Oliphant NP, Ray N, Bensaid K, Ouedraogo A, Gali AY, Habi O, et al. Optimising geographical accessibility to primary health care: a geospatial analysis of community health posts and community health workers in Niger. BMJ Glob Health. 2021 Jun 1;6(6):e005238. pmid:34099482
  54. 54.
    Oliphant NP, Sy Z, Koné B, Berthé M, Beebe M, Samake M, et al. Improving the efficiency of scale-up and deployment of community health workers in Mali: A geospatial analysis. PLOS Glob Public Health. 2022 Oct 19;2(10):e0000626.
  55. 55.
    Whidden C, Thwing J, Gutman J, Wohl E, Leyrat C, Kayentao K, et al. Proactive case detection of common childhood illnesses by community health workers: a systematic review. BMJ Glob Health. 2019;4(6):e001799. pmid:31908858
  56. 56.
    Johnson AD, Thiero O, Whidden C, Poudiougou B, Diakité D, Traoré F, et al. Proactive community case management and child survival in periurban Mali. BMJ Glob Health [Internet]. 2018;3(2). Available from: https://gh.bmj.com/content/3/2/e000634 pmid:29607100
  57. 57.
    Perry HB, Rassekh BM, Gupta S, Freeman PA. Comprehensive review of the evidence regarding the effectiveness of community-based primary health care in improving maternal, neonatal and child health: 7. shared characteristics of projects with evidence of long-term mortality impact. J Glob Health. 2017 Jun;7(1):010907. pmid:28685045
  58. 58.
    Brunie A, MacCarthy J, Mulligan B, Ribaira Y, Rabemanantsoa A, Rahantanirina L, et al. Practical Implications of Policy Guidelines: A GIS Model of the Deployment of Community Health Volunteers in Madagascar. Glob Health Sci Pract. 2020 Sep 30;8(3):466–77. pmid:33008858
  59. 59.
    Yates R. Universal health care and the removal of user fees. The Lancet. 2009 Jun 13;373(9680):2078–81. pmid:19362359
  60. 60.
    Dhatt R, Theobald S, Buzuzi S, Ros B, Vong S, Muraya K, et al. The role of women’s leadership and gender equity in leadership and health system strengthening. Glob Health Epidemiol Genomics. 2017;2. pmid:29868219
  61. 61.
    Feldhaus I, Silverman M, LeFevre AE, Mpembeni R, Mosha I, Chitama D, et al. Equally able, but unequally accepted: Gender differentials and experiences of community health volunteers promoting maternal, newborn, and child health in Morogoro Region, Tanzania. Int J Equity Health. 2015 Aug 25;14:70. pmid:26303909
  62. 62.
    Steege R, Taegtmeyer M, McCollum R, Hawkins K, Ormel H, Kok M, et al. How do gender relations affect the working lives of close to community health service providers? Empirical research, a review and conceptual framework. Soc Sci Med. 2018 Jul 1;209:1–13. pmid:29777956
  63. 63.
    Raven J, Wurie H, Baba A, Bah AJ, Dean L, Hawkins K, et al. Supporting community health workers in fragile settings from a gender perspective: a qualitative study. BMJ Open. 2022 Feb 1;12(2):e052577. pmid:35121601
  64. 64.
    Xia T, Song X, Zhang H, Song X, Kanasugi H, Shibasaki R. Measuring spatio-temporal accessibility to emergency medical services through big GPS data. Health Place. 2019 Mar 1;56:53–62. pmid:30703630
  65. 65.
    Maina J, Ouma PO, Macharia PM, Alegana VA, Mitto B, Fall IS, et al. A spatial database of health facilities managed by the public health sector in sub Saharan Africa. Sci Data. 2019 Jul 25;6(1):1–8.
  66. 66.
    Apparicio P, Gelb J, Dubé AS, Kingham S, Gauvin L, Robitaille É. The approaches to measuring the potential spatial access to urban health services revisited: distance types and aggregation-error issues. Int J Health Geogr. 2017 Aug 23;16(1):32. pmid:28830461
  67. 67.
    Faruqui N, Bernays S, Martiniuk A, Abimbola S, Arora R, Lowe J, et al. Access to care for childhood cancers in India: perspectives of health care providers and the implications for universal health coverage. BMC Public Health. 2020 Nov 3;20(1):1641. pmid:33143668
  68. 68.
    Viswanathan K, Hansen PM, Rahman MH, Steinhardt L, Edward A, Arwal SH, et al. Can community health workers increase coverage of reproductive health services? J Epidemiol Community Health. 2012 Oct 1;66(10):894–900. pmid:22068027



Source link

Leave a Reply

Your email address will not be published. Required fields are marked *