220 likes | 455 Views
Dynamic Resource Allocation in Conservation Planning. Daniel Golovin. Andreas Krause. Beth Gardner Sarah Converse Steve Morey. California Institute of Technology Center for the Mathematics of Information. 1. Ecological Reserve Design. How should we select land for conservation
E N D
Dynamic Resource Allocation in Conservation Planning Daniel Golovin Andreas Krause Beth Gardner Sarah Converse Steve Morey California Institute of Technology Center for the Mathematics of Information 1
Ecological Reserve Design How should we select land for conservation to protect rare & endangered species? Case Study: Planned Reserve in Washington State Mazamapocket gopher streaked horned lark Taylor’s checkerspot
Problem Ingredients • Land parcel details • About 5,300 parcels • soil types, vegetation, slope • conservation cost
Problem Ingredients • Land parcel details • Geography: Roads, Rivers, etc
Problem Ingredients • Land parcel details • Geography: Roads, Rivers, etc • Model of Species’ Population Dynamics • Reproduction, Colonization, Predation, Disease, Famine, Harsh Weather, …
Population Dynamics Modeled using a Dynamic Bayesian Network Our Choices Protected Parcels Time t Time t+1 Environmental Conditions (Markovian)
Population Dynamics Modeled using a Dynamic Bayesian Network Our Choices Protected Parcels Time t Time t+1 Environmental Conditions (Markovian)
Model Paramters • From the ecology literature, or • Elicited from panels of domain experts Annual Patch Survival Probability Patch Size (Acres)
From Parcels to Patches Patch 1 Patch 2 Most parcels are too small to sustain a gopher family So we group parcels into larger patches. We assume no colonization between patches, and model only colonization within patches. We optimize over (sets of) patches.
The Objective Function Pr[alive after 50yrs] 0.8 • In practice, use sample average approximation 0.7 0.5 f(R)= 2.0 (Expected # alive) Selected patches R • Choose R to maximize species persistence
“Static” Conservation Planning • Select a reserve of maximum utility, subject to budget constraint NP-hard But f is submodular We can find a near-optimal solution
Structure in Reserve Design A B • Diminishing returns: helps more in case A than in case B • Utility function f is submodular:
Solving the “Static” Conservation Planning Problem • More efficient algorithm with slightly weaker guarantees [Leskovec et al. ‘07] Theorem [Sviridenko‘04]: We can efficiently obtain reserve R such that
Selected patches are very diverse
Results: “Static” Planning • Can get large gain through optimization
Dynamic Conservation Planning • Build up reserve over time • At each time step t, the budget Bt and the setVtof available parcels may change • Need to dynamically allocate budget tomaximize value of final reserve Time t+1 Time t
Opportunistic Allocation forDynamic Conservation In each time step: • Available parcels and budget appear • Opportunistically choose near-optimal allocation Time t=2 Time t=1 Theorem: We get at least 38.7% of the value of the best clairvoyant algorithm* * Even under adversarial selection of available parcels & budgets.
Results: Dynamic Planning • Large gain from optimization & dynamic selection
Dynamic Planning w/Failures • Parcel selection may fail • Purchase recommendations unsuccessful • Patches may turn out to be uninhabitable • Can adaptively replan, based on observations • Opportunistic allocation still near-optimal • Proof uses adaptive submodularity[Golovin & Krause ‘10]
Dynamic Planning w/Failures • Failures increase the benefit of adaptivity 50% failure rate
Related Work • Existing software • Marxan [Ball, Possingham & Watts ‘09] • Zonation [Moilanen and Kujala‘08] • General purpose software • No population dynamics modeling, no guarantees • Sheldon et al. ‘10 • Models non-submodular population dynamics • Only considers static problem • Relies on mixed integer programming
Conclusions • Reserve design: prototypical optimization problem in CompSustAI • Large scale, partial observability, uncertainty, long-term planning, … • Exploit structure near-optimal solutions • General competitiveness result about opportunistic allocation with submodularity