Thermal Habitat Index of Many Northwest Atlantic Temperate Species Stays Neutral under Warming Projected for 2030 but Changes Radically by 2060

Abstract

Global scale forecasts of range shifts in response to global warming have provided vital insight into predicted species redistribution. We build on that insight by examining whether local warming will affect habitat on spatiotemporal scales relevant to regional agencies. We used generalized additive models to quantify the realized habitat of 46 temperate/boreal marine species using 41+ years of survey data from 35°N–48°N in the Northwest Atlantic. We then estimated change in a “realized thermal habitat index” under short-term (2030) and long-term (2060) warming scenarios. Under the 2030 scenario, ∼10% of species will lose realized thermal habitat at the national scale (USA and Canada) but planktivores are expected to lose significantly in both countries which may result in indirect changes in their predators’ distribution. In contrast, by 2060 in Canada, the realized habitat of 76% of species will change (55% will lose, 21% will gain) while in the USA, the realized habitat of 85% of species will change (65% will lose, 20% will gain). If all else were held constant, the ecosystem is projected to change radically based on thermal habitat alone. The magnitude of the 2060 warming projection (∼1.5–3°C) was observed in 2012 affirming that research is needed on effects of extreme “weather” in addition to increasing mean temperature. Our approach can be used to aggregate at smaller spatial scales where temperate/boreal species are hypothesized to have a greater loss at ∼40°N. The uncertainty associated with climate change forecasts is large, yet resource management agencies still have to address climate change. How? Since many fishery agencies do not plan beyond 5 years, a logical way forward is to incorporate a “realized thermal habitat index” into the stock assessment process. Over time, decisions would be influenced by the amount of suitable thermal habitat, in concert with gradual or extreme warming.

About the Authors:

Nancy L. Shackell

* E-mail: [email protected]

Affiliation: Oceans and Ecosystem Science Division, Department of Fisheries and Oceans, Bedford Institute of Oceanography, Dartmouth, Nova Scotia, Canada

Daniel Ricard

Affiliation: Biology Centre AS CR v.v.i., Institute of Hydrobiology, České Budějovice, Czech Republic

Christine Stortini

Affiliation: School for Resource and Environmental Studies, Dalhousie University Halifax, Nova Scotia, Canada

Introduction

Global scale meta-analyses of biological responses to climate change have provided vital insight into changes in distribution, phenology and species interactions over the last 2–3 decades [1]. Such information is necessary for policy makers working on general climate change adaptation strategies, but not sufficient for marine regulatory agencies where commercial species catches are highly regulated at smaller spatial and temporal scales. At regional and local scales (∼10–100 km2), the uncertainty in climate change projections increases [2] yet marine populations are highly responsive to regional temperature variation. It is indeed a challenge to provide regional scale information on how and when marine species will respond to climate change so that ocean managers can modify governance structures accordingly [3].

The biological response to climate variation can be complex but a population’s first acclimative response is typically a shift in spatial distribution in response to temperature change [1], [4]. In the ocean, global model projections predict that fish will gradually migrate poleward from 2005–2050 [5]. Recently, it has been shown that change is not necessarily poleward; there can be a variety of directional responses to local changes in temperature [6]. Range shifts on regional scales are already evident in the North Sea [7] and off the northeastern United States [8]. For a few commercial species on the US Northeastern seaboard, fisheries management strategies have lagged behind their northerly shifts [9]. In other areas the rate of temperature increase is slower [10] or well within the range of natural variability [11]. It has long been understood that temperature is a primary determinant of species distribution and we need to understand how species habitat availability will be affected by warming trends.

The Species Distribution Model (SDM) has been used in many studies to explore and forecast range shifts in response to warming. The SDM is based on the concept of an ecological “niche” where distribution is first and foremost determined by physiological constraints imposed by environmental variables (i.e. the “Fundamental” niche, reviewed in [12]). But an animal’s distribution is not just determined by environment; it is further constrained by species interactions such as competition and predation. This further constrained distribution has been referred to as the “Realized” niche. As such, one of the greatest weaknesses of most SDMs is that predictions do not account for species interactions, micro-evolutionary changes or dispersal abilities [13]. Most SDMs use a correlative approach whereas the much more difficult mechanistic approach can provide more insight and improve predictability [14], [15]. The uncertainty of correlative SDMs to predict species distribution is exacerbated by the uncertainty in climate change projections and the influence of other anthropogenic factors such as fishing and pollution [16]–[18]. Still, correlative SDMs can provide insight in appropriate situations and are constantly being modified to improve their application [18], [19].

In this study, we provide a broad regional-scale overview of changes in a realized thermal habitat index to projected short (2030) and long-term (2060) warming. We used a correlative SDM approach to quantify the realized habitat for each of 46 species using 41+ years of survey data from Cape Hatteras, North Carolina to Cape Breton Island, Nova Scotia, Canada (Figure 1). We estimated changes (and confidence limits) in the amount of each species’ realized thermal habitat index under short-term (2030) and long-term (2060) warming scenarios. Our goals are to provide an initial triage to gauge the severity of ecological/economic impacts and to flag vulnerable or expanding species in the USA and Canada. With such information, policy-makers and ocean managers could evaluate the risks of management action (or inaction) and use our approach to examine a subset of species at smaller spatial scales within the USA and Canada.

[Figure omitted. See PDF.]

Figure 1. Northeast Continental Shelf (United States) and Scotian Shelf (Canada) in the Northwest Atlantic.

The red line demarcates each country’s jurisdiction, separated by the Hague line. Most of the species studied are more common north of the Mid-Atlantic Bight.

https://doi.org/10.1371/journal.pone.0090662.g001

Materials and Methods

Data

Data used herein are publicly available from the Ocean Biogeographic Information System (OBIS) under the dataset names “DFO Maritimes Research Vessel Trawl Surveys Fish Observations” (OBIS-Canada) and the “Northeast Fisheries Science Center Bottom Trawl Survey Data” (OBIS-USA) [20], [21] (http://iobis.org/). Presence/absence data of 46 species (Table 1) from 35°N to 48°N from 1963–2012 were extracted (Figure 1). The species selected are common throughout the Gulf of Maine, Georges Bank and on the Scotian Shelf and form a temperate/boreal complex. We also obtained the environmental covariates of interest (date, depth, bottom temperature, latitude and longitude) from the OBIS datasets. Bottom temperature was recorded for each survey tow whether fish were present or not; these data are ideal for examining how local bottom temperature variation affects species distribution (Table S1).

[Figure omitted. See PDF.]

Table 1. Generalized Additive Model (GAM) result summaries for each species using Full dataset and 10 subsamples.

https://doi.org/10.1371/journal.pone.0090662.t001

Model Approach

There are a variety of SDM approaches to explore possible shifts in species distribution and an associated branch of research that has evaluated the assumptions and methodologies among approaches [22]–[24]. We used these evaluations to select the most appropriate approach, given our data. For this broad exercise of examining multiple species, we used a correlative approach [14]. We opted to use a binomial Generalized Additive Model (GAM) using a logit link function [25] as the SDM. Our decision was based on reviews that showed GAMs were superior to other SDM techniques under circumstances similar to ours [24], [26]. For each species, presence/absence data were modeled as a function of the following covariates: location (latitude, longitude), year, bottom temperature, and depth. The fitted values were estimates of the probability of occurrence at location for each species and reflect “realized habitat”. In preliminary analyses of several common species, we fit models using only bottom temperature in efforts to estimate “Potential” thermal habitat, but concluded that the model over-predicted occurrences for the purpose of estimating change under warming scenarios. Therefore we conservatively defined model fit as a “Realized Thermal Habitat Index” to reflect the correlative nature of our SDM [14]. All analyses were conducted using the R statistical language [27] using various packages including the “mgcv” package [28]. Further details are provided in the Supporting Information.

Model Evaluation

We evaluated each model’s accuracy using a method, similar in intent to the Wilcoxon Rank Sum test, called the “Area Under the Receiver Operating Characteristic Curve (AUC) available in R CRAN package “caTools”. This method evaluates how well the final model predicts a true positive and a true negative [29].

Regional Warming Scenarios

Global climate models (GCM), developed from general circulation models, or earth system models often have limited application at smaller spatial scales, yet organisms respond to smaller-scale variation. Methods to down-scale from large-scale GCM coarse resolution (grid size∼100–300 km2) to regional/local scale resolution (grid size:∼5–100 km2) are still evolving rapidly [2] and are not yet fully available for the Canadian portion of our study domain. Accuracy of SDMs predictions are confounded by uncertainties in GCM climate forecasts and uncertain “down-scaling” [22]. As described earlier, the field of predicting species distribution from SDMs is evolving as well [30]. All SDM climate change analyses suffer from these compounding sources of uncertainty. We used the most recent information available [11], [31] to create 2 likely scenarios in which to evaluate projected change in realized habitat indices.

The two scenarios are derived from regional syntheses of sea surface temperature (SST) projections [11] and trends [31] for the Scotian Shelf and Gulf of Maine. These trends and projections are for SST. Reliable bottom temperature projections are not yet available but we recognize that deeper waters warm more slowly than surface water. Therefore, the first scenario is derived from a long-term (50–60years) mid-range projection (3°C) from an ensemble of AR5 ESM August SST projections for the region [11]. Using each species dataset, we added 3°C to all bottom temperature data collected at depths equal to or less than 100 m, and 1.5°C to depths more than 100 m. The second scenario is based on SST empirical decadal trends in the region since 1985 [31] adjusted for 20 years (0.7°C added to bottom temperature at depths equal to or less than 100 m, and 0.35°C added to depths more than 100 m). We refer to Scenario 1 as Y2060 and Scenario 2 as Y2030.

Estimating Gain/Loss of Realized Thermal Habitat Index

We estimated area of habitat where probability of occurrence from the original model output was more than 0.54. Note that we could have selected any habitat >0.5 as our goal was to be conservative and estimate the change in the most likely realized thermal habitat index. We did the same for each of the projected scenarios’ model output (e.g. Y2060 and Y2030 realized thermal habitat index). Net change (gain/loss) was the percentage of current realized thermal habitat index less the future, divided by the current realized thermal habitat index.

The common method of sample-splitting data (using one set to calibrate a model and the other to test) does not address whether a model is suitable for extrapolation [24]. Further, Araújo et al. [32] argued that splitting the data into 2 sets does not provide an independent test as both sets are derived from the same source. However, we did wish to gauge the sample error in the gain/loss estimate. The range of estimates for each species is an index of uncertainty, which in turn discerns whether our method is appropriate for a given species. A large range would indicate that this method was not appropriate for the species, and other research avenues would have to be pursued. Therefore, for each species, we fit GAMs to 10 random sub-samples (0.6 of full dataset). This provided a range of all model output but more importantly a range of the estimated gain/loss as calculated based on probability estimates.

Results

Model Evaluation

Accuracy of predictions among models for 46 species, ranged from 0.77–0.99 (Median = 0.89) indicating strong accuracy for all species (Figure 2, Table 1). The power is derived from 41+ years of observations over a large geographic scale (∼35°N–48°N) and spanning a wide temperature range (∼–1.5°C to 25°C). While the models accurately reproduce the probability of occurrence for all species, it is evident that they are more accurate, and a higher amount of deviance is explained, for those species that are restricted by warmer temperatures throughout the study domain (Table 1). Snow crab, capelin, turbot, and northern shrimp would all be considered to be at the southern limit of their distribution within this study domain(∼41°N), indicating that the many zero values in the southern part of the domain (<41°N) contributed to increased accuracy.

[Figure omitted. See PDF.]

Figure 2. Area Under Receiving Operating Curve (AUC).

Expected Accuracy (of classification) of GAM model fit for each species, compared to the observations. An AUC value of 1 would indicate that predicted values resulted in a completely accurate classification of observations. For all but one species, the model was able to accurately classify >80% probability of occurrence.

https://doi.org/10.1371/journal.pone.0090662.g002

Between Countries and Scenarios

The median net change in the realized thermal habitat index is similar between countries and scenarios. Under the long-term (2060) warming scenario (Scenario 1), the majority of species in Canada and in the USA will lose realized thermal habitat (Figure 3). Under the short term (2030), Scenario 2, the median loss of both nations is neutral (within +/–10%). There are a few outliers indicating that there are extreme winners and losers, regardless of whether the majority stay relatively neutral (Figure 3). The majority of species are neutral as a result of the northern temperate/boreal nature of the species examined, and the scale of aggregation at the national level. It is important to note that there will be a more pronounced response for some species at smaller spatial scales. When the domain is divided into sub-regions, there can be a stronger negative response in the USA and on south-west Scotian shelf for species that are at the southern limit of their range (See Figure S1 for an example of cod’s response at sub-regional scales).

[Figure omitted. See PDF.]

Figure 3. Change in Realized Thermal Habitat Index among species within Canada and USA.

Boxplots represent distribution of net change among species within nation under the long-term (Y2060) and short-term (2030) scenarios.

https://doi.org/10.1371/journal.pone.0090662.g003

Ecological perspective

Trophic balance among functional groups contributes to ecosystem stability. Overfishing of one functional group can result in an explosion of their prey and competitors and lead to trophic imbalance [33]. Colder less diverse systems succumb more quickly to overfishing, and take longer to recover due to trophic imbalance than warmer, species rich regions [34]. From an ecosystem perspective, we examined whether one functional group is forecast to lose more realized thermal habitat than another, which would exacerbate any imbalance. Functional group membership for each species is listed in Table S2 and are derived from Shackell et al [35].

From an ecological perspective, in 2060 (Scenario 1) in Canada, functional groups differed significantly, according to a rank-sum test (Kruskal-Wallis chi-squared  =  15.3, df  =  6, p =  0.02) but not in the USA (Kruskal-Wallis chi-squared  =  4, df  =  6,p  =  0.7) (Figure 4A). Both nations follow a similar pattern, yet the difference is that small benthivores are expected to lose more realized thermal habitat in Canada than in the USA. The net change in planktivores (sandlance, herring, capelin) is more negative in both nations than piscivores; significantly so in Canada (Kruskal-Wallis chi-squared  =  6.8, df  = 1, p<0.01) but not in the USA (Kruskal-Wallis chi-squared  =  1.2, df  =  1, p  =  0.3).

[Figure omitted. See PDF.]

Figure 4. Change in Realized Thermal Habitat Index by functional group in Canada and USA.

Panels depict results under the long-term (Y2060) (A) and short-term (Y2030) (B) scenarios. Functional group member species are listed in Table S1. Piscivore refers to piscivore. Zoopisciv refers to zoopiscivore, Lg.Benthi refers to large benthivores, Md.Benthi refers to Medium-sized benthivores, Sm.Benthi refers to small-sized benthivores. Decapod B refers to decapod benthivores, Planktivo refers to planktivores, FilterFee refers to filter feeders (scallop).

https://doi.org/10.1371/journal.pone.0090662.g004

Differences among functional groups are not significant in Scenario 2 in Canada (Figure 4B) (Kruskal-Wallis chi-squared  =  12.12, df  =  6, p  =  0.06), nor in the USA: (Kruskal-Wallis chi-squared  =  3.1, df  =  6, p-value  =  0.8) (Figure 4B). Under a weak warming scenario, planktivores uniformly lose more realized thermal habitat than other functional groups, but this does not amount to a statistically significant difference.

Economic perspective

How will the change in the realized thermal habitat index affect non-commercial and commercial species? As an exercise, we weighted net change of the realized thermal habitat index by the 2011 value of landings in categories of non-commercial, low, medium and highly commercial species separately for Canada and the USA (Table S1). By 2060, currently highly commercial species in Canada would gain realized thermal habitat, while USA species of the same category would lose (Figure 5). In 2030, results are similar but dampened. These patterns are driven by lobster which was an order of magnitude of greater value than other highly valued species in Canada (scallop, snow crab) in 2011 (Table S1). Non-commercial species will lose realized thermal habitat while species of low commercial value will gain in the short and long-term scenarios in both countries. There is very little change in the medium commercial category in the long term, while in the short term, these species will gain realized thermal habitat in the USA.

[Figure omitted. See PDF.]

Figure 5. Change in Realized Thermal Habitat Index by commercial category in Canada and USA.

Panels depict results under the long-term (Y2060) (A) and short-term (Y2030) (B) scenarios. Commercial categories were defined by quantiles of 2011commercial value estimated separately for Canada and USA. The y-axis represents the weighted mean where the “Change in realized thermal habitat” index was weighted by species commercial value within categories. Values for each country were derived separately from 2011 commercial landing and value statistics for USA( http://www.st.nmfs.noaa.gov/commercial-fisheries/commercial-landings/monthly-landings/index) and Canada (http://www.dfo-mpo.gc.ca/stats/commercial/land-debarq/sea-maritimes/s2011av-eng.htm) and used to calculate weighted mean and confidence intervals.

https://doi.org/10.1371/journal.pone.0090662.g005

Species Perspective

Which species will gain or lose significant realized thermal habitat? The national, economic and ecological perspectives above provide overviews of the system whereas species-level information might arguably be used most by managers to set priorities. To lessen uncertainty in the forecast, we estimated net change in the realized thermal habitat index based solely on the most likely habitat, as defined by the probability estimate from the original model of >0.54 likelihood. In addition, we ran the model 10 times on subsamples (60%) to gauge the error surrounding net change in the realized thermal habitat index. In this manner, decisions can be based on level of certainty of results. By our definition, the original full model did not estimate highly probable habitat for the following species in the following regions: radiated shanny, cunner, and dory in Canada and USA, summer flounder in Canada, capelin, halibut, moustache sculpin and wolfish in the USA. Consequently, these species were not included in habitat gain/loss predictions for the relevant regions.

By 2060, 55% of the species examined will lose realized thermal habitat (<– 10%), 21% will gain (>10%), while 24% will stay neutral in Canada (Figure 6A). In the USA, 65% of the species will lose realized thermal habitat (<–10%), 20% will gain (>10%), while 15% will stay neutral (Figure 6B). By 2030, 12% of the species examined will lose realized thermal habitat (<– 10%), 14% will gain (>10%), while 74% will stay neutral in Canada (Figure 6C). By 2030, in Canada, losers with reasonable error bars include: moustache sculpin, sandlance, capelin, windowpane, ocean pout, while winners include: hagfish, lobster, jonah crab, shortfin squid, and monkfish. By 2030, in the USA, 8% of the species will lose realized thermal habitat (<– 10%), 10% will gain (>10%) while 83% will stay neutral (Figure 6D). By 2030, in the USA, losers include: sandlance, red crab and smooth skate while winners include: summer flounder, white hake and shortfin squid.

[Figure omitted. See PDF.]

Figure 6. Change in Realized Thermal Habitat Index of individual species for long-term (2060) and short-term (2030) scenarios

. Each boxplot represents the distribution of net change calculated from 10 models fit to subsamples of each species. Dashed lines are at 10%. Net change is considered neutral if median net change among 10 subsamples is between dashed lines. Positive outliers >100 are not shown. Change in habitat under scenario 1 (Y2060) in Canada (A) and the USA (B) and under scenario 2 (Y2030) in Canada (C) and USA (D). Abbreviated species names are listed beside full common names in Table S2.

https://doi.org/10.1371/journal.pone.0090662.g006

Discussion

We have provided a broad overview of how projected warming may affect species realized thermal habitat indices on the US Northeast Continental and Scotian Shelves. We narrowed down the most likely scenarios on a regional-scale and then provided an overview of the expected response to warming from national, ecological, economic and species perspectives. We have built on previous global-scale studies by taking the necessary step of analyzing finer-scale biological and oceanographic data [6]. On a regional scale (USA and Canada), most species realized thermal habitat will withstand warming projected for 2030, but not for 2060. Importantly, the magnitude of the 2060 warming projection (∼1.5–3C) was observed in 2012 across the northern part of the study domain [36]. We did not focus on the effect of current extremely anomalous years but it would be useful to pursue how species respond to extreme events in that it may affect distribution and productivity [37], [38] and ultimately local economies.

Populations may respond initially through a change in distribution, but over the long term, changes in productivity, predation, competition, species adaptations and fishing pressure will undoubtedly play much larger roles than warming [5], [16], [39], [40]. The long-term Scenario 1 (Y2060) results enabled us to gauge model performance, the extent that the community will change over time, and have raised the question of how current extreme events will affect productivity, and whether we should expect more anomalously warm years under climate change. Our broad overview on a regional scale showed that, under the more realistic short-term (Scenario 2: Y2030), only a few species are expected to lose realized habitat in both Canada and the USA, but this should be considered a starting point. As northern temperate/boreal species lose habitat, other species not considered in this analysis, including currently more southerly distributed species, will fill empty niches. This effect will be most pronounced in the Gulf of Maine/Georges Bank area where many temperate/boreal species are at the southern limit of their range. Initial analyses show that on smaller sub-regional scales, there is considerably more spatial variation in response where southerly distributed cod populations lose more realized thermal habitat than in the northeasterly subregions (Figure S1). Our approach can easily be used to aggregate at smaller spatial scales and this is an area of future research.

Planktivores, a key part of this region’s food chain, are expected to lose significant realized thermal habitat in both countries; this has implications for all species that eat planktivores. That is, changes in distribution and productivity of prey could result in indirect changes in predator distribution and productivity. However, planktivores, such as herring, spend only part of their time near the bottom and they are not caught as easily as groundfish by the surveys used in this analysis. To assess the potential direct and indirect impacts of warming on distribution and trophic balance, a next step, which is beyond the scope of this study, could be to investigate planktivore catchability. This would, of course, require accurate catchability indices for all species for each survey.

SDM as a Decision-making Tool: Caveats

Regarding those species whose expected net change is not neutral (more or less than 10%), it is important that any interpretation accounts for two issues, 1) the confidence of the net change estimate, and 2) the current extent of the species distribution. First, the range of net change among up to 10 models per species is an indication of response prediction uncertainty. When the range of values is large, definitive statements cannot be made and managers would have to seek further corollary information on the species of interest. In this analysis, those species would include: hagfish, jonah crab, halibut, barndoor skate, red crab and cusk in either USA or Canada. The first step towards improving these species estimates would be to aggregate results for each species on smaller spatial scales, and to seek and examine other environmental covariates that may be important to the species. Cusk, for example, are associated with terrain complexity and may not even be catchable on highly complex untrawlable bottom [41]. Hare et al. [41] modelled cusk habitat and included a covariate important to cusk: terrain ruggedness. While their model is not directly comparable to ours because they defined a broader window of “potential” habitat, their study illustrates the benefits of focussing on a single species.

Even when the sample error suggests confidence in our forecast, corollary information on the species of interest should be pursued. Sandlance is an expectant loser common across both scenarios and in both countries (Figure 6). Even though confidence is high, further research should be pursued to corroborate or qualify results as sandlance is an important forage fish.

Second, there is a positive relationship between abundance and distribution [42]. Widespread animals are generally more adaptable and can withstand a greater range of environmental conditions. There was no relationship between change in realized thermal habitat and current extent of distribution on the short-term, but when the thermal envelope is shifted enough, as in the long-term scenario, the change was only neutral for widespread species, whose extent of distribution is more than ∼22% of the total area (Figure 7). In effect, widespread species have a greater thermal window in the study domain. We might also expect that less widespread species will have a greater negative response in anomalously warm years. The relationship between net change and distribution (Figure 7) also provides context to the impact of a large net change; the realized thermal habitat index for a species may be projected to decrease by 75%, but the current total area can be quite small.

[Figure omitted. See PDF.]

Figure 7. Change in Realized Thermal Habitat Index habitat and current occupancy.

Lower rates of change are expected for widespread animals but the effect is only noticeable under Y2060 higher average temperatures or extreme events, as occurred in 2012.

https://doi.org/10.1371/journal.pone.0090662.g007

Our approach is flexible for managers. We arbitrarily designated +/–10% net change in the realized thermal habitat index as a window of neutrality. A fishery manager focused on depleted, threatened or endangered species may consider that window too narrow whereas those focussed on an increasingly abundant species would consider the window too broad. In either case, auxiliary information should always be used to provide context. Several species considered herein are currently depleted but were historically important ecological and commercial species. Some of these species have not fully recovered their size structure or former biomasses [43], [44] which means they are even more susceptible to climate variability [45]. It is probable that in 20 years, depleted species on the southern limit of their range will be extirpated from the USA as has been observed for terrestrial species [46]. In other words, the forecast of net change must be interpreted in the context of the species; extirpation will likely occur earlier for unhealthy populations.

SDM as a Decision-making Tool: Direction

In our view, forecast uncertainty is not unique to SDMs in climate change research and that uncertainty will not diminish fast enough for practical application. Natural resource management agencies still have to tackle how climate change will affect their ability to safeguard resources. Twenty and 50 year forecasts derived from SDMs are insightful, but are they useful to fishery agencies who plan on much shorter time scales? We echo the sentiment of others [47] by suggesting that an imminent and promising use of SDMs is to incorporate a thermal habitat indicator into current stock assessments. We also suggest this would be relatively easy to implement. There has been considerable discussion on, and progress towards, “Ecosystem-based Management” over the last 25 years (reviewed in [48]). Many of the “Ecosystem Indicators” are climate–related, such as bottom temperature and sea surface temperature, and are routinely monitored by many nations. It follows that oceanographic information gathered for ecosystem-based management can be used to incorporate climate change directly into stock assessments by, for example, using the realized thermal habitat index. Climate projections are uncertain partly because natural variability often swamps the signal of gradual change. A practical approach would be to adjust quotas based on amount of thermal habitat either solely or in a production model [49]. Over time, quotas would be adjusted up or down, in concert with gradual climate change. In fact, a “realized thermal habitat index” should be incorporated into stock assessments if only to monitor the effect of extreme events on subsequent productivity.

It is common to use a risk management approach when dealing with climate adaptation issues; priority issues are identified and auxiliary information is sought when uncertainty is high. SDMs from ecological, economic and species level perspectives as presented here can be used as a starting point to develop fuller assessments and adaptation strategies to counteract climate change impacts [50]. Habitat availability may not be the climate change “bottleneck” for common marine benthic species in the Northwest Atlantic in the next 20 years but the effect of current extremely warm years is unknown. To that end, we can now ask harder questions: will extreme warming events become more frequent? How will change in the base of the food chain affect overfished populations? Forecasts derived from SDMs may be only a starting point, but are necessary and practical as has been demonstrated for other marine species [51]. SDMs can provide context for a larger vulnerability assessment that would include some estimate of productivity response, diet availability, other important climate drivers such as dissolved oxygen and pH [52] and/or the inclusion of current stressors that would affect or weaken the ecosystem’s ability to respond [53].

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Percent change in realized habitat index for cod (Gadus morhua) in 6 subregions. Subregions are defined in Shackell et al 2012 [1] and describe ecoregions from Cape Breton Canada to Cape Hatteras, USA. MAB = Mid-Atlantic Bight, SNE = Southern New England, GB = Georges Bank, BoF = Bay of Fundy, WSS  = Western Scotian shelf, ESS = eastern Scotian Shelf (see Figure 1 in[1]). No estimates were made for MAB where amount of cod most likely habitat was too low.

https://doi.org/10.1371/journal.pone.0090662.s001

(DOCX)

Table S1.

Summary of bottom temperature records. Mean, standard deviation (SD), minimum (min) and maximum (max) bottom temperature and number of sets within month, year and region using cod (Gadus morhua) sub-file as an example.

https://doi.org/10.1371/journal.pone.0090662.s002

(DOCX)

Table S2.

Summary of Species, Functional Groups, Commercial Categories, Name abbreviations. Net change in Realized Thermal Habitat Index in Scenario 2 is also shown. NA refers to not enough probable habitat (insufficient probability estimates >0.5).

https://doi.org/10.1371/journal.pone.0090662.s003

(DOCX)

File S1.

Supplementary Information on Surveys and Analyses.

https://doi.org/10.1371/journal.pone.0090662.s004

(DOCX)

Acknowledgments

We greatly appreciate the thoughtful reviews of MJ Fogarty and B Helmuth. We also thank WK Li, CG Hannah, KT Frank and P Greyson for discussion.

Author Contributions

Conceived and designed the experiments: NLS. Performed the experiments: NLS DR CS. Analyzed the data: NLS. Contributed reagents/materials/analysis tools: NLS DR CS. Wrote the paper: NLS DR CS.

Citation: Shackell NL, Ricard D, Stortini C (2014) Thermal Habitat Index of Many Northwest Atlantic Temperate Species Stays Neutral under Warming Projected for 2030 but Changes Radically by 2060. PLoS ONE 9(3): e90662. https://doi.org/10.1371/journal.pone.0090662

References

1. Parmesan C, Yohe G (2003) A globally coherent fingerprint of climate change impacts across natural systems. Nature 421: 37–42

2. Chen D, Achberger C, Raisanen J, Hellstrom C (2006) Using statistical downscaling to quantify the GCM-related uncertainty in regional climate change scenarios: a case study of Swedish precipitation. Adv Atmos Sci 23: 54–60.

3. McCay BJ (2012) Anthropology: Shifts in fishing grounds. Nat Clim Chang 2: 840–841

4. Cheung WWL, Lam VWY, Sarmiento JL, Kearney K, Watson R, et al. (2009) Projecting global marine biodiversity impacts under climate change scenarios. Fish Fish 10: 235–251

5. Cheung W, Lam V, Pauly D (2008) Modelling present and climate-shifted distribution of marine fishes and invertebrates. Fish Cent Res Reports 16: 1–72.

6. Pinsky ML, Worm B, Fogarty MJ, Sarmiento JL, Levin SA, et al. (2013) Marine taxa track local climate velocities. Science 341: 1239–1242

7. Perry AL, Low PJ, Ellis JR, Reynolds JD (2005) Climate change and distribution shifts in marine fishes. Science 308: 1912–1915

8. Nye JA, Joyce TM, Kwon Y-O, Link JS (2011) Silver hake tracks changes in Northwest Atlantic circulation. Nat Commun 2: 412

9. Pinsky ML, Fogarty M (2012) Lagged social-ecological responses to climate and range shifts in fisheries. Clim Change 115: 883–891

10. Walther G, Post E, Convey P, Menzel A, Parmesan C, et al.. (2002) Ecological responses to recent climate change. Nature: 389–395.

11. Loder J, van der Baaren A (2013) Climate change projections for the Northwest Atlantic from six CMIP5 Earth System Models. Can Tech Rep Hydrogr Ocean Sci 286: xiv + 112 p.

12. Wake DB, Hadly EA, Ackerly DD (2009) Biogeography, changing climates, and niche evolution. Proc Natl Acad Sci 106: 19631–19636.

13. Pearson R, Dawson T (2003) Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful? Glob Ecol Biogeogr 12: 361–371.

14. Kearney M (2006) Habitat, environment and niche: what are we modelling? Oikos 115: 186–191

15. Helmuth B, Kingsolver JG, Carrington E (2005) Biophysics, physiological ecology, and climate change: does mechanism matter? Annu Rev Physiol 67: 177–201

16. Hale R, Calosi P, McNeil L, Mieszkowska N, Widdicombe S (2011) Predicted levels of future ocean acidification and temperature rise could alter community structure and biodiversity in marine benthic communities. Oikos 120: 661–674.

17. Beale C, Lennon J (2012) Incorporating uncertainty in predictive species distribution modelling. Philos Trans R Soc B Biol Sci 367: 247–258.

18. Thuiller W (2004) Patterns and uncertainties of species’ range shifts under climate change. Glob Chang Biol 10: 2020–2027.

19. Guisan A, Thuiller W (2005) Predicting species distribution: offering more than simple habitat models. Ecol Lett 8: 993–1009.

20. UNESCO/IOC (2012) Northeast Fisheries Science Center Bottom Trawl Survey Data (OBIS-USA). Available: iobis.org/mapper/?resource_id = 237. Accessed 1 July 2012.

21. UNESCO/IOC (2012) DFO Maritimes Research Vessel Trawl Surveys Fish Observations (OBIS Canada). Available: iobis.org/mapper/?resource_id = 237. Accessed 1 July 2012.

22. Seo C, Thorne J, Hannah L, Thuiller W (2009) Scale effects in species distribution models: implications for conservation planning under climate change. Biol Lett 5: 39–43.

23. Diniz-Filho J, Bini M, Rangel F, Loyola R, Hof C, et al. (2009) Partitioning and mapping uncertainties in ensembles of forecasts of species turnover under climate change. Ecography (Cop) 32: 897–906

24. Fronzek S, Carter TR, Luoto M (2011) Evaluating sources of uncertainty in modelling the impact of probabilistic climate change on sub-arctic palsa mires. Nat Hazards Earth Syst Sci 11: 2981–2995

25. Hastie T, Tibshirani R (1990) Generalized Additive Models. London: Chapman and Hall.

26. Araujo MB, Cabeza M, Thuiller W, Hannah L, Williams PH (2004) Would climate change drive species out of reserves? An assessment of existing reserve-selection methods. Glob Chang Biol 10: 1618–1626

27. R Core Team (2012) R: A Language and Environment for Statistical Computing. R Found Stat Comput Vienna, Austria. Available: http://www.r-project.org/.

28. Wood S (2006) Generalized Additive Models: An Introduction. Boca Raton, Florida.: R. Chapman and Hall/CRC.

29. Tuszynski J (2012) caTools: Tools: moving window statistics, GIF, Base64, ROC AUC, etc. R Packag version 114. Available: http://CRAN.R-project.org/package=caTools.

30. Gritti ES, Gaucherel C, Crespo-Perez M-V, Chuine I (2013) How can model comparison help improving species distribution models? PLoS One 8: e68823

31. Loder J, Wang Z, van der Baaren A, Pettipas R (2013) Trends and variability of sea surface temperature (SST) in the North Atlantic from the HadISST, ERSST and COBE datasets. Can Tech Rep Hydrogr Ocean Sci 292: viii + 36 p.

32. Araújo M, Pearson R, Thuiller W, Erhard M (2005) Validation of species–climate impact models under climate change. Glob Chang Biol 11: 1504–1513

33. Frank K, Petrie B, Choi J, Leggett W (2005) Trophic cascades in a formerly cod-dominated ecosystem. Science 308: 1621–1623.

34. Frank K, Petrie B, Shackell N, Choi J (2006) Reconciling differences in trophic control in mid-latitude marine ecosystems. Ecol Lett 9: 1096–1105.

35. Shackell NL, Bundy A, Nye JA, Link JS (2012) Common large-scale responses to climate and fishing across Northwest Atlantic ecosystems. ICES J Mar Sci 69: 151–162

36. Hebert D, Pettipas R, Brickman D, Dever M (2013) Meteorological, Sea Ice and Physical Oceanographic Conditions on the Scotian Shelf and in the Gulf of Maine during 2012. DFO Can Sci Advis Sec Res Doc 2013/058: v + 46 p.

37. Wernberg T, Smale D, Thomsen M (2012) A decade of climate change experiments on marine organisms: procedures, patterns and problems. Glob Chang Biol 18: 1491–1498.

38. Jentsch A, Kreyling J, Beierkuhnlein C (2007) A new generation of climate change experiments: Events, not trends. Front Ecol Environ 5: 65–374

39. Brander K (2007) Global fish production and climate change. Proc Natl Acad Sci USA 104: 19709–19714.

40. Simpson S, Jennings S, Johnson M, Blanchard J, Schon P-J, et al. (2011) Continental shelf-wide response of a fish assemblage to rapid warming of the sea. Curr Biol 21: 1565–1570.

41. Hare JA, Manderson JP, Nye JA, Alexander MA, Auster PJ, et al. (2012) Cusk (Brosme brosme) and climate change: assessing the threat to a candidate marine fish species under the US Endangered Species Act. ICES J Mar Sci 69: 1753–1768

42. Gaston K, Blackburn T (2008) Pattern and process in macroecology. Wiley.com.

43. Shackell N, Frank K (2007) Compensation in exploited marine fish communities on the Scotian Shelf, Canada. Mar Ecol Prog Ser 336: 235–247

44. Shackell NL, Frank KT, Fisher JAD, Petrie B, Leggett WC (2010) Decline in top predator body size and changing climate alter trophic structure in an oceanic ecosystem. Proc R Soc Biol Sci 277: 1353–1360

45. Planque B, Fromentin J-M, Cury P, Drinkwater K, Jennings S, et al. (2010) How does fishing alter marine populations and ecosystems sensitivity to climate? J Mar Syst 79: 403–417.

46. Franco A, Hill J, Kitschke C, Collingham Y, Roy D, et al. (2006) Impacts of climate warming and habitat loss on extinctions at species’ low-latitude range boundaries. Glob Chang Biol 12: 1545–1553.

47. Link JS, Nye JA, Hare JA (2011) Guidelines for incorporating fish distribution shifts into a fisheries management context. Fish Fish 12: 461–469

48. Sherman K (n.d.) Adaptive management institutions at the regional level: The case of Large Marine Ecosystems. Ocean Coast Manag. doi:http://dx.doi.org/10.1016/j.ocecoaman.2013.06.008.

49. Choi J, Zisserson B, Cameron B (2012) Assessment of Scotian Shelf Snow Crab in 2011. DFO Can Sci Advis Sec Res Doc 2012/024: iv + 95 p.

50. Hulme PE (2005) Adapting to climate change: is there scope for ecological management in the face of a global threat? J Appl Ecol 42: 784–794

51. Russell B, Connell S, Mellin C, Brook B, Burnell O, et al. (2012) Predicting the Distribution of Commercially Important Invertebrate Stocks under Future Climate. PLoS One 7: e46554

52. Cheung WWL, Dunne J, Sarmiento JL, Pauly D (2011) Integrating ecophysiology and plankton dynamics into projected maximum fisheries catch potential under climate change in the Northeast Atlantic. ICES J Mar Sci 68: 1008–1018

53. Niiranen S, Yletyinen J, Tomczak M, Blenckner T, Hjerne O, et al. (2013) Combined effects of global climate change and regional ecosystem drivers on an exploited marine food web. Glob Chang Biol 19: 3327–3342.