a Regional downscaling where three global climate models driven by the IPCC AR5 CMIP5 emission scenarios determine boundary conditions of the coupled ROMSNPZ high resolution oceanographic model for the Bering Sea, AK. b Biological downscaling of annual indices from the ROMSNPZ were used to drive thermal parameters in the CEATTLE model (i.e., weight-at-age and predation) as well as climate-enhanced spawner-recruitment relationships. c Annual harvest recommendations (ABC) from the assessment model which were translated into annual catch using the ATTACH social-economic model of the effect of EBFM policies on harvest

Scenario 1: Fixed exploitation rate. Exploitation rates = 0.05-0.4 in increments of 0.05. Floor assessed at the complex level. If the floor* of any complex is breached any further catch of this complex is considered a discard.

Scenario 2: The same as Scenario 1 with one addition. Floors assessed as the species level. If the floor* of any species is breached any further catch of this species is considered a discard.

Scenario 3: Variable exploitation rate. Starting exploitation rates = 0.05-0.4 in increments of 0.05. Each species complex is associated with a single fishing fleet. Species are still caught by multiple fleets, but when a species/complex becomes depleted, management actions only occur on the associated fleet. Each species complex is associated with the fleet that is considered the largest exploiter of the complex. The fishing fleets impact (exploitation) on the complex is adjusted through time5 (as depicted in HCR figure) when the complex biomass falls below 40% unfished biomass.

Scenario 4: The same as Scenario 3 with one addition. The fishing fleets impact on the complex is adjusted through time (as depicted in Figure 4) as the biomass of an individual species within the complex falls below 40% unfished biomass (50% for Elasmobranchs). This can be thought of as additional species protection