The scoring methodology translates merged species distribution models and extinction risk weights into cell-level sensitivity metrics, which are then rescaled and aggregated to management zones for decision support.
Cell-Level Scoring
For each grid cell, sensitivity is computed as the extinction risk weighted sum of species presence within each species category (Equation 7.1):
\[
\text{score}_{c,g} = \sum_{s=1}^{S_g} \frac{\text{er\_score}_s \times v_{s,c}}{100}
\tag{7.1}\]
where:
- \(\text{score}_{c,g}\) = sensitivity score for cell \(c\) and species category \(g\)
- \(\text{er\_score}_s\) = extinction risk score for species \(s\) (1–100; see Chapter 5)
- \(v_{s,c}\) = merged model value for species \(s\) in cell \(c\) (0–100%; see Chapter 6)
- \(S_g\) = number of valid species in category \(g\)
Division by 100 converts the product from a 0–10,000 range back to a 0–100 scale.
In plain terms, for each cell in the ocean the MST adds up the sensitivity contributions of all species found there. If a cell has many species that are both likely to be present and at high risk of extinction, it gets a higher sensitivity score. This helps identify places where rare or threatened species are concentrated.
Metric Keys
Cell-level scores are stored as named metrics in the DuckDB database:
Ecoregional Rescaling
Raw cell scores vary naturally across ecoregions due to differences in species richness, oceanographic conditions, and biogeographic patterns. To enable meaningful cross-region comparison, each metric is rescaled to a [0–100] range within each BOEM ecoregion:
\[
\text{score}'_{c} = \frac{\text{score}_{c} - \text{score}_{min}}{\text{score}_{max} - \text{score}_{min}} \times 100
\tag{7.2}\]
where \(\text{score}_{min}\) and \(\text{score}_{max}\) are the minimum and maximum cell values for that metric within the ecoregion. This ensures that a score of 100 represents the most sensitive cell within each ecoregion, not globally.
Primary Productivity Scoring
Net primary productivity (NPP) from the VGPM satellite model is included as a separate metric (primprod) that does not undergo extinction risk weighting. Instead, the raw NPP values (metric tons C km-2 yr-1) are directly rescaled within each ecoregion using the same [0–100] normalization.
Zone Aggregation
Rescaled cell scores are aggregated to three spatial levels using area-weighted averaging:
- Program Areas: the 20 BOEM planning regions used for offshore energy management
- Subregions: intermediate spatial groupings (e.g., Atlantic, Gulf of Mexico, Pacific, Alaska)
- Ecoregions: BOEM-defined ecological regions used as the rescaling baseline
For each zone, the aggregated score is:
\[
\text{score}_{zone} = \frac{\sum_{c \in zone} \text{score}'_c \times A_c}{\sum_{c \in zone} A_c}
\tag{7.3}\]
where \(A_c\) is the area of cell \(c\) (accounting for latitude-dependent cell size in the 0.05° grid).
Visualization
Flower Plot
The flower plot provides an intuitive summary of sensitivity scores by species category, inspired by the Ocean Health Index (Halpern et al. 2012) visualization approach:
- Petal length represents the rescaled sensitivity score (0–100) for each species category — longer petals indicate higher sensitivity
- Petal width represents the weight of each category contributing to the overall score — wider petals have more influence on the center score
- Center value shows the weighted mean across all categories
Figure 7.1 shows flower plots for all 20 BOEM Program Areas, generated from the v3 scoring run. Each flower summarizes the sensitivity profile of one program area with petals for Bird, Coral, Fish, Invertebrate, Mammal, Other, Turtle, and Primary Production. Alaska program areas (e.g., HOP, COK, GOA) show dominant Fish and Mammal petals, while Gulf areas (GAB, CEC, MAT) show stronger Invertebrate and Coral contributions.
The flower plot lets decision-makers quickly spot which ecological elements drive sensitivity in an area, aiding better planning and impact assessment.
The flower plot is implemented as an interactive ggiraph visualization using polar coordinates in R (ggplot2 + coord_polar()), with tooltips showing detailed score breakdowns on hover. The plot_flower() function in calc_scores.qmd generates static versions for reports, while the mapgl web application implements an interactive JavaScript version for responsive rendering (see Figure A.2).
Treemap
The treemap provides a detailed, hierarchical view of species contributions to sensitivity scores (Figure 7.2):
- Audience: scientists and analysts
- Shows the contribution of each species to the overall score within a pixel or management zone
- Elements are hierarchical: first by extinction risk category, then by taxonomic classification
- Interactive zooming allows drilling down to individual species contributions
Zooming into a specific category reveals individual species contributions:
Interactive Maps
The primary visualization tools are two Shiny web applications:
- mapgl: general sensitivity mapping app showing composite scores by program area with flower plots and treemaps (see Appendix: mapgl)
- mapsp: species distribution viewer showing individual and merged models with regulatory status information (see Appendix: mapsp)
Halpern, Benjamin S., Catherine Longo, Darren Hardy, Karen L. McLeod, Jameal F. Samhouri, Steven K. Katona, Kristin Kleisner, et al. 2012.
“An Index to Assess the Health and Benefits of the Global Ocean.” Nature.
https://doi.org/10.1038/nature11397.