In the first comprehensive survey of post-breeding shorebirds conducted along the remote deltaic mudflat habitats on the coastline of the Arctic National Wildlife Refuge, six species of shorebirds were documented aggregating prior to southbound migration. Energy reserves gained while foraging in these areas may be critical for southbound migration of long-distance migrant shorebirds, but these habitats are vulnerable to potential effects of oil and gas development and climate change. The study objective was to assess the relative density and species composition of shorebirds . Surveys were conducted at 13 major river deltas on the coast between late July and mid-August each year from 2006 to 2010. Double-observer methods were used in 2010 to estimate the detection rate in surveys of randomly-selected transect sections. Shorebird density varied significantly between years and among river deltas. Peak relative density estimates at three deltas, the Jago (247.8 birds/km2), the Kongakut (100.6 birds/km2) and the Hulahula (49.5 birds/km2), were significantly higher than the estimate for the Canning (16.0 birds/km2). Because shorebird density and abundance vary significantly among sites and years, and individuals likely move among multiple sites within a given year, shorebird conservation strategies for these habitats should consider them to be spatially and temporally interconnected.
In arctic Alaska, little is known about shorebird use of coastal areas during the post-breeding period, prior to migration (Taylor et al. 2010). Preliminary work indicated that shorebirds depend on the food resources found in specific coastal areas to acquire fat necessary for southward migration (Connors et al. 1979). Shorebirds move from inland tundra breeding areas (Brown et al. 2007) to coastal mudflats and other littoral habitats as the breeding season progresses, and have been found in higher densities at coastal sites after the breeding season (Connors et al. 1979; Connors et al. 1981; Smith and Connors 1993). Post-breeding shorebirds using the Arctic Refuge coastline are likely sequestering fuel reserves critical for southbound migration (Connors 1984; Lyons and Haig 1995).
Many species of shorebirds are of conservation concern due to ongoing population declines (Brown et al. 2001; Morrison et al. 2001; International Wader Study Group 2003; Bart et al. 2007). Several species of particular conservation concern use stopover sites within the Arctic Refuge, including four species (Whimbrel (Numenius phaeopus), Bar-tailed Godwit (Limosa lapponica), Dunlin (Calidris alpina and Buff-breasted Sandpiper (Tryngites subruficollis)) on the USFWS Birds of Conservation Concern list for Bird Conservation Region 3, Arctic Plains and Mountains (U.S. Fish and Wildlife Service 2008), and 16 species listed by the Alaska Shorebird Conservation Plan as moderate or high concern (Alaska Shorebird Group 2008). Placement of birds on these lists is intended to stimulate collaborative proactive conservation actions among federal, state, and private partners. Our study implemented such a partnership, and will help develop measures to conserve shorebirds by increasing our understanding of their use of vulnerable habitats during a critical phase in their life cycle.
The primary goal of this study was to assess the relative density and species composition of shorebirds using coastal areas of the Arctic Refuge prior to fall migration, which is critical information so that impacts of future habitat changes can be measured and potentially managed or mitigated. Understanding the factors affecting the relative density of shorebirds at these coastal habitats is an important step in managing populations. There is particular concern for the coastal resources of the Arctic Refuge because these areas are vulnerable to potential effects of offshore oil development in the eastern Beaufort Sea and to changing sea and weather conditions associated with climate change (Kendall et al. 2011; Nolan et al. 2011). A major objective of the Alaska Shorebird Conservation Plan (Alaska Shorebird Group 2008) is to ensure that adequate quantity and quality of habitat is identified and maintained to support shorebirds. Understanding habitat use is a critical first step in meeting this objective.
The relative importance for shorebirds of foraging sites along the coastline of the Arctic National Wildlife Refuge, Alaska, has never been measured. Aggregations occur in mudflats and associated wetlands, but most of these areas have never been systematically surveyed for shorebirds using the same methodology within or between years. These coastal mudflats are extremely remote and difficult to access because they are isolated by large expanses of shallow water and barrier islands, so no ground surveys have previously been conducted for shorebirds at most of these sites.
Pre-migration sites have commonly been eferred to as staging areas, but Warnock (2010) suggested the use of this term be restricted to sites where birds prepare for long flights requiring large energy reserves. Since the length of stay and accumulation of fat were not determined in this study, and the migration strategies used likely vary among the many species using these sites, we refer to the mudflat habitats simply as foraging areas, and determine the relative abundance at the many sites within the study area.
METHODS
Study Area and Field Methods
Our study area included the entire coastline of the Arctic National Wildlife Refuge, from the Staines River on the western boundary of the Refuge, to Demarcation Bay on the border with Canada (Fig. 1). We surveyed the mudflats of all river deltas as they enter the lagoons behind the barrier islands. The mudflats are extensive areas of very low gradient alluvial deposits from northward flowing arctic rivers, almost completely bare of vegetation. We did not survey the salt marsh habitats along the coast between river deltas, nor the barrier islands that form the outer edge of the coast over much of the area.
Figure 1.
The study area along the coastline of the Arctic National Wildlife Refuge, Alaska, showing 13 major river deltas where shorebird surveys were conducted between 2006 and 2010.

We established transects at all river deltas along the coastline. We used satellite imagery (Digital Globe 2003) to identify river deltas with mudflats, and explored all deltas during the first year of the survey in 2006 to determine if mudflat habitat was actually present. We identified a total of 13 river deltas with mudflats, and surveyed all of them in each year from 2006 to 2010. No mudflat habitat was located at the Turner River or in the rest of Demarcation Bay, so these were excluded from future surveys.
Within each river delta, we randomly selected from one to seven starting points for 0.5 km transects depending on the size of the delta, oriented generally east-west. We selected one transect segment for each 1 km of river delta width measured east to west, and laid out waypoints every 100m. Transect starting points were random with respect to longitude, but were constrained to be approximately 100m from the edge of the water. During the surveys, we allowed the transect waypoints to vary north or south to maintain the 100m distance from the edge of the water because the water level varied due to a small lunar tide and larger effects of wind. The total length of the survey transect on each delta varied between years due to changing water levels, so we recorded the total survey distance each year to use in density calculations. Transects were often moved inland due to high water levels at the time of the surveys, shortening the unsurveyed distance between them. The survey width was truncated to 300 m on each side of the observer, and birds detected beyond this distance were excluded from the analysis. The area surveyed on each transect, therefore, was the survey length multiplied by the survey width (600 m), and included the very shallow water at the ocean edge of the mudflat and the moist mudflat habitat inland, which were both used by foraging birds.
Three river deltas were divided into multiple segments because they had separate outflows that were not crossable on foot and had to be surveyed separately. These included the Jago, Aichilik and Kongakut. The segments were surveyed separately, but combined for the overall analysis by river delta. The total number of transects at each delta were as follows: Achilik 4, Canning 2, Egaksrak 2, Hulahula 6, Jago 6, Katakturuk 2, Kongakut 7, Niguanak 1, Okpilak 3, Sadlerochit 2 and Staines 1. We surveyed all transects once in each of the five study years, from 2006 to 2010, and all surveys were well after the peak of hatching in early July. Survey periods were as follows: 8/10/06 to 8/21/06; 7/25/07 to 8/5/07; 7/20/08 to 8/6/08; 7/28/09 to 8/5/09 and 7/28/10 to 8/4/10.
We usually worked in pairs, with one person observing and a second recording data and waypoints. On some surveys where access was difficult or time consuming, surveyors worked individually and observations were made and recorded by one observer. We recorded the species and counted the number of birds observed singly or in groups foraging on the ground, and their perpendicular distance from the survey line. Birds that flew over the transect and did not land were recorded as fly-overs, but not included in the analysis of density. We recorded age as hatch-year (HY) or after hatch-year (AHY) if it was possible to determine from plumage.
Standard distance-sampling methods (Buckland et al. 2001) were initially considered, but were determined to be inappropriate for the study because birds were not distributed uniformly across the landscape with respect to the survey transects, and therefore did not conform to the assumption of declining detectability with increasing distance from the transect. Because the birds were in open habitat and sometimes concentrated along the waterline, and therefore relatively easy to count, we instead analyzed the data using strip transect methods (simply dividing our count of all birds in the 600m wide transect by the area surveyed).
To estimate detection rate within the entire strip, we used standard double observer methods in 2010 to calculate an overall rate for all surveys. We randomly selected two waypoints in each transect to conduct double observer counts on a 50 m long section of the transect. If there were no birds present at the first location, we used the second location. The primary observer used a spotting scope from a vantage point outside of the transect to count all the birds within the marked-off 50 m section. The secondary observer then surveyed the area with binoculars using the normal techniques, and similar survey effort to the overall survey. The primary observer continued to watch the strip during the secondary observer's survey, and reduced the count if birds left the area before the secondary observer reached their location, because these birds would not be available for the secondary observer to count. We calculated detection rate as birds detected on the normal survey by the secondary observer / birds detected on the perpendicular count by the primary observer. We expected these rates to be high because the habitat is open and the birds are clearly visible.
Table 1.
Model selection results for candidate models comparing ln(shorebird density) on Arctic Refuge coastal mudflats, 2006–2010, as a function of explanatory variables, including AIC values, models compared, and number of parameters in each model. Minimum AIC for this model set was 227.6.

Statistical Analysis
We used a one-way analysis of covariance (ANCOVA) to analyze differences in shorebird densities among river deltas. We transformed the response (shorebird density expressed as birds/km2), by taking the natural logarithm (In), to improve normality. Because our primary interest was whether shorebird densities differed among deltas, our one-way ANCOVA had delta as the factor effect and included three covariates known to affect shorebird counts: Year, Julian Date and Julian Date2. We included the squared term to account for nonlinearity that occurs during shorebird occupancy, with numbers peaking and then falling off as the season progresses. When the second-order polynomial was included in the model, we made certain the lowerorder term “Julian Date” was also included (Kutner et al. 2005). We considered several models (Table 1) and selected the “best” model using Akaike's Information Criterion (AIC; (Akaike 1974; Burnham and Anderson 2002)) where the model with the smallest AIC value and largest model weight was considered the best model. Our base model included only the “River Delta” factor and the full model included River Delta, Julian Date, Julian Date2, Year and an interaction between Year and Julian Date2. Standard diagnostics for non-constant variance, normality of error terms and evaluation of outliers were conducted on the final model.
Our “best” model was used to estimate In (density) for each delta adjusted for covariates (i.e. the factorlevel means from the ANCOVA). These least-square means were computed by setting all covariate effects to their mean values. We evaluated confidence intervals for density estimates and used a Tukey-Kramer multiple comparison procedure to examine among-delta differences in shorebird densities (SAS Institute Inc. 2011). We back-transformed In (density estimates) for each river delta to aid in interpretation.
RESULTS
We documented that six species of shorebirds aggregate on the coastal plain of the Arctic National Wildlife Refuge after their breeding season and prior to fall migration, and another 13 species occur in smaller numbers. The number of shorebirds detected on transects at each river delta varied widely among years, as did the total number of birds observed across all deltas combined. We observed a total of 629 shorebirds in 2006; 4,469 in 2007; 8,984 in 2008; 5,277 in 2009 and 5,556 in 2010. The highest densities occurred in 2006 at the Katakturuk delta, in 2007 at the Kongakut delta, in 2008 at the Jago and Hulahula deltas, and in 2009 and 2010 at the Jago delta.
The most common species was Semipalmated Sandpiper, which accounted for more than 80% of the individuals observed in all years except 2006 (43%), and 83% of individuals detected across all years. Most Semipalmated Sandpipers observed were juveniles (78% of birds with age identified). The next most common species were Red-necked Phalarope (Phalaropus lobatus, 6%), Black-bellied Plover (Pluvialis squatarola, 4%), Dunlin (2%), Stilt Sandpiper (Calidris himantopus, 2%) and Pectoral Sandpiper (Calidris melanotos, 2%). All of these species except for Pectoral Sandpiper are listed as either moderate or high concern in the Alaska Shorebird Conservation Plan. All other species observed in our surveys each accounted for less than 1% of the total individuals recorded. These rarer species, listed in decreasing order of abundance, included: Sanderling (Calidris alba), American Golden-Plover (Pluvialis dominica), Ruddy Turnstone (Arenaria interpres), Long-billed Dowitcher (Limnodromus scolopaceus), Baird's Sandpiper (Calidris bairdii), Western Sandpiper (Calidris mauri), Semipalmated Plover (Charadrius semipalmatus), Buff-breasted Sandpiper, Red Phalarope (Phalaropus fulicarius), Bar-tailed Godwit, White-rumped Sandpiper (Calidris fuscicollis), Least Sandpiper (Calidris minutilla) and Whimbrel.
We conducted 47 comparisons of primary and secondary observer counts along randomly- selected 50m transect sections. Our mean detection rate for all species combined was 0.986 (SE 0.021). Because the detection rate was close to 1, we did not adjust our counts for detectability.
The top-ranked model for shorebird density chosen via AIC included a main effect of River Delta, and Year, Julian Date, and Julian Date2 as covariates (Table 1). While the best model included River Delta (F18,198, P = 0.0041, Table 2) and also a year effect, individual parameter estimates for each year were not statistically different from zero. We tested for differences among pairs of river deltas using Tukey-Kramer multiple comparisons and found three significant comparisons (Jago vs. Canning, P = 0.0006, Hulahula vs. Canning, P = 0.0271, Kongakut vs. Canning, P = 0.0452; Fig. 2). The higher densities in the Hulahula, Jago and Kongakut compared with the Canning were also supported by non-overlapping confidence intervals. Density estimates for each delta were back-transformed from the log scale for ease of interpretation (Table 3).
DISCUSSION
Our results show that relative densities of shorebirds differ significantly among river deltas on the Arctic Refuge coastline. The Jago River delta had the highest average density of any river delta, and the Jago, Kongakut and Hulahula deltas had significantly higher densities than the Canning River delta. These comparisons were likely significant because of the relatively low density at the Canning compared to all other deltas. The low density at the Canning may be a result of low availability of food resources at this site, and we are exploring the densities of invertebrates in related work that is currently underway. Some of the very small deltas, like the Staines, Angun and Sadlerochit, had relatively high densities of birds, although not significantly higher relative densities compared to other deltas, and small total numbers of birds because of their small size. Conversely, some other large deltas, like the Okpilak, Aichilik and Egaksrak, did not have significantly higher densities than other deltas, but still had high total numbers of birds because of their large area.
Table 2.
Parameter estimates from ANCOVA models used to estimate ln(shorebird density) from data collected in the Arctic National Wildlife Refuge, Alaska (2006–2010). River Delta was modeled as the factor effect in the ANCOVA model and Year, Julian Date and Julian Date2 were modeled as covariates.

While some deltas had significantly higher densities than others, all of the river delta mudflats supported shorebirds, and there was high variability in maximum counts among sites and years at all sites. Taylor et al. (2011) found that shorebirds on the North Slope did not necessarily move across the landscape systematically in the direction of their normal fall migration, thus it seems likely that birds may be using different sites as resource availability or environmental conditions change. Taylor et al. (2011) also showed that shorebirds move widely among north slope coastal sites within years, and found that Semipalmated Sandpipers moved more frequently than several other species.
Figure 2.
In(density estimates) and associated 95% confidence intervals for river deltas on the coast of the Arctic National Wildlife Refuge, Alaska, from surveys conducted between 2006 and 2010. Deltas with different letters had significantly different density estimates in the Tukey-Kramer multiple comparisons test.

Table 3.
Relative shorebird density estimates (least-squares means back-transformed from the log scale) for single annual surveys at coastal mudflats in 13 river deltas along the Arctic National Wildlife Refuge coastline, 2006–2010.

We cannot compare variability at sites within years in our data set, because the extreme remoteness of these sites and the time needed to travel between them limited us to only a single visit to each delta within each study year. However, the large variability within sites among years suggests that deltas are more or less valuable for foraging shorebirds at different times, and may be part of a complex that supports birds moving among sites as habitat quality varies. Although the overall fit of the model improved by including a year effect, the large standard errors of the year parameter estimates indicate it is not as useful for estimating relative shorebird densities compared with the other variables included in our best model.
The density values presented here are conservative. The model estimates densities averaged for the date of the single survey, and future work should model stopover time and estimate total numbers of birds using each site during the entire season (Frederiksen et al. 2001; Cohen et al. 2009). However, these density values are the first data ever presented that allow detailed comparisons among river deltas within the study area. Comparing the deltas to each other is an important first step toward development of management plans for these deltas, because this comparison helps determine the relative importance of each site within the Arctic Refuge. Additional work currently underway will help determine if the abundance patterns reported here persist with repeated measurements throughout the migration season.
There are very few other studies that measure densities of foraging shorebirds on the Arctic Refuge coastline or at nearby sites. The high abundance of Semipalmated Sandpipers on coastal mudflats in our surveys corresponds with results reported by Smith and Connors (1993) who found that the species most commonly used mudflat habitats near Barrow. Semipalmated Sandpipers are thought to have declined (Hitchcock and Gratto-Trevor 1997) as measured by significantly reduced counts on the wintering grounds (Andres et al. 2012; Morrison et al. 2012), so identification of key migration foraging areas is important. Andres (1994) reported overall densities for all shorebirds of 149.9/km2 in repeated surveys on the Colville River Delta, higher than all of the sites reported here for the Arctic Refuge except for the Jago. The predominant species at the Colville was Dunlin, which may be underrepresented in our sample due to survey timing being earlier than the peak of Dunlin migration. Dunlin are known to migrate later in the season (Warnock and Gill 1996), and preliminary data from our ongoing research at the sites with repeated measurements shows a late pulse of Dunlin moving through the area that we likely missed due to the timing of our coastal surveys. Dunlin were the most abundant species in the surveys conducted later in the season on the Colville river delta (Andres 1994), and adding Dunlin missed in our surveys as well as conducting repeated counts at each site would likely make the densities more comparable.
Taylor et al. (2010) found the same general pattern of high single counts at the Jago and the Kongakut river deltas in their surveys, and also found highly variable counts between years. Their results were based on different methods, including aerial surveys and counts that were not corrected for detectability, so only general comparisons are possible. Our maximum raw count at the Kongakut in 2007 was over two times the maximum count reported in their surveys, and our Jago count in 2008 was over nine times their maximum. Comparing our maximum counts to others from the North Slope in their much wider surveys, the Jago and the Kongakut would rank as the largest and second-largest concentration areas, followed by the Sagavanirktok in the Central Beaufort and the Ikpikpuk in the West Beaufort, which suggests that these Arctic Refuge foraging areas attract concentrations of birds that are large on the scale of the entire North Slope.
Shorebirds are known to depend on stopover sites during migration to regain depleted fat reserves, and survival has been shown to be correlated with body condition before long migratory flights (Wilson 1990; Baker et al. 2004; Atkinson et al. 2007; Morrison et al. 2007). Understanding the relative densities of birds using coastal sites to prepare for these southbound migrations will aid in planning conservation priorities in the face of changing environmental conditions. A major management concern at the coastal mudflat sites in this study is that reduced habitat quality could impact shorebird populations. Potential impacts could arise through effects of human-caused or natural events, including disturbance from increasing oil and gas development (National Academy of Sciences 2003), which could cause displacement of foraging shorebirds, or impacts from spills.
The habitats surveyed in this study are also likely vulnerable to effects of a changing climate. Climate change has been more pronounced in the Arctic than at lower latitudes, and the Arctic has experienced twice the warming of the global average (McCarty 2001; ACIA 2005; Hinzman et al. 2005; Serreze and Francis 2006; IPCC 2007). Warming has been particularly pronounced in northern Alaska (Martin et al. 2009). Climate change could affect shorebird habitats on the northern coast of Alaska in many ways, including trophic mismatch with invertebrate food sources (Visser and Both 2005; Tulp and Schekkerman 2008; Van Der Jeugd et al. 2009), habitat loss or shifts (Maclean et al. 2008), and extreme weather events or changes in inland hydrology or sea level (ACIA 2005; Martin et al. 2009). Determining the relative abundance of birds at coastal sites will help support conservation efforts to mitigate these impacts.
Sea level rise is predicted to accelerate, and could have significant impacts on foraging areas for shorebirds on the Arctic coastline. Similar concerns were reported by Galbraith et al. (2002) for shorebird stopover sites in temperate zones. Accessibility of appropriate habitats, which is controlled by water depth, is critical to their value for foraging shorebirds (Collazo et al. 2002). Because arctic coastal mudflats have very low gradients, small changes in sea level from increased ice melting and thermal expansion (IPCC 2007) can inundate large areas, making them accessible for shorter periods of time or inaccessible during the migration window. In some storms during our surveys, large mudflats like the Jago river delta were completely inundated, and shorebirds were displaced from a significant foraging area for approximately eight days during the period when they would normally be preparing for migration. Increases in sea level or increases in the severity or frequency of storms could have negative impacts on the ability of many species to utilize critical foraging habitats.
Our results document the relative density of shorebirds across a very large area of coastal mudflat habitat at a critical stage in their life cycle: their preparation for southbound migration. In general, our results suggest two overall conclusions relevant to management of coastal habitats. First, the large sites with consistently high average densities may attract a larger proportion of the post-breeding shorebird population than sites with lower densities, and should be managed to protect shorebird foraging habitat. Future work should be designed to confirm whether the total abundances of birds using each site follow the patterns found here for single counts. Second, because all sites were used heavily in at least one year of the study, management of the Arctic Refuge coastal areas should consider these sites as part of a habitat complex within which each site may be important at various times in preparing shorebirds for southbound migration. Taylor et al. (2010) suggested the same approach for the broader coastal environment of the northern Alaska coast, of which our study area was a small part. Managing these sites as part of an interconnected resource for foraging shorebirds will improve their ability to support shorebird abundance and diversity at a critical time of year, and provide a more resilient habitat in the face of predicted changes in climate and sea level.
ACKNOWLEDGMENTS
We are grateful to the ornithologists who participated in these surveys. In addition to the authors, surveys were conducted by R. Burner, H. Craig, L. DeCicco and B. Winn. J. Reynolds helped design the surveys. We acknowledge the Kaktovik Inupiat Corporation for permission to conduct surveys on their lands. We staged our survey out of Kaktovik, AK, and appreciate the support of the people of Kaktovik. Major funding for this study was provided by donors of Manomet Center for Conservation Sciences, by the U.S. Fish and Wildlife Service Challenge Cost Share Program, and by the Arctic National Wildlife Refuge.