Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Food insecurity in the Eastern Indo-Gangetic plain: Taking a closer look

  • Saumyadipta Pyne ,

    Contributed equally to this work with: Saumyadipta Pyne, Saurav Guha

    Roles Conceptualization, Investigation, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing

    spyne@ucsb.edu

    Affiliations Department of Statistics and Applied Probability, University of California Santa Barbara, Santa Barbara, California, United States of America, Health Analytics Network, Pittsburgh, Pennsylvania, United States of America

  • Saurav Guha ,

    Contributed equally to this work with: Saumyadipta Pyne, Saurav Guha

    Roles Data curation, Formal analysis, Investigation, Writing – original draft

    Affiliations Health Analytics Network, Pittsburgh, Pennsylvania, United States of America, ICAR-Indian Agricultural Statistics Research Institute, New Delhi, India

  • Sumonkanti Das,

    Roles Data curation, Project administration, Writing – review & editing

    Affiliation School of Demography, Australian National University, Canberra, Australia

  • Meghana Ray,

    Roles Data curation, Resources, Writing – original draft, Writing – review & editing

    Affiliation Health Analytics Network, Pittsburgh, Pennsylvania, United States of America

  • Hukum Chandra †

    † Deceased.

    Roles Conceptualization, Formal analysis, Investigation, Writing – original draft

    Affiliations Health Analytics Network, Pittsburgh, Pennsylvania, United States of America, ICAR-Indian Agricultural Statistics Research Institute, New Delhi, India

Abstract

Objective

Food security is an important policy issue in India. As India recently ranked 107th out of 121 countries in the 2022 Global Hunger Index, there is an urgent need to dissect, and gain insights into, such a major decline at the national level. However, the existing surveys, due to small sample sizes, cannot be used directly to produce reliable estimates at local administrative levels such as districts.

Design

The latest round of available data from the Household Consumer Expenditure Survey (HCES 2011–12) done by the National Sample Survey Office of India used stratified multi-stage random sampling with districts as strata, villages as first stage and households as second stage units.

Setting

Our Small Area Estimation approach estimated food insecurity prevalence, gap, and severity of each rural district of the Eastern Indo-Gangetic Plain (EIGP) region by modeling the HCES data, guided by local covariates from the 2011 Indian Population Census.

Participants

In HCES, 5915 (34429), 3310 (17534) and 3566 (15223) households (persons) were surveyed from the 71, 38 and 18 districts of the EIGP states of Uttar Pradesh, Bihar and West Bengal respectively.

Results

We estimated the district-specific food insecurity indicators, and mapped their local disparities over the EIGP region. By comparing food insecurity with indicators of climate vulnerability, poverty and crop diversity, we shortlisted the vulnerable districts in EIGP.

Conclusions

Our district-level estimates and maps can be effective for informed policy-making to build local resiliency and address systemic vulnerabilities where they matter most in the post-pandemic era.

Advances

Our study computed, for the Indian states in the EIGP region, the first area-level small area estimates of food insecurity as well as poverty over the past decade, and generated a ranked list of districts upon combining these data with measures of crop diversity and climatic vulnerability.

1. Introduction

Food security exists when all people, at all times, have physical and economic access to sufficient, safe and nutritious food that meets their dietary needs and food preferences for an active and healthy life, according to Food and Agriculture Organization (FAO) [1]. Conversely, food insecurity (FI) refers to a situation wherein members of society have “limited access to safe and healthy food” [2]. In the absence of food security, nutritional deficits may lead to stunting, wasting, and undernourishment that bring forth a sequalae of developmental, socioeconomic, and health consequences that impact societies in both short and long terms. In fact, it is believed that there can be no environmental security and sustainability without food security [3].

Over the past 50 years, gains in agricultural productivity from the Green Revolution have translated to increased food availability [4]. With about six-fold increase in food grain production from 50 million tonnes in 1950–1951 to about 291.95 million tonnes in 2019–2020, India has moved away from dependence on food aid to become a net food exporter country. Yet, it faces a complex challenge from multiple burdens of undernutrition, micronutrient deficiency, as well as obesity [5]. Recently, India ranked 107th among 121 countries in the 2022 Global Hunger Index (GHI) report (although this was disputed by its Government) [6]. Among Indian children of age 5 years and younger, 34.7% were stunted, 17.3% wasted, and 33.4% underweight; and among children of age 5–9 years, 21.9% were stunted, and 23% moderately-to-severely thin for their age [7]. More than half the women of reproductive age (15–49 years) in India are anemic while studies on long-term effects of early-life and pre-natal hunger among Indian populations are underway [8].

Research on climate change over the past decade has revealed the importance of climate pollutants (both short- and long-term) on monsoon patterns, precipitation changes, and increases in temperature [9]. A 15–30% decline in production of most cereals has been projected for Africa and South Asia [10]. In India, rice and wheat contribute to three-quarters of the country’s cereal production and their yields are susceptible to climate changes. Rainfed agriculture supports about 40% of India’s population and estimates of the past data show that changes in monsoon characteristics led to reduced rice yields in India by 1.7% during 1960–2002 [9]. For every 1-degree Celsius increase in temperature, loss of 3.7%-14.5% in India’s wheat yields was estimated. For rice, such estimates from multiple methods predict even larger temperature impact with an average reduction of 6.6 ± 3.8% per degree Celsius [11].

To study the complex interplay among socioeconomic conditions, agroecology, and climate change, and their combined effects on food security, few regions are as crucially important as India’s “breadbasket”, the Indo-Gangetic plain. The region comprises of a 2.5 million km2 fertile plain encompassing the northern regions of the Indian subcontinent. In particular, the Eastern Indo-Gangetic Plains (EIGP) region includes the states of Uttar Pradesh (UP), Bihar, and West Bengal (WB) in India, as well as parts of Nepal and Bangladesh. In India, EIGP consists of 39.27 million hectares and is home to 395.19 million people, i.e., 32.64% of India’s total population (2011 census). It is among the most densely populated (700–1200 persons/km2) regions in the world, and has high socioeconomic vulnerability.

EIGP is characterized by fertile soils with ample monsoon rainfall, continuous supply of surface and groundwater and a largely favorable climate that supports a predominantly rice-wheat cropping pattern. While UP and Bihar account for 32% and 5.76% of the country’s total wheat production respectively, WB, Bihar, and UP contribute respectively 13.26%, 7%, and 11.75% of its total rice production [12]. However, EIGP’s food security is potentially vulnerable to adverse effects of both anthropogenic and environmental factors. In 2012, the percentage of populations living below the national poverty line for UP, Bihar, and WB was 29.43, 33.74, and 19.98 respectively [13]. The rates of unemployment and seasonal out-migration are relatively high among its rural populations while land holdings are generally small.

EIGP’s environmental concerns include rising temperatures, high inter-annual variability of precipitation and frequent occurrence of adverse climatic events such as droughts, floods, and increasing cyclonic activity [14, 15]. Annual buildups of atmospheric pollutants in intensively farmed areas may have resulted in relative yield changes of -15% or greater in this region between 2006 and 2010 [16]. The Ganga basin also has severe groundwater contamination of Arsenic, which enters the food chain, especially through cultivation of rice [17]. In fact, disproportionately large contribution of rice production to resource use, greenhouse gases, and climate sensitivity (relative to its share of monsoon cereal calorie production) was observed in India [18]. Studies have noted the importance of addressing the impacts of climate change through solutions such as diet and crop diversification, improved farming technology, issues of governance, and other adaptations [19, 20].

Given the region’s diverse agroecological and biophysical conditions, availability of precise and timely disaggregate level statistics is essential for designing targeted policies to ensure food security in EIGP. In developing countries, however, the scarcity of reliable quantitative data can present a major challenge to policy-makers and researchers [21]. Such data, even when they exist, are often reported only at regional or state levels, and may be poorly correlated with local surveys [22]. In particular, FI presents a complex systemic challenge to policy planners, researchers, and government agencies at different administrative levels. To understand it closely, it is essential to disaggregate nationally-representative FI data for small geographic areas (e.g., districts) or demographic groups (e.g., gender-wise, social groups). However, in large-scale survey data, the sample sizes of such small areas may be very small or even zero. Small area estimation (SAE) provides us with an efficient approach to address this challenge with precision and rigor [23].

While there is no single universally accepted measure of food insecurity (FI), some classes of FI indicators are common in policy literature [24, 25]. A comparative performance analysis of the FAO estimates with anthropometric measures for nutritional assessment and estimates from household consumption expenditure surveys (HCESs) showed that the latter apparently has the most benefits among FI indicators [24]. HCESs can offer a less expensive yet more regularly-conducted means of detailed dietary intake assessment, especially for some developing countries. HCESs usually capture data on a detailed breakdown of household food expenditures, including food quantities and associated monetary transfers. This allows disaggregation of FI by geographic areas or socioeconomic groups, thus revealing any inequalities of food consumption within a country. A nationally representative HCES can reveal the share of households in the sample falling under a fixed nutritional threshold, and thus, provide an estimate of the prevalence of undernourishment in the country [24]. We captured this with the indicator FI Prevalence (FIP). Further, some indicators are more sensitive at certain severity levels [25] of FI, which are captured by the indicators FI Gap and Severity (FIG and FIS) and used in the present as well as past studies [26].

In rural India, FI is based on the average calorie intake of less than 2400 kcal per capita per day, a threshold based on the direct calorie intake method [27]. While small area estimates of FIP, FIG and FIS were recently computed for Bangladesh [26], no such FI data exist for the much larger, more populous and diverse EIGP region that lies within India. The HCES data are collected by the National Sample Survey Office (NSSO), Ministry of Statistics and Program Implementation, Government of India. However, researchers have noted the lack of publication of such data by the government after the 68th round of HCES in 2011–2012 [28]. However, while such nationwide surveys provide reliable state and national level estimates, they cannot directly produce reliable estimates at a lower administrative level, such as a district, due to small sample sizes. Therefore, they cannot reveal the existing local disparities in food security among the different populations living in hundreds of districts, especially in rural India.

The present study has several distinctive and innovative features. In this study, we used SAE to compute model-based district level estimates of FI prevalence (FIP), gap (FIG) and severity (FIS) in EIGP region covering the rural parts of the Indian states of UP, Bihar, and WB based on HCES 2011–2012 data and informed by local covariates. Previous small area studies [29, 30] have either estimated different outcome variables such as nutritional intake (protein, fat, etc.) or used area level models. In contrast, the present study is based on a more detailed small area model that allows us to estimate the measures of FI using individual level survey data by incorporating household level survey weight. Moreover, while the past studies [29, 30] have taken an isolated view of only one state, we have considered EIGP as a connected region known for its unique demographics, agroecology, and hydroclimatic risks. Thus, we included key regional parameters such as climatic vulnerability, crop diversity, poverty, and the disparities of FI across the Gangetic and non-Gangetic districts. This makes the present study distinctive in terms of its broader geographic scope, different model assumptions and outcome variables, and estimation of multiple additional district-specific indices.

The data used for SAE, and the SAE methodology of this study, are described in the next section. In Section 3, SAE results are discussed in detail while different maps shed light on the district-level disparities of food security in the EIGP region. Further, by comparison of FI with other indicators such as climate vulnerability, poverty, and crop diversity, we provide a shortlist of districts that could be useful for future studies. We end with concluding remarks in Section 4.

2. Materials and methods

2.1 SAE data

We used the latest round of available HCES in India, i.e., HCES 2011–2012, for rural areas of EIGP region, including UP, Bihar and WB. Usually, over the duration of a year, these surveys are conducted through interviews of a representative sample of households. The sampling design used is stratified multi-stage random sampling with districts as strata, villages as first stage units and households as second stage units. In this survey, a total of 12,791 households and 67,186 persons were surveyed from 127 districts of the EIGP region. In particular, 5915 (34429), 3310 (17534) and 3566 (15223) households (persons) were surveyed from 71, 38 and 18 districts of UP, Bihar and WB respectively. This study analyzed publicly available, deidentified and anonymized data. Data for analyses are within the manuscript and its Supporting Information file. The de-identified survey data are available from the National Sample Survey Office, Government of India (Website: http://mospi.nic.in/).

The 2011–2012 HCES provides information on the quantity and value of more than 142 food items with a reference period of last 30 days for a few food items and last 7 days for the remaining food items for each state separately for rural and urban areas. The state-level estimates of average calorie intake for UP, Bihar, WB and EIGP are 2200, 2242, 2199 and 2212 kcal per capita per day respectively. The quantities of food recorded as consumed by the household were converted into equivalent amounts of energy based on a nutrition chart, that provides the energy per unit of different foods in the Indian diet [31]. Note, our analyses and outputs cover all the rural districts surveyed in this HCES, some of which have since been split to form newer districts. Kolkata, an urban district and the capital city of WB, was excluded from the rural districts of EIGP.

The auxiliary variables taken from the 2011 Population Census of India [32] are available as counts at district level. We fit models with the direct survey estimates of FI indicators as the response variables and a selected set of 5 auxiliary variables as potential covariates: (1) Number of households (HH), (2) Proportion of scheduled caste population (SC), (3) Proportion of scheduled tribe population (ST), (4) Literacy rate (LR), and (5) Proportion of working population (WP). For our analysis, we used the average dietary intake in rural India to be 2400 kcal per capita per day, as recommended by the Ministry of Health and Family Welfare, Government of India [33].

2.2 District level poverty, crop diversity, and climate vulnerability

We used the same HCES data for estimating district-level measure of poverty given by a poverty index (PI) for UP, Bihar and WB as described in S1 File. We computed a district-specific entropy-based crop diversity index (CDI) based on land use statistics from the Ministry of Agriculture and Farmers’ Welfare, Government of India [34]. The CDI calculation is described in the S1 File. Recently, a report on the detailed national-scale assessment of climate vulnerability was released by the Department of Science & Technology (DST), Government of India; and the Swiss Agency for Development & Cooperation (SDC). The DST-SDC report [14] lists all districts in India ranked by their Climate Vulnerability Index (CVI), which we used for our analysis.

2.3 Measures of food insecurity

The target variable Y at the unit (household) level in the 2011–12 HCES survey data file is binary, corresponding to whether a household is “food insecure” (i.e., consumes less than 2400 Kcal per capita per day in rural India) or not. Let U denote the finite population of interest of size N partitioned into D disjoint small areas (here, districts), a sample s of size n is drawn from this population with a given survey design. Let the set of population units i.e., households in area i(i = 1, …, D) be denoted by Ui with known size Ni (i.e., the number of households in district i) such that and . Let si denote the collection of units making up the sample in area i such that and . Let yij denote the value of the variable of interest for unit j(j = 1, …, Ni) in area i. In this study, yij is the daily per capita calorie intake (averaged at the household level) for household j in district i.

2.4 FIP, FIG and FIS

The quantity that we are interested to estimate for area (i.e., district) i is the food insecurity indicator Fαi [35] defined as (1) where Fαij = (1 − yij/z)α I(yijz) and z is a preset threshold value of food insecurity (z = 2400 kcal in this study), I(yijz) indicator variable for a food insecure household j in district i and α is “sensitivity” parameter. The FI indicators FIP, FIG and FIS are calculated by setting α = 0, 1 and 2 respectively in [1].

2.5 SAE methodology

When unit level covariates are available, unit level SAE methodology dedicated to poverty estimation [3638] can be utilized to estimate food insecurity indicators. Since district-specific covariates extracted from full census data are available, area level SAE methods were implemented in this study. Fay–Herriot model [39] is a widely used area level method in SAE that assumes the availability of area-specific survey estimates, and that these follow a linear mixed model with area random effects. We used the Fay–Herriot model [39] to develop district level estimates of FIG and FIS. To compute model-based district level estimates of FIP, we modeled survey weighted area-specific proportions under a logistic linear mixed model (LLMM) [40]. In this study, a district (of EIGP) is considered as the small area of interest.

2.6 The direct estimator

Without using any auxiliary data, this estimator, denoted by Direct, of Fαi is , where wij denotes the survey weights associated with unit j in district i [21]. For α = 0, defines the survey weighted direct estimator of the proportion of food insecure households in district i.

The variance of the direct estimator is approximated by . Note, when a district i’s sample size is small, the sampling variance ψαi is large and hence the direct estimator could be unstable due to large standard error.

2.7 SAE of FIG and FIS

Let be the observed direct estimate of unobservable population parameter Fαi of variable of interest y for district i. Let xi be the k-vector of district level auxiliary variables, often obtained from various administrative and census records, related to the population parameter Fαi. The simple district-specific 2-stage Fay and Herriot model [39] is described as and . Alternatively, we can express it as (2)

Here β is a k-vector of unknown fixed effects parameters, uis are independently and identically distributed normal random errors with E(ui) = 0 and , and eαi’s are independent sampling errors normally distributed with E(eαi|Fαi) = 0, Var(eαi|Fαi) = ψαi. The two errors are independent of each other within and across districts. Usually, the sampling variance ψαi are assumed to closely approximate the unknown sampling variances, and in some cases, they are smoothed when they are too noisy due to very small sample size. The is unknown and needs to be estimated. Methods of estimating include maximum likelihood (ML) and restricted maximum likelihood (REML) under normality assumptions. Let denote the resulting estimator of and the empirical best linear unbiased estimator of β. Then under (2), the empirical best linear unbiased predictor (EBLUP) estimate of Fαi is (3)

Here, , where defines the shrinkage effect for district i. The mean squared error (MSE) estimation of EBLUP follows from Rao and Molina (2015) [23].

2.8 SAE of FIP

Estimation of small area proportions using such linear mixed models for continuous data, and derived with EBLUP, might be inconsistent in the sense that the prediction might not be within the [0,1] interval and a suitable transformation (for example, logarithm, square-root, arcsin) is needed for making the prediction within the [0–1] interval with the cost of bias correction for back transformation by utilizing the sampling error variances (noisy variances may lead to worse bias correction) [41, 42]. Thus, the small area proportions are assumed to follow a LLMM approach [23, 40] to produce model-based district level estimates of FIP. Let denote the sample count (e.g., number of food insecure households in a sample) in district i. If the sampling design is ignored and the number of food insecure households are independently and identically distributed within a district, then the sample count in district i can be assumed to follow a Binomial distribution with parameters ni and πi, i.e., Isi|vi ~ Bin(ni, πi), where πi is success probability [40]. Following Johnson et al. (2010) [43] and Chandra et al. (2011) [44], the model linking the probability πi with the covariates xi is the LLMM of form (4) i.e., . Here vi is the area-specific random effect that capture the area dissimilarities with and β is the k-vector of regression coefficients, often known as fixed effects parameters. Under (3), a plug-in empirical predictor of proportion (EPP) (e.g., FIP) F0i in district i is (5)

We used an iterative procedure that combines the Penalized Quasi-Likelihood estimation of β and v = (v1, …, vD)T with REML estimation of to estimate unknown parameters in (4) The MSE estimation of EPP is adopted from Chandra et al. (2011) [44]. Here, the survey weighted probability estimate for a district is modelled as a binomial proportion, with an “effective sample size” that equates the resulting binomial variance to the actual sampling variance of the survey weighted direct estimate for a district [40]. In model (4), the “actual sample size” and “actual sample count” were replaced with the “effective sample size” and “effective sample count” respectively to incorporate the sampling information. For details on model fitting, see S1 File.

2.9 Coefficient of variation

The coefficient of variation (CV) of a small area estimate is computed as , where is the estimate of standard error of and is the estimate of . It provides a measure of relative error and gives an indication of the precision of the model-based estimates when contrasted with the direct survey estimates. Therefore, estimates are more precise if they have smaller CVs. Conversely, the estimates with larger CVs are considered unreliable and unstable.

2.10 SAE validation and benchmarking

In SAE literature, two types of diagnostics are used: the model diagnostics, and the diagnostics for the small area estimates [45]. The former are used to verify the assumptions of the underlying model, i.e., how well the working model performed when it is fitted to data. The latter diagnostics are used to indicate validity and reliability of the small area estimates. We applied 3 diagnostics to our model: (a) the bias diagnostic, (b) the percent coefficient of variation (CV) diagnostic, and (c) the 95% confidence interval diagnostic. Moreover, we benchmarked our small area estimates against the direct estimates at a higher (i.e., state) level. See S1 File for details on diagnostics and benchmarks.

3. Results

A generalized linear model for FIP, and linear models for FIG and FIS, were initially fit with 5 auxiliary variables (see Methods). Akaike Information Criterion (AIC) was used for optimal model selection. We identified 2 of these, the proportion of Literacy rate (LR) and the proportion of working population (WP), which contribute significantly to district-specific FI in EIGP. The selected optimal model for FIP as outcome had significant (p<0.01) negative influences from LR and WP with regression coefficients -0.6800 and -2.9125 respectively. The selected optimal models for FIG and FIS as outcomes also had significant (p<0.05) influences from exactly the same 2 auxiliary variables. LR and WP had their respective regression coefficients as -0.0640 and -0.2706 for FIG; and -0.0195 and -0.0803 for FIS. Thus, LR and WP were used as covariates in our model-based SAE analysis. As noted in the Methods, we used empirical predictor of proportion (EPP) for SAE of FIP, and empirical best linear unbiased predictor (EBLUP) for SAE of FIG and FIS. Using a Wald statistic (W), we tested the model’s Goodness of Fit [45]. The small values of W = 5.68, 15.42, and 25.83 for FIP, FIG, and FIS respectively are significant at 5% level which indicates the consistency between the direct and the model-based SAE results. See S1 File for further details.

As is common for SAE studies, we conducted two types of diagnostics: (a) the model diagnostics, and (b) the diagnostics for the small area estimates [45]. The model diagnostics in S1 Fig show that the normality assumptions were satisfied. The bias diagnostics in S2 Fig indicate that the model-based estimates of FIP, FIG and FIS are less extreme when compared to the corresponding direct estimates, thus demonstrating that the typical SAE outcome shrinks the extreme values towards the average (see S1 File).

The district-specific estimates of the FIP, FIG and FIS along with their CVs and 95% CIs were computed by the direct and SAE methods (see Table 1). The S1 Table lists FIP, FIG and FIS estimates for all 127 rural districts of EIGP. We note in Table 1 that FIP of Bihar (0.62) and UP (0.62) are slightly lower than WB (0.70). Further, the value of FIS (0.03) is the same for each state whereas in WB, FIG (0.13) is marginally higher than FIG (0.11) in both UP and Bihar. The state level estimates shown in Table 1 are reasonably accurate in terms of the CVs. The S2 Table and S3 Fig summarize the district level distribution of sample size (i.e., number of households in sample), the number of food insecure households in sample (sample count) and the sampling fraction in 2011–2012 HCES.

thumbnail
Table 1. The estimate and percent coefficient of variation (CV) of FIP, FIG and FIS in EIGP region and the 3 states.

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

In Fig 1, we mapped the spatial distributions of FIP, FIG and FIS in UP, Bihar and WB based on our model-based small area estimates. The darker areas on the maps show districts with higher FI. Fig 2 presents the district-wise values of percentage CV for model-based estimates and direct estimates of (a) FIP, (b) FIG and (c) FIS in increasing order of sample sizes. In most districts, the CVs of the model-based estimates are significantly smaller than those of direct survey estimates, thus demonstrating higher precision of the former. Further, such improvement CV is higher for districts with smaller sample sizes as compared to larger sample sizes.

thumbnail
Fig 1. District-wise maps showing the spatial distribution of food insecurity prevalence (left panel), gap (center) and severity (right) generated by SAE method for the states Uttar Pradesh (top), Bihar (middle), and West Bengal (bottom) in the EIGP region.

Higher values of food insecurity indicators are shown in darker colors. (Base map source: https://gadm.org/download_country.html).

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

thumbnail
Fig 2. District-wise coefficient of variation (%CV) for the model-based small area estimates (solid line) and the direct estimates (dashed line) of (a) FIP, (b) FIG and (c) FIS.

Districts are arranged from left to right in increasing order of sample size.

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

Moreover, a series of diagnostic measures clearly demonstrate that the model-based estimates of FI indicators are more precise than the direct estimates (see S1 File). The district-wise maps of the 95% confidence intervals (CIs) for FIP, FIG and FIS produced by the direct and SAE methods are displayed in S3 Fig. Clearly, the 95% CIs for the direct estimates are wider than the 95% CIs for the model-based estimates. We note that the 95% CI for direct estimates are invalid (for example, values greater than 1.0 for FIP and negative for FIS, etc.) in many districts due to large standard errors.

A set of summary statistics for the direct and model-based estimates along with associated standard errors and CV of the FI indicators (FIP, FIG and FIS) for 127 districts are reported in Table 2. As expected, the averages of the model-based estimates of FI indicators are almost identical to those of the direct estimates but with lower variation. For example, the maximum standard error of FIP estimates generated by the direct and the SAE methods are 0.152 and 0.105 respectively. Moreover, we calculated the ratio of the CVs of the direct and the corresponding model-based estimates for each of FIP, FIG and FIS. The boxplots of these ratios in S4 Fig show that the model-based estimates have CVs that are smaller than those of the direct estimates, and hence are relatively more precise.

thumbnail
Table 2. Summary statistics of direct and model-based estimates along with their standard error (SE) and percent coefficient of variation (CV) of FIP, FIG and FIS.

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

For benchmarking the reliability of the small area estimates, we aggregated and compared them with estimates at higher spatial levels. Table 3 shows that aggregation of the district-level estimates closely matches the EIGP region and state level estimates of the same FI indicators. The internal benchmarking performance indicated that the model-based estimates do not need to be validated further. An overall geographical distribution of prevalence of FI in the EIGP region, including data from the neighboring country of Bangladesh [26], is interpolated using the estimates and represented as a contour map in S6 Fig.

thumbnail
Table 3. Aggregated level estimates of FIP generated by direct and model-based SAE methods in different states and EIGP region.

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

The SAE results clearly show unequal distribution of FI across the EIGP region (Fig 1). In particular, FIP, FIG and FIS values range from 26.4–90.3% (average 63.5%), 3.6–20.2% (average 11.5%) and 0.5–5.5% (average 2.9%) respectively across the districts of the region. In fact, the intra-state disparities across districts are also noteworthy. FIP across districts in UP, Bihar and WB vary from 47.6–85.7% (62.8%), 27–90.3% (61.2%) and 46.7–84.6% (71%) respectively, where the average is shown in parentheses. The corresponding ranges and averages for FIG are 6.2–20.2% (11.6%), 4.5–16% (10.6%) and 6.8–19.2% (13.2%); and for FIS 1.1–5.5% (2.9%), 0.9–4.3% (2.5%) and 1.6–5.1% (3.4%) in UP, Bihar and WB respectively. We note that the range of disparities in FIP is the highest in Bihar while those for FIG and FIS are the highest in UP.

Interestingly, we observed certain spatial patterns when we focused on the map of each state. In UP, the results indicate an east-west divide in the distribution of FIP, FIG and FIS. For example, in the western part of UP there are many districts with low levels of FIP, FIG and FIS. In contrast, the eastern part of UP has several districts with high rate of FIP, FIG and FIS. The districts with higher percentage of FIP (> 75%) are Mainpuri, Basti, Pratapgarh, Gonda, Kaushambi, Siddharth Nagar, Agra, Sant Ravidas Nagar (currently called Bhadohi), Deoria, Pilibhit, Sant Kabir Nagar, Kannauj, Ghaziabad, and Gorakhpur. In general, FIG and FIS status of these districts are also high except in few cases. For example, although Pilibhit, Ghaziabad and Deoria have high FIP but values of FIG and FIS index of these districts are nominal between 0.14–0.15 and 0.031–0.035 respectively. On the other hand, the UP districts with lower percentage of FIP (< 40%) are Jhansi, Lalitpur, Etawah, Mahoba and Meerut.

In Bihar, except for Sheikhpura (27.5%), FIP exceeds 40% in all other districts. However, only 3 districts—Siwan, Katihar, and Sheohar—have FIP > 80%. Further, high level of FIP exists in many districts from its northern part bordering with Nepal and its western part sharing border with eastern UP, e.g., West Champaran, Sitamarhi and Madhubani in northern Bihar; and Siwan, Buxar and Gopalganj in western Bihar. In WB, as in Bihar, FIP is greater than 50% in all districts except for Howrah (46.7%), which is across the river Ganga from the capital city of Kolkata. Districts with the highest FIP (>80%) are Jalpaiguri, Murshidabad, Nadia, and Hooghly. Notably, Jalpaiguri, in northern WB, which borders Bhutan and Bangladesh in the north and south respectively, and Murshidabad and Nadia in eastern WB, which border Bangladesh, have the highest values of FIP, FIG and FIS in that state. The least vulnerable districts are Howrah, Bankura and Bardhaman (currently split into Paschim Bardhaman and Purba Bardhaman districts). In WB, FI is generally concentrated in the northern and the eastern parts of the state.

Upon estimation of the district-specific FI indicators, we compared them with the CVI and PI of the same districts. As FIP is relatively independent of CVI and PI (Spearman’s ρ = 0.1 and -0.16 respectively), we depicted FIP using district-wise “bivariate” maps (Fig 3) for comparison with CVI (Fig 3a) and PI (Fig 3b). The red (high FIP) and blue (high CVI) districts exhibited spatial autocorrelation in Fig 3a. Further, we distinguished the districts into Gangetic (i.e., those located on the banks of the river Ganga as well as its major distributaries, Bhagirathi and Hooghly, in WB) and non-Gangetic groups, and compared the distributions of FIP across these 2 groups (Fig 4). While the distributions are similar in the overall region, UP, and WB, interestingly, they vary in the state of Bihar where the median FIP is higher for the Gangetic districts than the non-Gangetic ones.

thumbnail
Fig 3. Maps of EIGP showing each district’s FIP against (a) CVI and (b) PI.

A bivariate color key is used to depict the range of values for each pair of indices. (Base map source: https://gadm.org/download_country.html).

https://doi.org/10.1371/journal.pone.0279414.g003

thumbnail
Fig 4. Violinplot for comparing the FIP distributions of the Gangetic (G) and the non-Gangetic (NG) districts of EIGP region and the 3 states.

The former (G) are shown in blue and the latter (NG) in orange.

https://doi.org/10.1371/journal.pone.0279414.g004

Notably, we identified several districts that have relatively high values of both food insecurity as well as CVI (shaded black in Fig 3a). In fact, some of these also have high PI or low CDI. The resulting shortlist of 22 districts, as shown in Table 4, could be used for further investigation as well as prioritizing the need to address any systemic challenges in such areas via targeted policies that are suited to those conditions (e.g., cropped area, poverty, Gangetic location or otherwise) and requirements of concerned populations.

thumbnail
Table 4. A list of 22 districts in EIGP and the state, sample size (Sample), Food Insecurity Prevalence (FIP), Climate Vulnerability Index (CVI), Poverty Index (PI), Cropped Area (in hectares), and Crop Diversity Index (CDI).

Group denotes whether a district is Gangetic (G) or non-Gangetic (NG). Other abbreviations: Uttar Pradesh (UP), WB (West Bengal).

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

4. Discussion

The U.S. Department of Agriculture describes FI as a situation of “limited or uncertain availability of nutritionally adequate and safe foods or limited or uncertain ability to acquire acceptable foods in socially acceptable ways” [46]. This definition draws our attention to the fact that FI is more than just hunger, as is recognized by the 2030 sustainable development goals (SDGs) of the United Nations, and, in particular, SDG2, which pertains to ending hunger, improving food security and nutrition, and promoting sustainable agriculture [47]. Globally, almost 2 billion people face some form of FI. As per the World Food Programme, 135 million people suffer from acute hunger, and the COVID-19 pandemic is expected to double that number [48].

India has the natural resource capacity to achieve nutritional security, reduce health risks and improve environmental sustainability [49]. While the Green Revolution has substantially increased India’s food production over the past half century, its other more mixed outcomes include homogenization of cereal production, unsustainable resource use especially where agroecological conditions are not well-suited, and rising vulnerability to climate variations [20]. The 2022 GHI report specifically notes that “India shows the importance of considering the subnational context when designing programs and policies to target child stunting.” [6]. Yet, in the presence of its inter-regional disparities, India currently lacks the needed disaggregate level data on food security. Towards this, our study used a model-based SAE approach to generate reliable and representative direct and model-based estimates. The resulting maps illustrate unequal spatial distributions of FI prevalence, gap and severity among the populations in different districts of UP, Bihar and WB in EIGP region.

Our models were informed by various district-specific covariates, among which the literacy rate (LR) and the proportion of working population (WP) had significant negative influences respectively on a district’s food insecurity. Interestingly, while the share of LR in the EIGP is lower than the national average, other recent, independent studies on malnutrition in India also provide evidence in the same direction [50]. In small area literature [23], two standard types of validation are conducted: internal and external. We have executed both types of procedures thoroughly, and validated the model outputs using several diagnostic measures. They revealed significant gains in precision by our small area approach in obtaining the district-level estimates of FI. Since the area level model is developed at the district level by incorporating the household level variation as the sampling error variance, the model may not capture unknown variability due to any mid-level hierarchical components (such as a village or a sub-district area) in between a household and a district. Such misspecification is an important issue for understated standard errors of model-based estimator in the literature of unit level SAE method for poverty estimation [41, 42].

EIGP has long been studied for its relatively low productivity, poor infrastructure, limited capacity for private investment, and climate sensitivity [51]. Multiple effects of climate change, including recurrent droughts, floods, and cyclones, are causes of serious concern in EIGP [52]. Notably, UP, Bihar and WB appeared among the 10 most climatically vulnerable states in a recently published report [14]. Additionally, UP, Bihar and WB also rank (5, 9, and 14 respectively) among the poorest states in India [13]. Thus, both climate vulnerability and poverty might act as stressors that can exacerbate the prevailing FI, and disparities therein, among different populations in these states. On the other hand, crop diversity, especially of an agricultural district, might confer it with local resiliency and adaptability against both climate change as well as FI. Hence, we also computed crop diversity and poverty indices for each district in EIGP. Taken together with our estimates of FI measures, they provide a nuanced understanding of the overall conditions in the region.

The COVID-19 pandemic has impacted food security and livelihoods in many parts of the world [53, 54]. This is especially true for India which reported the second highest number of COVID-19 cases of any country in the world (>32 million) by mid-2021, with UP, Bihar and WB together contributing to almost 4 million cases [55]. Recommendations to address the ensuing FI range from expansion of governmental support for existing welfare schemes to investment in targeted social protection strategies and safety nets [56, 57]. Our SAE can provide disaggregated data for identifying the districts in EIGP with quantifiable impacts on their food security, as the diagnostic measures clearly demonstrate that our model-based estimates of FI indicators are more precise than direct estimates. For instance, relatively higher FIP in the Gangetic districts, as we observed in Bihar and WB, may benefit from systems that are resilient to FI while being adaptive to the stresses of climate change. In fact, while the number of Gangetic districts may be less than non-Gangetic ones, the policies to mitigate risks in the two groups should preferably be distinctive—with due consideration to their respective extreme event scenarios, and focused on adaptation strategies that are suited to each group and sustained implementation [19].

The Gangetic basin is considered to be under extreme pressure and one of the most vulnerable regions in terms of climate change [52]. Recurrent floods and droughts, extreme heat and seismic events, river bank erosion and sediment dynamics, urbanization and pollution, livelihoods and internal migration, land use and crop management, etc.–a variety of determinants would need to be considered for any integrated data-driven policy framework to ensure food security in EIGP. Towards this, data such as from the present study could provide spatially detailed insights to existing mechanisms, e.g., the Member of Parliament Local Area Development Scheme in India. For instance, post-pandemic assessment of the impact on population-specific food insecurity may lead to creation of durable community assets as well as support for micro, small and medium enterprises (or MSMEs [58]) in the vulnerable districts to jointly address local unemployment, poverty and food insecurity. The SAE approach allows us to take a closer look in that strategic direction as we plan to extend it to studies of the other states and regions in our future work. The inclusion of data for both urban and rural districts collected from multiple rounds of survey can lead to the development of a powerful and robust spatio-temporal area level model and, thus, better understanding of this complex phenomenon.

5. Conclusion

We estimated the district-specific food insecurity indicators, and mapped their local disparities over the EIGP region. By comparing food insecurity with indicators of climate vulnerability, poverty and crop diversity, we shortlisted the vulnerable districts in EIGP. These estimates and maps can be effective for informed policy-making to build local resiliency and address systemic vulnerabilities where they matter most in the post-pandemic era.

Supporting information

S1 Table. District wise estimates of FIP, FIG, FIS, and other variables.

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

(DOCX)

S2 Table. District level household summary in 2011–2012 HCES data.

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

(DOCX)

S1 Fig. District level residuals, histograms and q-q plots.

https://doi.org/10.1371/journal.pone.0279414.s003

(DOCX)

S2 Fig. Bias diagnostic plot of small area versus direct estimates.

https://doi.org/10.1371/journal.pone.0279414.s004

(DOCX)

S3 Fig. District-wise sampling fractions for the 3 EIGP states.

https://doi.org/10.1371/journal.pone.0279414.s005

(DOCX)

S4 Fig. District-wise confidence interval plots.

https://doi.org/10.1371/journal.pone.0279414.s006

(DOCX)

S5 Fig. Coefficient of variation boxplot of small area versus direct estimates.

https://doi.org/10.1371/journal.pone.0279414.s007

(DOCX)

S6 Fig. Geographical distribution of FIP in EIGP region.

https://doi.org/10.1371/journal.pone.0279414.s008

(DOCX)

Acknowledgments

The co-authors dedicate the article to the memory of Dr. Hukum Chandra who passed away in 2021 due to Covid-19.

References

  1. 1. World Food Summit—Final Report. In: Report of the World Food Summit [Internet]. Food and Agriculture Organization of the United Nations; 1996 [cited 2021 Jun 1]. http://www.fao.org/3/w3548e/w3548e00.htm
  2. 2. Hunger and food insecurity [Internet]. Food and Agriculture Organization of the United Nations. [cited 2021 Jun 1]. http://www.fao.org/hunger/en/
  3. 3. Poppy GM, Jepson PC, Pickett JA, Birkett MA. Achieving food and environmental security: New approaches to close the gap. Philos Trans R Soc B Biol Sci. 2014;369(1639). pmid:24535384
  4. 4. SOFI 2017—The State of Food Security and Nutrition in the World [Internet]. [cited 2021 Jan 13]. http://www.fao.org/state-of-food-security-nutrition/en/
  5. 5. Mendu V, Perugu S, Manne M, Bolla R, Aruru M, Pyne S. Nutrition Atlas of ICMR-National Institute of Nutrition: an informatics platform on nutrition in India. Stat Appl. 2018;17:209–19.
  6. 6. von Grebmer, K., J. Bernstein, D. Resnick, M. Wiemers, L. Reiner, M. Bachmeier, et al. 2022. 2022 Global Hunger Index: Food Systems Transformation and Local Governance. Bonn: Welthungerhilfe; and Dublin: Concern Worldwide. https://www.globalhungerindex.org/
  7. 7. Ministry of Health and Family Welfare (MoHFW), Government of India U, Council. and P. Comprehensive National Nutrition Survey 2016–2018 [Internet]. 2019. https://www.popcouncil.org/uploads/pdfs/2019RH_CNNSreport.pdf
  8. 8. Kumaran K, Krishnaveni G V., Suryanarayana KG, Prasad MP, Belavendra A, Atkinson S, et al. Protocol for a cluster randomised trial evaluating a multifaceted intervention starting preconceptionally-Early Interventions to Support Trajectories for Healthy Life in India (EINSTEIN): A Healthy Life Trajectories Initiative (HeLTI) Study. BMJ Open. 2021;11(2). pmid:33593789
  9. 9. Auffhammer M, Ramanathan V, Vincent JR. Climate change, the monsoon, and rice yield in India. Clim Change [Internet]. 2012 Mar 2 [cited 2021 Jun 1];111(2):411–24. Available from: https://link.springer.com/article/10.1007/s10584-011-0208-4
  10. 10. Knox J, Hess T, Daccache A, Wheeler T. Climate change impacts on crop productivity in Africa and South Asia. Environ Res Lett [Internet]. 2012 Sep 14 [cited 2021 Jun 1];7(3):34032. Available from: https://iopscience.iop.org/article/10.1088/1748-9326/7/3/034032
  11. 11. Zhao C, Liu B, Piao S, Wang X, Lobell DB, Huang Y, et al. Temperature increase reduces global yields of major crops in four independent estimates. Proc Natl Acad Sci U S A [Internet]. 2017 Aug 29 [cited 2021 Jun 1];114(35):9326–31. Available from: https://www.pnas.org/content/114/35/9326 pmid:28811375
  12. 12. Area Production Statistics. National Horticulture Board. Ministry of Agriculture and Farmer’s Welfare. Government of India. [Internet]. [cited 2021 Jun 2]. http://nhb.gov.in/Statistics.aspx?enc=WkegdyuHokljEtehnJoq0KWLU79sOQCy+W4MfOk01GFOWQSEvtp9tNHHoiv3p49g
  13. 13. Reserve Bank of India—Publications [Internet]. [cited 2021 Jun 1]. https://www.rbi.org.in/scripts/PublicationsView.aspx?id=18621
  14. 14. Department of Science & Technology G of I. Climate Vulnerability Assessment for Adaptation Planning in India Using a Common Framework. [Internet]. 2021. https://dst.gov.in/sites/default/files/Full Report %281%29.pdf
  15. 15. Joint Global Change Research Institute and Battelle Memorial Institute PND. India: Impact of Climate Change to 2030. A Commissioned Research Report [Internet]. Washington, D.C.; 2009. https://www.dni.gov/files/documents/climate2030_india.pdf
  16. 16. Burney J, Ramanathan V. Recent climate and air pollution impacts on indian agriculture. Proc Natl Acad Sci U S A [Internet]. 2014 Nov 18 [cited 2021 Jun 1];111(46):16319–24. Available from: www.pnas.org/cgi/doi/10.1073/pnas.1317275111 pmid:25368149
  17. 17. Chakraborti D, Singh SK, Rahman MM, Dutta RN, Mukherjee SC, Pati S, et al. Groundwater arsenic contamination in the ganga river basin: A future health danger [Internet]. Vol. 15, International Journal of Environmental Research and Public Health. MDPI AG; 2018 [cited 2021 May 31]. p. 180. www.mdpi.com/journal/ijerph
  18. 18. Davis KF, Chhatre A, Rao ND, Singh D, Ghosh-Jerath S, Mridul A, et al. Assessing the sustainability of post-Green Revolution cereals in India. Proc Natl Acad Sci U S A [Internet]. 2019 [cited 2021 Jun 1];116(50):25034–41. Available from: https://pubmed.ncbi.nlm.nih.gov/31754037/ pmid:31754037
  19. 19. Aryal JP, Sapkota TB, Khurana R, Khatri-Chhetri A, Rahut DB, Jat ML. Climate change and agriculture in South Asia: adaptation options in smallholder production systems. Environ Dev Sustain 2019 226 [Internet]. 2019 Jul 9 [cited 2021 Aug 21];22(6):5045–75. Available from: https://link.springer.com/article/10.1007/s10668-019-00414-4
  20. 20. Pingali PL. Green revolution: Impacts, limits, and the path ahead. Proc Natl Acad Sci U S A. 2012;109(31):12302–8. pmid:22826253
  21. 21. Chandra H, Aditya K, Sud UC. Localised estimates and spatial mapping of poverty incidence in the state of Bihar in India—An application of small area estimation techniques. PLoS One [Internet]. 2018 Jun 1 [cited 2021 Jun 2];13(6):e0198502. Available from: pmid:29879202
  22. 22. Aleksandrowicz L, Tak M, Green R, Kinra S, Haines A. Comparison of food consumption in Indian adults between national and sub-national dietary data sources. Br J Nutr [Internet]. 2017 Apr 2 [cited 2017 Nov 21];117(07):1013–9. Available from: https://www.cambridge.org/core/product/identifier/S0007114517000563/type/journal_article pmid:28462737
  23. 23. Rao JNK, Molina I. Small Area Estimation: Second Edition [Internet]. Small Area Estimation: Second Edition. Wiley; 2015.
  24. 24. Jones AD, Ngure FM, Pelto G, Young SL. What are we assessing when we measure food security? A compendium and review of current metrics. Adv Nutr. 2013;4(5):481–505. pmid:24038241
  25. 25. Maxwell D, Vaitla B, Coates J. How do indicators of household food insecurity measure up? An empirical comparison from Ethiopia. Food Policy [Internet]. 2014;7:107–16. Available from:
  26. 26. Hossain MJ, Das S, Chandra H, Islam MA. Disaggregate level estimates and spatial mapping of food insecurity in Bangladesh by linking survey and census data [Internet]. Vol. 15, PLoS ONE. Public Library of Science; 2020 [cited 2021 Jun 2]. p. e0230906. Available from: pmid:32275683
  27. 27. Madhavan Nair KP, Augustine LF. Country-specific nutrient requirements & recommended dietary allowances for Indians: Current status & future directions [Internet]. Vol. 148, Indian Journal of Medical Research. Wolters Kluwer Medknow Publications; 2018 [cited 2021 Jun 3]. p. 522–30.
  28. 28. Chancel L. Indian Inequality Updates (2015–2019). Technical Note. [Internet]. 2020. https://wid.world/document/indian-inequality-updates-2015-2019-world-inequality-lab-technical-note-2020-09/
  29. 29. Guha S, Chandra H. Measuring and Mapping Disaggregate Level Disparities in Food Consumption and Nutritional Status via Multivariate Small Area Modelling. Soc Indic Res [Internet]. 2021 Jan 3 [cited 2021 Dec 16];154(2):623–46. Available from: https://link.springer.com/article/10.1007/s11205-020-02573-8
  30. 30. Guha S, Chandra H. Measuring disaggregate level food insecurity via multivariate small area modelling: evidence from rural districts of Uttar Pradesh, India. Food Secur. 2021;13(3):597–615.
  31. 31. Gopalan C, Sastri, Rama B., Balasubramaniam S. Nutritive values of Indian foods. 2004.
  32. 32. CensusInfo India 2011–Dashboards, Data Query, Houselisting and Housing data, Population and Education [Internet]. [cited 2021 Jun 2]. http://www.dataforall.org/dashboard/censusinfo/
  33. 33. Indian Council of Medical Research. Nutrient Requirements and Recommended Dietary Allowances for Indians: A Report of the Expert Group of the Indian Council of Medical Research. Hyderabad; 2010.
  34. 34. Web Based Land Use Statistics Information System [Internet]. [cited 2021 Aug 22]. https://aps.dac.gov.in/LUS/Public/Reports.aspx
  35. 35. Foster J, Greer J, Thorbecke E. A class of decomposable poverty measures. Econometrica [Internet]. 1984;52(3):761–6. Available from:
  36. 36. Elbers C, Lanjouw JO, Lanjouw P. Micro-level estimation of poverty and inequality. Econometrica. 2003; 71: 355–364. [cited 2022 Oct 5]
  37. 37. Molina I, Rao JNK. Small area estimation of poverty indicators. The Canadian Journal of Statistics. 2010; 38: 369–385. [cited 2022 Oct 5]
  38. 38. Tzavidis N, Salvati N, Pratesi M, Chambers R. M-quantile models with application to poverty mapping. Statistical Methods and Applications. 2008; 17(3): 393–411. [cited 2022 Oct 5]
  39. 39. Fay RE, Herriot RA. Estimates of Income for Small Places: An Application of James-Stein Procedures to Census Data. J Am Stat Assoc [Internet]. 1979 Jun [cited 2021 Jun 2];74(366a):269–77. Available from: https://www.tandfonline.com/doi/abs/10.1080/01621459.1979.10482505
  40. 40. Chandra H, Chambers R, Salvati N. Small area estimation of survey weighted counts under aggregated level spatial model. Surv Methodol Stat Canada, Cat No 12-001-X [Internet]. 2019 [cited 2021 Jun 2];45(1):31–59. Available from: https://www150.statcan.gc.ca/n1/pub/12-001-x/2019001/article/00006-eng.htm
  41. 41. Das S, Chambers R. Robust mean-squared error estimation for poverty estimates based on the method of Elbers, Lanjouw and Lanjouw. Journal of the Royal Statistical Society: Series A (Statistics in Society). 2017, 180 (4), 1137–1161. [cited 2022 Oct 5]
  42. 42. Casas-Cordero C., Encina J., & Lahiri P. 2016. Poverty Mapping for the Chilean Comunas. Pages 379–403 of: Analysis of Poverty by Small Area Estimation. John Wiley & Sons. [cited 2022 Oct 5]
  43. 43. Johnson FA, Chandra H, Brown JJ, Padmadas SS. District-level estimates of institutional births in ghana: Application of small area estimation technique using census and DHS data. J Off Stat. 2010;26(2):341–59.
  44. 44. Chandra H, Salvati N, Sud UC. Disaggregate-level estimates of indebtedness in the state of Uttar Pradesh in India: an application of small-area estimation technique. [Internet]. 2011 Nov [cited 2021 Aug 22];38(11):2413–32. Available from: https://www.tandfonline.com/doi/abs/10.1080/02664763.2011.559202
  45. 45. Brown G, Chambers R, Heady P, Heasman D. Evaluation of small area estimation methods—an application to unemployment estimates from the UK LFS. Proceedings of Statistics Canada Symposium 2001. Achieving Data Quality in a Statistical Agency: a Methodological Perspective. 2001. [cited 2022 Oct 5]
  46. 46. Bickel G, Nord M, Price C, Hamilton W, Cook J. Guide to Measuring Household Food Security, Revised 2000 [Internet]. 2000 [cited 2021 Jun 2]. https://www.scirp.org/(S(i43dyn45teexjx455qlt3d2q))/reference/referencespapers.aspx?referenceid=1418378
  47. 47. The state of food security and nutrition in the world 2019. Safeguarding against economic slowdowns and downturns. [Internet]. Rome; 2019. http://www.fao.org/3/ca5162en/ca5162en.pdf
  48. 48. COVID-19 will double number of people facing food crises unless swift action is taken [Internet]. Food and Agriculture Organization of the United Nations. 2020 [cited 2017 Dec 3]. http://www.fao.org/about/what-we-do/en/
  49. 49. Damerau K, Davis KF, Godde C, Herrero M, Springmann M, Bhupathiraju SN, et al. India has natural resource capacity to achieve nutrition security, reduce health risks and improve environmental sustainability. Nat Food [Internet]. 2020;1(10):631–9. Available from: http://dx.doi.org/10.1038/s43016-020-00157-w
  50. 50. Choudhary K, Das S, Ghodajkar P. Trends of adult height in India from 1998 to 2015: Evidence from the National Family and Health Survey. PLoS One. 2021;16(9):e0255676. pmid:34534216
  51. 51. Aggarwal PK, Joshi PK, Ingram JSI, Gupta RK. Adapting food systems of the Indo-Gangetic plains to global environmental change: Key information needs to improve policy formulation. Environ Sci Policy. 2004;7(6):487–98.
  52. 52. Sharma B, Amarasinghe U, Xueliang C, de Condappa D, Shah T, Mukherji A, et al. The Indus and the Ganges: river basins under extreme pressure. In: Fisher M, Cook S, editors. Water, food and poverty in river basins: defining the limits. Routledge; 2012. p. 30–58.
  53. 53. Das S, Rasul MG, Hossain MS, Khan AR, Alam MA, Ahmed T, et al. Acute food insecurity and short-term coping strategies of urban and rural households of Bangladesh during the lockdown period of COVID-19 pandemic of 2020: report of a cross-sectional survey. BMJ Open [Internet]. 2020 Dec 1 [cited 2021 Dec 16];10(12):e043365. Available from: https://bmjopen.bmj.com/content/10/12/e043365 pmid:33310813
  54. 54. Kang Y, Baidya A, Aaron A, Wang J, Chan C, Wetzler E. Differences in the early impact of COVID-19 on food security and livelihoods in rural and urban areas in the Asia Pacific Region. Glob Food Sec. 2021 Dec 1;31:100580. pmid:36570721
  55. 55. COVID-19 Dashboard, Government of India [Internet]. [cited 2021 Aug 1]. https://www.mygov.in/covid-19
  56. 56. Nguyen PH, Kachwaha S, Pant A, Tran LM, Ghosh S, Sharma PK, et al. Impact of COVID-19 on household food insecurity and interlinkages with child feeding practices and coping strategies in Uttar Pradesh, India: a longitudinal community-based study. BMJ Open [Internet]. 2021 Apr 1 [cited 2021 Dec 16];11(4):e048738. Available from: https://bmjopen.bmj.com/content/11/4/e048738 pmid:33883156
  57. 57. Sinha D. Hunger and food security in the times of Covid-19. J Soc Econ Dev 2021 232 [Internet]. 2021 Feb 15 [cited 2021 Dec 16];23(2):320–31. Available from: https://link.springer.com/article/10.1007/s40847-020-00124-y pmid:34720482
  58. 58. Mishra K, Rampal J. The COVID-19 pandemic and food insecurity: A viewpoint on India. World Dev [Internet]. 2020 Nov 1 [cited 2021 Jun 1];135:105068. Available from: https://www.sciencedirect.com/science/article/abs/pii/S0305750X20301947?via%3Dihub