spatial production allocation - esri€¦ · 2015 esri international user conference --...

32
Not your usual SPAM but the Spatial Production Allocation Model Ulrike Wood-Sichra IFPRI, Washington DC ESRI User Conference San Diego, July 20-24, 2015

Upload: others

Post on 10-Aug-2020

6 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

Not your usual SPAM but the

Spatial Production Allocation Model

Ulrike Wood-SichraIFPRI, Washington DC

ESRI User ConferenceSan Diego, July 20-24, 2015

Page 2: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

SPAM – How it works and challenges encountered

1. What is SPAM2. The SPAM process3. Crop list and statistical data coverage4. Visualization of downscaling5. SPAM results6. Equations behind SPAM7. SPAM on the Web8. Challenges & Issues9. Validation10. Users

Page 3: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

1. SPAM ?

Page 4: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

…. continue

1. Spatial Production Allocation Model (SPAM)

Drawing on a variety of inputs • SPAM uses an entropy-based, data-fusion

approach to• “plausibly” assess cropping system distribution

and performance at a• “meso-gridded” scale: 5-minute globally• 30-seconds at country level (if data is available)

Page 5: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

1. Spatial Production Allocation Model (SPAM)Major input layers for SPAM• SPAM uses an entropy-based, data-fusion approach to• National and subnational production statistics: area,

yield (production)• Production system• Rural population density• Cropland extent• Irrigation map• Crop suitability• Existing crop distribution maps• Crop prices

…. continue

Page 6: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

Sub-national area (crop)

irrigated_ land

crop_land

suitable_area (crop)

crop distribution (crop)

SPAM

other data: production systems, ... (crop)

pixelized area (crop)

PRIOR irrigated

PRIOR subsistence

PRIOR rainfed_high

PRIOR rainfed_low

2. The SPAM process

42 crops simultaneously

…. continue

Page 7: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

…. continue2. The SPAM process

Behind the scene adjustments/calculations

Sub-national area (crops)

irrigated_ land

crop_land

suitable_area (crop)

FIXED!!

Sum >=

pix = max(ag,irr)

pix = f(suit,ag,irr)

Pre-process

irrigatedsubsistence

rainfed-highrainfed-low

PRIORS

Sum >=

Sum >=

Page 8: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

3. SPAM2005 Crop List and ...1 wheat 14 dry beans 28 sugarcane2 rice 15 chickpea 29 sugar beet3 maize 16 cowpea4 barley 17 pigeon pea 30 cotton5 pearl millet 18 lentils 31 other fibres6 finger/small millets 19 other pulses 32 arabica coffee7 sorghum 33 robusta coffee8 other cereals 20 soybeans 34 cocoa

21 groundnuts 35 tea9 potato 22 coconuts 36 tobacco

10 sweet potato 23 oil palm 37 banana11 yam 24 sunflower 38 plantains12 cassava 25 rapeseed 39 tropical fruit13 other roots & tub. 26 sesame seed 40 temperate fruit

27 other oil crops 41 vegetables42 all the rest

Page 9: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

3. and ... Statistical Area Coverage – subnational 1

Page 10: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

3. and ... Statistical Area Coverage – subnational 2

Page 11: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

4. Visualization of Downscalingexample: Ghana harvested area rice per region

Page 12: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

layer 1: Northern Region - harvested area rice 4. Visualization of Downscaling ... continued

Page 13: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

layer 2: Northern Region – irrigated area

irrigated area

harvested area rice

4. Visualization of Downscaling ... continued

Page 14: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

layer 3: Northern Region – agricultural area

cropland area

harvested area rice

irrigated area

4. Visualization of Downscaling ... continued

Page 15: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

layer 4: Northern Region – rural population

rural population

harvested area rice

irrigated area

cropland area

4. Visualization of Downscaling ... continued

Page 16: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

layer 5: Northern Region – suitable irrigated area rice

suitable irrigated area rice

harvested area rice

irrigated area

cropland area

rural population

4. Visualization of Downscaling ... continued

Page 17: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

layer 6: Northern Region – suitable rainfed high area rice

suitable rainfed high area rice

harvested area riceirrigated area

cropland area

rural population

suitable irrigated area rice

4. Visualization of Downscaling ... continued

Page 18: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

layer 7: Northern Region – suitable rainfed low area rice

suitable rainfed low area rice

harvested area riceirrigated area

cropland area

rural population

suitable irrigated area rice

suitable rainfed high area rice

4. Visualization of Downscaling ... continued

Page 19: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

some more data

result 1: Northern Region – allocated irrigated area rice4. Visualization of Downscaling ... continued

Page 20: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

some more data

result 2: Northern Region – allocated rainfed high area rice4. Visualization of Downscaling ... continued

Page 21: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

some more data

result 3: Northern Region – allocated rainfed low area rice4. Visualization of Downscaling ... continued

Page 22: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

some more data

result 4: Northern Region – allocated subsistence area rice4. Visualization of Downscaling ... continued

Page 23: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

some more data

result 5: Northern Region – allocated rice area total4. Visualization of Downscaling ... continued

scaled to FAO Country Totals

Page 24: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

5. SPAM ResultsCrop Distribution – Vegetables 2005

Page 25: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

5. SPAM ResultsCrop Distribution – Maize 2005

Page 26: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

6. Equations behind SPAM Minimize difference between prior and allocated area share for all pixels, crops and production systems in a cross entropy equation

subject to constraints (limits) dictated by existing• cropland area• irrigated area• suitable area• crop area statistics

and solve in GAMS

But first calculate the priors for each pixel, crop and productionsystem as a function of potential revenue, irrigated area and agricultural land,

where potential revenue is a function of percentage of crop area, crop price, rural population density and potential yield

Page 27: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

7. SPAM on the Web: MapSpam.info

Page 28: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

8. Challenges & …

• Different sources -> ‘contradictory’ information• Raster data not at same scale• Sub-national data complete, at least level1, better level2• Conform national crops -> FAO/SPAM crops• Consistencies between layers – constraints met

crop_land >= stats, irr >= crop_irr, suit_land >= crop_land >= stats

• Cropping intensities & production systems shares not consistent with data and model

• Validation of results (lack of established validation data, methods, and protocols)

Page 29: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

8. … and Issues

• Unreliable, often carelessly processed/validated statistics• Statistics do not match admin. area shape files • Unreliable cropland extent/area intensity estimates• Lack of data on cropping patterns and systems (e.g. cropping

intensity – converting harvested to physical land footprint)• Unsatisfactory data on location-specific biophysical conditions

(e.g., soils) and economic behaviour (e.g., prices and risks)• Lack of established validation data, methods, and protocols• Scientific peer review does not imply data are “fit-for-purpose”• Unspecified reliability of results• Consistency of approach has potential tradeoff with reliability

(e.g. patchwork of best national data vs consistent regional data)

Page 30: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

9. Validation

• Validation process by other CGIAR centers (e.g. IRRI, CIAT, ILRI, CIP, CYMMT). Each focuses on their mandate crops.

• Crop map view ‘parties’ attended by local experts and agronomists

• Crowd-sourcing on a dedicated website (MapSPAM.info)

Page 31: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

10. Users and Applications

• CGIAR centers such as IRRI, CYMMT, CIP, CIAT, ILRI. FAO, World Bank, and universities.

• HarvestChoice, Agricultural Water Management, AgFutures

• Fill the gaps between micro-macro linkage, between biophysical models and economic models

• Widely applied in country strategy work within IFPRI, regional priority settings such as ASARECA, CORAF, and in ReSAKSS, CAADP, AGRA.

Page 32: Spatial Production Allocation - Esri€¦ · 2015 Esri International User Conference -- Presentation Keywords: 2015 Esri International User Conference -- Presentation, EsriUC, Not

Thank you !