Micronutrient deficiencies (MNDs) remain widespread among people in sub-Saharan Africa1,2,3,4,5, where access to sufficient food from plant and animal sources that is rich in micronutrients (vitamins and minerals) is limited due to socioeconomic and geographical reasons4,5,6. Here we report the micronutrient composition (calcium, iron, selenium and zinc) of staple cereal grains for most of the cereal production areas in Ethiopia and Malawi. We show that there is geospatial variation in the composition of micronutrients that is nutritionally important at subnational scales. Soil and environmental covariates of grain micronutrient concentrations included soil pH, soil organic matter, temperature, rainfall and topography, which were specific to micronutrient and crop type. For rural households consuming locally sourced food—including many smallholder farming communities—the location of residence can be the largest influencing factor in determining the dietary intake of micronutrients from cereals. Positive relationships between the concentration of selenium in grain and biomarkers of selenium dietary status occur in both countries. Surveillance of MNDs on the basis of biomarkers of status and dietary intakes from national- and regional-scale food-composition data1,2,3,4,5,6,7 could be improved using subnational data on the composition of grain micronutrients. Beyond dietary diversification, interventions to alleviate MNDs, such as food fortification8,9 and biofortification to increase the micronutrient concentrations in crops10,11, should account for geographical effects that can be larger in magnitude than intervention outcomes.
Globally, more than two billion people are affected by one or more MNDs and the risks of deficiency are greater in sub-Saharan Africa (SSA) than in most other regions1,2,3,12. These MNDs, which are also known as ‘hidden hunger’, remain a major challenge for achieving the United Nations’ Sustainable Development Goal 2 (zero hunger) by 203012. Causes of MNDs include the inadequate intake of micronutrients —in particular, calcium (Ca), iron (Fe), iodine (I), selenium (Se), zinc (Zn) and vitamin A —especially in regions in which diets are dominated by cereals and where access to foods from plant and animal sources that are richer in nutrients is limited11. Most cereal grains have inherently small micronutrient concentrations, especially once bran and embryo fractions are removed during milling10. Cereal grains also contain large concentrations of anti-nutritional compounds such as phytates (inositol phosphate compounds), which inhibit the absorption of Ca, copper (Cu), Fe, magnesium (Mg) and Zn in the human gut2,3,5,10.
The prevalence of MNDs must be estimated if policy responses are to be developed in the context of socioeconomic and environmental drivers of food system change6. The prevalence of MNDs can be determined from multiple sources of evidence. Biomarkers of status, including micronutrient concentrations or enzyme activities in blood and other tissues, are often used to assess population status7,12,13. However, establishing thresholds of sufficiency for biomarkers can be challenging due to variation in the ranges considered to be ‘healthy’ between demographic groups, physiological buffering and the influence of infection and inflammation, which can have short-term effects on circulating concentrations of micronutrients in the human body7,14. Biomarker studies also impose burdens on participants and technical challenges for the collection, storage and analysis of samples, especially in low-income settings. Complementary methods to estimate MND risks include measuring micronutrient intake from composite dietary analyses15 or, more commonly, from estimates of intake from dietary recall5, household food consumption and expenditure data5,16 and food balance sheets1,2,3.
Estimates of the prevalence of MNDs from food supply and intake require reliable data on the micronutrient composition of food; national or international food-composition tables are routinely used in SSA2,5. However, data are typically limited to single values for each food–nutrient combination and are therefore of limited value for exploring spatial drivers of micronutrient intake. This is problematic given the widespread consumption of locally produced staple foods such as cereal grains among smallholder farming communities, and given that food micronutrient composition varies within countries5,17.
Subnational food composition data have not yet been used to inform estimates of the prevalence of MNDs in SSA. This is despite a long history of studies in the fields of medical geology and veterinary sciences, which have linked human and livestock nutritional status to soil and landscape factors. These include studies on the micronutrients I18, Se4,18,19, molybdenum (Mo)20 and cobalt (Co)21. Substantial variation in the Ca, Fe, Se and Zn concentrations of staple crops was recently reported from nine districts in Uganda and some of this variation was associated with soil characteristics17. The case for using subnational food composition data to estimate the prevalence of MNDs is strengthened by recent evidence of long-range geospatial variation in grain Se concentration being linked to biomarkers of Se dietary status. In parts of the Amhara region of Ethiopia, nutritionally important variation in the grain Se concentration of teff (Eragrostis tef (Zucc.) Trotter) and wheat (Triticum aestivum L.) was associated with soil and landscape covariates at distances extending further than 100 km, which is similar to patterns seen among human biomarkers of Se dietary status22,23. In Malawi, human Se dietary status was associated with variation in the grain Se concentration in maize (Zea mays L.) at localized scales4,5,15,24,25,26. Grain sampling for micronutrient quality has not yet been conducted systematically at wider geographical scales in SSA.
Grain surveys of cultivated land in Ethiopia and Malawi
Grain micronutrient concentrations in cereal crops are reported here from 1,389 locations in Ethiopia, based on a spatially balanced sample in the Amhara, Oromia and Tigray regions, during the late-2017 and late-2018 harvest seasons. The sampling frame represents a subset of around 354,000 km2 of cultivated land—representative of most of the cereal production area in Ethiopia (Fig. 1)—which was constrained to accessible locations that were less than 2.5 km from a known road. At each location, a grain sample and a co-located composite soil sample were taken with the informed consent of the farmer. Grain samples reported from Ethiopia included teff (n = 373), wheat (n = 328), maize (n = 302), sorghum (Sorghum bicolor (L.) Moench; n = 138), barley (Hordeum vulgare L.; n = 181) and finger millet (Eleusine coracana (L.) Gaertn.; n = 39), with a smaller number of triticale (× Triticosecale Wittm. ex A. Camus; n = 20) and rice (Oryza sativa L.; n = 8) samples. In Malawi, grain micronutrient concentrations are reported from 1,812 locations, which were sampled during the April–June 2018 harvest season. The Malawi sampling frame represents around 66,000 km2 of cultivated land (Fig. 1). Cereal production in Malawi is less diverse than in Ethiopia and grain samples comprised mostly maize (n = 1,608) samples, together with sorghum (n = 117), rice (n = 54), pearl millet (Pennisetum glaucum (L.) R. Br.; n = 32) samples and a single finger millet sample. Sampling designs, grain analyses and geostatistical methods expanded on those described for part of the Amhara region22 (Extended Data Figs. 1–9 and Extended Data Tables 1–3).
Variation between and within crops
Grain Ca, Fe, Se and Zn concentrations varied substantially between and within crop species (Fig. 2 and Extended Data Table 1). Maize typically had the lowest concentration of all four micronutrients; people relying on maize-based diets are therefore likely to have the lowest micronutrient intakes. Finger millet is a good potential source of Ca, with a median grain Ca concentration (4,574 mg kg–1; range, 3,203–6,264 mg kg–1) that is almost two orders of magnitude greater than in maize (median, 64.5 mg kg–1) in Ethiopia. The single finger millet sample from Malawi had a grain Ca concentration of 3,564 mg kg–1, which is consistent with the large Ca concentrations that have been reported for finger millet flour from markets in Malawi26 (range, 2,900–4,700 mg kg–1). A screen of finger millet varieties in India27 showed a grain Ca concentration range of 1,350–3,120 mg kg–1 (n = 49 genotypes), indicating that the trait is conserved across the species. Although the bioavailability of Ca in finger millet will be limited by phytate, as it is for teff, which also had a high grain Ca concentration in Ethiopia (median, 1,473 mg kg–1; range, 46.8–7,925 mg kg–1), fermenting flour to produce injera (a thin flat bread that is commonly consumed in Ethiopia) increases the bioavailability of Ca and other mineral micronutrients by stimulating endogenous phytase enzymes28.
Within-species variation in grain Ca, Fe, Se and Zn concentrations will arise due to spatial variation in soil and landscape factors, and the effects of extrinsic soil dust. The potential effects of soil dust are pronounced for Fe; soil total Fe concentrations (medians, 92,744 and 28,804 mg kg–1 for Ethiopia and Malawi, respectively) are more than three orders of magnitude larger than median grain Fe concentrations of 20.3 and 21.3 mg kg–1 for maize grain in Ethiopia and Malawi, respectively. Grain concentrations of Se and Zn are less sensitive to soil dust due to the much lower total concentrations of these elements in soils (medians, 0.35 and 0.32 mg kg–1 for Se and 100.1 and 33.7 mg kg–1 for Zn, in Ethiopia and Malawi, respectively).
Geospatial mapping and dietary contributions
This study focused on mineral concentrations in grain of teff and wheat for Ethiopia and maize for Malawi. These crops were chosen because they comprise a large proportion of the energy intake in national diets and have good spatial coverage in the survey. Grain concentration maps were based on ordinary kriging. Kriging variances—the expected squared error of the predictions—quantify the uncertainties in the maps (Extended Data Figs. 3, 4). The dietary contribution for each crop–nutrient combination was then mapped as a percentage of dietary requirements and visualized on quartile scales from yellow (small) to dark red (large). These calculations used food balance sheets from the Food and Agriculture Organization29 and assumed estimated average requirement (EAR) thresholds for a representative demographic group—adult women aged 18–24 years eating an unrefined (that is, high phytate) diet30.
There is spatially dependent variation in grain Ca, Fe, Se and Zn concentrations over large distances in Ethiopia and Malawi (Figs. 3, 4). These observations are likely to be of nutritional importance given that most cereals are grown, milled and consumed locally in Ethiopia31 and Malawi5, because they show that the dietary supply and intake of these nutrients varies substantially from one location to another. In Ethiopia, spatial dependencies were seen over distances from 100 to 200 km for the concentration of Ca in teff, for the concentrations of Fe in teff and wheat, and for the concentration of Se in wheat (Extended Data Fig. 2 and Extended Data Table 3). For grain Se and Zn concentrations in teff, and grain Ca and Zn concentrations in wheat, longer-range spatial variation extends beyond 250 km. In Malawi, the spatial dependence of the variation in grain Ca, Fe and Se concentration of maize occurs at distances of 50–80 km. For the concentration of Zn in maize grain in Malawi, more of the variation was attributable to differences over distances that were too short to be resolved by our sampling frame, although the value of the variogram still increased at distances up to 100 km.
In Ethiopia, the Ca concentrations in teff and wheat grain were generally greater in crops sampled from north, northeast and east Tigray region, north and northwest Amhara region and from areas of the Rift Valley than crops from most of the Amhara region (Fig. 3). A trend of increasing grain Se concentrations in teff and wheat22 and increasing Se dietary status of children31 from the west to the east of the Amhara region has been reported previously. The largest grain Fe and Zn concentrations in teff and wheat were in north Amhara, Tigray and northwest and southwest Oromia regions. In Malawi, maize grain Ca, Se, Zn and—to a lesser extent—Fe concentrations are typically greater in the Shire Valley of the southern region, than in the northern and central regions (Fig. 4). These observations are consistent with observations that crops growing on localized soil types (for example, Vertisols) have greater grain micronutrient concentrations than crops growing on more weathered, acidic soils4,5,25,26. As in Ethiopia, geographical areas with the largest maize grain Se concentration in Malawi co-located with areas of Se sufficiency that have been reported in both cross-sectional studies15 and representative surveillance analyses24 of human biomarkers of Se dietary status.
Links between grain micronutrient concentration and a dietary outcome were observed for Se, for which there is a reliable biomarker of dietary status7,13. Biomarker data for Se for 321 enumeration areas in Ethiopia23 and 101 enumeration areas in Malawi24 are available from recent National Micronutrient Surveys. For each enumeration area, the nearest grain sample site for any crop type was identified. In both countries, there is a statistically significant positive relationship between the Se concentration in grain and the Se concentration in serum (Ethiopia) and plasma (Malawi) (Extended Data Fig. 5). Direct evidence of links between the grain concentration of other micronutrients, biomarkers of dietary status and health outcomes remains a major research challenge.
Soil and environmental factors that influence variation in grain micronutrients—for teff, wheat and maize in Ethiopia, and for maize in Malawi—were observed (Extended Data Figs. 6–9 and Extended Data Table 3). For example, soil pH was statistically significantly correlated with grain Se concentration for all crops in both countries (positive relationships). Grain Zn concentration was significantly correlated with soil pH for teff in Ethiopia (negative relationship) and for maize in both countries (positive relationships). Further studies of contrasting responses to soil pH, observed in teff and maize, are needed. However, the generally weak predictive value of soil pH on grain Zn concentration (Extended Data Figs. 7, 9) is consistent with a survey of staple crops in Uganda17. Soil organic carbon was significantly correlated with grain Se concentration in wheat in Ethiopia (negative relationship), grain Zn concentration in wheat in Ethiopia (positive relationship) and grain Zn concentration in maize in both countries (positive relationships). Grain micronutrient composition was also significantly correlated with at least one environmental covariate for each micronutrient, crop and country. In Ethiopia, grain Se concentration was negatively correlated with mean annual precipitation for each crop and was positively correlated with mean annual temperature for teff and wheat in Ethiopia and maize in Malawi. Mean annual temperature and topographic index—a measure of soil wetness at a site—were positively correlated with grain Zn concentration for maize in Malawi. Multivariate spatial statistical modelling could be informative about soil and environmental factors that jointly influence the variation of grain micronutrient concentration; however, this requires additional statistical assumptions and is beyond the scope of this study.
In Ethiopia, the proportion of dietary Ca requirements that can be met by the consumption of teff and wheat is much greater than for maize in Malawi, however, it is still likely to be less than 25% of the required amount for most of the population (Fig. 3). In Malawi, maize intake provides less than 3% of dietary Ca requirements for many, despite providing more than 50% of dietary energy requirements (Fig. 4). Cereals provide a greater proportion of the dietary Fe, Se and Zn requirements than Ca in both countries. For example, for some rural households in Malawi, more than 100% of dietary Fe, Se and Zn requirements will be met by eating locally produced maize alone. However, most households will receive less than 25% of Se, less than 50% of Fe and less than 75% of their Zn requirements from a typical consumption pattern for maize.
Surveillance and interventions
Surveillance of MNDs and food system foresighting activities currently rely on national or international food composition data to estimate micronutrient supply as a proxy for intake6,7. Subnational food composition data are not yet used, probably because of the logistical and conceptual challenges in generating data in the required forms. Challenges include sampling crops across large areas during short harvest periods, analysing large numbers of samples and associated data in areas in which access to suitable laboratories and trained personnel is often lacking, and communicating sparse data and associated uncertainties. Spatial dependencies in the grain micronutrient concentrations found in this study can be communicated simply using localized administrative level boundaries. Mean grain Ca and Zn concentrations for teff and wheat show spatial trends when mapped to woreda administrative levels in the Amhara region of Ethiopia using geostatistical models and block kriging (Fig. 5).
Ideally, improved access to diverse diets would alleviate many MNDs1,2,3,4,5,6. However, this is unlikely to be feasible in the shorter term for many people for socioeconomic reasons6. A crop micronutrient survey within a country could inform shorter-term interventions to alleviate MNDs. Such interventions include supplementation32, food fortification8,9 and biofortification of staple crops through breeding and agronomic approaches10,11. For example, Zn-biofortified wheat varieties, which were recently released in India and Pakistan33,34, have been bred with a target to increase Zn concentrations by 8–12 mg kg–1 above a notional baseline grain Zn concentration of 25 mg kg–1. Similarly, increases in grain Zn concentration in new Zn-biofortified hybrid maize varieties of 15% in Guatemala (ICTA HB-18 and ICTA B-15) and 36% in Colombia (BIO-MZN01), with target levels of around 30 mg kg–1, have been reported35. Here, geographical differences in grain Zn concentration, in both wheat (Ethiopia) and maize (Malawi), can exceed these breeding targets. Subnational spatial sources of variation could support priority areas for the release of new crop varieties or micronutrient fertilization strategies and improve impact evaluations of biofortification interventions. The persistence of MNDs is clearly more complex than micronutrient supply in cereal grains36. For example, livestock are important components of diverse diets and income; identifying areas of lower concentrations of micronutrients in forage crops could inform strategies to improve livestock health and production20,21.
This study did not explore the effects of crop variety (genotype) or farmer management strategies (management), which will contribute substantially to the variation in yield and grain micronutrient concentrations, even over short distances within and between fields of the same farm37,38. For example, the preferential use of locally sourced organic materials by smallholder farmers on certain fields can improve the quality and yield of grain micronutrients in maize-based38 and wheat-based39 systems in SSA. Temporal variation in environmental (environment) factors, such as the projected decreases in cereal grain micronutrient concentrations due to increased atmospheric CO2—of 6% for Fe and 9% for Zn in wheat by the mid-twenty-first century40—and increased leaching of soil Se under higher rainfall41, should also be considered. However, rising temperatures may compensate for some of these effects in terms of grain micronutrient quality42. A better understanding of how the complex interactions between genotype, management and environment drive crop micronutrient quality within diverse, climate-smart farming systems is essential for a more-sustainable global food system.
Grain and soil samplings from farmers’ fields were completed in November–December 2017 (for most of Amhara region) and November 2018–January 2019 (Amhara, Oromia and Tigray regions) in Ethiopia, and in April–June 2018 in Malawi. The work involved sampling of grain and soil from farmers’ fields and grain stores with the informed consent of the farmers. The work was conducted under ethical approvals from the University of Nottingham, School of Sociology and Social Policy Research Ethics Committee (REC); BIO-1718-0004 and BIO-1819-001 for Ethiopia and Malawi, respectively. These REC approvals were recognized formally by the Directors of Research at Addis Ababa University (Ethiopia) and Lilongwe University of Agriculture and Natural Resources (Malawi), who also reviewed the study protocols.
The objective of this study was to support the spatial mapping of grain Ca, Fe, Se and Zn concentration of cereal crops. We sought reasonable spatial coverage over the target sample frame and used ‘main-site’ and ‘close-pair’ sampling to support the estimation of parameters of the spatial linear mixed model (LMM)22,43.
In Ethiopia (Amhara, Oromia and Tigray regions), target sample frames were constrained to locations at which the probability of the land being under crop production was ≥0.9 based on predictions produced on a 250-m grid. These predictions were derived from the interpretation of high-resolution satellite imagery by trained observers and by machine learning methods applied to multiple covariates derived from remote sensing data and digital elevation models44,45. The sample frame was further constrained to include only those locations from a 250-m grid that fell within 2.5 km of a known road. A map indicating nodes on a 250-m grid (with the same origin as the agricultural land-use grid) that met this requirement was prepared. These constraints may introduce possible biases into predictions made at locations outside the designed sample frame, however, it would not otherwise have been possible to visit all of the sample locations in the time available. Information on the distribution of roads was taken from OpenStreetMap (www.openstreetmap.org). Of a total land area of around 1.1 million km2 in Ethiopia, the total cropland mask represented 354,325 km2, of which 220,467 km2 was within 2.5 km of a known road. In Malawi, the cropland area was determined from the European Space Agency Climate Change Initiative46. The agricultural area used was defined as including all raster cells that included the category of ‘cropland’ in their description. In Malawi, where road access to cropped areas is generally better than in Ethiopia, no constraint to road distance was imposed on sample locations. The mapped cropland areas are shown in Fig. 1.
In Ethiopia, a total of 1,825 primary sample locations were selected a priori, with each 250-m grid node within the sampling frame allocated an equal prior inclusion probability. This was done using the lcube function from the BalancedSampling library for the R platform47,48. The lcube function implements the cube method, to enable random sampling according to specified inclusion probabilities while aiming for balance and spread with respect to the specified covariates49. Here, inclusion probabilities were uniform across the sample frame and sample locations were selected for spatial balance, which entails that the mean coordinates of sample locations are close to the mean coordinates of all points in the sample frame, and spatial spread, which ensures that the observations are spread out rather than clustered with respect to the spatial coordinates50. A subset of 175 of these locations were selected as close-pair sites at which an additional nearby sample was taken to support the estimation of parameters of the spatial LMM43.
In Malawi, a different sampling method was used to achieve good spatial coverage of a total of 1,710 main-site locations. These included 820 fixed sample points from the 2015/16 Demographic and Health Survey of Malawi24,51. The stratify function in the spcosa library for the R platform52 divides a sampling domain into Delaunay polygons centred on a set of fixed points and with the remaining polygon centroids selected to partition the domain into approximately equal-area regions. The centroids of the polygons were selected as sample points. An additional 890 sample points were found in addition to the 820 fixed ones, with the stratify function. Once these were obtained, a further 190 locations were selected at random as close-pair sites for an additional nearby sample, as in Ethiopia.
Sampling was conducted by teams who were trained in standard procedures and risk assessments. Each team planned to visit five main-site locations per day. Main-site locations were loaded onto a computer tablet and printed on paper maps for each team. A team would navigate to the target location, using a handheld GPS device for the last few kilometres. At each sample location, the team would identify the nearest field with a mature cereal crop within a 1-km radius, and sample grain and soil, subject to farmer consent. If a field with a standing mature cereal crop was not apparent, that is, the crop had been harvested, or a non-cereal crop had been grown, the team would ask the farmer to identify a field from which a cereal crop had recently been harvested and stored, and from which a sample could be obtained. If sampling was not possible, then the team would either look beyond a 1-km radius for an alternative location, or abandon the location. At designated close-pair locations, a second field was identified ideally within around 500 m (range, 100–1,000 m) of the main-site location. If a close-pair location could not be found, then a close-pair location would be selected at the next sample location that was not already earmarked for a close-pair sample.
Within a selected field, samples were taken from a 100 m2 (0.01 ha) circular plot. This was centred as close as practical to the middle of the field unless this area was unrepresentative due to disease or crop damage. Five subsample points were located (Extended Data Fig. 1). The first point was at the centre of the plot. Two subsample points were then selected at locations on a line through the plot centre along the crop rows, and two more points on a line orthogonal to the first through the plot centre. Where possible, the central sampling location was fixed between crop rows, and the long axis of the sample array (with sample locations at 5.64 and 4.89 m) was oriented in the direction of crop rows with the short axis perpendicular to the crop rows. A single soil subsample was collected at each of the five subsample points with a Dutch auger with a flight of length of 150 mm and diameter of 50 mm. The auger was inserted vertically to the depth of one flight and the five subsamples were stored in a single bag. Where a mature or ripe crop was still standing in the field, grain samples were taken close to each augering position by a different operator, to minimize further contamination by dust and soil. For maize, a single cob was taken at each of the five points. Maize kernels were stripped from around 50% of each cob lengthways and composited into a single sample envelope for each location. For smaller-grained crops, sufficient stalks were taken so that approximately 20–50% of the sample envelope was filled (dimensions 15 cm × 22 cm), with samples placed grain‐first into the sample bag and the stalks were twisted off the grain heads and discarded. If a crop was in field stacks, then a subsample, comprising five cobs for maize, or a representative sample for other crops was taken from each available stack, taking material from inside the stack to minimize contamination by dust and soil (Extended Data Fig. 1). If a crop was in a farmers’ store, that is, already averaged from across the field, then a representative sample was taken, while avoiding grain from the store floor if grain was loosely stored and avoiding grain with visible soil or dust contamination.
Photographs at sample locations and of sample bags were recorded for quality assurance along with site GPS locations. In Ethiopia, 1,385 of the 1,389 locations from where grain data are reported had positional uncertainties of ≤8 m as recorded by the GPS. The other four locations had positional uncertainties of 9–16 m. In Malawi, 1,790 of the 1,812 locations had positional uncertainties of ≤9 m. A further 16 locations had positional uncertainties of 10–17 m, and 6 locations had positional uncertainties of 2,900–5,000 m. We took a decision not to exclude any data based on positional uncertainties for this study. We used robust estimators of the variograms (Extended Data Fig. 2), which are resistant to effects of spatial outliers due to a small number of points being in the wrong position53 and these models were validated. Any effect of position error on the broad mapped pattern of long-range spatial variation at these scales will therefore be very limited and localized.
Whole-grain samples were air-dried in their sample bags. Each sample was then ground in a domestic stainless-steel coffee grinder, which was wiped clean before use and after each sample with a non-abrasive cloth. Whole grains are consumed in many settings, although more-refined fractions—with concomitant losses of more micronutrient-rich bran and endosperm fractions—are often consumed. All preparation was done away from sources of contamination by soil or by dust. A 20-g subsample of the ground material was then shipped to the University of Nottingham. Soil samples were oven-dried at 40 °C for 24–48 h depending on the moisture content of the soil. Preparation took place in a soil laboratory to avoid cross-contamination with grain samples. Plant material was removed from each soil sample, which was then disaggregated and sieved to pass 2 mm. This material was then coned and quartered to produce subsample splits. A 150-g subsample of soil was poured into a self-seal bag, labelled and shipped to the UK for analysis in the laboratories at the University of Nottingham and Rothamsted Research. Soil pH(water) and soil organic carbon (SOC) content were measured using standard methods22.
Grain micronutrient analyses
Grain micronutrient analysis methods followed standard methods54. Approximately 0.2 g of each ground sample was weighed and digested using a microwave system. For samples collected in the Amhara region in 2017, a Multiwave 3000 48-vessel MF50 rotor (Anton Paar) was used; digestion vessels were perfluoroalkoxy tubes in polyethylethylketone pressure jackets (Anton Paar). Samples were digested in 2 ml 70% trace-analysis-grade HNO3, 1 ml Milli-Q water (18.2 MΩ cm; Fisher Scientific) and 1 ml H2O2. Settings were: 1,400 W, 140 °C, 2 MPa, for 45 min. For samples collected in 2018–2019, we used a Multiwave Prom 41HVT56 rotor and pressure-activated venting vessels made of modified polytetrafluoroethylene (56-ml ‘SMART VENT’, Anton Paar). Samples were digested in 6 ml of 70% trace-analysis-grade HNO3. Settings were: 1,500 W, 10 min heating to 140 °C, 20 min holding at 140 °C, and 15 min cooling to 55 °C. Two operational blanks were typically included in each digestion run. Duplicate samples of a certified reference material (Wheat flour SRM 1567b, NIST) were included in approximately every fourth digestion run. Following digestion, each tube was made up to a final volume of 15 ml by adding 11 ml Milli-Q water, then transferred to a 25-ml universal tube (Sarstedt) and stored at room temperature. Samples were further diluted 1:5 with Milli-Q water into 13-ml tubes (Sarstedt) before analysis.
Multi-elemental analysis of grain (Ag, Al, As, B, Ba, Be, Ca, Cd, Cr, Co, Cs, Cu, Fe, K, Li, Mg, Mn, Mo, Na, Ni, P, Pb, Rb, S, Sr, Ti, Tl, U, V and Zn) was undertaken using inductively coupled plasma mass spectrometry (ICP–MS; iCAPQ, Thermo Fisher Scientific). The instrument used a helium collision cell with kinetic energy discrimination to reduce polyatomic interferences. Samples were introduced from an autosampler incorporating an ASXpress rapid uptake module (Cetac ASX-520, Teledyne Technologies) through a perfluoroalkoxy Microflow PFA-ST nebuliser (Thermo Fisher Scientific). Internal standards were introduced to the sample stream on a separate line via the ASXpress unit and included Sc (20 μg l–1), Rh (10 μg l–1), Ge (10 μg l–1) and Ir (5 μg l–1) in 2% HNO3 (Primar Plus grade; Fisher Scientific). An external multi-element calibration standard (Claritas-PPT grade CLMS-2; SPEX Certiprep) was used to calibrate Ag, Al, As, B, Ba, Be, Cd, Ca, Co, Cr, Cs, Cu, Fe, K, Li, Mg, Mn, Mo, Na, Ni, P, Pb, Rb, S, Se, Sr, Tl, U, V and Zn, in the range of 0–100 μg l–1 (0, 20, 40, 100 μg l–1). A bespoke external multi-element calibration solution (PlasmaCAL, SCP Science) was used to create Ca, K, Mg and Na standards in the range of 0–30 mg l–1. B, P and S calibration used in-house standard solutions (KH2PO4, K2SO4 and H3BO3); Ti was determined semiquantitatively. Sample processing was undertaken using Qtegra software (Thermo Fisher Scientific) with external cross-calibration between pulse-counting and analogue detector modes when required. Se was determined separately using a triple quadrupole ICP–MS (iCAP TQ; Thermo Fisher Scientific) using an oxygen cell to mass shift the isotope 80Se to m/z 96 (80Se16O) to reduce interference from the 40Ar dimer. Drift correction was achieved using Rh as an internal standard; calibration used the CLMS-2 multi-element standard (Certiprep).
Analyses were conducted in batches of around 240 samples per run on the ICP–MS instrument (Extended Data Table 1). Individual grain concentration data were corrected for run-specific operational blanks and then converted to concentration on a dry-matter basis. Element-specific limits of detection were reported as 3× the standard deviation of the operational blank concentrations, assuming a notional starting dry weight of 0.2 g of sample material (Extended Data Table 1). For samples for which the grain element concentration was less than the limit of detection, data were removed before the statistical analyses. No adjustment was made for potential contamination of grain samples, for example, with soil dust from the field or store using typical markers (for example, Fe, V). Two sorghum samples, taken from a grain store in Malawi, were excluded from the data analysis based on high concentrations of Ca, Mg and other elements that were considered unlikely to have arisen from soil contamination.
Summary statistics were computed for concentrations of all elements in grain and histograms were examined. It is common for geochemical variables to be log-normally distributed, and the coefficient of skewness of the data was examined using octile skewness as a robust measure of asymmetry of the distribution55. The data were analysed on the original scales of measurement (mg kg–1) if the octile skewness was within the interval [−0.2, 0, 2] as described previously56. If the octile skewness fell outside this range (with a positive value), and the absolute value of the octile skewness after loge transformation was smaller than on the original scale, then the data were analysed on the loge scale (Extended Data Tables 1, 2).
Variograms were estimated for each variable using three estimators (Extended Data Table 2): the standard estimator57 and two alternatives58,59. The two alternative estimators are more robust than the standard estimator to outlying data, be these marginal outliers (apparent on the histogram) or spatial outliers (apparent in local spatial context). The variogram estimates were formed from all pairwise comparisons among observations up to a maximum lag distance (330 km in Ethiopia and 100 km in Malawi), and with lag bins with a width of 10 km. The distance between any two locations (specified by latitude and longitude) was computed as the great circle distance using the distVincentySphere function from the geosphere library of the R platform47,60. Exponential variogram functions were then fitted to the estimates by weighted least squares61. The exponential variogram was selected because it ensures positive definite covariance matrices for distances on the sphere62. Each model was then tested by cross-validation. The selected models are shown in Extended Data Fig. 2.
In cross-validation, each observation was removed and predicted from the remaining observations by ordinary kriging. This was done using each variogram model fitted for each variable to the sets of point estimates computed by the three estimators. The standardized squared prediction error (SSPE) was computed for each observation as the squared difference between the observed and predicted values divided by the prediction error variance (kriging variance) as computed in the standard ordinary kriging equations53. The median value of the SSPE has an expected value of 0.455 in the case of a valid underlying variogram model with normally distributed kriging errors53. The standard estimator due to Matheron57 is more statistically efficient than the robust alternatives, so if the model fitted to these estimates appeared correct from the cross-validation results (the median SSPE is within the 95% confidence interval) then the alternatives were not considered. If the SSPE suggests that the model fitted to Matheron estimates are affected by outliers, then the models fitted to robust estimates were also cross-validated, and one was selected on the cross-validation results53.
Once a variogram model was selected it was used to compute predictions of the grain concentration (with or without transformation) on nodes of a fine square grid. This was done by ordinary kriging61. In ordinary kriging, it is assumed that the mean value of the variable of interest is locally constant, but unknown. An estimate is found that is a linear combination of the observations such that the expected squared error of the prediction (the ordinary kriging variance) is minimized. The kriging variance is also computed as a measure of the uncertainty of the prediction. For a given variogram, the kriging variance reflects the proximity of the location at which a prediction is made to observations of the variable of interest. In the case of variables that had been transformed to logarithms, the prediction at each location was computed by exponentiation of the prediction on the loge scale. The resulting back-transformed estimate is a median-unbiased predictor, which is appropriate for variables with a skewed distribution63. Note, the kriging variance cannot be meaningfully back-transformed. However, it still indicates how the uncertainty of the predictions varies in space and so it is mapped here on the loge scale for transformed variables. The kriging variance maps are shown in Extended Data Figs. 3, 4. In Ethiopia, where the sample sites are irregularly distributed over the mapped area, the kriging variance differs markedly depending on the local sample density.
In Ethiopia, we mapped the percentage of dietary requirement potentially met by the intake of wheat and teff (Fig. 3); similarly, in Malawi for maize (Fig. 4). In Ethiopia, 97.6 g per capita per day of wheat and 89.3 g per capita per day for teff (based on the Food and Agriculture Organization database item, ‘other cereals’) were used as reference intakes29; in Malawi, 342.8 g per capita per day was used as a reference maize intake29. Estimated average requirement (EAR) thresholds of 860, 22.4, 45 and 10.2 mg per capita per day for Ca, Fe, Se and Zn, respectively, were chosen as a representative threshold, based on an adult woman aged 18–24 years eating an unrefined (that is, high phytate) diet30. These thresholds are similar to other demographic groups. Because a comparable measure of uncertainty to the kriging variance for this derived variable is not available, we defined a mask for Ethiopia, where the sparsity of sampling leads to larger uncertainty compared with Malawi. We considered the Zn concentration in teff, a variable with spatial dependence over long distances, and identified those areas in which the kriging variance for this variable exceeded 75% of the variance of the variable itself because the sampling was sparse. These areas defined the mask used for the maps of the percentage of dietary requirement for all variables and is shown in grey in Fig. 3.
The maps shown in Figs. 3, 4 are made by point kriging—that is, they are predictions of a measurement at an unsampled site. Ordinary kriging can also be used to predict the mean value of a variable across a region or ‘block’61. To illustrate the communication value of this approach, the mean concentrations of Zn and of Ca in grain (teff and wheat) are shown at the woreda level in the Amhara region (Fig. 5). These were obtained by block kriging of woreda means from the same sample data, and with the same variogram models as used to produce the point kriging predictions.
Determining the links between a dietary biomarker of status and grain micronutrient concentration focused here on Se, for which there is a reliable biomarker of dietary status13; this is not the case for Zn7 and the other micronutrients analysed in this study. Data on the concentration of Se in the blood of women of reproductive age were available from micronutrient surveys in Ethiopia (serum23) and Malawi (plasma24). Mean concentrations were computed for each enumeration area available—321 in Ethiopia and 101 in Malawi—for which the latitude and longitude for each enumeration area centroid were available. For each enumeration area, the nearest grain sample site (regardless of crop) was found. In Ethiopia, the median distance to the nearest grain sample site over all enumeration areas was 17 km. In Malawi, the median distance was much shorter (1.4 km), which is attributable both to the denser crop sampling in Malawi and to the fact that enumeration areas used for the micronutrient survey were targeted for sampling.
The enumeration areas do not comprise a simple independent random sample, so the relation between the serum or plasma Se concentration and the Se concentration in the nearest grain sample could not be quantified by a statistic such as the correlation coefficient. It was therefore studied with an appropriate LMM incorporating a spatially correlated random effect, modelled with a Matérn correlation function64. Exploratory analysis indicated that a loge transformation of all serum or plasma Se concentrations was necessary for the assumption of normal random effects to be plausible. The observed serum or plasma Se concentration was modelled as a linear function of the concentration in grain at the nearest sample location. Residual maximum likelihood was used to estimate the random effects parameters. The fixed effects were then estimated by weighted least squares65 along with their standard errors. The evidence against the null hypothesis that the regression coefficient for serum or plasma Se on grain Se concentrations was zero was tested by a log-likelihood ratio test.
Plots of the serum or plasma Se concentration in women of reproductive age, in Ethiopia and Malawi, respectively, show a positive correlation between the variables (Extended Data Fig. 5), which is supported by the statistical models. For Ethiopia, the estimated regression coefficient (log[ng serum Se per ml] and log[mg of grain Se per kg]) was 0.08 with a standard error of 0.02. The null hypothesis that the coefficient was zero is rejected on the grounds of the log-likelihood ratio statistic (L = 14.48, P = 1.4 × 10−4). If there was no relationship between the biomarker and the grain Se concentrations, the probability of obtaining an L-statistic this large or larger would be very small. Similarly, for Malawi, the estimated regression coefficient was 0.09 with a standard error of 0.03 (L = 11.56, P = 6.7 × 10−4).
Data on the concentration of Se and of Zn in grain (maize in Malawi; teff, wheat and maize in Ethiopia) were extracted along with the corresponding data on soil pH(water) and SOC. Observations for three environmental covariates were also extracted for these same locations: (1) mean annual temperature; (2) mean annual precipitation values from the CHELSA datasets66,67, which are downscaled to a spatial resolution of 30 s; (3) topographic index values from the 30-s resolution MERIT Digital Elevation Model68. The topographic index—which is sometimes called the topographic wetness index—is a measure of the tendency for drainage to accumulate at a site and is widely used as a predictive measure for soil properties. Following exploratory analysis, grain Se concentration data were loge-transformed to make the assumption of normal random effects plausible; this was not necessary for grain Zn. Measurements of SOC from Malawi showed a marked positive skew and were therefore loge-transformed.
Data were analysed using a LMM in which a regression-type function of environmental covariates was considered as a fixed effect along with a spatially autocorrelated random effect, as described for the biomarker data. Random-effect parameters were estimated by maximum likelihood, and then the fixed-effect parameters by weighted least squares. The first model was fitted with a constant mean as the only fixed effect. The second model was fitted with mean annual precipitation as a fixed effect. The evidence for this predictor was evaluated by a log-likelihood ratio test of the second model against the first. The predictor was retained initially if the probability of obtaining a value of L as large or larger than observed under the null hypothesis of no effect of the predictor was <0.05. Additional predictors were then considered in the order of mean annual temperature then topographic index. Finally, the P values were evaluated as the outcome of multiple tests, controlling the false-discovery rate (FDR)69 at 0.05. The FDR is the expected proportion of incorrectly rejected null hypotheses among all rejected ones. For each micronutrient in each crop and country, the predictors were identified that could be regarded as significant with FDR control at 0.05. The same procedure was followed to produce a comparable model based on the two soil properties, considering first soil pH and then SOC.
Plots of grain Se and Zn concentrations against the environmental covariates and soil properties are shown in Extended Data Figs. 6–9. The evidence for environmental covariates or soil properties as predictors of micronutrient content in the grain is summarized in Extended Data Table 3. Extended Data Table 3 also shows the estimated coefficients, and their standard errors, in separate models for each micronutrient; for maize, teff or wheat in Ethiopia and for maize in Malawi. The predictors in these models are only those that were selected with FDR control.
Further information on research design is available in the Nature Research Reporting Summary linked to this paper.
All data are freely available from the corresponding author and available online at https://github.com/rmlark/GeoNutrition.
All code is freely available from the corresponding author and available online at https://github.com/rmlark/GeoNutrition.
Wessells, K. R. & Brown, K. H. Estimating the global prevalence of zinc deficiency: results based on zinc availability in national food supplies and the prevalence of stunting. PLoS ONE 7, e50568 (2012).
Joy, E. J. M. et al. Dietary mineral supplies in Africa. Physiol. Plant. 151, 208–229 (2014).
Kumssa, D. B. et al. Dietary calcium and zinc deficiency risks are decreasing but remain prevalent. Sci. Rep. 5, 10974 (2015).
Ligowe, I. S. et al. Selenium deficiency risks in sub-Saharan African food systems and their geospatial linkages. Proc. Nutr. Soc. 79, 457–467 (2020).
Joy, E. J. M. et al. Dietary mineral supplies in Malawi: spatial and socioeconomic assessment. BMC Nutrition 1, 42 (2015).
Nelson, G. et al. Income growth and climate change effects on global nutrition security to mid-century. Nat. Sustain. 1, 773–781 (2018).
King, J. C. et al. Biomarkers of nutrition for development (BOND)—zinc review. J. Nutr. 146, 858S–885S (2015).
Das, J. K., Salam, R. A., Kumar, R. & Bhutta, Z. A. Micronutrient fortification of food and its impact on woman and child health: a systematic review. Syst. Rev. 2, 67 (2013).
Marks, K. J. et al. Review of grain fortification legislation, standards, and monitoring documents. Glob. Health Sci. Pract. 6, 356–371 (2018).
White, P. J. & Broadley, M. R. Biofortification of crops with seven mineral elements often lacking in human diets—iron, zinc, copper, calcium, magnesium, selenium and iodine. New Phytol. 182, 49–84 (2009).
Bouis, H. E. & Saltzman, A. Improving nutrition through biofortification: a review of evidence from HarvestPlus, 2003 through 2016. Glob. Food Sec. 12, 49–58 (2017).
Gödecke, T., Stein, A. J. & Qaim, M. The global burden of chronic and hidden hunger: trends and determinants. Glob. Food Sec. 17, 21–29 (2018).
Fairweather-Tait, S. J. et al. Selenium in human health and disease. Antioxid. Redox Signal. 14, 1337–1383 (2011).
Caulfeld, L. E., Richard, S. A., Rivera, J. A., Musgrove, P. & Black, R. E. in Disease Control Priorities in Developing Countries 2nd edn (eds Jamison, D. T. et al.) 551–567 (The International Bank for Reconstruction and Development/The World Bank and Oxford Univ. Press, 2006).
Hurst, R. et al. Soil-type influences human selenium status and underlies widespread selenium deficiency risks in Malawi. Sci. Rep. 3, 1425 (2013).
Bermudez, O. I., Lividini, K., Smitz, M. F. & Fiedler, J. L. Estimating micronutrient intakes from Household Consumption and Expenditures Surveys (HCES): an example from Bangladesh. Food Nutr. Bull. 33 (Suppl), S208–S213 (2012).
Bevis, L. E. & Hestrin, R. Widespread heterogeneity in staple crop mineral concentration in Uganda partially driven by soil characteristics. Environ. Geochem. Health 43, 1867–1889 (2021).
Vanderpas, J. B. et al. Iodine and selenium deficiency associated with cretinism in northern Zaire. Am. J. Clin. Nutr. 52, 1087–1093 (1990).
Fordyce, F. in Essentials of Medical Geology (eds Selinus, O. et al.) 373–415 (Elsevier, 2005).
Faye, B., Grillet, C., Tessema, A. & Kamil, M. Copper deficiency in ruminants in the rift valley of east Africa. Trop. Anim. Health Prod. 23, 172–180 (1991).
Lark, R. M. et al. Mapping trace element deficiency by cokriging from regional geochemical soil data: a case study on cobalt for grazing sheep in Ireland. Geoderma 226–227, 64–78 (2014).
Gashu, D. et al. Spatial prediction of the concentration of selenium (Se) in grain across part of Amhara Region, Ethiopia. Sci. Tot. Environ. 733, 139231 (2020).
Belay, A. et al. Selenium deficiency is widespread and spatially dependent in Ethiopia. Nutrients 12, 1565 (2020).
Phiri, F. P. et al. The risk of selenium deficiency in Malawi is large and varies over multiple spatial scales. Sci. Rep. 9, 6566 (2019).
Chilimba, A. D. C. et al. Maize grain and soil surveys reveal suboptimal dietary selenium intake is widespread in Malawi. Sci. Rep. 1, 72 (2011).
Joy, E. J. M. et al. Soil type influences crop mineral composition in Malawi. Sci. Total Environ. 505, 587–595 (2015).
Badigannavar, A. & Ganapathi, T. R. Genetic variability for mineral nutrients in indigenous germplasm lines of finger millet (Eleusine coracana Gaertn.). J. Cereal Sci. 84, 1–6 (2018).
Gabaza, M., Shumoy, H., Muchuweti, M., Vandamme, P. & Raes, K. Enzymatic degradation of mineral binders in cereals: impact on iron and zinc bioaccessibility. J. Cereal Sci. 82, 223–229 (2018).
FAO. FAOSTAT Database 2017 http://faostat.fao.org/ (FAO, accessed 30 May 2020).
Allen, L. H., Carriquiry, A. L. & Murphy, S. P. Perspective: proposed harmonized nutrient reference values for populations. Adv. Nutr. 11, 469–483 (2020).
Gashu, D., Marquis, G. S., Bougma, K. & Stoecker, B. J. Spatial variation of human selenium in Ethiopia. Biol. Trace Elem. Res. 189, 354–360 (2019).
Allen, L. H. Iron supplements: scientific issues concerning efficacy and implications for research and programs. J. Nutr. 132 (Suppl), 813S–819S (2002).
Velu, G. et al. Reaching out to farmers with high zinc wheat varieties through public-private partnerships – an experience from Eastern-Gangetic Plains of India. Adv. Food Technol. Nutr. Sci. Open J. 1, 73–75 (2015).
Zia, M. H. et al. Site-specific factors influence the field performance of a Zn-biofortified wheat variety. Front. Sustain. Food Syst. 4, 135 (2020).
Maqbool, A. M. & Beshir, A. Zinc biofortification of maize (Zea mays L.): status and challenges. Plant Breed. 138, 1–28 (2019).
Barrett, C. B. & Bevis, L. E. The self-reinforcing feedback between low soil fertility and chronic poverty. Nat. Geosci. 8, 907–912 (2015).
Bevis, L. E. M. & Barrett, C. B. Close to the edge: high productivity at plot peripheries and the inverse size-productivity relationship. J. Dev. Econ. 143, 102377 (2020).
Manzeke, M. G. et al. Fertilizer management and soil type influence grain zinc and iron concentration under contrasting smallholder cropping systems in Zimbabwe. Sci. Rep. 9, 6445 (2019).
Wood, S. A., Tirfessa, D. & Baudron, F. Soil organic matter underlies crop nutritional quality and productivity in smallholder agriculture. Agric. Ecosyst. Environ. 266, 100–108 (2018).
Smith, M. R. & Myers, S. S. Impact of anthropogenic CO2 emissions on global human nutrition. Nat. Clim. Change 8, 834–839 (2018).
Jones, G. D. et al. Selenium deficiency risk predicted to increase under future climate change. Proc. Natl Acad. Sci. USA 114, 2848–2853 (2017).
Köhler, I. H., Huber, S. C., Bernacchi, C. J. & Baxter, I. R. Increased temperatures may safeguard the nutritional quality of crops under future elevated CO2 concentrations. Plant J. 97, 872–886 (2019).
Lark, R. M. & Marchant, B. P. How should a spatial-coverage sample design for a geostatistical soil survey be supplemented to support estimation of spatial covariance parameters? Geoderma 319, 89–99 (2018).
AfSIS. New cropland and rural settlement maps of Africa. http://africasoils.net/2015/06/07/new-cropland-and-rural-settlement-maps-of-africa/ (2015; accessed 10 January 2020).
Walsh, M. & Wu, W. GeoSurvey data prediction workflows. Open Science Framework (OSF) Repository https://osf.io/vxc97/ (2020).
ESA. Land Cover CCI Data version 2.0.7, land cover data for 2015 http://maps.elie.ucl.ac.be/CCI/viewer/download.php (European Space Agency Climate Change Initiative, accessed 21 November 2017).
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2017).
Grafström, A. & Lisic, J. BalancedSampling: balanced and spatially balanced sampling. R package version 1.5.2 https://CRAN.R-project.org/package=BalancedSampling (2016).
Deville, J. C. & Tillé, Y. Efficient balanced sampling: the cube method. Biometrika 91, 893–912 (2004).
Grafström, A. & Schelin, L. How to select representative samples. Scand. J. Stat. 41, 277–290 (2014).
NSO Malawi & ICF. Malawi Demographic and Health Survey 2015–16 (NSO and ICF, 2017).
Walvoort, D. J. J., Brus, D. J. & de Gruijter, J. J. An R package for spatial coverage sampling and random sampling from compact geographical strata by k-means. Comput. Geosci. 36, 1261–1267 (2010).
Lark, R. M. A comparison of some robust estimators of the variogram for use in soil survey. Eur. J. Soil Sci. 51, 137–157 (2000).
Khokhar, J. S. et al. Variation in grain Zn concentration, and the grain ionome, in field-grown Indian wheat. PLoS ONE 13, e0192026 (2018).
Brys, G., Hubert, M. & Struyf, A. A robust measure of skewness. J. Comput. Graph. Stat. 13, 996–1017 (2004).
Rawlins, B. G., Lark, R. M., O’Donnell, K. E., Tye, A. M. & Lister, T. R. The assessment of point and diffuse metal pollution of soils from an urban geochemical survey of Sheffield, England. Soil Use Manage. 21, 353–362 (2005).
Matheron, G. Traité de Géostatistique Appliqué, Tome 1. Memoires du Bureau de Recherches Géologiques et Minières (Technip, 1962).
Cressie, N. & Hawkins, D. M. Robust estimation of the variogram: I. Math. Geol. 12, 115–125 (1980).
Dowd, P. A. in: Geostatistics for Natural Resources Characterization (eds Verly, G. et al.) 91–106 (Reidel, 1984).
Hijmans, R. J. Geosphere: spherical trigonometry. R package version 1.5-7 https://CRAN.R-project.org/package=geosphere (2017).
Webster, R. & Oliver, M. A. Geostatistics for Environmental Scientists 2nd edn (John Wiley & Sons, 2007).
Gneiting, T. Strictly and non-strictly positive definite functions on spheres. Bernoulli 19, 1327–1349 (2013).
Pawlowsky-Glahn, V. & Olea, R. A. Geostatistics for Computational Data (Oxford Univ. Press, 2004).
Diggle, P. & Ribeiro, P. J. Model-Based Geostatistics (Springer, 2007).
Lark, R. M., Cullis, B. R. & Welham, S. J. On spatial prediction of soil properties in the presence of a spatial trend: the empirical best linear unbiased predictor (E-BLUP) with REML. Eur. J. Soil Sci. 57, 787–799 (2006).
Karger, D. N. et al. Climatologies at high resolution for the Earth’s land surface areas. Sci. Data 4, 170122 (2017).
Karger, D. N. et al. Data from: climatologies at high resolution for the Earth’s land surface areas. Dryad https://doi.org/10.5061/dryad.kd1d4 (2017).
Yamazaki, D. et al. A high accuracy map of global terrain elevations. Geophys. Res. Lett. 44, 5844–5853 (2017).
Lark, R. M. Controlling the marginal false discovery rate in inferences from a soil dataset with α‐investment. Eur. J. Soil Sci. 68, 221–234 (2017).
This work was supported by GeoNutrition projects funded by the Bill & Melinda Gates Foundation (INV-009129) and the UKRI Biotechnology and Biological Sciences Research Council (BBSRC)/Global Challenges Research Fund (GCRF) (BB/P023126/1). The funders were not involved in the study design; the collection, management, analysis, and interpretation of data; the writing of the paper or the decision to submit the paper for publication. The boundaries, denominations and any other information shown on these maps do not imply any judgment about the legal status of any territory or constitute any official endorsement or acceptance of any boundaries on the part of any Government. We acknowledge the contributions made to this research by the participating farmers and field sampling teams. In Ethiopia, field sampling teams were from the Amhara, Oromia and Tigray Regional Bureau of Agriculture. In Malawi, field sampling teams were from the Department of Agricultural Research Services and Lilongwe University of Agriculture and Natural Resources. Training support for field activities in Malawi was facilitated by funding from the Royal Society-UK Foreign, Commonwealth & Development Office (FCDO), under project AQ140000, “Strengthening African capacity in soil geochemistry for agriculture and health”. Mineral analytical support was provided by B. Broadley, K. Davis, P. Muleya, S. Vasquez Reina, S. Dunham, J. Carter and J. Hernandez with support from BBSRC Institute Strategic Project Soil to Nutrition (BBS/E/C/000I0310). E.L.A.’s contribution is published with the permission of the Director of the British Geological Survey (UKRI). The CHELSA project is acknowledged for making the downscaled climate data available from https://climatedataguide.ucar.edu/.
The authors declare no competing interests.
Peer review information Nature thanks Christopher Barrett, Ken Giller, Bernard Vanlauwe and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Target layout of the five sample points (black circles) for a standing crop (left) or for a crop harvested in the field (right).
a, Teff in Ethiopia. b, Wheat in Ethiopia. c, Maize in Malawi.
a–h, Teff (a–d) and wheat (e–h) grain concentration of Ca (a, e), Fe (loge) (b, f), Se (loge) (c, g) and Zn (d, h).
a, Ca (loge). b, Fe (loge). c, Se (loge). d, Zn.
Extended Data Fig. 5 Relationships between Se concentration in blood fractions and grain Se concentrations.
a, Ethiopia (serum). b, Malawi (plasma).
Extended Data Fig. 6 Relationships between grain Se concentration and environmental and soil covariates in Ethiopia.
Data for teff (top), wheat (middle) and maize (bottom) are shown.
Extended Data Fig. 7 Relationships between grain Zn concentration and environmental and soil covariates in Ethiopia.
Data for Teff (top), wheat (middle) and maize (bottom) are shown.
Extended Data Fig. 8 Relationships between grain Se concentration and environmental and soil covariates in Malawi.
Data for maize are shown.
Extended Data Fig. 9 Relationships between grain Zn concentration and environmental and soil covariates in Malawi.
Data for maize are shown.
About this article
Cite this article
Gashu, D., Nalivata, P.C., Amede, T. et al. The nutritional quality of cereals varies geospatially in Ethiopia and Malawi. Nature 594, 71–76 (2021). https://doi.org/10.1038/s41586-021-03559-3
Population and Environment (2021)
Huangjing (Polygonati rhizoma) is an emerging crop with great potential to fight chronic and hidden hunger
Science China Life Sciences (2021)