Translator Disclaimer
30 November 2020 Elevation and Child Linear Growth in Nepal
Author Affiliations +
Abstract

The relationship between elevation of residence and a child's linear growth was studied using data for 8824 children below the age of 5 years born between 2001 and 2016 at elevations ranging from 50 to 3200 m above sea level in Nepal. Multiple regression was used to measure the role of a variety of household and community factors in explaining the observed elevation effect. A negative association was found between elevation and linear growth that varied substantially across the sample but retained a significant marginal effect across model specifications. Controlling for household wealth, access to markets, indoor air quality, and a range of other factors associated with elevation, for each 1000 m gain in elevation, height for age z score (HAZ) declined by between 0.10 and 0.20 points for an average child, and by between 0.35 and 0.42 points for a child with the characteristics of those living at the highest elevations. Results underscore the potential developmental risks for children living at high elevations and call attention to factors that help to mitigate these risks.

Introduction

Background and motivation

Understanding the role of the physical environment in a child's nutrition, growth, and development is important for targeting public interventions and for assessing trade-offs and synergies associated with potential interventions. In Nepal, where a large percentage of the population lives a kilometer or more above sea level, the topography creates substantial challenges to providing infrastructure, delivering basic services, and improving agriculture. One empirical consequence of these development challenges is a strong observed relationship between elevation and rates of child stunting. High rates of stunting in the mountains is often interpreted as an artifact of factors associated with elevation. In this study, we investigated the relationship between elevation and linear growth directly, using height for age z scores (HAZs) for 8824 children below the age of 5 years born between 2001 and 2016 at elevations ranging from 50 to 3200 m above sea level (masl). We focused on HAZ because impaired growth resulting from long-term malnutrition undermines an individual's lifetime health and human capital accumulation, with implications for lifetime earnings and a country's long-term economic development (Almond et al 2018).

We used 3 waves of the Nepal Demographic and Health Survey (DHS). This allowed us to analyze the stature of a large number of children observed between 0 and 5 years of age, a particularly important period in a child's development. Leveraging the richness of the DHS data, we modeled a child's height at a point in time as a function of location, multiple interacting components of socioeconomic status, and health inputs, including exposure to indoor air pollution. This provided us the opportunity to limit the extent to which unobserved socioeconomic status confounds estimates of the magnitude of the relationship between elevation and linear growth. The impact of hypoxia on linear growth is a well-established question of importance in high-elevation studies (see, for example, Frisancho and Baker 1970; Hoff 1974; Frisancho et al 1975; Pawson 1977; Mueller et al 1978). We added to this and more recent literature in 5 important ways. First, although several studies have focused on birth outcomes, the first year of life, and growth among school-aged children, evidence regarding growth between birth and 5 years of age, a truly critical period of child development, remains sparse. Second, relatively few studies of Nepali children have been conducted. Third, a large and detailed database allowed us to perform a more comprehensive analysis than has been previously conducted. Fourth, we explored the impact of hypoxia in conjunction with information on indoor air pollution, a suspected confounder. Fifth, we examined multiple aspects of socioeconomic status, allowing us to more accurately remove many of these effects than has been possible in previous studies.

Related research

Although the terms altitude and elevation are used somewhat interchangeably in the published literature, in this paper we follow McVicar and Körner (2013) and use the term elevation to refer to the vertical distance between mean sea level and the point of household residence, rather than altitude, which denotes the vertical distance between the land surface and an object not in direct contact with that surface (eg an airplane). In past work, elevation was found to be a significant predictor of linear growth at birth, 6 months, and 1 year for Bolivian children living at 3600 masl, when compared to their counterparts living at 400 masl. Similar findings were later reported for a different Bolivian sample, as well as for children in the US mountain west and Tibet (Frisancho and Baker 1970; Haas et al 1982; Greksa et al 1985; Dang et al 2007; Zahran et al 2014; Wang et al 2016). Recent research has reinforced these findings in studies similar to ours. Studies using data from high-elevation areas in Tibet (Argnani et al 2008), Argentina (Roman et al 2015), and Ethiopia (Mohammed et al 2020) have shown a negative and statistically significant association between elevation and linear growth, with children living above 2000 masl facing increased odds of stunting of 28–29%, and children living at higher elevations (3000 masl or more) facing even greater risk.

Evidence suggests a possible biological component to these results, including slowed growth of the fetus during pregnancy due to hypoxia (Moore et al 2011). However, some research has shown no growth differences between affluent high-elevation Bolivian children and affluent low-elevation Guatemalan children, suggesting socioeconomic factors could account for the elevation–stunting relationship. Mortality rates are higher in high-elevation ecological zones of Nepal than at lower elevations, and among rural and low-wealth populations (Chin et al 2011), patterns that may be related to stunting and driven by better access to maternal and health services and food markets (Hotchkiss 2001; Chikhungu et al 2014; Mulmi et al 2016).

Among household-level factors that strongly vary with elevation, maternal education stands out as a particularly important predictor of linear growth (Bicego and Boerma 1993; Boyle et al 2006; Semba et al 2008). Similarly, community-level socioeconomic status modifies the effects of household-level socioeconomic status, which suggests that improvements in access to basic community resources may contribute to individual child health (Stinson 1982; Fotso and Kuate-Defo 2005). Recent work from Nepal supports this conjecture. As a number of studies have noted that while nutrition data, measured using HAZ, stunting rates, or other measures have improved for children below age 5 since 2000, gains have been largest for children from wealthier and better educated households (Dorsey et al 2017; Nepali et al 2019; Budhathoki et al 2020; Hanley-Cook et al 2020). In a recent study from Nepal, characteristics of children and households were found to explain most of the variance in height for age and weight for height, with relatively smaller but statistically significant contributions from community-level factors (Smith and Shively 2019). One confounding factor in the elevation–status–stunting chain may be indoor air quality. Evidence from Bangladesh and India indicates that the use of low-quality heating and cooking fuels is strongly and significantly associated with risk of stunting (Hong et al 2006; Mishra and Retherford 2007), and children living at higher elevations may face greater exposure to smoke and respiratory illness due to both household poverty and heating needs.

Methods

Data source and key measurements

Our data came from the 2006, 2011, and 2016 waves of the Nepal DHS (MOH et al 2007, 2012, 2017). The DHS is a nationally representative household survey, and the dataset combining these 3 waves included 9937 children below age 5. After removing children with incomplete sets of information needed for the analysis (under the assumption that data are missing at random), our final sample included 8824 children residing in locations bounded by 26°37′80.9″N; 80°05′81″E and 28°90′04″N; 88°15′68″E. The DHS sampling methodology for Nepal was consistent over this period. The DHS includes anthropometric data, information on household health, sanitation, and demographic conditions, district-level geographic references, and a community-level measure of elevation. Our observations constituted a pooled cross section of children under 5 years of age measured in a specific survey at a specific point in time. Our outcome variable was height for age z score (HAZ), a common measure of long-term nutritional status in children that is calculated by comparing a child's linear growth attainment at a particular age to that of the World Health Organization (WHO) reference population (WHO 1995). We expected the relationship between elevation and linear growth to be potentially nonlinear and mediated through other continuous variables, and therefore we used the continuous measure of HAZ as our dependent variable.

We relied on DHS-reported measurements of elevation, which were available for each survey cluster. This provided a reliable estimate of the elevation a child experienced because these clusters (ie small enumeration units roughly comparable to a community) are sufficiently small such that it is unlikely that a household's elevation would deviate markedly from the community value. We also accounted for several additional household-level factors. We used year fixed effects (for 2006 and 2011, holding 2016 as the reference year) in all regressions. We also included fixed effects for the 3 officially designated agroecological zones of Nepal. The designation of these 3 distinct zones—the Terai, the hills, and the mountains—is based on biophysical and climate characteristics. These vary from tropical in the Terai (17% of land area), through subtropical and cool, but snow-free, in the hills (68% of area), to alpine and seasonally cold in the mountains (15% of land area). These agroecological zone fixed effects are correlated but not collinear with elevation. For example, some hill communities are at higher elevations than low-elevation mountain communities. The designation also reflects economic, cultural, and political circumstances distinct from elevation. The Terai served as the reference zone. Coefficient estimates for these fixed effects are listed in the tables of results and noted as such.

We also included demographic characteristics, including indicators for the gender of the head of the household, child sex, child age with a squared term, urban residency, and several variables capturing aspects of healthcare utilization and sanitation. The DHS provides measures of maternal body mass index (BMI) and education, and a household wealth index constructed by DHS analysts by applying principal components analysis to a set of household assets and characteristics (Filmer and Prichett 2001; see  https://dhsprogram.com/topics/wealth-index/Wealth-Index-Construction.cfm). Maternal education can affect childhood health and nutrition directly. It also reflects parental background and earning potential. Similarly, a mother's BMI is a direct measure of her health. In addition, to the extent BMI reflects household food availability, dietary diversity, sanitation, or disease burden, inclusion of maternal BMI in the child's regression controls for these important factors that we cannot observe directly.

Finally, to separate the influence of elevation from that of isolation, we measured remoteness with 2 variables that we added to the DHS dataset: Euclidean distance (in km) to Kathmandu from the centroid of the district in which each child resided, and a quality-weighted measure of road density (in km/km2) based on Thapa and Shively (2018). These variables are particularly important because accounting for remoteness and isolation independent of elevation in our regressions reduces potential bias in our estimate of the elevation effect that would be introduced by omission of these separate but closely related factors. We have strong reason to believe that infrastructure development supports food security and linear growth for children, particularly in Nepal. For example, Shively (2017) found that additional infrastructure reduces children's vulnerability to negative rainfall shocks, and Shively and Thapa (2016) found that additional infrastructure moderates both the level and volatility of food prices, which matter for food security and dietary quality.

Empirical strategy

We measured the association between HAZ and elevation using unweighted ordinary least squares (OLS) regression models. We did not weight our regressions to match the national population because we were explicitly interested in the heterogeneity of elevation effects between individual children. Because relatively few children live at very high elevations, and because these children represent only about 8% of our sample, the relational importance of high elevation on stunting would be measured less precisely in a weighted regression that implicitly placed more emphasis on low-elevation populations.

While our regressions captured the extent to which the elevation association can be explained by competing hypotheses related to variables correlated with elevation in Nepal, there was one possible shortcoming in our analysis. To the extent family migration decisions systematically coincide with birth timing, our data on growth outcomes could reflect hidden bias due to the spatial sorting of children. In Nepal, internal migration flows tend to run from high to low elevations. As a result, bias in our estimates of elevation effects could arise if migration decisions depend on factors other than those we observed, and if children born at high elevations to nonmigrating mothers are systematically advantaged or disadvantaged compared to children born at lower elevations. According to census figures, the former increased from 9% in 1981 to 15% in 2011, and the latter increased from 7% in 1981 to 10% in 2011 (CBS 2003, 2014). In addition, migrants are primarily male. Furthermore, citing census data, Thapa et al (2019) reported that internal migration in Nepal largely occurs from the hills to the Terai. As a result, we do not expect our regression findings to suffer from a failure to account for endogenous internal migration. However, we acknowledge that international migration, which is also important (Thapa et al 2019), could play a role that we did not account for here. Unfortunately, available data do not allow us to study the migration issue.

Results

Descriptive statistics

We used a pooled dataset combining children from all available years. HAZ values have risen over time in Nepal, but each survey round contains substantial variation in HAZ and elevation. We summarized this variation in Figure 1, in which we plotted a locally weighted scatterplot smoothing (LOWESS) regression line of HAZ against elevation. Figure 1 illustrates the strong negative relationship between HAZ and elevation, although the wide confidence interval at high elevations indicates both a shortage of observations at the highest elevations and the fact that not all high-elevation children have low HAZ. We summarized all variables used in our regressions in Table 1. Table 2 contains a summary and synthesis of the regressions and the main findings. Primary regression results are presented in Table 3, and auxiliary regression results are presented in  Tables S1 (mred-40-03-07_s01.pdf) and  S2 (mred-40-03-07_s01.pdf) (Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1).

FIGURE 1

Height for age z score (HAZ) of children below 5 years of age by elevation. Locally weighted scatterplot smoothing (LOWESS) regression line is shown with bands indicating 95% confidence interval. Line plots a running mean least squares regression based on a standard tricube weighting function with a bandwidth of 0.8.

img-z3-2_R11.jpg

TABLE 1

Descriptive statistics for the sample.

img-z4-2_R11.gif

TABLE 2

Synthesis of findings.

img-z5-2_R11.gif

TABLE 3

Regression results for main models.

img-z6-2_R11.gif

Regression results

We began by fitting a set of baseline models in which HAZ was the dependent variable and elevation was the independent variable of interest. We also included several sets of variables, described below, that may mediate the HAZ–elevation relationship, and a set of control variables, that is, variables that would be expected to correlate with HAZ, and the omission of which from the regression would potentially result in a biased estimate of the marginal effect of elevation. In models 1–3 (Table 3), our independent variables were elevation and our control variables, with no interaction terms. In models 4–8 ( Table S1 (mred-40-03-07_s01.pdf), Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1), we report regressions that allowed interactions between elevation and controls. In models 9–10 ( Table S2 (mred-40-03-07_s01.pdf), Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1), we allowed the marginal effect of elevation on HAZ to vary across values of elevation, but we included no interactions. Unless otherwise specified, each control variable entered the model linearly and additively. Our control variables included indicators for year, urban residency, and agroecological zone, and household demographic variables, including individual sex and age (as a level and a square), an indicator for residence in a female-headed household, asset-based wealth index, maternal educational attainment, and mother's age at the time of birth. We also accounted for the presence of health practices that could affect HAZ, including whether the child was currently being breastfed, whether the child had ever received vitamin A supplementation or a vaccination, whether the household had access to clean water, and whether the household used biomass fuel indoors, as well as maternal BMI, as a proxy for short-term calorie availability, diet, and health environment. We also controlled for aspects of remoteness, including the distance to Kathmandu and road density. We treated educational attainment, wealth, maternal BMI, smoke fuel, and the remoteness variables as mediators, and we allowed them to interact with elevation in auxiliary regression specifications. The mediators also entered individually in the models with interactions.

Across all specifications, the estimated coefficient on elevation was negative and statistically different from zero. In terms of magnitude, a 1 km increase in elevation approximately offset the HAZ gain associated with having a mother with secondary education rather than primary education. For example, the coefficient on elevation was –0.219 in model 3 (Table 3), whereas the difference between the coefficients on secondary and primary educational attainments, 0.329 and 0.082, was 0.247. The robustness of the finding on the elevation coefficient to the inclusion of many household controls in model 2, and to the addition of remoteness controls in model 3, was the primary finding of interest in this series of models.

The results in Table 3 establish the role of elevation, as its estimated coefficient remained negative, large, and significant across specifications. Because the models summarized in  Table S1 (mred-40-03-07_s01.pdf) (Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1) included elevation interaction terms, we cannot directly interpret the coefficient on elevation in these regressions as the marginal effect, as the magnitude of elevation's effect on HAZ depended on additional coefficients. Instead, we present in Figures 24 the marginal effects of elevation on HAZ, evaluated at different values of the interacting variables. The marginal effect of elevation, its sign, and significance were robust to these changes. In many cases, the marginal effect of elevation on HAZ remained large in absolute value compared to other marginal effects in  Table S1 (mred-40-03-07_s01.pdf) (Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1). Figures 2 and 3 illustrate that the marginal effect of elevation dropped substantially in absolute value as household wealth and mother's BMI increased. Only at the very highest values of these mediator variables did the marginal effect of elevation lose significance. We saw a similar result regarding road density and distance to Kathmandu: High road density substantially lowered the magnitude of the marginal effect of elevation, although it remained significant throughout almost the entire range of road density (see Figure 4). We interpret these results as evidence that, while elevation effects do vary, in some cases quite substantially, over commonly proposed mediators, they remain robust to the combination of additional controls and interaction terms, and the relationship with HAZ remains strong and statistically significant.

FIGURE 2

Average marginal effect of elevation on height for age z score (HAZ) of children below 5 years of age, by mother's body mass index (BMI). Vertical bars indicate 95% confidence intervals for marginal effects.

img-z7-2_R11.jpg

FIGURE 3

Average marginal effect of elevation on height for age z score (HAZ) of children below 5 years of age, by household wealth. Vertical bars indicate 95% confidence intervals for marginal effects.

img-z7-4_R11.jpg

FIGURE 4

Average marginal effect of elevation on height for age z score (HAZ) of children below 5 years of age, by road density. Vertical bars indicate 95% confidence interval.

img-z7-6_R11.jpg

In  Table S2 (mred-40-03-07_s01.pdf) (Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1), we summarize 2 models in which we allowed the strength of the HAZ–elevation relationship to vary across elevation levels. In the first specification, we included a cubic polynomial on elevation, and in the second, we interacted elevation with binary variables taking values of 1 if elevation fell between 1000 and 2000, 2000 and 3000, or above 3000 masl. Results show that higher elevations are disproportionately associated with low HAZ values, beyond the linear trend we imposed in the previous models. However, these model results did not differ significantly from one another or from those summarized in  Table S1 (mred-40-03-07_s01.pdf) (ie F-tests failed to reject the null hypothesis of equal variance when comparing these models; see Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1). This means that, although the higher-order coefficients are individually significant, adding them jointly does not improve the model fit, indicating that a model that allows a closer nonlinear fit in the relationship between HAZ and elevation does not outperform statistically the linear form assumed in the models of  Table S1 (mred-40-03-07_s01.pdf) (Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1).

Discussion

Elevation maintained a significant and negative marginal effect across our specifications. In specifications with interaction terms, the average marginal effect maintained a relatively large absolute value. Our investigation focused on the extent to which various mediators or confounders can explain away this relationship: If household and community factors can reduce or remove the elevation effect, one might conclude that elevation does not matter per se. As the coefficients on the interaction effects in  Table S1 (mred-40-03-07_s01.pdf) (Supplemental material,  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1) suggest, and as Figures 24 show more clearly, household wealth, maternal BMI (a proxy for household food security and dietary diversity), and local infrastructure do alter this relationship substantially. At the very highest values of each of these mediators, the marginal effect of elevation loses its significance, but these figures show that at low (and routinely observed) levels of maternal BMI, wealth, and infrastructure, elevation has a significant and more negative marginal effect on HAZ than at higher values of the mediators. The indoor air pollution coefficients, including the interaction coefficient, were insignificant across all specifications, and therefore for the sake of space, we did not produce figures analogous to Figures 24 for those variables. Although previous studies from Nepal have found evidence of indoor pollution affecting stunting (eg Paudel et al 2013), it could be that exposure to indoor air pollution is so ubiquitous across the DHS sample, and such small proportions of households use nonbiomass fuels, that we were unable to detect a pollution–HAZ association. The effect of indoor air pollution may also be captured by wealth or education.

These results underscore the importance of maternal education and socioeconomic status in improving nutritional outcomes. This is consistent with recent work on nutrition in Nepal that included direct measures of healthcare access (Dorsey et al 2017; Nepali et al 2019; Budhathoki et al 2020; Hanley-Cook et al 2020). These studies, which did not account for elevation, nevertheless found strong negative associations between stunting and both maternal education and socioeconomic status. The results on vaccination and breastfeeding are unexpected and at odds with the study of Hanley-Cook et al (2020). The negative correlation observed here between vaccination and linear growth could be a selection problem, wherein children who consume more healthcare do so because they are sicker, but we do not have any evidence from field-level institutional details to say this with confidence. The negative sign for “still breastfeeding” could be an artifact of breastfeeding older children for whom it may not be nutritionally helpful, but we lack information to test this conjecture.

To further aid in the interpretation of our findings, we quantified the predicted effects with a thought experiment designed to isolate the elevation effect from other factors. Consider the possibility of moving a hypothetical child from a high-elevation to low-elevation setting while maintaining the household characteristics representative of a typical high-elevation household. To do this, we multiplied a hypothetical change in elevation by marginal effects of elevation in the interaction models summarized in Figures 24, evaluated at the mean values of mother's BMI, wealth index, and road density for children in the top and bottom deciles of the elevation distribution. These calculations suggest that, holding all else constant, the HAZ of a child at the 95th percentile of elevation would improve by 0.35–0.42 if relocated to the 5th percentile of elevation while maintaining the characteristics of the high-elevation household. In other words, the effect of a decrease in elevation of approximately 2000 m would result in an improvement in linear growth slightly greater than the improvement associated with having a mother with secondary education versus one with no education. For a somewhat less dramatic example, moving a child from the 90th percentile of elevation to the 50th percentile of elevation, a change of approximately 1200 m, would improve HAZ by 0.22–0.27, an improvement slightly smaller than the gains associated with having a mother with primary education rather than no education. Mean values of wealth and road density are much higher at low elevations than at high elevations, so the implied effects of a move from low to high elevations, maintaining original mediator values, are not symmetric. The wealth and road density interaction models suggest that moving a child to higher elevation—from the 5th percentile to the 95th percentile—would reduce the child's HAZ by 0.29 and 0.18, respectively, and that moving a representative child from the 50th to the 95th percentile of elevation would reduce HAZ by 0.18 and 0.11. The thought experiment of moving from high to low elevation underscores the costs high-elevation children face, but the asymmetry between hypothetical moves highlights the degree to which higher values of mediator variables weaken the elevation effect. This is important because it suggests that improvements in basic indicators of development in mountain settings, for example, in infrastructure and household wealth, could substantially mitigate the high-elevation growth drag for young children.

Elevation is strongly associated with restricted linear growth in particularly remote areas. However, even in places with average values of wealth, maternal BMI, and road density, the marginal effect of elevation (approximately –0.20) is statistically significant and nontrivial in magnitude. This suggests elevation matters on its own, because even the most convincing of the confounders we can add to the model does not eliminate the elevation effect entirely. Although we cannot fully rule out differences in agricultural production and diets as explanations of differences in HAZ across elevations, we believe the inclusion of maternal BMI in our models captures a substantial component of underlying household food access. For unobserved variation in agricultural conditions and food access to negate the elevation finding, the following would be required: (1) wealthier, high maternal BMI, and high-infrastructure households have better access to sufficient food and nutrition; (2) these improvements weaken the penalty imposed by living at high elevations; and (3) this effect absorbs any effect of the mediators it confounds. While our data do not allow testing of these suppositions directly, the first likely holds. We are more skeptical about the second and third, however: Diet enhancements can clearly improve HAZ, but it is not obvious that diets should matter more at high elevations than at low elevations, or that hidden variation in diets could completely account for mediators. Future research for Nepal and elsewhere might attempt to account for dietary factors directly. For example, in a recent study of Ethiopian children living at high elevation, Mohammed et al (2020) found a negative association between height and elevation while controlling for household dietary diversity, measured using 24 hour recall data on 7 food groups. This suggests that the elevation relationship is robust to the inclusion of dietary measures.

Micronutrient deficiencies arising from the micronutrient content of soils at high elevations could also contribute to this relationship, as zinc and iron deficiencies are correlated with stunting (Prasad 1991; Gibson et al 2008; Bevis 2015). Zinc deficiencies in diets can arise due to zinc deficiencies in soils, particularly in places where people rely on their own production and consume limited amounts of animal protein because the zinc content of plants comes directly from the soil in which they grow (Mayer et al 2007; Singh 2009). There is also a growing body of evidence showing that zinc deficiencies can hinder the absorption of iron, so that a zinc deficiency can contribute to iron deficiency, which can also slow growth (Graham et al 2012). In Nepal, soils at high elevations are generally deficient in zinc and iron (Andersen 2007). If soils in the Terai have more iron and zinc than those at higher elevations, children in subsistence households at higher elevations may have higher rates of stunting due to lack of iron and zinc in diets. When households at high elevation in Nepal purchase foods, they are likely purchasing foods grown at lower elevations, and to the extent that soils at lower elevations contain sufficient iron and zinc to support adequate levels in diets, this access to market-provided food could reduce the risk of stunting. Bevis et al (2019) reported evidence that low-zinc soils are associated with elevated rates of stunting in Nepal's Terai. If high-elevation soils are disproportionately prone to zinc deficiencies, this could be playing a role in our findings related to elevation. There is also evidence that boron toxicity, if present, could contribute to slow linear growth for some children living at high elevations (Hjelm et al 2019).

We cannot fully rule out differences in agricultural production and diets as explanations of the differences in HAZ across the range of elevation studied here. As we noted above, maternal BMI may capture a component of household diets, as it is at least partially correlated with underlying food access in the short to medium term. Our inability to measure these factors directly is a limitation. A larger sample of children from the mountain zone, which includes the highest elevations in Nepal, would also improve our study, because these children comprised only 8% of our sample. Many hill-zone children live at elevations similar to mountain-zone children, somewhat reducing this shortcoming. However, having more observations of children in the mountain zone would improve confidence in our findings.

Last, growth and development under hypoxic conditions could not be observed directly in this study; however, hypoxia has been shown to affect growth in animal models and should be considered as a factor in the linear growth differences observed at elevation. Farahani et al (2007) argued that chronic constant hypoxia (CCH) and chronic intermittent hypoxia (CIH) produce different body and organ growth patterns. In CCH, changes to liver and kidney size are proportionate to body size, while heart, lungs, and brain are spared or increase in size relative to the body. In animals, measurable changes due to hypoxia appear as quickly as 7 days after exposure; individual mice exposed to CIH in the study were able to catch up from the loss of growth, but those under CCH were not (Farahani et al 2007). Additionally, elevation may play a role in hypoxia inducible factor (HIF) signaling, impairing normal development in oxygen delivery tissues (Harrison et al 2015). This suggests that hypoxia-mediated regulation of growth may be an important aspect of normoxic development in all animals, including humans.

Conclusion

Arguments in favor of improving conditions in remote and very high-elevation mountain areas are often dismissed as cost-ineffective due to low population density and difficult access. Further, development funding in Nepal, and elsewhere, often focuses on metrics such as the number of individuals reached per unit of development investment, which can leave the remotest areas further marginalized. Cost-effectiveness is clearly an important consideration whenever budgets are constrained and trade-offs must be made, but often the full spectrum of payoffs from mountain development is not fully appreciated or accounted for. In this paper, we examined the relationship between child growth and elevation in Nepal, with an emphasis on evaluating the extent to which other factors, particularly remoteness, household resources, and indoor air pollution, mediate or explain the observed negative relationship.

We found a pernicious effect of elevation on child growth that was very robust to the inclusion of many variables and consisted of more than mere remoteness. This suggests that elevation is strongly associated with child growth in Nepal, and that some forms of development could help to ameliorate this effect, although defining and validating an exact mechanism that would let us make a credible causal case are outside the scope of this study. In this sample, elevation effects were mediated by household wealth, maternal BMI, and infrastructure. These are all amenable to public interventions. While elevation maintained a negative and significant sign at all but the highest values of road density, the magnitude of this marginal effect dropped substantially as density rose. Improvements in infrastructure provide a strong source of mediation, as even a moderate change in road density resulted in a large change in the marginal effect of elevation. Figure 4 suggests that if all children in our sample, regardless of their elevation, experienced the average road density observed among children in the bottom elevation decile, then those children living at the 90th percentile of elevation (approximately 2000 masl) would exhibit HAZ values approximately 0.25 points higher than those observed at lower elevations given their current levels of road density. At the highest levels of road density, the elevation effect would disappear entirely. The robustness of this HAZ–elevation relationship suggests a role for nutrition interventions specifically targeting children at high elevations, as they face more restricted linear growth than do their lower-elevation cohorts. Importantly, differences did not disappear when we accounted for key potential confounders, including household wealth, remoteness, and indoor air quality. We expect that researchers studying children living in similarly remote communities (as has been documented in, for example, Tibet, Peru, and Bolivia) would find similar patterns, both because economic circumstances in these similarly remote areas are sufficiently comparable, and because our socioeconomic controls reduced the extent to which socioeconomic differences between these locations could drive regression results. Our findings suggest a moderate reduction in linear growth for children at high elevations, and they are consistent with findings reported from other locations, including Tibet (Argnani et al 2008), Argentina (Roman et al 2015), and Ethiopia (Mohammed et al 2020). Future studies that combine samples would strengthen our confidence that elevation, and not unobserved country-specific factors, is driving this pattern.

ACKNOWLEDGMENTS

Financial support for this research was provided by the Feed the Future Innovation Lab for Nutrition, which is funded by the United States Agency for International Development (USAID). USAID provides financial support for the DHS Program, but it had no role in the collection, analysis, or interpretation of the data, and it played no role in the decision to submit the paper for publication. The opinions expressed herein are those of the authors and do not necessarily reflect the views of the sponsoring agency.

REFERENCES

1.

Almond D, Currie J, Duque V. 2018. Childhood circumstances and adult outcomes: Act II. Journal of Economic Literature 56:1360–1446. Google Scholar

2.

Andersen PA. 2007. Review of micronutrient problems in the cultivated soil of Nepal. Mountain Research and Development 27:331–335. Google Scholar

3.

Argnani L, Cogo A, Gualdi-Russo E. 2008. Growth and nutritional status of Tibetan children at high altitude. Collegium Anthropologicum 32:807–812. Google Scholar

4.

Bevis L. 2015. Soil-to-human mineral transmission with an emphasis on zinc, selenium, and iodine. Springer Science Reviews 3:77–96. Google Scholar

5.

Bevis L, Kim K, Guerena D. 2019. Soils and South Asian Stunting: Low Soil Zinc Availability Drives Child Stunting in Nepal. Selected paper prepared for presentation at the 2019 Agricultural & Applied Economics Association Annual Meeting, Atlanta, GA, 21–23 July 2019.  https://ageconsearch.umn.edu/record/290956/files/Abstracts_19_05_15_14_24_21_61__70_62_42_116_0.pdf; accessed on 30 December 2020. Google Scholar

6.

Bicego GT, Boerma JT. 1993. Maternal education and child survival: A comparative study of survey data from 17 countries. Social Science and Medicine 36:1207–1227. Google Scholar

7.

Boyle MH, Racine Y, Georgiades K, Snelling D, Hong S, Omariba W, Hurley P, Rao-Melacini P. 2006. The influence of economic development level, household wealth and maternal education on child health in developing world. Social Science and Medicine 63:2242–2254. Google Scholar

8.

Budhathoki SS, Bhandari A, Gurung R, Gurung A, Ashish KC. 2020. Stunting among under 5-year-olds in Nepal: Trends and risk factors. Maternal and Child Health Journal 24:39–47.  https://doi.org/10.1007/s10995-019-02817-1Google Scholar

9.

CBS [Central Bureau of Statistics]. 2003. Population Monograph of Nepal: Volume 1 (Population Dynamics). Kathmandu, Nepal: Government of Nepal, National Planning Commission Secretariat. Google Scholar

10.

CBS [Central Bureau of Statistics]. 2014. Population Monograph of Nepal: Volume 1 (Population Dynamics). Kathmandu, Nepal: Government of Nepal, National Planning Commission Secretariat. Google Scholar

11.

Chikhungu LC, Madise NJ, Padmada SS. 2014. How important are community characteristics in influencing children's nutritional status? Evidence from Malawi population-based household and community surveys. Health and Place 30:187–195. Google Scholar

12.

Chin B, Montana L, Basagana X. 2011. Spatial modeling of geographic inequalities in child mortality across Nepal. Health and Place 17:929–936. Google Scholar

13.

Dang S, Yan H, Yamamoto S. 2007. High altitude and early childhood growth retardation: New evidence from Tibet. European Journal of Clinical Nutrition 62:342–348. Google Scholar

14.

Dorsey JL, Manohar S, Neupane S, Shrestha B, Klemm RDW, West KP. 2017. Individual, household, and community level risk factors of stunting in children younger than 5 years: Findings from a national surveillance system in Nepal. Maternal and Child Nutrition 14:e12434.  https://doi.org/10.1111/mcn.12434Google Scholar

15.

Farahani R, Kanaan A, Gavrialov, A, Brunnert S, Douglas RM, Morcillo P, Haddad GG. 2007. Differential effects of chronic intermittent and chronic constant hypoxia on postnatal growth and development. Pediatric Pulmonology 43:20–28. Google Scholar

16.

Filmer D, Prichett L. 2001. Estimating wealth effects without expenditure data—or tears: An application to educational enrollments in states of India. Demography 38:115–132. Google Scholar

17.

Fotso JC, Kuate-Defo B. 2005. Socioeconomic inequalities in early childhood malnutrition and morbidity: Modification of the household-level effects by the community SES. Health and Place 11:205–225. Google Scholar

18.

Frisancho AR, Baker PT. 1970. Altitude and growth: A study of the patterns of physical growth of a high altitude Peruvian Quechua population. American Journal of Physical Anthropology 32:279–292. Google Scholar

19.

Frisancho AR, Velasquez T, Sanchez J. 1975. Possible adaptive significance of small body size in the attainment of aerobic capacity among high altitude Quechua natives. In: Watts E, Johnston FE, Lasker G, editors. Biosocial Interrelations in Population Adaptation. The Hague, The Netherlands: Mouton Publishers, pp 55–64. Google Scholar

20.

Gibson RS, Hess SY, Hotz C, Brown KH. 2008. Indicators of zinc status at the population level: A review of the evidence. British Journal of Nutrition 99(Supplement3):S14–S23. Google Scholar

21.

Graham RD, Knez M, Welch RM. 2012. How much nutritional iron deficiency in humans globally is due to an underlying zinc deficiency? Advances in Agronomy 115:1–40. Google Scholar

22.

Greksa LP, Spielvogel H, Caceres E. 1985. Effect of altitude on the physical growth of upper-class children of European ancestry. Annals of Human Biology 12:225–232. Google Scholar

23.

Haas JD, Moreno-Black G, Frongillo EA, Pabon J, Pareja G, Ybarnegaray J, Hurtado L. 1982. Altitude and infant growth in Bolivia: A longitudinal study. American Journal of Physical Anthropology 59:251–262. Google Scholar

24.

Hanley-Cook G, Alemayehu A, Dahal P, Chitekwe S, Rijal S, Bichha RP, Parajuli KR, Kolsteren P. 2020. Elucidating the sustained decline in under-three child linear growth faltering in Nepal, 1996–2016. Maternal and Child Nutrition 2020:e12982.  https://doi.org/10.1111/mcn.12982Google Scholar

25.

Harrison JF, Shingleton AW, Callier V. 2015. Stunted by developing in hypoxia: Linking comparative and model organism studies. Physiological and Biochemical Zoology 88:455–470. Google Scholar

26.

Hjelm C, Harari F, Vahter M. 2019. Pre- and postnatal environmental boron exposure and infant growth: Results from a mother–child cohort in northern Argentina. Environmental Research 171:60–68. Google Scholar

27.

Hoff C. 1974. Altitudinal variations in the physical growth and development of Peruvian Quechua. Homo 24:87–99. Google Scholar

28.

Hong R, Banta JE, Betancourt JA. 2006. Relationship between household wealth inequality and chronic childhood under-nutrition in Bangladesh. International Journal of Equity in Health 5:15.  https://doi.org/10.1186/1475-9276-5-15Google Scholar

29.

Hotchkiss DR. 2001. Expansion of rural health care and the use of maternal services in Nepal. Health and Place 7:39–45. Google Scholar

30.

Mayer AM, Latham M, Duxbury J, Hassan N, Frongillo E. 2007. A food-based approach to improving zinc nutrition through increasing the zinc content of rice in Bangladesh. Journal of Hunger and Environmental Nutrition 2(1):19–39. Google Scholar

31.

McVicar TR, Körner C. 2013. On the use of elevation, altitude, and height in the ecological and climatological literature. Oecologia 171:335–337.  https://doi.org/10.1007/s00442-012-2416-7Google Scholar

32.

Mishra V, Retherford R. 2007. Does biofuel smoke contribute to anaemia and stunting in early childhood? International Journal of Epidemiology 36:117–129. Google Scholar

33.

MOH [Ministry of Health, Nepal], New ERA, ICF. 2012. Nepal Demographic and Health Survey 2011: Key Indicators. Kathmandu, Nepal, and Calverton, MD: MOH, New ERA, and ICF International. Google Scholar

34.

MOH [Ministry of Health, Nepal], New ERA, ICF. 2017. Nepal Demographic and Health Survey 2016: Key Indicators. Kathmandu, Nepal, and Calverton, MD: MOH, New ERA, and ICF International. Google Scholar

35.

MOH [Ministry of Health, Nepal], New ERA, Macro International. 2007. Nepal Demographic and Health Survey 2006: Key Indicators. Kathmandu, Nepal: MOH, New ERA, and Macro International. Google Scholar

36.

Mohammed SH, Habtewold TD, Abdi DD, Alizadeh S, Larijani B, Esmaillzadeh A. 2020. The relationship between residential altitude and stunting: Evidence from >26,000 children living in highlands and lowlands of Ethiopia. British Journal of Nutrition 123(8):934–941.  https://doi.org/10.1017/S0007114519003453Google Scholar

37.

Moore LG, Charles SM, Julian CG. 2011. Humans at high altitude: Hypoxia and fetal growth. Respiratory Physiology and Neurobiology 178:181–190. Google Scholar

38.

Mueller WH, Schull VN, Schull WJ, Soto P, Rothhammer F. 1978. A multinational Andean genetic and health program: Growth and development in an hypoxic environment. Annals of Human Biology 5(4):329–352. Google Scholar

39.

Mulmi P, Block S, Shively S, Masters W. 2016. Climatic conditions and child height: Sex-specific vulnerability and the protective effects of sanitation and food markets in Nepal. Economics and Human Biology 23:63–75. Google Scholar

40.

Nepali S, Simkhada P, Davies I. 2019. Trends and inequalities in stunting in Nepal: A secondary data analysis of four Nepal demographic health surveys from 2001 to 2016. BMC Nutrition 5(19):1–10.  https://doi.org/10.1186/s40795-019-0283-xGoogle Scholar

41.

Paudel R, Pradhan B, Wagle R, Pahar D, Onta S. 2013. Risk factors for stunting among children: A community based case control study in Nepal. Kathmandu University Medical Journal 10:18–24. Google Scholar

42.

Pawson IG. 1977. The growth characteristics of populations of Tibetan origin in Nepal. American Journal of Physical Anthropology 47(3):357–501. Google Scholar

43.

Prasad A. 1991. Discovery of human zinc deficiency and studies in an experimental human model. American Journal of Clinical Nutrition 53:403–412. Google Scholar

44.

Roman EM, Bejarano IF, Alfaro EL, Abdo G, Dipierri JE. 2015. Geographical altitude, size, mass and body surface area in children (1–4 years) in the Province of Jujuy (Argentina). Annals of Human Biology 42:431–438. Google Scholar

45.

Semba RD, de Pee S, Sun K, Sari M, Akhter N, Bloem MW. 2008. Effect of parental formal education on risk of child stunting in Indonesia and Bangladesh: A cross-sectional study. The Lancet 371:322–328. Google Scholar

46.

Shively G. 2017. Infrastructure mitigates the sensitivity of child growth to local agriculture and rainfall in Nepal and Uganda. Proceedings of the National Academy of Science of the United States of America 114:903–908. Google Scholar

47.

Shively G, Thapa G. 2016. Markets, transportation infrastructure and food prices in Nepal. American Journal of Agricultural Economics 99:660–682. Google Scholar

48.

Singh MV. 2009. Micronutrient nutritional problems in soils of India and improvement for human and animal health. Indian Journal of Fertility 5(4):11–16. Google Scholar

49.

Smith T, Shively G. 2019. Multilevel analysis of individual, household, and community factors influencing child growth in Nepal. BMC Pediatrics 19:91–105. Google Scholar

50.

Stinson S. 1982. The effect of high altitude on the growth of children of high socioeconomic status in Bolivia. American Journal of Physical Anthropology 59:61–71. Google Scholar

51.

Thapa G, Shively G. 2018. A dose–response model of road development and child nutrition in Nepal. Research in Transportation Economics 70:112–124. Google Scholar

52.

Thapa NR, Adikhari S, Budhathoki PK. 2019. Influence of internal migration on the use of reproductive and maternal health services in Nepal: An analysis of the Nepal Demographic and Health Survey 2016. PLOS One 14(5):e0216587. Google Scholar

53.

Wang W, Liu F, Zhang Z, Zhang Y, Fan X, Liu R, Dang S. 2016. The growth pattern of Tibetan infants at high altitudes: A cohort study in rural Tibet region. Nature Scientific Reports 6:1–8. Google Scholar

54.

WHO [World Health Organization]. 1995. Physical Status: The Use and Interpretation of Anthropometry. Report of a WHO Expert Committee. Technical Report Series No. 854. Geneva, Switzerland: WHO. Google Scholar

55.

Zahran S, Breunig I, Link B, Snodgrass J, Wiler S. 2014. A quasi-experimental analysis of maternal altitude exposure and infant birth weight. American Journal of Public Health 104:S166–S174. Google Scholar

Appendices

Supplemental material

 TABLE S1 (mred-40-03-07_s01.pdf) Regression results for models 4–8 (with variable interactions).

 TABLE S2 (mred-40-03-07_s01.pdf) Regression results for nonlinear elevation specifications.

Found at:  https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1.S1.

© 2020 Shively et al. This open access article is licensed under a Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/). Please credit the authors and the full source.
Gerald Shively, Timothy Smith, and Megan Paskey "Elevation and Child Linear Growth in Nepal," Mountain Research and Development 40(3), R11-R20, (30 November 2020). https://doi.org/10.1659/MRD-JOURNAL-D-19-00063.1
Received: 1 April 2020; Accepted: 1 September 2020; Published: 30 November 2020
JOURNAL ARTICLE
PAGES


SHARE
ARTICLE IMPACT
Back to Top