Next Article in Journal
Selenium Nanoparticles as Potential Antioxidants to Improve Semen Quality in Boars
Next Article in Special Issue
To Flee or Not to Flee: How Age, Reproductive Phase, and Mate Presence Affect White Stork Flight Decisions
Previous Article in Journal
Sheep Face Detection Based on an Improved RetinaFace Algorithm
Previous Article in Special Issue
Variations in the Reproductive Strategies of Different Charadrius alexandrinus Populations in Xinjiang, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effects of Climate Change on the Nesting Phenology of Three Shorebird Species in the United States

1
Department of Biological Sciences, Clemson University, Clemson, SC 29631, USA
2
GAI Consultants, Homestead, PA 15120, USA
3
Department of Geology and Environmental Geosciences, University of Charleston, SC at the College of Charleston, Charleston, SC 29424, USA
4
Graduate School of Education & Human Development, George Washington University; Washington, DC 20052, USA
5
School of Medicine Greenville, University of South Carolina, Greenville, SC 29605, USA
6
East Cooper OB/GYN, Mount Pleasant, SC 29464, USA
*
Author to whom correspondence should be addressed.
Animals 2023, 13(15), 2459; https://doi.org/10.3390/ani13152459
Submission received: 19 June 2023 / Revised: 20 July 2023 / Accepted: 27 July 2023 / Published: 29 July 2023

Abstract

:

Simple Summary

Previous studies have shown how climate change can lead some bird species to alter the timing of their breeding and migration. To understand this better in shorebirds, we investigated if the nesting times of three shorebird species in the United States have fluctuated due to increases in temperature and precipitation at their breeding sites. We found that the laying dates of our species have become earlier and temperatures at our nest sites have significantly increased over time. Precipitation did not show a clear trend over time in our study, though previous studies indicate precipitation has been increasing in the U.S. While we were unable to find evidence that precipitation affects when our species begin laying eggs, one species, the Willet, had significantly earlier laying dates as temperatures during the breeding months increased. While another species showed this same trend, it was not significant, and the third species showed later laying dates as temperatures increased. Our results demonstrate that changes in climate affect species differently, warranting further study in this field.

Abstract

Previous research suggests that a frequent response of organisms to the ongoing climate crisis is the adjustment of their reproductive timing or breeding phenology. Shorebirds may be especially vulnerable to increasing temperatures and precipitation, as many are migratory and depend on coastal habitats for wintering and breeding. These particular habitats could be at risk due to changes in climate, and nesting times often depend on food availability, which is often directly influenced by temperature. We investigated if clutch initiation dates (CID) for three shorebird species in the United States have become earlier over time with increasing temperatures and precipitation. We used nest records from Cornell’s NestWatch program and various museum databases and weather station data from the National Oceanic and Atmospheric Administration. We found evidence that CIDs have become earlier over time, though this was only a significant factor for one species. While temperature in our study areas has increased significantly over time, precipitation changes were more variable and not always significantly predicted by time. We found evidence that one species may be responding to increasing temperatures by nesting earlier, but there was no support for our hypothesis that CID has changed due to changes in precipitation for any species. Results varied for each species, indicating the importance of further studies on shorebirds as the effects of climate change on their nesting phenology may not be fully realized and will likely depend on the species’ biology and distribution.

1. Introduction

Increasing temperatures across the globe [1] have led to observed changes in the environment that can affect migratory bird species, such as the timing of snow melt in the Arctic [2,3], the emergence of spring vegetation [4], and essential food sources for breeding adults and their young [5,6]. Migratory species may be more vulnerable to these types of changes than resident species because they rely on different wintering and breeding habitats, as well as stopover sites during migration. While some migratory bird species might be able to adjust their timing of migration and nesting based on these recent environmental changes [4,7], others may lag behind or arrive too early to the breeding grounds before conditions are optimal for nesting [8,9,10]. Any shifts or mismatches in nesting phenology with respect to the availability of critical resources, such as nesting habitat and food, could lead to decreased reproductive success and population decline for these species [3,9,11,12].
Over the last few decades, several studies have found that some bird species have begun nesting earlier (e.g., [13,14,15,16,17]), which might be a response to increasing temperatures and the earlier emergence of spring vegetation and food sources, as was shown for the Willet (Tringa semipalmata semipalmata) by Smith et al. [4], and for several passerine species in the U.S. by Mayor et al. [10]. In contrast, Kwon et al. [11] found that three shorebird species that breed in the Arctic have delayed clutch initiation and shortened their incubation period over two decades, potentially in response to a cooling climate trend at the start of their breeding season. Kwon et al. [3] found evidence that phenological mismatches in Arctic shorebirds were likely related to changes in temperature and timing of snowmelt, which may be one cause that has led to population decline in the species nesting in the eastern part of the breeding range. These studies suggest that changes in global temperatures can either lead to changes in the timing of reproductive behaviors in birds or increase phenological mismatches, which may have negative consequences on their populations.
Increasing temperatures can also influence other aspects of climate, such as precipitation rates. As surface temperatures increase due to increasing greenhouse gas emissions, the rate of evaporation of surface moisture increases, which can lead to greater amounts of rainfall as reviewed by [18]. Some sources indicate that the precipitation rate in the U.S. has been increasing since the 1970s [19,20]. Further research predicts rainfall events will increase dramatically in intensity yet occur at a less frequent rate due to temperature increases exceeding evaporation rates [18,21]. These changes in precipitation could have direct effects on the habitat and reproductive success of certain bird species. For example, shorebirds nesting along coastlines could be subjected to more flooding and extreme weather events, such as tropical storms or hurricanes, as well as a loss of habitat due to sea level rise [22,23,24,25,26].
In this study, we investigated the hypothesis that rising temperatures and levels of precipitation have led to earlier laying dates in migratory bird species. Specifically, we focused on three shorebird species that breed in the United States of America: the willet (family Scolopacidae, hereafter “willet”), Wilson’s plover (Charadrius wilsonia, family Charadriidae, hereafter “plover”), and the black-necked stilt (Himantopus mexicanus, family Recurvirostridae, hereafter “stilt”). These species were chosen because they all may be especially vulnerable to a warming climate [22], and none are closely related, reducing phylogenetic effects on our results. All species in our study have been found to nest in or close to vegetation [4,27,28], the emergence of which could be directly affected by spring temperatures. These species also often nest along coastlines, putting them at a higher risk of exposure to extreme weather events and changes in habitat due to factors caused by climate change [22,23,24,25,26,29]. Shorebirds also have a fixed clutch size and typically produce only one brood a year [30,31].
We used laying dates as a measure of nest phenology change and asked three primary questions: (1) Have laying dates changed over year, latitude, or longitude? (2) Has temperature and precipitation changed over year, latitude, or longitude in our study areas? (3) Does temperature or precipitation predict laying date? We predicted that laying dates will have become earlier over time due to overall increases in surface temperature [1] and increased amounts of precipitation in the U.S. [18].

2. Materials and Methods

2.1. Study Species and Sites

In this study, all stilt nests were located in California, while all plover and willet nests were located along the east coast in the following states: Florida, Georgia, South Carolina, North Carolina, and Virginia (Figure 1). California is incredibly variable in climate and habitats [32]. The stilt nests were found in several types of ecoregions within this state, the top three being the Central California Valley, the Sonoran Basin and Range, and the Southern California/Northern Baja Coast (see [32] and Figure 1A). The Central California Valley is largely flat, agricultural land characterized by long, hot summers with little rainfall [32]. Likewise, the Sonoran Basin and Range is a hot, dry, desert habitat, while the Southern California/Northern Baja Coast is characterized by coastal sage scrub and chaparral with both low-lying areas by the coast and small hills. Many of the nests in the eastern states were found within a region called the North American Coastal Plain (NACP) (see [33] and Figure 1B,C). This region borders the Atlantic Ocean and consists of habitats such as closed forests, pine savannas, grasslands, and wetlands, including swamps, marshes, and bogs [29,33]. Several factors influenced by climate change, including rising sea level, increased episodes of flooding or drought, hurricanes that are more intense, and increase levels of saline in ground and surface waters, are threatening the habitats in this region [29].
Stilts have an incubation period of around 25–26 days and lay an average of four eggs per clutch over 4–5 days [34]. In our study, we assumed that stilts laid one egg each day with an incubation period of 25 days based on data collected in the Tulare Basin (average and median incubation period was 25 days, n = 43) [34], which is geographically close to the majority of the nests used in this study (Figure 1A; Supplementary Table S1). The timing of incubation for stilts appears to vary depending on ambient temperatures, but we assumed that stilts began incubating after laying the final egg. Plovers typically lay three eggs per clutch over 4–6 days. We assumed an average of five days for the laying period for a clutch of three eggs [35]. Regular incubation begins after the last egg is laid and also appears to be around 25 days, though this is based on limited observations. Willets typically lay four eggs per clutch, usually over six days. We assumed incubation began after the last egg was laid, but this has not been well-studied [36]. The incubation period is around 25–26 days, with an average of 25.2 days or a median of 25.5 days. Based on this, we assumed an incubation period of 26 days (Table 1 summarizes the metrics used to help us estimate laying date).

2.2. Estimating Clutch Initiation Dates

To determine if the laying date or clutch initiation date (hereafter referred to as “CID”) for the willet, plover, and stilt have changed over time and with changing climate, we collected nesting information for these species from 1852 to 1983 using several museum databases, as well as nest observation cards from 1966 to 1989 provided by the Cornell NestWatch program, giving us information on a total of 994 nests (Table 2 and Figure 2A). The database website Arctos (https://arctosdb.org, accessed in October 2021) was used to retrieve records from some smaller museums (Table 2). Once CIDs were estimated for all clutches, they were converted into a Julian calendar date for use in the statistical analyses (Figure 2B and Supplementary Table S1).
From museum records, we were usually able to estimate CID using the clutch size and the collector’s description of incubation stage at the time the clutch was collected, following similar methods to McNair [43] (see also [44]). If a collector wrote down a specific number of days or a range of days in which they thought the eggs had been incubated upon collection, we used that number (or the median of the range rounded up to the nearest whole number) to estimate the length of incubation time at collection. However, many collectors described the incubation stage in categorical terms (e.g., about 1/2 incubated, advanced, etc.). Different collectors also used different terms that likely had similar meanings. For consistency, we combined similar incubation stage terms collectors used into our own incubation stage categories (e.g., fresh, slight incubation, 25% incubation, etc.) that were then translated into a specific number of days based on the species’ known incubation period (Table 1 and Table 3). Clutch size was used to determine the laying period based on what is currently known about each species’ laying behavior (Table 1). We assumed stilts laid one egg per day, so their laying period equaled their clutch size (for a clutch of four eggs, the laying period would be four days) [34]. We assumed the plover and willet laying periods were always the clutch size plus two days (for a clutch of four eggs, the laying period would be six days) [35,36]. We assumed all birds began incubating on the final day of the laying period. To account for this, we added the laying period-1 to the estimated incubation period. This total was the number of days we counted backward from the date the clutch was collected to reach an estimated CID (Supplementary Table S1) (see also [43]).
If a clutch size was two or more less than a typical full clutch, and if no incubation stage term was given by the collector (or if incubation stage was listed as “unknown”), we assumed the bird was still in the laying period, so the incubation period was 0. For example, according to Robinson et al. [34], a typical clutch size for the stilt is four. Any nests with 1–2 eggs listed as the clutch size were treated as if they were still in the laying period when they were collected if no incubation stage was recorded. Thus, a stilt nest with one egg as the clutch size collected on 4 May was given an estimated CID of 4 May. A stilt nest with two eggs as the clutch size collected on 4 May was given an estimated CID of 3 May. If the incubation stage was not recorded or labeled “unknown” and the clutch appeared to be completed upon collection, we did not attempt to estimate CID and these nests were not included in analyses where an estimation of CID was required (Supplementary Table S1).
NestWatch nest observation cards had to be treated differently from the museum records, as most of these included observations of an active nest, and observers typically did not collect and remove egg contents for long-term storage. Therefore, they could not determine what incubation stage the nest was in upon discovery simply by looking at the eggs. We used the observer’s description of the bird’s behavior and changes in nest contents over a few days to estimate CID. For example, if an observer recorded one egg on day 1 for a stilt nest, and two eggs on day 2, we assumed the nest was found during the laying period and day 1 was the CID. If an observer said the egg was pipping, but the actual hatching date was not recorded, we assumed hatching occurred the following day and counted backward from that date based on the typical incubation time and laying period for that species (incorporating clutch size as stated before). If the observer recorded that chicks were still wet, we assumed hatching occurred that day. If hatching occurred between visits within a 1–5 day range, we always assumed it occurred on the median day of the range (always rounding up to the nearest whole number). If the observer estimated the age of the young to be a range of days, we always used the median age of the range, rounding up to the nearest whole number. We did not include any clutches for the following situations:
  • Only one visit was recorded, making it impossible to know the true stage of the nest;
  • The nest was only found during incubation, with no indication of when the eggs were laid or when they started hatching;
  • Hatching occurred between visits with more than a 5-day range with no age estimate of the young;
  • Only young were observed, but there was no clear age estimate, or the number of young seen was less than the typical clutch size for that species;
  • One plover nest was not used because it may have been a later nesting attempt of the same pair that had a failed nest earlier that year in that area.

2.3. Nest Location and Climate Data

If latitude and longitude were provided by the collector or observer, these were considered to be the true coordinates of the nest. If only a location description was provided, we estimated latitude and longitude coordinates from the most specific locality recorded using the website epsg.io from MapTiler Team [45]. The coordinates were uploaded as a layer into ArcGIS Pro 2.9 by ESRI (Figure 1). We downloaded a list of weather stations, along with their location coordinates, for each state within our study from the National Oceanic and Atmospheric Administration (NOAA), which was also uploaded as a layer into ArcGIS. In ArcGIS, we used a tool called “Generate Near Table” to produce a list of the ten closest weather stations to each clutch. The global monthly climate summaries for each of these stations from 1851–2022 were downloaded from NOAA. From this, we used the average precipitation (cm) and temperature (Celsius) during the first four months of the breeding season for our species (March–June). We excluded July as only stilts initiated nests during the first-half of this month (Figure 2B). Unfortunately, not all stations had both precipitation and temperature data for all years the station was active, and most stations were only active a portion of the years from which we sampled data. In total, we retrieved climate data for 671 weather stations with 239 from California close to the stilt nests and 432 in the east (Florida, Georgia, North Carolina, South Carolina, Virginia) close to the plover and willet nests (Supplementary File S1).

2.4. Have CIDs Changed Based on Year, Latitude or Longitude for Our Three Species?

To answer this question, we ran a full linear regression model (LM) to see if CIDs for all species combined have changed over time or space. Our predictor variables included the year of the nest, nest latitude, nest longitude, species, and all interactions between species and the other predictor variables. We used a total of 733 nests in this analysis where CID and nest coordinates were able to be reasonably estimated or had been provided (stilt = 339, willet = 198, plover = 196). A second analysis was run with each species separately, allowing us to understand which factors were the most important for each species. Similar methods were used by Torti and Dunn [16].

2.5. Has Temperature and Precipitation Changed over Year, Latitude or Longitude in Our Study Areas?

In this analysis, we treated each weather station as a separate data point and controlled for station ID since most had been sampled multiple times over many years (Supplementary File S1). Torti and Dunn [16] separated their temperature data into four regions based on latitude (latitudes greater than vs. less than 40 degrees) and longitude (longitudes greater than vs. less than 110 degrees). In our study, all stilt nests had longitudes greater than 110 degrees, while all plover and willet nests had longitudes less than 110 degrees (Figure 1). Additionally, all nest latitudes were below 40 degrees. Based on this, we separated temperature and precipitation data into eastern and western regions and ran mixed linear models (MLM) with station ID as a random effect and year, station latitude, and station longitude as fixed effects. The response variable was either east temperature, west temperature, east precipitation, or west precipitation for each month of the main part of the breeding season (March–June). In total, we ran 16 MLMs.

2.6. Does Temperature or Precipitation Predict CID?

For this analysis, we used similar methods described by Torti and Dunn [16]. We estimated the temperature and precipitation at a nest for each of the first four months of the breeding season by averaging the average monthly temperature and precipitation of the ten closest weather stations to a given nest that were active during the same year the nest was active. However, in many cases, there may have been fewer than ten stations (or no stations) active during a particular year a nest was active. Therefore, our sample size for this analysis was much smaller and, in some years, may have only included data from a single weather station for a particular nest (Supplementary File S2). Next, we determined which average monthly temperature and precipitation best predicted CID for each of our three species [16]. We performed LM analyses for each species separately, using CID as our response variable and average monthly temperature or precipitation as the predictor variable for March, April, May, and June, giving us a total of four LMs for temperature for each species and four LMs for precipitation for each species. The best model was considered to be the one with the highest R2.
For all analyses, any interactions or individual predictor variables that were not significant (p > 0.05) were removed in a stepwise fashion (with the most insignificant factors being removed first) from the final model that is reported. When comparing models in the first and third analyses, we used either R2 (when only one predictor variable was included) or adjusted R2 (when more than one predictor variable was included) to determine which model was best. We considered the best model to be the one with the highest R2 or adjusted R2. In the second analysis, we report the models with the lowest Akaike Information Criterion (AIC) score because in several of these models, the adjusted R2 was minimally lower when some significant effects were removed, making it more difficult to interpret. In some cases, the best model may have included both significant and insignificant effects, though this was rare. All statistical tests in this study were performed in JMP Pro 16.0.0 with an alpha of 0.05.

3. Results

According to our estimated CIDs, willets began nesting first, with the earliest CID being on 24 March (83 Julian date) and the latest being on 23 June (174 Julian date) (Figure 2B). The average willet CID was 5 May (125.3 Julian date). Plovers began nesting as early as 3 April (93 Julian date), with the latest nest initiated on 21 June (172 Julian date) and the average CID on 12 May (132.5 Julian date). Stilt CIDs ranged from 12 April (102 Julian date) to 11 July (192 Julian date), with an average date of 22 May (142.5 Julian date).

3.1. Have CIDs Changed Based on Year, Latitude or Longitude for Our Three Species?

In the model where all species were combined, all predictor variables were significant, as was the interaction between species and nest longitude (LM, adjusted R2 = 0.24, F7, 725 = 33.9, p < 0.001; Table 4). Stilts began nesting later than plovers and willets and were the only species in our study to start nests in July (Figure 2B). Though the majority of nests for all species were started in May, stilts had significantly later CIDs, on average, than the other two species (LM: Tukey HSD, p < 0.05; Figure 3). All species showed later CIDs as latitude increased, likely due to the later arrival of spring and warmer temperatures in more northern latitudes (Latitude: estimate = 1.75, p = 0.008). Species showed different trends in their CIDs with relation to longitude, which was explored more in the second analysis where separate LMs were run for each species. Overall, species’ laying dates have become earlier over time (Year: estimate = −0.09, p = 0.003), but year was not the most important predictor of CID for all species (Figure 4).
In the LM where only stilt nests were included (n = 339), the only significant predictor of stilt CID was longitude (LM: adjusted R2 = 0.04, F1, 337 = 16.5, p < 0.001; Longitude: estimate = 2.26, p < 0.001). The nests that were further inland (higher longitude values) had later CIDs than nests closer to the coast. However, there was a lot of variation in this dataset, with many nests being found at the same longitudes but showing a wide range of CIDs (Figure 4A). Plover CIDs (n = 196) were only significantly affected by latitude (LM: adjusted R2 = 0.10, F1, 194 = 22.01, p < 0.001; Latitude: estimate = 2.22, p < 0.001; Figure 4B), becoming later at more northern latitudes. Willet CIDs (n = 198) were significantly affected by year, becoming earlier over time (LM: adjusted R2 = 0.28, F2, 195 = 39.9, p < 0.001; Year: estimate = −0.13, p = 0.002), and latitude, becoming later at more northern latitudes (Latitude: estimate = 3.03, p < 0.001; Figure 4C).

3.2. Has Temperature and Precipitation Changed over Year, Latitude, or Longitude in Our Study Areas?

For this analysis, we only report results from the best models with the lowest AIC (Table 5). In our study areas, average monthly temperature has increased significantly over time for all four months (March-June) in both the west (California) and east, and these models all showed high adjusted R2 values (MLM: adjusted R2 = 0.69–0.92, p < 0.03 for all models; Table 5(A)). In the east, and in the west during March, temperature significantly decreased as latitude increased (colder temperatures occurred at more northern latitudes; p < 0.001 for all models). In June in the west, temperature increased significantly at more northern latitudes (p = 0.008). When longitude had a significant effect on temperature, we found for both the east and the west, temperature significantly increased further from the coast (longitude and temperature showed a positive relationship in the west but a negative relationship in the east, p < 0.002 for all models). In our dataset, longitude values are negative because all are in the western hemisphere. Therefore, larger values in the west indicate more inland locations, while larger values in the east indicate locations on the coast (Figure 1).
All precipitation models showed much lower adjusted R2 values (0.14–0.28 for the west and 0.07–0.20 for the east; Table 5(B)), which may indicate precipitation is more variable than temperature or may not be as predictable using a linear model over time or space. In the west, the year never had a significant effect on precipitation. In the east, year was always a significant effect (p < 0.001); however, the relationship between year and precipitation differed depending on the month. The best model for predicting precipitation changes in the west was for March (MLM: adjusted R2 = 0.28), which showed a significant decrease in precipitation at more northern latitudes (p = 0.006) and more inland longitudes (p < 0.001), but this relationship was reversed in May and June. The best model in the east was for June (MLM: adjusted R2 = 0.20), which showed an increase in precipitation over time (p < 0.001) and a decrease in precipitation at more northern latitudes (p < 0.001).

3.3. Question 3: Does Temperature or Precipitation Predict CID?

Each species had different results for the temperature analysis (Table 6). For stilts, CIDs always became later as temperature increased (positive relationship), which was contrary to our prediction. The average June temperature best predicted CID for stilts, and this was a significant effect (LM: R2 = 0.04, F1, 219 = 8.52, estimate = 0.86, p = 0.004; Figure 5A). The average May temperature was also a significant predictor of stilt CID (p = 0.03). For willets, CIDs always became earlier as temperature increased (negative relationship), and temperature was always a significant effect for each month (p < 0.005 for all models), as predicted (Table 6(A)). The average May temperature best predicted CID for willets (LM: R2 = 0.15, F1, 110 = 20.11, estimate = −4.07, p < 0.001; Figure 5B). Temperature was never a significant predictor of plover CID for any month (Table 6(A)), but the average April temperature was the best predictor (LM: R2 = 0.02, F1, 77 = 1.83, estimate = −1.17, p = 0.18; Figure 5C). Plover CIDs became earlier with increasing temperature in March and April, but later with increasing temperature in May and June.
Precipitation never significantly predicted CID for any species for any month of the breeding period (LM: R2 = 0–0.03 and p > 0.07 for all models; Table 6(B)). Still, out of the four breeding months, March precipitation was the best predictor of stilt CID (LM: R2 = 0.01, F1, 221 = 3.12, estimate = −0.07, p = 0.08), June precipitation was the best predictor of willet CID (LM: R2 = 0.02, F1, 125 = 2.35, estimate = −0.03, p = 0.13), and May was the best predictor of plover CID (LM: R2 = 0.03, F1, 96 = 2.75, estimate = −0.05, p = 0.10). Even though precipitation was not a significant effect, CID became earlier as precipitation increased (negative relationship) for the best models for each species, which matched our predictions.

4. Discussion

4.1. Main Findings Regarding Temperature

We found partial support for our hypothesis that increasing temperatures have caused shorebird clutch initiation dates to become earlier over time. While our first model showed that CIDs became earlier when all species were combined, year was not a significant predictor of CID for the stilt or the plover when the models were run with each species separately. However, willet CIDs have become significantly earlier over time (Figure 4C) and with increasing temperatures (Figure 5B). Even though temperature was not a significant predictor of CID for plovers, the trend was in the predicted direction (Figure 5C). Unexpectedly, stilt CIDs became significantly later with increasing temperatures (Figure 5A).
Temperatures have been increasing in our study areas over time (Table 5(A)); thus, we did find support that the timing of nest initiation for willets has likely become earlier over time, due, at least in part, to increasing temperatures. These results for the willet match well with those found by Smith et al. [4], who investigated factors affecting the timing of breeding events in this same species. In their study, Smith et al. [4] found temperatures at their sites from Georgia to Maine had increased over time, leading to earlier springs, and that spring temperature was a strong predictor of willet nest initiation, along with salt marsh biomass. However, they also noted that the correlation between temperature and nest initiation was strongest in the more northern latitudes (above 36°). This pattern makes sense as locations at higher latitudes have been warming more quickly than more southern locations [46] (p. 688), [47]. A study on multiple types of animal and plant species from across the globe demonstrated that species at latitudes above 50° have exhibited stronger phenological shifts compared to those below 50° [48] (see also [49]). All of the nests in our study were below 40°. This might explain why our models looking at how well the temperature in our areas predicted CIDs had lower R2 values than what Smith et al. [4] found, or were insignificant in the case of the plover.

4.2. Main Findings Regarding Precipitation

Our hypothesis that increasing precipitation has caused shorebird CIDs to become earlier over time was unsupported. We found no evidence that precipitation was changing over time around our western nests, and whether it was increasing or decreasing over time near the eastern nests depended on the month (Table 5(B)). We found no evidence that precipitation changes could reliably predict changes in CID for any of our species (Table 6(B)). Although the best models showed the predicted trend of CIDs becoming earlier as precipitation increased, this was never significant, and R2 values in all models were low (0–0.03). Studies have indicated that obvious increases in precipitation in the U.S. only recently started occurring (since the 1970s) [19,20,50]. Only 66 out of 994 of the nest records we found were active from 1970 onwards, and even fewer of these nests were used in the final analysis comparing CID to precipitation near nests. This is one possible explanation for why we were unable to detect any significant trends using precipitation data. Additionally, it may be harder to use precipitation data to predict CID for birds simply because changes in precipitation, unlike changes in surface temperatures, appear to be more variable across locations due to many factors [18]. The average monthly amount of precipitation may not be as important a predictor of changes in nest phenology as changes in the frequency, duration, and intensity of rainfall events [18,51]. Further, precipitation changes in different parts of the country may be different; Trenberth [52] found rainfall occurred for 30% of an hour in the Northwest, butonly occurred 2% of an hour in California, where all of our western nests were located.
In the west, latitude and longitude significantly predicted precipitation, and the two best models showed a decrease in precipitation at more inland locations in March and April (Table 5(B)). This was interesting because stilt laying dates became later with increasing temperatures and at locations further inland (which also had warmer temperatures than locations on the coast). Stilts also did not begin nesting, in general, until mid-April. This could suggest that stilts breeding further from the coast may be delayed by higher temperatures and lower amounts of precipitation. Multiple studies on other species in North America and Europe indicate the importance of successful nesting on food availability (e.g., [4,5,6]). As stilts largely forage on macroinvertebrates in wetlands [34], both temperature and precipitation could affect their food supply. However, we could not find previous studies that specifically focused on factors influencing the timing of stilt nesting behavior. One study found that weather patterns, such as air temperature, precipitation, and storms, all seemed to influence nest initiation in stilts at a constructed wetland, but the authors did not provide evidence for this claim in the form of data or statistical analysis [53].

4.3. Different Responses among Species

Similar to our results, multiple studies in North America and Europe, even those that investigated many more species and nests over wider geographic areas, have found species-specific differences in response to climate changes (e.g., [7,10,13,14,16,17]). While both of our eastern species showed a trend towards earlier CID as temperatures increased, the western species showed the opposite trend. Additionally, temperatures and precipitation averages in different months predicted CID differently for each species. Mayor et al. [10] also found regional differences in passerine arrival dates in western forests vs. eastern forests and cautioned against applying patterns found in one species or one region to other species or regions not included in the study. Differences in how species respond to changes in climate can also occur at more local scales [7]. Species-specific factors, including migration distance [6,22], reproductive strategy [7], and trophic level [14], can also account for mixed responses to climate change. Additionally, changes in temperature, precipitation, or other factors related to climate change on wintering grounds or stopover sites used during migration may affect nesting phenology in different ways depending on the region and species [54]. Thus, it is imperative that more studies on various species across wide geographic ranges and time be conducted to increase our understanding of how individual species in specific regions have been affected by climate change and how flexible these species are at responding and adapting to these rapid changes in their environments. This is especially important for our study species, as Galbraith et al. [22] indicated that all three shorebird species are more vulnerable to extinction—in part due to climate change—than currently considered based on the U.S. Shorebird Conservation Plan Risk Categories. In particular, the plover was listed as “critical” by Galbraith et al. [22], the highest risk category (see also [55,56]), yet this species has been sorely understudied [27].

4.4. The Use of Museum Records to Study Nest Phenology

As a final point, we encourage future studies to increase their use of museum egg records to investigate nesting phenological changes in avian species (see also [17]). Museum collections are a valuable resource for scientific research due to their worldwide sampling across temporal, geographical, and taxonomic species, and museum records are typically freely available and easily accessible online [57,58]. Broad sampling can be impossible to collect while performing fieldwork, and long-term data collection is vital for phenological studies. However, the use of museum collections for phenological investigations has been criticized in the past for a possible lack of accuracy in the data from the collectors and collection biases [59]. Collection bias may lead to over-representation of certain egg traits, such as color and size [60], preference for parasitized clutches, and for clutches laid much earlier in the season than the natural variation seen in the environment [61]. However, concern about these biases may be exaggerated, as accuracy issues with a particular collector or collection are often obvious and can be accounted for if used with caution [58,62]. Further, while the information on a clutch card may be too broad or missing data, this and other studies [17,43,44,62] have demonstrated ways to reasonably estimate clutch initiation dates from the data recorded on clutch cards. Gathering enough nest records can also help mitigate issues with bias and missing information. Ultimately, museum collections are beneficial for long-term phenological analyses because researchers can obtain and utilize valuable information from the collections to answer important scientific questions over a long span of time and space [44,57,58].

5. Conclusions

Our findings demonstrate that different species are affected in varying ways by their environments. Stilt laying dates became later at higher temperatures and at more inland locations. Willet laying dates became earlier at higher temperatures, as predicted, and while this trend was also present for plovers, latitude was the most important factor for this species. The variability we and previous studies have found among species’ responses to changes in climate makes it even more important for additional studies on shorebird nesting phenology to be conducted. To determine the best conservation practices that should be implemented for shorebirds, we need a better understanding of exactly how these species might be influenced by changes in climate, especially considering how vulnerable they likely are to such changes [22].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani13152459/s1, Table S1: Records of all Nests with Justification for Estimated CID; File S1: Temperature and Precipitation Tables for the East and West Regions; File S2: CID vs. Temperature and Precipitation Tables for all Three Species.

Author Contributions

Conceptualization, V.E.A., A.G., A.B., M.B., E.B., H.W., E.C. and G.P.; methodology, V.E.A., A.G., A.B., M.B., E.B., H.W., E.C. and G.P.; data curation, V.E.A., A.G., A.B., M.B., E.B., H.W., E.C., G.P. and F.H.; investigation, V.E.A., A.G., A.B., M.B., E.B., H.W., E.C., G.P. and F.H; software, V.E.A.; formal analysis, V.E.A., A.G. and E.B.; writing—original draft preparation, V.E.A., A.G., A.B., M.B. and F.H.; writing—review and editing, V.E.A., A.B. and M.B.; visualization, V.E.A., A.G. and E.B.; supervision, V.E.A.; project administration, V.E.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding for publication from Clemson University, specifically, the Creative Inquiry program, the Libraries Open Access Publishing Fund and the Department of Biological Sciences.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials. The datasets used in this study are available in Supplementary Table S1, File S1 and File S2.

Acknowledgments

This study was made possible through the Clemson University Creative Inquiry program. We thank the director of this program, Cora Allard, and the chair of the Biological Sciences Department, Saara DeWalt, for all of their assistance in the development of the Creative Inquiry course that provided an opportunity for undergraduate students to conduct this research, as well as in funding the publication of this work. We acknowledge D. Russell at the Natural History Museum in Tring, UK, for providing the initial inspiration for this project. We thank Melissa Fuentes who gave us access to the egg collection records at the Clemson University Bob and Betsy Campbell Natural History Museum. We thank Robyn Bailey and the other individuals with the NestWatch program at the Cornell Lab of Ornithology, who provided us access to the data cards from NestWatch and were so helpful in sending us resources to explain how to interpret the cards. We thank the volunteers who collected the nesting data for this study through The Cornell Laboratory of Ornithology’s Nest Box Network, The Birdhouse Network, NestWatch, the Smithsonian Migratory Bird Research Center’s Neighborhood NestWatch, and various state nest monitoring projects that have contributed their data. We appreciate the statistical advice for this manuscript provided voluntarily by Matt Koski, and we thank Jeff Robel, who gave us instructions on how to access and download data from NOAA. Finally, we are indebted to Jeffrey Adams, who assisted us with organizing and extracting information from the large amounts of data used in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hegerl, G.C.; Broennimann, S.; Cowan, T.; Friedman, A.R.; Hawkins, E.; Iles, C.E.; Mueller, W.; Schurer, A.; Undorf, S. Causes of climate change over the historical record. Environ. Res. Lett. 2019, 14, 123006. [Google Scholar] [CrossRef]
  2. Liebezeit, J.R.; Gurney, K.E.B.; Budde, M.; Zack, S.; Ward, D. Phenological advancement in arctic bird species: Relative importance of snow melt and ecological factors. Polar Biol. 2014, 37, 1309–1320. [Google Scholar] [CrossRef]
  3. Kwon, E.; Weiser, E.L.; Lanctot, R.B.; Brown, S.C.; Gates, H.R.; Gilchrist, H.G.; Kendall, S.J.; Lank, D.B.; Liebezeit, J.R.; McKinnon, L.; et al. Geographic variation in the intensity of warming and phenological mismatch between Arctic shorebirds and invertebrates. Ecol. Monogr. 2019, 89, e01383. [Google Scholar] [CrossRef]
  4. Smith, J.A.M.; Regan, K.; Cooper, N.W.; Johnson, L.; Olson, E.; Green, A.; Tash, J.; Evers, D.C.; Marra, P.P. A green wave of saltmarsh productivity predicts the timing of the annual cycle in a long-distance migratory shorebird. Sci. Rep. 2020, 10, 20658. [Google Scholar] [CrossRef]
  5. Van Noordwijk, A.J.; McCleery, R.H.; Perrins, C.M. Selection for the timing of great tit breeding in relation to caterpillar growth and temperature. J. Anim. Ecol. 1995, 64, 451–458. [Google Scholar]
  6. Both, C.; Visser, M.E. The effect of climate change on the correlation between avian life-history traits. Glob. Chang. Biol. 2005, 11, 1606–1613. [Google Scholar] [CrossRef] [Green Version]
  7. McGuire, R.L.; Lanctot, R.B.; Saalfeld, S.T.; Ruthrauff, D.R.; Liebezeit, J.R. Shorebird reproductive response to exceptionally early and late springs varies across sites in arctic alaska. Front. Ecol. Evol. 2020, 8, 577652. [Google Scholar] [CrossRef]
  8. Both, C.; Visser, M.E. Adjustment to climate change is constrained by arrival date in a long-distance migrant bird. Nature 2001, 411, 296–298. [Google Scholar] [CrossRef] [Green Version]
  9. Carey, C. The impacts of climate change on the annual cycles of birds. Philos. Trans. R. Soc. B Biol. Sci. 2009, 364, 3321–3330. [Google Scholar] [CrossRef]
  10. Mayor, S.J.; Guralnick, R.P.; Tingley, M.W.; Otegui, J.; Withey, J.C.; Elmendorf, S.C.; Andrew, M.E.; Leyk, S.; Pearse, I.S.; Schneider, D.C. Increasing phenological asynchrony between spring green-up and arrival of migratory birds. Sci. Rep. 2017, 7, 1902. [Google Scholar] [CrossRef] [Green Version]
  11. Kwon, E.; English, W.B.; Weiser, E.L.; Franks, S.E.; Hodkinson, D.J.; Lank, D.B.; Sandercock, B.K. Delayed egg-laying and shortened incubation duration of Arctic-breeding shorebirds coincide with climate cooling. Ecol. Evol. 2017, 8, 1339–1351. [Google Scholar] [CrossRef] [PubMed]
  12. Halupka, L.; Czyż, B.; Dominguez, C.M.M. The effect of climate change on laying dates, clutch size and productivity of Eurasian Coots Fulica atra. Int. J. Biometeorol. 2020, 64, 1857–1863. [Google Scholar] [CrossRef]
  13. Crick, H.Q.P.; Sparks, T.H. Climate change related to egg-laying trends. Nature 1999, 399, 423. [Google Scholar] [CrossRef]
  14. Both, C.; van Asch, M.; Bijlsma, R.G.; Burg, A.B.v.D.; Visser, M.E. Climate change and unequal phenological changes across four trophic levels: Constraints or adaptations? J. Anim. Ecol. 2008, 78, 73–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Watts, H.E.; Jimenez, D.; Pacheco, V.; Vilgalys, T.P. Temperature-correlated shifts in the timing of egg-laying in House Finches Haemorhous mexicanus. Ibis 2018, 161, 428–434. [Google Scholar] [CrossRef]
  16. Torti, V.M.; Dunn, P.O. Variable effects of climate change on six species of North American birds. Oecologia 2005, 145, 486–495. [Google Scholar] [CrossRef]
  17. Bates, J.M.; Fidino, M.; Nowak-Boyd, L.; Strausberger, B.M.; Schmidt, K.A.; Whelan, C.J. Climate change affects bird nesting phenology: Comparing contemporary field and historical museum nesting records. J. Anim. Ecol. 2022, 92, 263–272. [Google Scholar] [CrossRef]
  18. Trenberth, K.E.; Dai, A.; Rasmussen, R.M.; Parsons, D.B. The changing character of precipitation. Bull. Am. Meteorol. Soc. 2003, 84, 1205–1218. [Google Scholar] [CrossRef]
  19. Ross, R.J.; Elliott, W.P. Tropospheric water vapor climatology and trends over North America: 1973–1993. J. Clim. 1996, 9, 3561–3574. [Google Scholar] [CrossRef]
  20. Ross, R.J.; Elliott, W.P. Radiosonde-based northern hemisphere tropospheric water vapor trends. J. Clim. 2001, 14, 1602–1612. [Google Scholar] [CrossRef]
  21. Hennessy, K.J.; Gregory, J.M.; Mitchell, J.F.B. Changes in daily precipitation under enhanced greenhouse condi-tions. Clim. Dyn. 1997, 13, 667–680. [Google Scholar] [CrossRef]
  22. Galbraith, H.; DesRochers, D.W.; Brown, S.; Reed, J.M. Predicting vulnerabilities of North American shorebirds to climate change. PLoS ONE 2014, 9, e108899. [Google Scholar] [CrossRef]
  23. Cohen, J.B.; Houghton, L.M.; Fraser, J.D. Nesting density and reproductive success of piping plovers in response to storm- and human-created habitat changes. Wildl. Monogr. 2009, 173, 1–24. [Google Scholar] [CrossRef]
  24. Convertino, M.; Bockelie, A.; Kiker, G.A.; Muñoz-Carpena, R.; Linkov, I. Shorebird patches as fingerprints of fractal coastline fluctuations due to climate change. Ecol. Process. 2012, 1, 9. [Google Scholar] [CrossRef] [Green Version]
  25. Cohen, M.C.L.; de Souza, A.V.; Liu, K.-B.; Rodrigues, E.; Yao, Q.; Pessenda, L.C.R.; Rossetti, D.; Ryu, J.; Dietz, M. Effects of Beach Nourishment Project on Coastal Geomorphology and Mangrove Dynamics in Southern Louisiana, USA. Remote Sens. 2021, 13, 2688. [Google Scholar] [CrossRef]
  26. Reed, K.A.; Wehner, M.F.; Zarzycki, C.M. Attribution of 2020 hurricane season extreme rainfall to human-induced climate change. Nat. Commun. 2022, 13, 1905. [Google Scholar] [CrossRef]
  27. DeRose-Wilson, A.; Fraser, J.D.; Karpanty, S.M.; Catlin, D.H. Nest-site selection and demography of Wilson’s Plovers on a North Carolina barrier island. J. Field Ornithol. 2013, 84, 329–344. [Google Scholar] [CrossRef]
  28. Ackerman, J.T.; Herzog, M.P.; Takekawa, J.Y.; Hartman, C.A. Comparative reproductive biology of sympatric species: Nest and chick survival of American avocets and black-necked stilts. J. Avian Biol. 2014, 45, 609–623. [Google Scholar] [CrossRef]
  29. White, E.E.; Ury, E.A.; Bernhardt, E.S.; Yang, X. Climate change driving widespread loss of coastal forested wetlands throughout the north american coastal plain. Ecosystems 2021, 25, 812–827. [Google Scholar] [CrossRef]
  30. Smith, P.A.; Gilchrist, H.G.; Smith, J.N.M. Effects of nest habitat, food, and parental behavior on shorebird nest success. Condor 2007, 109, 15–31. [Google Scholar] [CrossRef]
  31. Arnold, T.W. What limits clutch size in waders? J. Avian Biol. 1999, 30, 216–220. [Google Scholar] [CrossRef]
  32. Griffith, G.E.; Omernik, J.M.; Smith, D.W.; Cook, T.D.; Tallyn, E.; Moseley, K.; Johnson, C.B. Ecoregions of California (Poster): U.S. Geological Survey Open-File Report 2016–1021, with Map, Scale 1:1,100,000; U.S. Geological Survey: Reston, VA, USA, 2016. [Google Scholar] [CrossRef]
  33. Noss, R.F.; Platt, W.J.; Sorrie, B.A.; Weakley, A.S.; Means, D.B.; Costanza, J.; Peet, R.K. How global biodiversity hotspots may go unrecognized: Lessons from the North American Coastal Plain. Divers. Distrib. 2014, 21, 236–244. [Google Scholar] [CrossRef]
  34. Robinson, J.A.; Reed, J.M.; Skorupa, J.P.; Oring, L.W. Black-necked Stilt (Himantopus mexicanus), version 1.0. In Birds of the World; Poole, A.F., Gill, F.B., Eds.; Cornell Lab of Ornithology: Ithaca, NY, USA, 2020. [Google Scholar] [CrossRef]
  35. Zdravkovic, M.G.; Corbat, C.A.; Bergstrom, P.W. Wilson’s Plover (Charadrius wilsonia), version 1.0. In Birds of the World; Rodewald, P.G., Ed.; Cornell Lab of Ornithology: Ithaca, NY, USA, 2020. [Google Scholar] [CrossRef]
  36. Lowther, P.E.; Iii, H.D.D.; Gratto-Trevor, C.L. Willet (Tringa semipalmata), version 1.0. In Birds of the World; Poole, A.F., Gill, F.B., Eds.; Cornell Lab of Ornithology: Ithaca, NY, USA, 2020. [Google Scholar] [CrossRef]
  37. American Museum of Natural History, New York, USA. Available online: https://data.nhm.ac.uk/dataset/collection-specimens/resource/05ff2255-c38a-40c9-b657-4ccb55ab2feb?filters=collectionCode%3Azoo (accessed on 18 June 2023).
  38. Denver Museum of Nature and Science, Colorado, USA. Available online: https://arctos.database.museum/collection/DMNS:Egg (accessed on 18 June 2023).
  39. Jurica-Suchy Nature Museum at Benedictine University, Illinois, USA. Available online: https://arctos.database.museum/collection/JSNM:Egg (accessed on 18 June 2023).
  40. Museum of Vertebrate Zoology (MVZ), University of California, Berkeley, California, USA. Available online: https://arctos.database.museum/collection/MVZ:Egg (accessed on 18 June 2023).
  41. Smithsonian National Museum of Natural History, Smithsonian Institution, Washington, DC, USA. Available online: https://siarchives.si.edu/ (accessed on 18 June 2023).
  42. Western Foundation of Vertebrate Zoology, California, USA. Available online: https://collections.wfvz.org/ (accessed on 18 June 2023).
  43. McNair, D.B. A Comparison of Oology and Nest Record Card Data in Evaluating the Reproductive Biology of Lark Sparrows, Chondestes grammacus. Southwest. Nat. 1985, 30, 213. [Google Scholar] [CrossRef]
  44. Scharlemann, J.P.W. Museum egg collections as stores of long-term phenological data. Int. J. Biometeorol. 2001, 45, 208–211. [Google Scholar] [CrossRef]
  45. MapTiler Team. Epsg.io. MapTiler. 2022. Available online: https://epsg.io/map#srs=4326&x=0.000000&y=0.000000&z=1&layer=streets (accessed on 18 June 2023).
  46. Intergovernmental Panel on Climate Change (IPCC). Impacts, adaptation, and vulnerability. Part B: Regional aspects. In Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Barros, V.R., Field, C.B., Dokken, D.J., Mastrandrea, M.D., Mach, K.J., Bilir, T.E., Chatterjee, M., Ebi, K.L., Estrada, Y.O., Genova, R.C., et al., Eds.; Cam-bridge University Press: Cambridge, UK, 2014; p. 688. [Google Scholar]
  47. AMAP. Arctic Climate Issues 2011: Changes in Arctic Snow, Water, Ice and Permafrost. SWIPA 2011; Overview Report; AMAP: Oslo, Norway, 2012; 112p. [Google Scholar]
  48. Root, T.L.; Price, J.T.; Hall, K.R.; Schneider, S.H.; Rosenzweig, C.; Pounds, J.A. Fingerprints of global warming on wild animals and plants. Nature 2003, 421, 57–60. [Google Scholar] [CrossRef]
  49. Shave, A.; Garroway, C.; Siegrist, J.; Fraser, K.C. Timing to temperature: Egg-laying dates respond to temperature and are under stronger selection at northern latitudes. Ecosphere 2019, 10, e02974. [Google Scholar] [CrossRef] [Green Version]
  50. Fischer, E.M.; Knutti, R. Observed heavy precipitation increase confirms theory and early models. Nat. Clim. Chang. 2016, 6, 986–991. [Google Scholar] [CrossRef]
  51. Berg, P.; Moseley, C.; Haerter, J.O. Strong increase in convective precipitation in response to higher temperatures. Nat. Geosci. 2013, 6, 181–185. [Google Scholar] [CrossRef]
  52. Trenberth, K.E. Atmospheric moisture residence Times AND cycling: Implications for rainfall rates and climate change. Clim. Chang. 1998, 39, 667–694. [Google Scholar] [CrossRef]
  53. Davis, D.E.; Hanson, C.H.; Hansen, R.B. Constructed wetland habitat for American avocet and black-necked stilt foraging and nesting. J. Wildl. Manag. 2008, 72, 143–151. [Google Scholar] [CrossRef]
  54. Raquel, A.J.; Devries, J.H.; Howerter, D.W.; Clark, R.G. Reproductive consequences of climate variability in migratory birds: Evidence for species-specific responses to spring phenology and cross-seasonal effects. Oecologia 2019, 191, 217–229. [Google Scholar] [CrossRef]
  55. Andres, B.A.; Smith, P.A.; Morrison, R.I.G.; Gratto-Trevor, C.L.; Brown, S.C.; Friis, C.A. Population estimates of North American shorebirds, 2012. Wader Study Group Bull. 2012, 119, 178–194. [Google Scholar]
  56. Zdravkovic, M.G. Conservation Plan for the Wilson’s Plover (Charadrius wilsonia); Version 1.0.; Manomet Center for Conservation Sciences: Manomet, MA, USA, 2013. [Google Scholar]
  57. Green, R.E.; Scharlemann, J.P.W. Egg and skin collections as a resource for long-term ecological studies. Bull. Br. Ornithol. Club 2003, 123, 165–176. [Google Scholar]
  58. Marini, M.; Hall, L.; Bates, J.; Steinheimer, F.D.; McGowan, R.; Silveira, L.F.; Lijtmaer, D.A.; Tubaro, P.L.; Córdoba-Córdoba, S.; Gamauf, A.; et al. The five million bird eggs in the world’s museum collections are an invaluable and underused resource. Ornithology 2020, 137, ukaa036. [Google Scholar] [CrossRef]
  59. Lack, D. Clutch and brood size in the robin. Br. Birds 1946, 39, 98–109, 130–135. [Google Scholar]
  60. Thompson, J.E.; Birkhead, T.R. Avian egg collections: Museum collection bias driven by shape and size. J. Avian Biol. 2020, 51, jav.02507. [Google Scholar] [CrossRef]
  61. Kiff, L.F. History, Present Status, and Future Prospects of Avian Eggshell Collections in North America. Auk 2005, 122, 994–999. [Google Scholar] [CrossRef]
  62. McNair, D.B. Egg data slips-are they useful for information on egg-laying dates and clutch size? Condor 1987, 89, 369–376. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Nest locations of the stilt in California (A), and the willet (B) and plover (C) on the east coast of the U.S. Circle size represents the relative number of nests in an area with the actual number included inside each circle. Map created by A.G. in ArcGIS Pro 2.9 by ESRI with further edits in Photoshop by E.B.
Figure 1. Nest locations of the stilt in California (A), and the willet (B) and plover (C) on the east coast of the U.S. Circle size represents the relative number of nests in an area with the actual number included inside each circle. Map created by A.G. in ArcGIS Pro 2.9 by ESRI with further edits in Photoshop by E.B.
Animals 13 02459 g001
Figure 2. Histograms showing (A) the distribution of nests of each species collected from 1852 to 1989 for all 944 nests (stilt = 412, willet = 287, plover = 295) and (B) the timing of clutch initiation during the breeding season for each species using only nests where clutch initiation date could be reasonably estimated (total nests = 733, stilt = 339, willet = 198, plover = 196). In (B), the x-axis shows calendar dates above with the corresponding Julian date below.
Figure 2. Histograms showing (A) the distribution of nests of each species collected from 1852 to 1989 for all 944 nests (stilt = 412, willet = 287, plover = 295) and (B) the timing of clutch initiation during the breeding season for each species using only nests where clutch initiation date could be reasonably estimated (total nests = 733, stilt = 339, willet = 198, plover = 196). In (B), the x-axis shows calendar dates above with the corresponding Julian date below.
Animals 13 02459 g002
Figure 3. Least square means of clutch initiation dates (CID) ± 95% CI for all three species. An asterisk (*) indicates a significant difference between groups (LM: Tukey HSD, p < 0.05). Total nests = 733, stilt = 339, willet = 198, and plover = 196 for nests from 1852–1989.
Figure 3. Least square means of clutch initiation dates (CID) ± 95% CI for all three species. An asterisk (*) indicates a significant difference between groups (LM: Tukey HSD, p < 0.05). Total nests = 733, stilt = 339, willet = 198, and plover = 196 for nests from 1852–1989.
Animals 13 02459 g003
Figure 4. Leverage plots showing the relationship between clutch initiation dates (CID) and nest year, nest latitude, and nest longitude for the stilt (A), plover (B), and willet (C) in a linear model. Results are reported for the best linear model for each species which only included significant fixed effects. The red solid line is a least squares line of best fit surrounded by 95% confidence bands. The blue horizontal line represents the hypothesis that the effect in question has zero influence on CID. Stilt nests = 339, willet nests = 198, and plover nests = 196 for nests from 1852–1989.
Figure 4. Leverage plots showing the relationship between clutch initiation dates (CID) and nest year, nest latitude, and nest longitude for the stilt (A), plover (B), and willet (C) in a linear model. Results are reported for the best linear model for each species which only included significant fixed effects. The red solid line is a least squares line of best fit surrounded by 95% confidence bands. The blue horizontal line represents the hypothesis that the effect in question has zero influence on CID. Stilt nests = 339, willet nests = 198, and plover nests = 196 for nests from 1852–1989.
Animals 13 02459 g004
Figure 5. Relationship between clutch initiation date (CID) and nest temperature for the month that best predicted CID for the stilt (A), n = 221, years from 1896–1975; willet (B), n = 112, years from 1896–1989; and plover (C), n = 79, years from 1897–1970. The p-value and R2 values have been provided on the graph for each model.
Figure 5. Relationship between clutch initiation date (CID) and nest temperature for the month that best predicted CID for the stilt (A), n = 221, years from 1896–1975; willet (B), n = 112, years from 1896–1989; and plover (C), n = 79, years from 1897–1970. The p-value and R2 values have been provided on the graph for each model.
Animals 13 02459 g005
Table 1. Summary of metrics used to determine laying date for each species.
Table 1. Summary of metrics used to determine laying date for each species.
SpeciesAverage Clutch SizeAverage Incubation PeriodAssumed Laying Period
Stilt4254
Willet4266
Plover3255
Table 2. Number of clutches for each species for each source used in this study. Records for sources with asterisks (*) were retrieved from the Arctos Database website (https://arctos.db.org, accessed on October 2021).
Table 2. Number of clutches for each species for each source used in this study. Records for sources with asterisks (*) were retrieved from the Arctos Database website (https://arctos.db.org, accessed on October 2021).
Source for ClutchStiltWilletPloverTotal
American Museum of Natural History [37]0202747
Clemson University Bob and Betsy Campbell Natural History Museum7362669
Denver Museum of Nature and Science * [38]155626
Jurica-Suchy Nature Museum at Benedictine University * [39]1146
Museum of Vertebrate Zoology * [40]272938
Smithsonian National Museum of Natural History [41]592741
Western Foundation of Vertebrate Zoology [42]314203181698
Cornell NestWatch Program43111569
Total412287295994
Table 3. Actual terms used by collectors to describe incubation stages of clutches at collection. Similar terms used by collectors were combined into our own incubation stage categories based on the estimated percentage of incubation the clutch had undergone at collection. This percentage was translated into an estimated number of days the eggs of a clutch had likely been incubated based on the typical incubation period for each species. The assumed incubation period for the stilt and plover was 25 days and 26 days for the willet.
Table 3. Actual terms used by collectors to describe incubation stages of clutches at collection. Similar terms used by collectors were combined into our own incubation stage categories based on the estimated percentage of incubation the clutch had undergone at collection. This percentage was translated into an estimated number of days the eggs of a clutch had likely been incubated based on the typical incubation period for each species. The assumed incubation period for the stilt and plover was 25 days and 26 days for the willet.
Collector Incubation Stage TermsIncubation Stage CategoryEstimated % of Incubation at CollectionEstimated Number of Days of Incubation
All fresh, Fresh, Nearly fresh, Perfectly freshFresh0%0
Fresh to slightMedian of “Fresh” and “Slight”8%Only applies to two willet clutches: 2
Slight, Slightly incubated, Very slight, Began, Begun, Begun in all, Just begun, Commenced, Barely commenced, Just commenced, Trace, Trace only, Small embryos, Started, Just started, Few daysSlight Incubation15–16%4
1/4 incubated, 25%, About 25%1/4 incubation25%Stilt and plover: 6, willet: 7
1/3, 1/3 incubated, 1/3 advanced, About 1/3, About 1/3 advanced, About one third1/3 incubation33%Stilt and plover: 8, willet: 9
2/5 incubated2/5 incubation40%Only applies to one stilt clutch: 10
1/2, 1/2 incubated, Halfway, 1/2 advanced, About 1/2, About 1/2 advanced, Half, About half, About one half, Medium embryos, Well along, Well begun, Well started, 50% incubated, Medium, Nearly 1/21/2 incubation50%13
Advanced, Advanced-far, Highly advanced, Large embryos, Heavy, Considerable, Far advanced, Very far advanced, Well advancedAdvanced Incubation80%Stilt and plover: 20, willet: 21
Ready to hatch, Nearly ready to hatch, Feathers ready to hatchTwo days left of incubation96%Stilt and plover: 24, willet: 25
Table 4. Summary of the linear model results for our first analysis investigating if CID for all species combined has changed over time or space. Results reported are only for the best model and only include significant effects and interactions. Total nests = 733, stilt = 339, willet = 198, and plover = 196 for nests from 1852–1989.
Table 4. Summary of the linear model results for our first analysis investigating if CID for all species combined has changed over time or space. Results reported are only for the best model and only include significant effects and interactions. Total nests = 733, stilt = 339, willet = 198, and plover = 196 for nests from 1852–1989.
Fixed EffectsEstimatep-Value
Nest Year−0.090.003
Nest Latitude1.750.008
Nest Longitude2.03<0.001
Species: Stilt77.34<0.001
Species: Willet−54.91<0.001
Species: Plover−22.430.03
Species (Stilt) × Longitude2.390.02
Species (Willet) × Longitude−0.570.37
Species (Plover) × Longitude−1.820.02
Table 5. Summary of the mixed linear model results for our analysis investigating if temperature (A) or precipitation (B) in our study areas in the west (for stilt nests) and east (for plover and willet nests) have changed over time and space for the first four months of the breeding season. In these models, the weather station ID where temperature and precipitation data were retrieved was always treated as a random effect to control for multiple samples taken over the years the station was active. Results reported are only for the best models (effects with dashed lines were insignificant and removed from the final model). Climate data were retrieved from 239 stations in California close to the stilt nests and 432 stations in the east (Florida, Georgia, North Carolina, South Carolina, Virginia) close to the plover and willet nests from 1851–2022.
Table 5. Summary of the mixed linear model results for our analysis investigating if temperature (A) or precipitation (B) in our study areas in the west (for stilt nests) and east (for plover and willet nests) have changed over time and space for the first four months of the breeding season. In these models, the weather station ID where temperature and precipitation data were retrieved was always treated as a random effect to control for multiple samples taken over the years the station was active. Results reported are only for the best models (effects with dashed lines were insignificant and removed from the final model). Climate data were retrieved from 239 stations in California close to the stilt nests and 432 stations in the east (Florida, Georgia, North Carolina, South Carolina, Virginia) close to the plover and willet nests from 1851–2022.
A West TemperatureEast Temperature
MonthFixed EffectsAdjusted R2Estimatep-ValueAdjusted R2Estimatep-Value
MarchYear0.830.02<0.0010.800.0040.02
Latitude−1.17<0.001−1.11<0.001
Longitude----
AprilYear0.840.01<0.0010.860.01<0.001
Latitude--−0.84<0.001
Longitude0.90<0.001--
MayYear0.870.01<0.0010.800.004<0.001
Latitude--−0.43<0.001
Longitude0.88<0.001−0.32<0.001
Year0.920.01<0.0010.690.01<0.001
JuneLatitude0.960.01−0.22<0.001
Longitude1.68<0.001−0.240.001
B West PrecipitationEast Precipitation
MonthFixed EffectsAdjusted R2Estimatep-valueAdjusted R2Estimatep-value
MarchYear0.28--0.09−0.14<0.001
Latitude−4.880.0066.76<0.001
Longitude−13.15<0.001−6.68<0.001
AprilYear0.23--0.090.21<0.001
Latitude−1.760.053 13.66<0.001
Longitude−4.57<0.001−4.08<0.001
MayYear0.14--0.07--
Latitude1.68<0.001−4.71<0.001
Longitude--7.02<0.001
Year0.15--0.200.28<0.001
JuneLatitude1.12<0.001−9.71<0.001
Longitude0.410.04--
1 This model had the lowest AIC score even though latitude is only near significant.
Table 6. Summary of the linear model results for our analysis investigating if nest temperature (A) or nest precipitation (B) could significantly predict clutch initiation date (CID) for the stilt, willet, and plover, and which monthly average was the best predictor. In this model, the only fixed effect was either nest temperature or nest precipitation.
Table 6. Summary of the linear model results for our analysis investigating if nest temperature (A) or nest precipitation (B) could significantly predict clutch initiation date (CID) for the stilt, willet, and plover, and which monthly average was the best predictor. In this model, the only fixed effect was either nest temperature or nest precipitation.
ASpeciesMonthR2F-Value with dfEstimatep-Valuen
StiltMarch0.012.28 1, 2100.850.13212
April0.013.00 1, 2110.930.08213
May0.025.02 1, 2120.810.03214
June 10.048.52 1, 2190.860.004221
WilletMarch0.089.22 1, 113−2.100.003115
April0.1316.10 1, 111−2.93<0.001113
May 10.1520.11 1, 110−4.07<0.001112
June0.078.71 1, 108−3.960.004110
PloverMarch0.010.84 1, 75−0.810.3677
April 10.02331.83 1, 77−1.170.1879
May0.02271.80 1, 771.950.1879
June0.0171.33 1, 752.090.2577
BSpeciesMonthR2F-value with dfEstimatep-valuen
StiltMarch 10.0143.12 1, 221−0.070.08223
April00.15 1, 215−0.020.70217
May00.01 1, 207−0.010.93209
June0.0112.32 1, 213−0.380.13215
WilletMarch00.15 1, 1240.010.70126
April0.0172.05 1, 1180.040.15120
May00.08 1, 1260.010.78128
June 10.0182.35 1, 125−0.030.13127
PloverMarch0.010.60 1, 950.020.4497
April00.03 1, 93−0.010.8795
May 10.032.75 1, 96−0.050.1098
June0.011.05 1, 96−0.020.3198
1 This model had the highest R2 and was considered to be the best predictor of CID for that species.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abernathy, V.E.; Good, A.; Blanchard, A.; Bongiovanni, M.; Bonds, E.; Warner, H.; Chaknis, E.; Pulsifer, G.; Huntley, F. The Effects of Climate Change on the Nesting Phenology of Three Shorebird Species in the United States. Animals 2023, 13, 2459. https://doi.org/10.3390/ani13152459

AMA Style

Abernathy VE, Good A, Blanchard A, Bongiovanni M, Bonds E, Warner H, Chaknis E, Pulsifer G, Huntley F. The Effects of Climate Change on the Nesting Phenology of Three Shorebird Species in the United States. Animals. 2023; 13(15):2459. https://doi.org/10.3390/ani13152459

Chicago/Turabian Style

Abernathy, Virginia E., Abby Good, Autum Blanchard, Marlisa Bongiovanni, Emily Bonds, Hampton Warner, Eleni Chaknis, Gabriella Pulsifer, and Faith Huntley. 2023. "The Effects of Climate Change on the Nesting Phenology of Three Shorebird Species in the United States" Animals 13, no. 15: 2459. https://doi.org/10.3390/ani13152459

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop