This study describes a method for modeling and predicting, from biological and physical variables, habitat use by a commercially harvested groundfish species. Models for eastern Bering Sea flathead sole Hippoglossoides elassodon were developed from 3 relationships describing the response of organism abundance along a resource continua. The model was parameterized for 1998 to 2000 trawl survey data and tested on 2001 and 2002 data. Catch per unit effort (CPUE) of flathead sole had a curvilinear relationship with depth, peaking at 140 m, a proportional relationship with bottom water temperature, a positive curvilinear relationship with potential cover (invertebrate sheltering organisms such as anemones, corals, sponges, etc.), a negative relationship with increasing mud:sand ratio in the sediment, and an asymptotic relationship with potential prey abundance. The predicted CPUE was highly correlated (r 2 = 0.63) to the observations (1998 to 2000) and the model accurately predicted CPUE (r 2 = 0.58) in the test data set (2001 and 2002). Because this method of developing habitat-based abundance models is founded on ecological relationships, it should be more robust for predicting fish distributions than statistically based models. Thus, the model can be used to examine the consequences of fishing activity (e.g. reduction in sheltering organisms), changes in temperature (e.g. climate effects) and interaction between variables, and can be modified to incorporate new variables as more information is collected about a species.
KEY WORDS: Fish habitat · Habitat model · Bering Sea · Soft sediment · Flatfish
Resale or republication not permitted without written consent of the publisherMar Ecol Prog Ser 290: 251-262, 2005al. 1992, O'Brian & Rago 1996, Stoner et al. 2001. General additive models (GAMs) have the advantage of exhibiting nonlinear response curves describing the relationship between fish distributions and habitat, but the shape of the response curves may be difficult to explain in terms of organism ecology. It is notable that existing statistical habitat models are rarely used for prediction (but see Stoner et al. 2001 for exception), because the statistical relationships describing existing distributions may not be applicable to future distributions (Beutel et al. 1999). Utilization of ecological relationships provides the advantages of nonlinearity in modeling distributions across habitats, a distribution response for each habitat variable that is justifiable based on organism ecology, and therefore a more robust prediction of future organism distribution across habitats than traditional statistical models.The characterization of habitat use is also dependent on the habitat variables used in the analysis. For example, Swartzman et al. (1992) modeled 5 eastern Bering Sea flatfish species (rock sole Lepidopsetta sp., Alaska plaice Pleuronectes quadrituberculatus, flathead sole Hippoglossoides elassodon, yellowfin sole Pleuronectes asper and Greenland turbot Reinhardtius hippoglossoides) using a GAM wit...