GCRI · Food Nexus

Spatial statistics for food security

Spatial statistics turns sparse surveys, soil samples and market prices into maps with stated uncertainty. In food systems it estimates food insecurity for small areas, predicts soil properties between sampled profiles and tests whether price anomalies or malnutrition cluster in space.

01Register

Spatial statistics for food security on the record.

Method recordSTEM method
  1. MethodSpatial statistics and geostatistics
  2. Applied inFood Nexus
  3. Canonical pageSTEM Commons
02Application

What spatial statistics and geostatistics does here.

Spatial statistics and geostatistics analyse data that vary across space, accounting for spatial dependence. They estimate values at unsampled locations and test for spatial patterns. In food and agriculture, trade, resources and supply chains and climate and hazards, spatial statistics produce continuous hazard and exposure surfaces with uncertainty. They reveal geographic clustering of risks and vulnerabilities across sectors. The canonical description is on the STEM Commons.

Small-area estimation of food insecurity

Household survey results are combined with satellite, census and market covariates to estimate food insecurity for districts the survey sampled too thinly, with an uncertainty for each estimate.

Digital soil mapping

Soil properties measured at sampled profiles are predicted across the landscape from terrain, climate and spectral covariates by regression kriging or machine learning, with prediction intervals mapped.

Clustering of malnutrition and price anomalies

Spatial autocorrelation and hotspot tests show whether high malnutrition or price anomalies cluster, before any cause is proposed.

03Geography

Growing conditions in the Black Sea region.

Crop greenness, root-zone soil moisture and daytime surface heat over the Black Sea region, from NASA GIBS. Each layer carries its own date.

The interactive map needs JavaScript; the layer list and sources are below.

Layers 3

Growing conditions, for this page

About this layer
Evidence
Observed (clouds leave gaps)
Source
MODIS Terra NDVI (NASA LANCE)
Provider
NASA LANCE
Licence
NASA open data
Vintage
daily; Satellite index; clouds leave gaps
Newest at build
2026-10-13
Resolution
about 250 m grid
Uncertainty
Satellite greenness index over the 8 days to the date shown; an indicator of crop and pasture condition, not a yield estimate or a food security warning.
Reached through
gibs.earthdata.nasa.gov, only after you allow the live map
About this layer
Evidence
Modelled
Source
SMAP Level 4 (NASA NSIDC DAAC)
Provider
NASA NSIDC DAAC
Licence
NASA open data
Vintage
daily; Model with satellite data assimilated
Newest at build
2026-10-11
Resolution
about 2 km grid
Uncertainty
Model with satellite data assimilated. Satellite or model raster as the provider publishes it. Clouds, swath gaps and retrieval error leave gaps and noise; values are indicators, not measurements at a place.
Reached through
gibs.earthdata.nasa.gov, only after you allow the live map
About this layer
Evidence
Observed (persistent cloud still leaves gaps)
Source
VIIRS NOAA-20 Level 3, 8-day (NASA LANCE)
Provider
NASA LANCE
Licence
NASA open data
Vintage
8-day; Satellite 8-day composite; persistent cloud still leaves gaps
Newest at build
2026-09-14
Resolution
about 1 km grid
Uncertainty
Satellite skin temperature of the ground over the 8 days ending on the date shown, a heat exposure indicator; not air temperature, heat index or a heat warning.
Reached through
gibs.earthdata.nasa.gov, only after you allow the live map
About this layer
Evidence
Observed
Source
VIIRS NOAA-20 corrected reflectance (NASA LANCE)
Provider
NASA LANCE
Licence
NASA open data
Vintage
daily; Satellite image
Newest at build
2026-10-14
Resolution
about 250 m grid
Uncertainty
Satellite image. Satellite or model raster as the provider publishes it. Clouds, swath gaps and retrieval error leave gaps and noise; values are indicators, not measurements at a place.
Reached through
gibs.earthdata.nasa.gov, only after you allow the live map
About this layer
Evidence
Modelled (latest year GIBS publishes)
Source
MODIS Terra and Aqua land cover (MCD12Q1) (NASA LP DAAC)
Provider
NASA LP DAAC
Licence
NASA open data
Vintage
yearly; Satellite classification; latest year GIBS publishes
Newest at build
2024-01-01
Resolution
about 500 m grid
Uncertainty
Satellite land cover classes for the year shown, including croplands; a classification with known errors, not a farm register.
Reached through
gibs.earthdata.nasa.gov, only after you allow the live map
About this layer
Evidence
Modelled (latest 8-day period GIBS publishes)
Source
MODIS Aqua GPP (MYD17) (NASA LP DAAC)
Provider
NASA LP DAAC
Licence
NASA open data
Vintage
8-day; Satellite-driven model; latest 8-day period GIBS publishes
Newest at build
2026-09-30
Resolution
about 500 m grid
Uncertainty
Model estimate of carbon taken up by plants over the 8-day period shown; an indicator of growth, not harvest or production.
Reached through
gibs.earthdata.nasa.gov, only after you allow the live map

Latest available

NASA GIBS / EOSDIS · ASTER GDEM (METI, NASA) · Blue Marble (NASA) · © OpenStreetMap contributors (ODbL) · Natural Earth

Loading the latest NASA imagery.

Drag to move. Ctrl or ⌘ and scroll, pinch or the buttons to zoom; hover for the values under the pointer; reset returns to the region.

Growing conditions, framed on the Black Sea region.Satellite and model products as NASA publishes them, each with its own date. They are indicators, not assessments of any place.

Data: NASA GIBS · ASTER GDEM · Blue Marble · © OpenStreetMap contributors · Natural Earth · NASA EONET · NASA JPL OPERA DIST-ALERT · NASA SMAPSources and licences

  • We acknowledge the use of imagery provided by services from NASA’s Global Imagery Browse Services (GIBS), part of NASA’s Earth Science Data and Information System (ESDIS). Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). NASA open data
  • We acknowledge the use of imagery provided by services from NASA’s Global Imagery Browse Services (GIBS), part of NASA’s Earth Science Data and Information System (ESDIS). Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). NASA open data
  • MODIS land products (vegetation indices, gross primary production, land cover), NASA LP DAAC and LANCE; imagery via NASA GIBS. Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). NASA open data
  • SMAP Level 4 soil moisture and snow mass, NASA Global Modeling and Assimilation Office; distributed by the NASA National Snow and Ice Data Center DAAC; imagery via NASA GIBS. Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). NASA open data
  • VIIRS near-real-time products (flood, land surface temperature, snow cover, corrected reflectance), NASA LANCE; imagery via NASA GIBS. Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). NASA open data
  • ASTER GDEM is a product of METI and NASA (ASTER Global Digital Elevation Model V003, NASA LP DAAC); imagery via NASA GIBS. Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). Cite as: NASA/METI/AIST/Japan Spacesystems and U.S./Japan ASTER Science Team (2019). ASTER Global Digital Elevation Model V003. doi:10.5067/ASTER/ASTGTM.003. NASA open data
  • Made with Natural Earth. Free vector and raster map data @ naturalearthdata.com. Public domain
  • © OpenStreetMap contributors (ODbL), via NASA GIBS (coastlines, boundaries, roads and place names). Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). ODbL
  • Blue Marble, NASA Earth Observatory; imagery via NASA GIBS. Colours redrawn from the published colour map; values averaged into hexagonal cells where shown (modified). NASA open data
  • Context on the globe and flat map: coastlines, boundaries, roads and place names from NASA GIBS reference layers (OpenStreetMap data, ODbL, © OpenStreetMap contributors); Blue Marble imagery as an optional base; 3D terrain and the elevation read-out from the ASTER GDEM elevation colour index (5 m to 200 m classes, exaggerated for legibility). Boundaries are drawn as published and imply no position on any territory.
  • Live events: NASA Earth Observatory Natural Event Tracker (EONET); events are indicative and keep their original source links.
  • Live: NASA EONET Events are indicative, not authority warnings.
  • Observed (EO): NASA JPL OPERA DIST-ALERT (contains modified Copernicus Sentinel data 2026). We acknowledge the use of imagery provided by services from NASA’s Global Imagery Browse Services (GIBS), part of NASA’s Earth Science Data and Information System (ESDIS).
  • Observed (EO): NASA SMAP. We acknowledge the use of imagery provided by services from NASA’s Global Imagery Browse Services (GIBS), part of NASA’s Earth Science Data and Information System (ESDIS).
  • NASA GIBS: We acknowledge the use of imagery provided by services from NASA’s Global Imagery Browse Services (GIBS), part of NASA’s Earth Science Data and Information System (ESDIS).
  • River discharge forecast, 30 days with range (GloFAS via Open-Meteo): Contains modified Copernicus Emergency Management Service information 2026; not a flood warning.
  • Past year, daily temperature (ERA5 via Open-Meteo): Contains modified Copernicus Climate Change Service information 2026.
  • Past year, daily temperature and sunshine (NASA POWER): These data were obtained from the NASA Langley Research Center (LaRC) POWER Project funded through the NASA Earth Science/Applied Science Program..
  • Weather data by Open-Meteo.com (https://open-meteo.com/), licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/).
  • Optional layers (13 sources in the layer list: live, forecast, observed and context) print their provider credit on the map and here when switched on; forecasts are model output and derived indices are estimates.

Layers and dates

  • Vegetation greenness (NDVI), rolling 8 days: MODIS Terra NDVI, NASA LANCE · 13 Oct 2026
  • Root-zone soil moisture, 0 to 100 cm: SMAP Level 4, NASA NSIDC DAAC · 11 Oct 2026
  • Land surface temperature, day, 8-day composite: VIIRS NOAA-20 Level 3, 8-day, NASA LANCE · 14 Sep 2026
  • True colour: VIIRS NOAA-20 corrected reflectance, NASA LANCE · 13 Oct 2026
  • Land cover type (IGBP), yearly: MODIS Terra and Aqua land cover (MCD12Q1), NASA LP DAAC · 2024
  • Gross primary production, 8-day total: MODIS Aqua GPP (MYD17), NASA LP DAAC · 30 Sep 2026

04Data

What the method consumes.

Datasets and records used by spatial statistics and geostatistics on this platform, with the role each plays and the label its values carry. Sources are listed in the source register.

DataRoleLabel
Household and nutrition surveysResponse variablesOBSERVED
Soil profile databasesPoint observationsOBSERVED
Satellite covariatesPredictorsMODELLED
Market price monitoringSpatial price fieldsPublished statistics
05Simulator

A shock through an input-output economy.

A loss of output in one sector propagates through the technical coefficients to its suppliers and customers, and recovers at the rate you set; the coefficients are teaching values, not national accounts. Teaching model with stated illustrative values, not a prediction. The map on this page shows The Black Sea region; the model itself is not tied to a place.

Controls

Sector hit first (position in the list)
0
Size of the shock
30% of its output
Duration of the shock
30 days
Recovery rate
0.08 per day
Days simulated
180 days

Readouts

Output lost over the run
0.65% of a year’s output
Total over direct (multiplier)
1.63 times
Most affected other sector
Energy
Its peak inoperability
3.8%
Numbers, methods and sources
SectorPeak inoperability (%)Output weight (%)
Farming30.4616.0
Fertiliser and inputs0.2711.2
Food processing0.3726.7
Transport0.5012.8
Wholesale and retail3.5416.6
Energy3.8016.7

Methods

  1. Input-output analysis Leontief (1936); Miller and Blair (2009), Input-Output Analysis, Cambridge University Press
  2. Inoperability input-output model Haimes and Jiang (2001), Journal of Infrastructure Systems 7
  3. Dynamic inoperability and recovery Lian and Haimes (2006), Systems Engineering 9
  4. Uncertainty fan (parameter Monte Carlo) Spiegelhalter, Pearson and Short (2011), Visualising uncertainty about the future, Science 333
  5. One-at-a-time sensitivity (tornado) Saltelli and others (2008), Global Sensitivity Analysis: The Primer, Wiley

Data

The model runs on the teaching values set by the controls. To calibrate it for a place, read it against these admitted layers and series, which this page’s map can draw:

  • Seaports Natural Earth
  • Maritime chokepoints Nexus (named straits and canals at their published coordinates)
  • Railways Natural Earth
  • Major roads Natural Earth
A shock through an input-output economy. Illustrative teaching model: it runs on the values set by the controls, not on measured data. Methods: Input-output analysis (Leontief (1936); Miller and Blair (2009), Input-Output Analysis, Cambridge University Press); Inoperability input-output model (Haimes and Jiang (2001), Journal of Infrastructure Systems 7); Dynamic inoperability and recovery (Lian and Haimes (2006), Systems Engineering 9).

Across the six pillars, led by Nexus Consortium

06Validation

How results are checked.

Predictions are tested using spatial cross-validation that respects dependence. Uncertainty maps are compared with independent observations.

CheckHow it is done
Spatial cross-validationBlocks of locations are withheld in turn so spatial dependence does not inflate apparent skill
Uncertainty mapsEvery published map has a companion map of its prediction error
Scale and zoning testsResults are re-run at different aggregation units to check that patterns are not artefacts of boundaries

Results depend on the scale and zoning of analysis. Sparse or clustered sampling limits prediction in unsampled areas.

07Twin and cascade simulator

Model spatial statistics for food security as a system. Run the cascade.

Assembled from the Nexus Standards area profile FAO and food security: the hazards it names, the systems and sectors that carry the loss, the controls that dampen it and the data that observe it. Choose the components, pick an initiating hazard and step through the cascade. The Map view places the run on an example geography with open data. Every coupling weight is a stated assumption you can change. A teaching model, not a forecast, an assessment or a certification. To model a place of your own choosing, open the full System Builder.

Hazards

  • Famine
  • Ice-jam flood
  • Tropical cyclone
  • Heatwave
  • HIV epidemic
  • Flash flood
  • Storm surge (set off downstream)
  • Intense rainfall (set off downstream)
  • River flood (set off downstream)
  • Dangerous surf (set off downstream)
  • Widespread power outage (set off downstream)
  • Mass displacement (set off downstream)
  • Coastal flood (set off downstream)
  • Extreme fire weather (set off downstream)
  • Wildfire (set off downstream)
  • Heat-driven grid stress (set off downstream)
  • Agricultural drought (set off downstream)
  • Flash drought (set off downstream)

Controls

  • NS-3 Risk record UNIs (grammar, legacy crosswalk, allocation authority)
  • NS-10 Policies and clauses (machine-readable obligations: definition, binding, evaluation and enforcement)
  • NRS-N3 Food and agriculture risk and resilience
  • NRS-C7 Early warning, triggers and anticipatory action
  • NRS-H61 Hazard family: Cascading failure (NHS-1 sys.cascade)
  • NRS-K40 Domain pack profile: Islands small states (DOMAIN-040)
  • NRS-S9 System profile: NXS-AAP (Anticipatory action and triggers)
  • NRS-H57 Hazard family: Humanitarian crises (NHS-1 soc.humanitarian)

Observation

  • FAOSTAT API (listed for reference)
  • GIEWS FPMA food price monitoring tool (listed for reference)
  • FEWS NET data (listed for reference)
  • Prevalence of undernourishment
  • Prevalence of moderate or severe food insecurity in the population, based on the Food Insecurity Experience Scale (FIES)
  • Indicator of food price anomalies

Example geography

  • The Ganges and Brahmaputra Basin: 178 hexagonal cells
  • Surface soil moisture (SMAP), NASA GIBS
  • Surface water through cloud, 30 m (OPERA DSWx-S1 radar), NASA GIBS
  • Rain rate, last half hour (IMERG), NASA GIBS
  • Surface air temperature, daytime (AIRS), NASA GIBS
  • Night lights tonight, outages (VIIRS Day/Night Band), NASA GIBS
09Sector signature

The sector signature of spatial statistics for food security.

The sectors spatial statistics for food security touches and how strongly they depend on each other: what a study has to observe and model together.

Select corners · hover lines6 / 15 relations

Signature

Food and its partners

Food gives nutrition to health (0.7), land conversion to biodiversity (0.6) and withdrawals and runoff to water (0.3).

W · E · I · F · H · B

Sectors
4 of 6
Relations
6 of 15
Coupling
0.46 share of total

Strongest relations in play

  • Food → Biodiversity · land conversion0.6
  • Water → Food · irrigation0.8
  • Water → Biodiversity · environmental flows0.6
  • Food → Health · nutrition0.7
Teaching values for a typical city-region, not measurements for any place. Coupling is the summed strength of the relations in play, as a share of all fifteen.
10Stewardship

Who maintains the method.

Space and Earth Observation stewards spatial statistics and geostatistics for the Consortium and reviews how it is applied across platforms.

11Methods and evidence

How spatial statistics for food security is modelled, measured and tested.

The model, its assumptions, its data and its uncertainty, in the form a scientist or engineer would check. Structures follow published methods. The values drawn are illustrative, not findings.

System model

The couplings between the systems in spatial statistics for food security: signs of influence, a reinforcing and a balancing loop, and cross-links where stress jumps.

++−++FoodTradeClimateCommunitiesEconomyR reinforcingB balancing
Causal loop diagram from the page’s systems; signs show the direction of influence. Illustrative structure for a system dynamics or network model, not estimated couplings.
FAIR guiding principles · GO FAIR and FORCE11W3C PROV · World Wide Web ConsortiumGuidance on consistent treatment of uncertainties · IPCCCEOS Analysis Ready Data · Committee on Earth Observation Satellites
13Questions

Questions and answers.

What are the limits of spatial statistics and geostatistics in this field?
Results depend on the scale and zoning of analysis. Sparse or clustered sampling limits prediction in unsampled areas.
How are results checked?
Predictions are tested using spatial cross-validation that respects dependence. Uncertainty maps are compared with independent observations.
Are the values on this page real?
The data names, roles and checks are real descriptions of practice. The canvas shows a design of the method’s use, not a run; no result is reported here.
Where else is this method applied?
The canonical method page lists every platform that applies it.
ReportsContribute

Contribute evidence from this page.

Reports is the Consortium’s quarterly publication. What you know about this subject can become evidence, a case or a correction in the next edition, credited to you and published under CC BY 4.0 where it is your own work.

GCRIFood NexusGlobal Sector Risk and ResilienceCountry Systems Risk and ResilienceEvidence Brief
Food Nexus

Bring a method, a dataset or a question.

Information and enquiries are open today. An enquiry receives a reference and names the next step; applications, payments, compute allocations and registrations open only when their team, terms and review capacity exist.

The Global Centre for Risk and Innovation

Science · evidence · measurement

Not for profit

Our articles, maps and twins are free to read and free to republish. Readers and partners keep it that way.

Partner with us

Who funds us

GCRI is a not-for-profit institution. To ask who funds a piece of work, use the funding questions form.

Independence

Editors and the expert reviewers named on each article decide what we publish. Funders and partners do not review or approve articles before publication.

Licence

Articles: CC BY 4.0. Figures, maps, twins and data: CC BY 4.0. Other uses: ask through the permissions form.