Optimizing Dam Removals in Lake Erie: An Interactive Exploration

Based on: Zheng, P. Q., B. F. Hobbs, and J. F. Koonce (2009), Optimizing multiple dam removals under multiple objectives: Linking tributary habitat and the Lake Erie ecosystem, Water Resour. Res., 45, W12417, doi:10.1029/2008WR007589. [Link to Paper]

Lake Erie, the shallowest and most biologically productive of the Great Lakes, faces numerous environmental challenges. One significant issue is the impact of thousands of dams built on its tributary rivers. Many of these dams are old, serve no current economic purpose, and block fish migration routes, harming native species like the popular walleye.

Removing these dams sounds like a straightforward solution, right? Restore the rivers, help the fish! However, it's not that simple. Dam removal is expensive, and paradoxically, it can also help invasive species like the destructive sea lamprey by giving them access to new spawning grounds. Furthermore, removing one dam can affect the situation at others upstream or downstream.

So, managers face a complex puzzle: With a limited budget, which dams should be removed to achieve the *best overall outcome* for the Lake Erie ecosystem, balancing the benefits for desirable fish like walleye against the costs and the risks of helping invasive species like sea lamprey?

This article explores the approach taken by Zheng, Hobbs, and Koonce (2009) to tackle this problem using multi-objective optimization. We'll break down their model and use interactive visualizations to understand the core concepts.

The Balancing Act: Benefits, Costs, and Trade-offs

At its heart, the decision involves trade-offs:

Imagine you have several dams you *could* remove. Each has a different cost, a different potential benefit for walleye, and a different risk for lamprey. How do you choose?

Quantifying the Factors

To make informed decisions, we need ways to measure these benefits and costs.

1. Measuring Ecosystem Health: The MVA Score

The researchers didn't just look at walleye. They considered the overall health of the Lake Erie ecosystem using multiple criteria, combined into a single score using Multi-Criteria Value Analysis (MVA). They identified key indicators representing ecological structure, function, and societal value (like fishing).

These indicators included:

Each indicator ($x_i$) is assigned a weight ($W_i$) reflecting its importance, and its current level is converted to a value ($V_i(x_i)$) on a 0-1 scale (0=worst, 1=best). The overall ecosystem health score ($V$) is the weighted sum:

$$ V(x_1, ..., x_8) = \sum_{i=1}^{8} W_i V_i(x_i) $$

Interactive MVA Score Calculator (Simplified)

Let's simulate how this works. Imagine we have just 4 criteria. Adjust the sliders to see how the levels of different criteria affect the overall score. (Weights are pre-set based roughly on the paper's revised weights, normalized). Assume linear value functions $V_i(x_i)$ for simplicity, mapping raw values to a 0-100 scale.

50
50
50
50
Overall Ecosystem Health Score: 50.0 / 100

The key challenge is that removing dams changes these underlying criteria ($x_i$). The researchers used a complex simulation model of Lake Erie (LEEM - Lake Erie Ecological Model) to predict how changes in tributary habitat (specifically, walleye young-of-year survival) would ripple through the food web and affect these lake-wide indicators.

2. Estimating Costs

The model also needs cost estimates:

The total cost ($x_9$) for a portfolio of removed dams is the sum of removal costs and necessary lamprey control costs for all selected dams.

The Optimization Puzzle: Selecting the Best Portfolio

Now we have the pieces: a way to measure overall ecosystem health (the MVA score) and ways to estimate the costs. The goal is to select a set of dams for removal that maximizes the ecosystem health score, without exceeding a given budget ($B$).

This is formulated as a Mixed Integer Linear Program (MILP). The core decision is binary for each potential dam $j$: either remove it ($d_j = 1$) or keep it ($d_j = 0$).

The objective is to Maximize $z_1$ (the MVA score):

$$ \text{Maximize } z_1 = \sum_{i=1}^{8} W_i V_i(x_i) $$

Subject to:

  1. Budget Constraint: The total cost ($x_9$) cannot exceed the budget ($B$).
  2. $$ x_9 = \sum_{j \in J} (C_j^{dam} + C_j^{lamprey}) d_j \leq B $$
  3. Ecosystem Linkage: The ecosystem criteria ($x_i$) depend on the baseline level ($X_i^{base}$) plus the cumulative effect of removed dams on walleye recruitment ($\Delta YOY^{walleye}$), modeled via impulse responses ($IR_{x_i}$).
  4. $$ x_i = X_i^{base} + IR_{x_i} \Delta YOY^{walleye} \quad \text{where} \quad \Delta YOY^{walleye} = \sum_{j \in J} (\Delta YOY_j^{walleye}) d_j $$
  5. Habitat Logic: A dam ($n$) cannot be removed if a dam ($j$) immediately downstream of it is *not* removed (ensures connectivity).
  6. $$ d_n \leq d_j \quad \text{for all } n \text{ directly upstream of } j $$
  7. Binary Decision: Each $d_j$ must be 0 or 1.
  8. $$ d_j \in \{0, 1\} \quad \text{for all } j \in J $$

Interactive Dam Portfolio Selector (Simplified Toy Example)

Let's try solving a mini-version of this puzzle. Below are 5 hypothetical dams. Each has a removal cost, an estimated positive impact on the overall MVA score (Ecosystem Benefit), and an estimated negative impact via lamprey cost/effects (Lamprey Penalty, also folded into the MVA score calculation implicitly). Your budget is shown below. Click dams to select them for removal. Can you maximize the Net Ecosystem Benefit without exceeding the budget?

Selected Dams: 0
Total Cost: $0 M
Total Ecosystem Benefit: 0 points
Status: Within Budget

Note: This is highly simplified. The real model considers complex ecosystem interactions (LEEM), specific habitat suitability (HSI), upstream/downstream logic, and many more dams (~139 candidates in the paper).

Seeing the Trade-offs: The Efficient Frontier

Solving the MILP for a single budget gives one optimal portfolio. But what if the budget changes? Or what if we want to see the *range* of possibilities?

By running the optimization model many times with different budget levels ($B$), the researchers traced out an efficient frontier (also called a Pareto frontier). This curve shows the maximum possible ecosystem health score achievable for any given level of spending.

Points on the curve represent "efficient" portfolios – you can't improve the ecosystem score further without spending more money, and you can't spend less money without accepting a lower ecosystem score.

Exploring the Efficient Frontier

The chart below shows a hypothetical efficient frontier based on the paper's findings (Figure 4a). Use the slider to set a budget. The corresponding point on the frontier shows the best ecosystem health achievable for that budget, and the tooltip (hover over points) indicates how many dams might be removed in that optimal portfolio.

For Budget $15 M: Max Ecosystem Health ~68% (relative to max possible), ~31 dams removed.

Observe the "diminishing returns": initial spending yields large gains in ecosystem health, but achieving the highest levels requires much larger budget increases for smaller improvements.

What if Priorities Change? Sensitivity Analysis

The MVA score depends on the weights ($W_i$) assigned to each criterion. What if stakeholders are particularly concerned about the sea lamprey invasion risk? The researchers explored this by adding an objective to *minimize* potential lamprey habitat, giving this objective a weight ($WL$).

Maximize: $(1 - WL) \times (\text{Ecosystem Health } z_1) - WL \times (\text{Lamprey Habitat } z_3)$

Increasing the weight ($WL$) on minimizing lamprey habitat fundamentally changes which dams are selected in the optimal portfolio, even for the same budget.

Impact of Lamprey Aversion Weight (WL)

Let's simulate how changing the priority (weight $WL$) placed on avoiding lamprey habitat affects the outcome, keeping the budget fixed at $15M (based roughly on Table 3).

Resulting Portfolio (at $15M Budget):
- Dams Removed: 31
- Ecosystem Health (z1 %): 68.3%
- Lamprey Habitat Created (z3 % of potential): 49.1%

As you increase $WL$, the model prioritizes removing dams that open up less (or zero) lamprey habitat. This leads to removing *more* dams (often smaller, cheaper ones with less lamprey risk), but the overall ecosystem health score ($z_1$, which heavily weights walleye) decreases significantly.

Conclusion

The work by Zheng, Hobbs, and Koonce provides a powerful framework for tackling the complex problem of multiple dam removals. Key takeaways include:

While powerful, these models are simplifications. They rely on data and sub-models (like LEEM and cost regressions) that have uncertainties. The paper emphasizes that the MILP should be used as a screening tool to identify promising portfolios of dam removals. Further site-specific studies are crucial before making final decisions on removing any particular dam.

This interactive exploration aimed to provide intuition for the core concepts. By making the trade-offs tangible and allowing exploration of parameters, we can better appreciate how mathematical modeling helps navigate complex environmental management decisions.