ohi logo
OHI Science | Citation policy 

## This file makes it easier to process data for the OHI global assessment
##  by creating the following objects:
## 
##  * dir_M = identifies correct file path to Mazu (internal server) based on your operating system
##  * mollCRS = the crs code for the mollweide coordinate reference system we use in the global assessment
##  * regions_shape() = function to load global shapefile for land/eez/high seas/antarctica regions
##  * ohi_rasters() = function to load two rasters: global eez regions and ocean region
##  * region_data() = function to load 2 dataframes describing global regions 
##  * rgn_syns() = function to load dataframe of region synonyms (used to convert country names to OHI regions)
##  * low_pop() = function to load dataframe of regions with low and no human population
##  * UNgeorgn = function to load dataframe of UN geopolitical designations used to gapfill missing data
## loads 2 dataframes: rgns_all and rgns_eez 
##  rgns_all = includes eez/high seas/antarctica regions, IDs correspond with region shapefile and raster
##  rgns_eez = includes only eez regions

1 Source

We use the methods and code from Casey O’Hara: https://github.com/oharac/spp_risk_dists to prepare the data for the species subgoal.

Reference:

2 Summary

Spatial species range data and extinction risk data from IUCN is used to generate regional scores for the Species subgoal (part of the Biodiversity goal) and resilience layers.

Mean risk status per cell: Species ranges are converted to a global spatial raster of 10 km resolution.

The mean extinction risk for each cell, \(\bar{R}_{cell}\), is calculated by averaging the IUCN extinction risk of the species with ranges overlapping the cell.

Risk is a scaled value representing the species extinction risk category: * ‘LC’ = 0.0, ‘NT’ = 0.2, ‘VU’ = 0.4, ‘EN’ = 0.6, ‘CR’ = 0.8, ‘EX’ = 1.0

\[\bar{R}_{cell} = \frac{\displaystyle\sum_{species}(Risk)}{n_{spp}}\]

Mean risk status per region: The mean extinction risk for a region, \(\bar{R}_{SPP}\), is estimated by averaging the risk values of the raster cells falling within each OHI region, with each cell’s contribution weighted by the number of species in the cell.

Species goal model

The regional risk values are converted to species status scores by subtracting the risk values from 1 and rescaling so a risk value of $0.75 $ receives a score of zero.

From Halpern et al (2012):

The target for the Species sub-goal is to have all species at a risk status of Least Concern. We scaled the lower end of the biodiversity goal to be 0 when 75% species are extinct, a level comparable to the five documented mass extinctions and would constitute a catastrophic loss of biodiversity.

\[X_{SPP} = \frac{((1 - \bar{R}_{SPP}) - 0.25)}{(1 - 0.25)} * 100%\]

where:

  • \(X_{SPP}\) is Species goal status
  • \(\bar{R}_{SPP}\) is mean extinction risk for a region

Trend is calculated using time series of risk categories based on current and past assessments.

3 Updates from previous assessment

Additional year of data

4 Data Sources

  • Reference:

  • Description: Shapefiles containing polygons of assessed species ranges; each shapefile represents all assessed species within a comprehensively-assessed (i.e. >90% assessed) taxonomic group.

  • Native data resolution: NA

  • Time range: NA

  • Format: Shapefile

  • Reference:
    *IUCN and UNEP-WCMC (2020), The World Database on Protected Areas (WDPA) [On-line], April 2020. Cambridge, UK: UNEP-WCMC.

    • Zipped shapefile available at: www.protectedplanet.net.
    • Downloaded: April 24, 2020

  • Description: Shapefile of World Database on Protected Areas

  • Time range: 1819 - 2020; some protected areas do not have an associated “status year” and are reported as year 0.

  • Format: Shapefile

*File location: Mazu:git-annex/globalprep/_raw_data/wdpa_mpa/d2020/WDPA_Apr2020-shapefile/

5 Methods

There are several steps that need to be taken to get to this point in the data prep.

Here is an overview of the organization of files and data that are run prior to this:

5.0.1 Code

5.0.1.1 Run this first! Setup directory: spp_risk_dists/_setup

In this directory are a sequence of files used to generate the bits and pieces that are later assembled into the rasters of biodiversity risk.

.Rmd files are sequenced with a prefix number (and letter) to indicate the order of operations. Briefly:

  1. Pull information from the IUCN Red List API to determine an overall species list, habitat information, and current risk (conservation status).
    • 1_set_up_iucn_habs_and_risk.Rmd
  2. Pull information from API on risk from regional assessments; also recode the regions according to Marine Ecoregions (Spalding et al, 2007) for later spatialization.
    • 2_set_up_iucn_risk_regional.Rmd
    • 2a_set_up_regional_codes.Rmd
  3. Pull historical assessment information from API for possible trend analysis.
  • 3_set_up_iucn_trends.Rmd
  1. Set up spatial layers in Gall-Peters, 100 km2 cells. Layers include:
    • cell ID (cells are sequentially numbered for combining with tabular data)
    • ocean area
    • marine protected area (classification, year of protection, proportion of protection)
    • Exclusive Economic Zones (EEZ) and FAO fishing regions
    • Marine Ecoregions of the World
    • bathymetry
    • NOTE: these layers are all saved in the spp_risk_dists/_spatial directory.
      • 4a_set_up_ocean_area_and_mpa_pct.Rmd
      • 4b_set_up_eez_meow_rasts.Rmd
  2. Convert species range maps to rasters. (skip 5a .Rmd because all that is doing is checking the non-comprehensive data, which is not used for OHI and was not downloaded for the subgoal.)
    • For maps provided directly by IUCN, aggregate into multispecies files based on family. There is some cleaning done at this stage to fix problematic extents and attributes.
    • From the list of all available maps, generate a master list of all mapped, assessed species for inclusion in the study.
    • Rasterize each species to a .csv that includes cell ID and presence. A .csv format was used for file size and ease of reading and binding into dataframes.
      • 5b_generate_spp_map_list.Rmd
      • 5c_rasterize_spp_shps.Rmd
  3. Aggregate individual species ranges into larger taxonomic groups, and summarize key variables (mean risk, variance of risk, number of species, etc) by group.
    • Technically this is not necessary but makes it easier to quality check the process along the way, and supports mapping at the level of taxonomic group rather than the entire species list level.
    • This process is done twice: once for uniform weighting and once for range-rarity weighting. Resulting files are saved separately.
      • 6a_aggregate_spp_ranges.Rmd
      • 6b_aggregate_rr_spp_ranges.Rmd

5.0.1.2 Then run this! Root directory: spp_risk_dists

At this level there are several scripts, prefixed 1x_biodiversity_maps, that collate the various taxonomic group level files (generated in setup part 6) and summarize to the global level. These need to be run before spp_data_prep.Rmd! * 1a_biodiversity_maps_comp_assessed.Rmd * 1c_biodiversity_maps_all_spp.Rmd * spp_data_prep.Rmd

  • Note each creates a specific aggregation - comprehensively assessed species vs all available species; uniform vs range-rarity weighting.
  • The rasters generated in these scripts are saved in the _output folder.

5.0.2 Data and output files

The spp_risk_dists/_data folder contains tabular data about IUCN species used throughout the processing of this analysis. These files are generated by scripts in the setup directory.

The spp_risk_dists/_spatial folder contains general spatial data generated and/or used in the setup scripts. These include:

  • rasters for cell ID, EEZ ID, marine ecoregion ID, ocean area, and bathymetry masks.
  • tabular data of region names and lookups for IUCN regional assessment to marine ecoregion.
  • tabular data of marine protected area level/year/coverage to cell ID.
  • shapefiles used for map plotting from Natural Earth.

The spp_risk_dists/_output folder contains the rasters of biodiversity risk, species richness, variance of risk, etc generated from the scripts in the base directory.

5.1 Compare all vs. comprehensively assessed species

In the past, we have used all species with IUCN risk assessments to calculate the species subgoal. However, some of Casey’s work suggests it is better to use the taxa groups that have been comprehensively assessed by IUCN (> 90% of species assessed). The general concern is that IUCN tends to oversample species in the Atlantic portion of the ocean, relative to other regions. This is indicated by the larger number of species with IUCN status in this region. However, the Atlantic falls in line with the other regions when looking at the comprehensively assessed species.

5.2 SPP: Status

5.2.2 Status: estimate for previous years

We use the trend data to estimate risk values for previous years (vs. using the same values for all assessment years). The change in species status across years is based on a linear model.

Trend is calculated using the same method as the risk calculation. For each cell, we multiply the average species trend by the number of species in order to weight each cell’s contribution by the number of species in the cell. We sum these values for each OHI region and calculate for each region: (average species trend * number species)/number of species

5.2.4 Status: Converting regional mean risk to status

We rescale the data so a risk factor of 0.75 is equal to zero.

5.3 SPP: Trend

5.3.2 Trend: Gapfilling missing data

Check there are data for every region. Region 232 (Bosnia) does not have a value which is not surprising because their coast is super small and results are erratic for this region. We estimate this using the mean of the 2 surrounding regions.

5.4 Trend: Final formatting for ohi-global

5.4.1 Trend: Compare to last year

We use a very different approach for calculating trend this year. Previously we used proxy data rather than actual change in IUCN status over time. Our previous method overestimated the magnitude of the trend. It is not surprising there is poor correlation with trend estimates in previous years, but it is reassuring the the values mainly fall in the same quadrant.

5.5 Resilience data

We use species condition data as a resilience measure as well. We also calculate species condition at 3nm of shoreline, because for some goals, nearshore species condition is the relevant metric.

5.5.3 Resilience: estimate for previous years

We use the trend data to estimate risk values for previous years (vs. using the same values for all assessment years). The change in species status across years is based on a linear model.

Trend is calculated using the same method as the risk calculation. For each cell, we multiply the average species trend by the number of species in order to weight each cell’s contribution by the number of species in the cell. We sum these values for each OHI region and calculate for each region: (average species trend * number species)/number of species

5.5.6 Resilience: Converting regional mean risk to status

We rescale the data so a risk factor of 0.75 is equal to zero.