Habitat selection by American mink during Summer is related to hotspots of crayfish prey

Habitat selection by mammalian carnivores may be driven by prey availability, physical characteristics of the habitat, and landscape context. However, the cues used by carnivores to select habitat are often unclear. We examined the seasonal diet of American mink Neovison vison and determined if the abundance of a primary prey, crayfish, was an important driver of habitat use during summer in an agricultural landscape in Illinois. We also evaluated effects of stream size, water depth, riparian buffer width, and urbanization on occupancy of stream segments by mink. We collected mink scats during three seasons and tested for seasonal differences in the percentage of occurrence and volume percentage of prey classes in the diet of mink. Crayfish remains were the dominant component of mink scats during summer. In summer 2012, we performed occupancy surveys for mink and concurrently measured crayfish densities and habitat features in 59 stream segments. Site occupancy by mink was related positively to presence of local areas with high crayfish concentrations (hotspots) instead of local habitat characteristics that might indicate high prey densities. Mink also were associated negatively with degree of urbanization and stream size. Our study highlights the effectiveness of integrating data on diets and occupancy modeling to obtain insights on cues used by carnivores to select habitat.

Habitat selection by mammalian carnivores may be driven by prey availability, physical characteristics of the habitat, and landscape context. However, the cues used by carnivores to select habitat are often unclear. Within a home range, carnivores should select habitat patches that maximize the probability of encountering prey. Th ey may choose locations for hunting based on habitat characteristics that should indicate high prey abundance (Irwin et al. 2007, Slauson et al. 2007, or based directly on prey abundance (O ' Donoghue et al. 1998, Fukui et al. 2006. Because prey habitat does not necessarily guarantee actual occurrence of prey at any given time, using direct measures of prey abundance could be more informative than using the amount of prey habitat as a surrogate for prey abundance (Keim et al. 2011). Moreover, habitat selection studies should be paired with mechanistic data such as diet metrics (Keim et al. 2011) to understand carnivore behavior more fully.
One such carnivore that may select habitat based on prey abundance is the American mink Neovison vison . Th e species is an opportunistic, generalist carnivore native to North America and invasive in Europe, South America and Asia (Macdonald andHarrington 2003, Bonesi andPalazon 2007). Th is semiaquatic mustelid is closely associated with watercourses, but is capable of spending extended periods of time away from water. For example, in central Illinois, USA, radio-marked mink were located Ͼ 100 m away from the water ' s edge 14% of the time (A. Ahlers pers. comm.). Mink have seasonally variable diets due to temporal changes in relative availabilities of prey (Gerell 1967, Dunstone and Birks 1987, Brzezi ń ski 2008. Th e diet of mink has been examined within its native range (Dearborn 1932, Sealander 1943, Korschgen 1958, Racey and Euler 1983, Arnold and Fritzell 1987, Ben-David et al. 1997, Hoff man et al. 2009, and some studies have examined habitat use by mink in its native range (Burgess and Bider 1980, Arnold and Fritzell 1990, Ben-David et al. 1995, Stevens et al. 1997, Loukmas and Halbrook 2001. However, studies linking diet with habitat selection behavior of mink are scarce. Crayfi sh often occur in the diet of American mink where crayfi sh are available both in the native (Dearborn 1932, Burgess andBider 1980) and invasive range (Gerell 1967, Day and Linn 1972, Ward et al. 1986, Melero et al. 2008, Fischer et al. 2009). In the invasive range, mink population density is correlated positively with the proportion of crayfi sh in the diet, and home range size is correlated negatively with the crayfi sh contribution (Melero et al. 2014). Spatial variability in abundance of this aquatic prey could explain habitat selection by mink during seasons when crayfi sh are plentiful.
Mink may select habitat in relation to crayfi sh abundance in two general ways. Mink could use habitat cues to decide where to hunt for crayfi sh, or they could respond directly to relative crayfi sh densities. More specifi cally, because crayfi sh are patchily distributed within streams (DiStefano et al. 2003), mink may respond to presence of prey hotspots instead of to average prey densities for sites. Studies across multiple taxa indicate prey hotspots can be important determinants of predator space use and foraging behavior (Th ompson et al. 2001, Davoren et al. 2003, Gende and Sigler 2006, L ó pez-Bao et al. 2011. Mink use some parts of their home range more intensively than others (Yamaguchi et al. 2003), which could refl ect prey distributions and suitable hunting places (Gerell 1970).
In addition to prey abundance, local habitat and landscape characteristics may be important correlates of habitat selection by mink. In the Midwestern United States, agricultural intensifi cation, including drainage for row crops, and urbanization have reduced the amount of suitable habitat for semiaquatic wildlife (Zucker and Brown 1998), and limited the terrestrial habitat available to mink surrounding the small streams and agricultural ditches in which they occur. Th ese riparian buff ers vary considerably in characteristics such as width and vegetation structure (Ahlers et al. 2010). Larger riparian buff ers provide mink with increased foraging space and terrestrial alternatives to the aquatic prey located within the stream channel. Large streams have a greater diversity of aquatic prey than do small streams (Sheldon 1968, Osborne andWiley 1992), and mink are associated positively with water depth of streams (Schooley et al. 2012). Additionally, mink might avoid areas of human development (Racey andEuler 1983, Brzezi ń ski et al. 2012), although information regarding eff ects of urbanization on mink is limited (Gehrt et al. 2010).
Climate change is expected to increase variability in precipitation in the Midwestern United States; climate models predict an increase in frequency of summer drought and spring fl ooding events (Wuebbles andHayhoe 2004, Trenberth 2011). Th ese changes will lead to increased temporal fl uctuations in water depths and fl ow regimes of streams that could have consequences for predator -prey interactions in riparian ecosystems. Th e severe drought of 2012 (Illinois Department of Natural Resources 2013) presented an ideal opportunity to examine predator -prey interactions during environmental conditions expected to occur more frequently under climate change.
We examined the seasonal diet of American mink and determined if a primary prey, crayfi sh, was an important driver of habitat selection by mink during summer of a severe drought year in an agricultural landscape. We hypothesized that diet of mink would refl ect seasonal availability of prey. We predicted crayfi sh would be most important in the diet during summer when they are most available, whereas consumption of crayfi sh would decrease during winter when mink shift to prey such as mammals and fi sh. We also hypothesized that habitat selection by mink during summer would be related directly to prey abundance. Hence, we expected a higher probability of occupancy for mink at sites with crayfi sh hotspots or higher mean densities. We also hypothesized that stream and landscape characteristics would infl uence habitat selection by mink. We predicted site occupancy by mink would be related positively to water depth and stream size because larger streams contain more aquatic resources. We also predicted a positive association between occupancy and riparian buff er width because wider buff ers provide more foraging opportunities. Finally, we expected a negative association between site occupancy by mink and degree of urbanization.

Study area
Our study was conducted in east -central Illinois, USA and was centered on Champaign-Urbana (40 ° 12 ′ N, 88 ° 26 ′ W). Fifty-nine sites were randomly selected from a previous stratifi ed random sample of 90 sites, 50% of which were located within a 2-km radius of incorporated towns ( Ͼ 2500 people), and 50% of which were located outside of this urban buff er (Cotner and Schooley 2011). Sampling for crayfi sh density and site occupancy by mink occurred at these 59 study sites. Each site was a 200-m stretch of wadeable stream, ranging from 1st to 5th order in size (average site wetted width ϭ 0 -23.3 m; average site depth ϭ 0 -0.97 m). Sites represented potential resource patches for mink so our measure of site occupancy corresponded to habitat use (Schooley et al. 2012). Th e median distance between one site and the nearest site was 3.4 km (range ϭ 0.5 -13.5 km). Sites were distributed across an urbanization gradient (proportion of impervious surface within a 500-m buff er around each site, Cotner and Schooley 2011), and had a wide range of riparian buff er widths (0 -466 m). We also collected mink scats at an additional 60 locations to increase our sample size.

Scat collection
To evaluate seasonal variation in composition of mink diets, we collected scats during fall, winter and summer. Mink scat was identifi ed by its unique, twisted appearance with tapered ends. Mink scat is distinguishable from other mammal scat in the area based on size and appearance (Rezendes 1999). However, if uncertainties in identifi cation existed, the scat was not collected. Scat samples from a single site were combined in plastic bags and stored at Ϫ 15 ° C. In fall 2011 (21 September -10 November), our 59 study sites were surveyed for mink scat by two trained searchers; scat was found at 18 sites. Th ere was uncertainty regarding the source of one collected scat sample, resulting in an analysis of 17 scats from 18 sites for fall. In winter 2012 (4 January -9 March), the 18 sites were revisited and searched, and all found scats were collected. To maximize effi ciency of our searching, we did not revisit the 41 sites at which we did not fi nd scat in fall. To increase our sample size, we instead searched an additional 60 locations within our study area. Th ese locations were underneath and within 50 m of bridges in areas that often included rock and other substrates where mink typically deposit scat. Th e collection locations were Ն 1.5 km apart. Th is separation distance was based on an estimate of average length of mink home ranges in the study area (A. Ahlers pers. comm.), and increased the likelihood of independence of scat samples. In summer 2012 (25 June -31 July), the 18 sites and 60 scat collection locations were revisited. Scat was also collected opportunistically within the same three seasons during a concurrent study of radio-marked mink.

Diet analysis
Each scat was soaked in warm water to facilitate separation of prey remains, washed through a sieve (0.8-mm mesh), and air-dried. We sorted remains into seven prey classes (crayfi sh, mammal, bird, fi sh, insect, reptile and unknown) under a dissecting microscope (10 ϫ ) based on hair, teeth, bones, feathers, scales and exoskeleton fragments. Th e unknown prey class was mostly comprised of unidentifi able bone fragments. Th ese general prey classes were adequate for testing our main prediction regarding seasonal shifts for crayfi sh in mink diets.
Diet composition was recorded for each of the three seasons (fall, winter, summer) using three metrics: percentage of occurrence (PO; the number of occurrences of each prey class divided by the total number of scat samples, times 100), relative frequency of occurrence (RFO; the number of occurrences of each prey class divided by the total number of occurrences of identifi ed prey), and volume percentage (VOL; visually estimated as the percentage of each prey class in each scat). PO and RFO were highly correlated positively for all prey classes (r Ͼ 0.80), so we used PO for analyses. Each metric comes with caveats. PO may overestimate less digestible prey and underestimate more digestible prey, and overestimate the contribution of prey taken regularly but in small amounts. VOL does not account for variation in scat sizes and also could underestimate the contribution of highly digestible prey (Klare et al. 2011). To minimize the infl uence of biases associated with a single metric, we used both PO and VOL methods in our analyses (Zabala and Zuberogoitia 2003). However, extrapolations from scat samples to the actual relative amounts of each prey type consumed should be made with caution. When making seasonal comparisons with PO data, it should be more diffi cult to detect diff erences among prey classes because the contribution of rare food items is exaggerated. VOL data are more likely to reveal diff erences in the relative amounts of prey types in scats, and thus detect seasonal specialization even if the overall variety of diet items per scat changes little. Because we expect biases due to digestibility of prey types to remain constant across time, and our analyses focus on relative diff erences among seasons, conclusions about seasonal variation in mink diet should be robust.
For the PO metric, we tested for variation in mink diet among seasons using χ 2 -tests and Fisher ' s exact tests when Ն 20% of the expected frequencies were Ͻ 5 (Zar 1984). For the VOL metric, we tested for between-season diff erences in diet composition using multi-response permutation procedures (MRPP - Mielke andBerry 2001, Roberts andTaylor 2008). We employed a Bonferroni correction when making multiple, pair-wise seasonal comparisons ( α ϭ 0.0167). An eff ect size ' A ' was calculated to measure the overall dietary agreement among scat samples within the same season (McCune and Grace 2002, Roberts and Taylor 2008). Within-season homogeneity of scat samples is greater than expected by chance when A Ͼ 0, equal when A ϭ 0, and less when A Ͻ 0 (Roberts and Taylor 2008). Th e reptile prey class occurred in only one scat sample during fall, and was excluded from analyses of VOL. Th e MRPP analysis was conducted in PC-ORD 6.0 (McCune and Meff ord 2011).
To quantify the diversity of prey found in the mink diet, we calculated the Shannon diversity index (Shannon 1948) and dietary evenness in each season using both PO and VOL metrics. We also calculated food niche breadth using Levins (1968) in which p i is the proportion of scats containing prey class i . Diff erences in Shannon diversity index values between seasons were assessed using Student ' s t-tests, with variances calculated according to Zar (1984). Statistical tests were performed in SAS 9.2 (SAS Inst.).

Occupancy surveys and crayfi sh sampling
Occupancy surveys for mink and crayfi sh sampling were conducted at the 59 study sites from 18 May -26 July 2012, during the core months of the severe drought of 2012, which included the second driest January to July period on record in Illinois (Illinois State Water Survey 2012a, Illinois Department of Natural Resources 2013). Two trained observers independently surveyed each site once. Each observer walked along both sides of the 200-m stream segment and searched for sign within 5 m of the water ' s edge (Schooley et al. 2012). A site was considered occupied by mink if scat or tracks were detected, but we dealt with imperfect detection in our modeling. Sign was likely left by mink with established home ranges because our sampling was not during the dispersal period (Larivi è re 1999).
Concurrently, crayfi sh ( Orconectes , Procambarus , Cambarus ) density was measured in each of the 59 sites by sampling 1-m 2 areas in 10-m segments of stream (20 samples per site). We divided the stream ' s wetted width in half, randomly selected either the left or right side to begin sampling, and alternated sides as we sampled upstream. A sample was taken at the fi rst ' high-quality ' crayfi sh habitat encountered on the selected area within each 10-m segment. In streams Յ 1 m wide, the entire stream width was sampled. High-quality habitat consisted of in-stream gravel, cobble, and rocks; anchored woody debris; or submerged vegetation. If there was no high-quality habitat within a segment, we sampled the fi rst 1-m 2 area at the downstream end of the segment. A seine net was placed perpendicular to the stream fl ow, and a 1-m 2 area of substrate upstream of the net was disturbed so that all crayfi sh were washed into the seine net (Mather and Stein 1993, Flinders and Magoulick 2003, Taylor and Soucek 2010. When necessary in low-fl ow areas, we dragged the seine net through the sampling area while disturbing the substrate to collect crayfi sh. Upon collection, crayfi sh were classifi ed as juvenile ( Ͻ 15 mm carapace length) or adult, and adults were identifi ed to species. All crayfi sh species were pooled for analyses.
We quantifi ed three habitat covariates that could be correlated with crayfi sh density (Riggert et al. 1999, DiStefano et al. 2003: substrate particle size, number of crayfi sh burrows, and number of woody debris accumulations. At each sampling location, we created a transect perpendicular to the stream fl ow. We dropped a metal rod at 1-m intervals for the entire wetted width of the stream and measured the substrate particle size that the rod landed on using a gravelometer. Th ese particle sizes were averaged to create one measure of substrate particle size per 10-m segment of stream. Th e number of active crayfi sh burrows was counted in each 10-m segment of stream. We also counted within-stream and landscape covariates. We did not include two covariates in the same model if they were correlated at r Ն 0.60. We ran this candidate set alone and with each of the three measures of crayfi sh density (48 total models), and ranked models using Akaike ' s information criterion (AIC; Burnham and Anderson 2002). Akaike weights (w i ) were summed for models that contained each of the three measures of crayfi sh density to determine (using across-model inference) which measure of prey abundance was the best predictor of occupancy by mink.

Seasonal diet
We analyzed 103 scat samples (fall ϭ 17, winter ϭ 43, summer ϭ 43). Crayfi sh and mammals were the most common diet items for mink throughout the year, occurring in Ͼ 76% of all scats and comprising Ͼ 83% of the diet by volume. However, percentage of occurrence of diet items diff ered among seasons ( χ 2 ϭ 23.7, DF ϭ 12, p ϭ 0.02). Crayfi sh occurred most frequently in summer (Fisher ' s exact test: winter -summer, p ϭ 0.003; fall -summer, p ϭ 0.016; Fig. 1), when PO of mammals was lowest (Fisher ' s exact test, fall -summer, p ϭ 0.025; Fig. 1). PO of mammals in the diet was greatest during fall and winter (Fig. 1). Fish occurred in 35.9% of scats, 8.1% by volume. PO of fi sh increased from fall to winter (Fisher ' s exact test p ϭ 0.076), and both PO (44.2%) and VOL (17.5%) were greatest in winter. Birds occurred in 22.3% of scats, but only comprised 3% of the total scat volume. Similarly, insects occurred commonly in the diet (PO ϭ 48.5%) but made up a small percentage of the dietary volume (3.6%). Th e unknown prey class was most likely the remains of mammals and birds but could have included other vertebrate prey. Based on PO, the contribution of the unknown prey class to the diet did not diff er between seasons (p Ͼ 0.72 for all seasonal comparisons).
MRPP results for the VOL metric showed that diet composition for summer diff ered from diet composition for fall and winter (fall -winter, A ϭ 0.010, p ϭ 0.15; wintersummer, A ϭ 0.098, p Ͻ 0.001; fall -summer, A ϭ 0.155, p Ͻ 0.001). MRPP results highlight the greater contribution the number of woody debris accumulations anchored in the streambed in each segment of stream.

Detection and occupancy covariates
We recorded covariates that could infl uence detection (p) of mink sign including observer, Julian date, number of days since rain, and recent rainfall -total rainfall (cm) for the seven days prior to each survey (Illinois State Water Survey, station no. 118740, Urbana, Illinois). Recent rainfall could wash away sign or raise water levels to hide sign (Schooley et al. 2012).
We measured covariates for occupancy including crayfi sh density, water depth (m), average riparian buff er width (m), degree of urbanization (0 -1), drainage area (km 2 ), wetted width (m), and stream order (1 -5). Water depth (thalweg) and riparian buff er width were measured at 50-m intervals (fi ve values) and averaged for each site (Cotner and Schooley 2011). Drainage area, wetted width, and stream order were correlated positively (all r Ͼ 0.47), so we used principal components analysis (PCA) to create orthogonal principal components (PC). Th e fi rst PC (sizePC) explained 75.9% of the variation and was correlated positively with all three variables (r ϭ 0.74 -0.94), so we used sizePC as a measure of stream size in our models (Cotner and Schooley 2011). We excluded water depth from the PCA because small streams have dynamic fl ow regimes tied to local precipitation events; they fl ood and subside faster than do large streams (Ahlers et al. 2010). Th us, small streams can have deep waters during a sampling period because water depth is infl uenced by more than stream size. In addition, water depth alone can be a predictor of site occupancy and colonization by mink (Schooley et al. 2012).
Crayfi sh were patchily distributed within sites, so we evaluated three measures of crayfi sh density. For each site, we calculated average density of adult crayfi sh (adults m Ϫ2 ), average density of total crayfi sh (adults ϩ juveniles m Ϫ2 ), and presence -absence of a crayfi sh ' hotspot ' . A site was considered to contain a prey hotspot if Ն 1 of the 20 kick-seine samples had a total crayfi sh density in the top 15% ( Ն 15 crayfi sh m Ϫ2 ; n ϭ 171) of all sampled crayfi sh densities (median ϭ 2 crayfi sh m Ϫ2 , n ϭ 1,180; Supplementary material Appendix 1 Fig. A1). Although mink are likely to prey only on larger, adult crayfi sh, we used total crayfi sh density to defi ne hotspots because mink could use the presence of juvenile crayfi sh as an indicator of where adults occur.

Occupancy modeling
We estimated the probability of site occupancy ( ψ ) by mink using single-season occupancy models that accounted for imperfect detection (MacKenzie et al. 2006) in Program PRESENCE 5.8 (Hines 2006). We fi rst selected the best model for detection (p), and then modeled occupancy. We used a maximum of two covariates for p and four covariates for ψ in a single model to avoid over-parameterization. We developed a candidate set of 12 occupancy models that included within-stream covariates (stream size, water depth) and landscape-level covariates (riparian buff er width, degree of urbanization). Th e candidate set contained models that tested each covariate separately, and also combinations of   of crayfi sh to the diet in summer than in fall and winter (Fig. 2). Repeating the analysis excluding the unknown prey class did not qualitatively aff ect MRPP results. Based on the PO metric, the diversity indices for mink diet did not diff er strongly among seasons (Table 1; Shannon diversity index: fall -winter, t ϭ 0.45, DF ϭ 86, p ϭ 0.65; winter -summer, t ϭ 1.03, DF ϭ 211, p ϭ 0.31; fall -summer, t ϭ 0.19, DF ϭ 70, p ϭ 0.85). Based on the VOL metric, summer diets of mink had the lowest evenness, narrowest niche breadth, and lowest diversity ( Table 1). Diet of mink was least diverse during summer when they focused on crayfi sh (Shannon diversity index: fall-winter, t ϭ 0.19, DF ϭ 91, p ϭ 0.85; winter -summer, t ϭ 4.00, DF ϭ 244, p Ͻ 0.001; fall -summer, t ϭ 3.23, DF ϭ 103, p Ͻ 0.01).
Mink sign was detected at 18 of 59 sites (na ï ve occupancy ϭ 0.305). We decided the best model for detection was the competitive model ranked second ( Δ AIC ϭ 0.14; Table 2), and we used this model for subsequent evaluation of occupancy covariates. Th e top detection model contained observer alone, but adding rainfall to that model improved the log likelihood substantially (Table 2). Detection (p) was related positively to rainfall ( β rainfall ϭ 0.652, SE ϭ 0.501).
Akaike weights (w i ) summed across occupancy models indicated that among the three measures of crayfi sh abundance, the best predictor of site occupancy by mink was presence of a hotspot (hotspot w i ϭ 0.821, total crayfi sh density w i ϭ 0.100, adult crayfi sh density w i ϭ 0.032). All competitive models ( Δ AIC Ͻ 2) contained crayfi sh hotspot as a covariate (Table 3). Occupancy probability was related positively to the presence of a crayfi sh hotspot at a site ( β hotspot ϭ 1.721, SE ϭ 0.625; Fig. 3). Estimated occupancy from the hotspot model was 0.562 for sites with crayfi sh hotspots, and 0.187 for sites without hotspots. Site occupancy by mink was related negatively to stream size ( β sizePC ϭ -0.598, SE ϭ 0.463; Fig. 3) and urbanization Table 3. Ranking of occupancy models for American mink in Illinois based on Akaike ' s information criterion (AIC). Detection covariates included observer and rainfall for the seven days prior to each survey (rainfall). Occupancy covariates included presence-absence of a crayfi sh hotspot (hotspot), stream size (sizePC), water depth, and degree of urbanization. Δ AIC ϭ AIC for a given model minus AIC for the best model. K ϭ number of model parameters, w i ϭ Akaike weights, and LL is the log-likelihood. Competitive models ( Δ AIC Ͻ 2) and the intercept-only model without occupancy covariates are presented.   (Table 3). Gray triangles indicate estimated occupancy at sites with crayfi sh hotspots. Gray diamonds indicate estimated occupancy at sites without crayfi sh hotspots. Black circles indicate na ï ve occupancy for sites (1 ϭ occupied, 0 ϭ unoccupied).
( β urban ϭ -1.793, SE ϭ 1.606; Fig. 3). Although water depth occurred in competitive occupancy models (Table 3), water depth did not substantially increase model fi t based on log likelihoods; the top two models were essentially unchanged by including depth as a covariate. In addition, the model with depth alone performed worse than the intercept-only model. Riparian buff er width also was not a good predictor of site occupancy.

Discussion
Th e spatial variation in abundance of a common prey, crayfi sh, was the primary predictor of habitat selection during summer by American mink within a humandominated landscape. Mink seemed to select sites based directly on high prey concentrations instead of habitat characteristics that might indicate high prey densities. Mink shifted strongly to feeding mainly on crayfi sh during summer, and mink were more likely to occupy stream segments that contained crayfi sh hotspots. Crayfi sh hotspots were a far better predictor of occupancy by mink than were average crayfi sh densities. Habitat occupancy by mink was also aff ected by stream and landscape characteristics; mink were associated negatively with both stream size and urbanization.
To increase the chance of encountering prey and maximizing energetic gains in a patchy environment (MacArthur andPianka 1966, Charnov 1976), predators may either select prey habitat or select locations most used by prey (Flaxman andLou 2009, Keim et al. 2011). Mink in our study appeared to select foraging habitat based on locations most used by prey -hotspots with high densities of adult and juvenile crayfi sh. Given that crayfi sh density was not highly correlated with any measured habitat variable, mink likely cue directly on crayfi sh. Th is tactic may be particularly eff ective in human-dominated landscapes in which stream habitat structure is altered due to channelization, habitat heterogeneity is low, and prey distribution might mostly refl ect spatial population dynamics (i.e. crayfi sh distributions may be patchy due to demographic processes instead of differences in habitat quality). In our study, habitat selection by mink also was infl uenced by urbanization. Mink had lower occupancy probabilities in urban areas (Fig. 3). Th e ability of mammalian carnivores to adapt to urban areas is infl uenced by characteristics such as body size, reproductive potential, diet, behavior, and habitat requirements (Gehrt et al. 2010). Urbanization fragments natural habitats (McKinney 2002), which can cause the decline or local extinction of carnivore species (Crooks 2002). Roads and human development act as barriers to dispersal (Forman andAlexander 1998, Riley et al. 2006) and increase mortality risk from vehicle collisions (Tigas et al. 2002, Dickson et al. 2005). However, little information exists regarding the eff ects of urbanization on mink (Gehrt et al. 2010). In Canada, cottage development around lakes reduced habitat heterogeneity, decreased mink activity, and altered diets of mink (Racey and Euler 1983). Radio-marked American mink around Polish lakes also avoided areas near human settlements (Brzezi ń ski et al. 2012). Th ese studies were not set in urban areas, but they demonstrate the negative impact of human disturbance typical of urban areas on the behavior of mink, and agree with our result.
Contrary to our prediction, site occupancy for mink was related negatively to stream size (Fig. 3). We expected mink to select larger streams because we assumed that larger streams had more available resources, and we do not have an explanation for this surprising result. In the invasive range, Zuberogoitia et al. (2006) found the use of small versus large streams by radiomarked mink was sex-dependent; females used small streams and males used large streams. Unfortunately, we cannot assess sex-dependent habitat selection with our occupancy surveys.
We did not detect a strong, positive relationship between occupancy probability for mink and water depth in contrast to previous research in our study area (Schooley et al. 2012 Schooley et al. (2012) noted that colonization of vacant sites was variable if water depths were Յ 0.4 m, but consistently high if water depths were Ͼ 0.4 m. Seventy-one precent (42 of 59) of sites during our study had water depths below 0.4 m. Th us, we acknowledge that our results from this severe drought year may not extend to years of high precipitation. However, because climate models for the Midwestern United States predict an increase in the frequency of summer drought and spring fl ooding events (Wuebbles andHayhoe 2004, Trenberth 2011), our study may off er a glimpse into occupancy dynamics under climate change. Mink and other semiaquatic mammals are likely to experience increased environmental stochasticity that could create temporal variability in factors that infl uence habitat selection.
We measured a range of habitat variables typically used for assessing crayfi sh abundance (Riggert et al. 1999, DiStefano et al. 2003, but mink could have keyed in on an unmeasured habitat factor. For example, submerged vegetation provides protective cover for crayfi sh (Kershner and Lodge 1995) and reduces predation risk from aquatic and terrestrial predators (Wolff 2013). We did not measure submerged vegetation cover at crayfi sh hotspots in this study. However, we performed a post hoc analysis using submerged vegetation data collected between July and September 2012 at our study sites (Wolff 2013) to examine the relationship between submerged vegetation cover, crayfi sh hotspots, and site occupancy (Supplementary material Appendix 2). Crayfi sh hotspots were associated positively with submerged vegetation cover, but submerged vegetation did not explain much variation (R 2 Ͻ 0.17). Moreover, submerged vegetation cover was a poor predictor of site occupancy (worse than intercept-only model) compared to crayfi sh hotspots (Supplementary material Appendix 2 Table A1). Th us, we found little evidence that mink recognize areas with high cover of submerged vegetation as good locations to hunt.
Th e diet of mink in our study refl ected the seasonal availability of prey. Although crayfi sh and mammals were important year-round food sources for mink, crayfi sh occurred most frequently and in the highest volume in the summer diet. Low diversity, evenness and niche breadth values indicate greater dietary specialization on crayfi sh during summer. Th is increased contribution of crayfi sh to the diet appears to drive the habitat selection behavior of mink. Some species of crayfi sh move into deeper water and severely reduce activity at low water temperatures (Aiken 1968). Th is behavior may make crayfi sh less available to terrestrial and semiaquatic predators during winter. Conversely, fi sh increased in frequency and volume in the diet from fall to winter. Other diet studies of mink have observed increased fi sh consumption during winter as a result of increased vulnerability of fi sh (Gerell 1967, Magnusdottir et al. 2012. Low water temperatures decrease the mobility of fi sh, making them more vulnerable to predation (Parsons andSmiley 2003, Brown et al. 2011). Our study area experienced an unseasonably mild winter during 2012. January and February temperatures in Illinois were 3.67 ° C and 2.61 ° C warmer than average, respectively (Illinois State Water Survey 2012b, 2012c). Th us, many streams did not freeze over, but were still cold enough to decrease the mobility of fi sh, thus allowing mink constant access to vulnerable fi sh prey.
Our results have implications for the effi cacy of habitat models for species management. For instance, Loukmas and Halbrook (2001) concluded that the poor performance of a habitat suitability model for mink was primarily due to lack of habitat variables for key prey. If mink cue in on their prey directly, then improving model performance will be diffi cult if habitat variables are used as surrogates for prey abundance. However, quantifi cation of prey abundance can be diffi cult, time-consuming, and expensive, especially for a generalist predator that consumes a variety of prey items. Our approach of using seasonal diet data to identify key prey and direct resource measurement could be the most effi cient strategy for developing predictive habitat models for predators.
Our combined results for diet and occupancy modeling indicate the spatial distribution of mink during summer is greatly infl uenced by the abundance patterns of their preferred prey. Our data support recent evidence from Europe suggesting that availability of crayfi sh, and an increased prevalence of crayfi sh in the diet of American mink, may alter the spatial distribution of mink and aid in expansion of mink across its invasive range (Melero et al. 2014, Rodrigues et al. 2014. Further examination of this predator -prey interaction in North America could lead to a better understanding of mink ecology in places where the co-occurrence of mink and crayfi sh is novel. More generally, the importance of prey hotspots merits further investigation to understand the mechanisms underlying habitat selection by carnivores.