1,599
Views
0
CrossRef citations to date
0
Altmetric
Bee management

A national survey of managed honey bee colony losses in the USA: Results from the Bee Informed Partnership for 2020–21 and 2021–22

ORCID Icon, ORCID Icon, , ORCID Icon & ORCID Icon
Pages 1-14 | Received 24 Mar 2023, Accepted 16 Aug 2023, Published online: 16 Oct 2023

Abstract

In recent years, high losses of Western honey bee (Apis mellifera) colonies have been reported in the United States. Because honey bees are important to agricultural systems, it is critical to document when high losses occur and to explore patterns of loss among beekeeper subgroups. We used a voluntary, retrospective questionnaire to ask U.S. beekeepers (backyard, sideline, and commercial) about their colony loss during the 2020–21 and 2021–22 beekeeping seasons. We found that U.S. beekeepers lost 50.8% (38.0–63.1; 95% b.CI) of their colonies in 2020–21, the highest annual (year-long) loss reported to date. During the following year, 2021–22, beekeepers lost 39.0% (31.5–47.9; 95% b.CI) of their colonies, which is close to the average of previously published annual loss rates (40.9%). In contrast to previous years, backyard beekeepers (managing 50 or fewer colonies) had an elevated summer loss across both years and commercial beekeepers (managing >500 colonies) reported “weather” as an important cause of colony loss in the summer of 2021. Our results show that severe colony losses still occur periodically. From these data on colony loss rates, it does not appear that the health of honey bee colonies in the U.S. has improved since this survey began in 2008. We suggest that research should continue to focus on stressors that beekeepers most often perceive to be leading causes of loss—Varroa destructor mites, queen issues, and weather.

Introduction

To document and explain managed Western honey bee (Apis mellifera L.) colony losses in the United States, the non-profit Bee Informed Partnership (BIP) launched its Loss and Management Survey in 2011 as a continuation of the national Winter Loss Survey which began in 2007 (vanEngelsdorp et al., Citation2007). The survey initially focused on colony losses in the winter season (1 October to 1 April), when colonies depend on the extended survival of individual honey bees during a time of limited resource availability and opportunities to raise young bees (van Dooremalen et al., Citation2012). However, summer colony loss (1 April to 1 October) and annual colony loss (1 April to 1 April) have been reported since the 2012–13 survey year (Steinhauer et al., Citation2014). This and subsequent studies showed that substantial colony losses occur not only during winter but also during summer (Bruckner, Wilson, et al., Citation2023; Kulhanek et al., Citation2017; Lee et al., Citation2015; Seitz et al., Citation2015; Steinhauer et al., Citation2014).

While loss rates fluctuate across years, winter losses are usually higher (22–38%; mean of published reports 29.0%) than those in summer (18–32%; mean 23.5%). Losses over the annual season have ranged from 34 to 45%, with a mean of 40.9% of colonies lost during the whole year. Backyard beekeeping operations (i.e., those managing 50 or fewer colonies) typically have higher colony loss rates than sideline and commercial beekeepers (i.e., those managing 51–500 and more than 500 colonies, respectively) (Bruckner, Wilson, et al., Citation2023; Kulhanek et al., Citation2017; Lee et al., Citation2015; Seitz et al., Citation2015; Spleen et al., Citation2013; Steinhauer et al., Citation2014; vanEngelsdorp et al., Citation2007, Citation2008, Citation2010, Citation2011, Citation2012). Even during years with comparatively low winter losses, beekeepers have reported losing more colonies than they consider acceptable (Kulhanek et al., Citation2017).

To counterbalance losses, beekeepers can increase their number of colonies from within their operations through splitting—the act of dividing colonies into two or more smaller colonies (Connor, Citation2014)—or by obtaining colonies from other beekeepers. However, high colony losses could potentially lead to declines in the total colony population if newly generated colonies are insufficient to compensate for losses (vanEngelsdorp & Meixner, Citation2010). Despite reports of high losses, estimates of the number of honey producing colonies in the U.S. have remained relatively stable since 2006, fluctuating between 2.3 and 2.8 million colonies (USDA-NASS, Citation2009, Citation2020). While it is important to have periodic information on the number of colonies in the U.S., it is also crucial to track colony loss rates because high colony losses can suggest vulnerabilities in the system that are not apparent when only considering annual reports of colony populations.

In past years of the survey, queen issues, Varroa destructor mites (Anderson and Trueman; hereafter “Varroa”), weather, starvation, and pesticides have been reported by beekeepers as important causes of colony losses (e.g., Bruckner, Wilson, et al., Citation2023). Researchers also view these as important causes of loss (Guzmán-Novoa et al., Citation2010). We now understand colony loss as a multifactorial process in which concurrent and interacting stressors can lead to colony losses (Potts et al., Citation2010; Steinhauer et al., Citation2018). The prevalence and impact of specific stressors, such as Varroa and viruses, can vary with time of year (Faurot-Daniels et al., Citation2020; Martin, Citation2001), and may change over time (Martin et al., Citation2012). Therefore, it is valuable to document what beekeepers perceive as the leading causes of loss at different times of the year and to continue this monitoring over a period of years.

It has been suggested that commercially managed colonies encounter stressful environments and beekeeping practices which could increase losses (Simone-Finstrom et al., Citation2016). Some specific stressors of concern are exposure to pesticides and acquisition of pests and pathogens during pollination events (Alger et al., Citation2018; Frazier et al., Citation2015), as well as stress during long-distance colony trucking (Ahn et al., Citation2012; Simone-Finstrom et al., Citation2016).

A majority of colonies in the U.S. are transported across state lines, including to California for almond pollination (Goodrich et al., Citation2019). In past years of the survey, no association has been found between almond pollination or migratory beekeeping and elevated winter losses, and occasionally these activities have been associated with reduced winter losses (Lee et al., Citation2015; vanEngelsdorp et al., Citation2010). However, colonies that are transported across state lines and used to pollinate California almonds constitute a large and economically important segment of the U.S. honey bee population (Goodrich et al., Citation2019). Therefore, we continue to explore the potential association of inter-state movement and almond pollination on winter losses to monitor for problems among colonies that are managed in these ways.

The purpose of this study was to continue the long-term documentation of managed honey bee colony losses in the United States by reporting the 2020–21 and 2021–22 summer, winter, and annual loss rates. We reported colony losses nationally, by operation type (backyard, sideline, and commercial), and by state. We also documented what beekeepers perceive as the causes of summer and winter loss, and what they consider an acceptable rate of winter loss. Finally, we investigated possible associations of winter loss with transportation of colonies across state lines, and with the use of colonies to pollinate California almonds.

Materials and methods

Survey

Following methods from previous years (Bruckner, Wilson, et al., Citation2023; Kulhanek et al., Citation2017; Lee et al., Citation2015; Seitz et al., Citation2015; Spleen et al., Citation2013; Steinhauer et al., Citation2014; vanEngelsdorp et al., Citation2007, Citation2008, Citation2010, Citation2011, Citation2012), the 2020–21 and 2021–22 BIP Loss and Management Survey was a retrospective annual survey of beekeepers’ colony losses and management practices. Participants were recruited directly via email (BIP mailing list, n = 26,404), and indirectly by asking well-connected partners (e.g., beekeeping journals, beekeeping clubs, state beekeeping associations, and national beekeeping organizations) to advertise the survey to beekeepers. Additionally, snowball recruitment was used, whereby respondents were invited to promote the survey among their contacts. All responses were accepted through an online portal: The 2020–21 survey was hosted on a third-party service (www.SelectSurvey.com), whereas the 2021–22 survey was hosted on a domain owned by the Bee Informed Partnership (www.beeinformed.org). In both years, the online portals accepted responses from 1 to 30 April. The survey consisted of two distinct sections: “Loss” and “Management.” Here we report on responses to the Loss section only. Other publications will report on beekeepers’ responses to the Management section.

As in previous publications, we categorized beekeepers into three operational types based on the number of colonies managed on 1 October of each survey year. “Backyard,” “sideline,” and “commercial” beekeepers were defined as those who managed 50 or fewer, 51–500, and more than 500 colonies, respectively. We refer to sideline and commercial beekeepers collectively as “large-scale” beekeepers. We also grouped respondents by whether they participated in specific beekeeping activities: moving colonies across state lines or using colonies for California almond pollination. Finally, we categorized respondents as single-state operations or multi-state operations depending on whether they managed colonies in one or several states, regardless of whether they reported moving colonies across state lines.

For each respondent, we determined the colony loss rate for three distinct seasons: summer, winter, and annual. Previously, the three seasons were defined as follows: summer (1 April to 1 October), winter (1 October to 1 April), and annual (1 April to 1 April) (e.g., Kulhanek et al., Citation2017). For backyard beekeepers, we retained these definitions because they appear to work well for the respondents. This also maintained consistency with past surveys. During focus sessions, large-scale beekeepers reported they often begin splitting colonies before the traditional 1 April start date for the survey’s summer season. Therefore, the 1 April cutoff date is not suited to their seasonal calendar. According to these beekeepers, the moment preceding splitting typically represents the annual minimum number of colonies for large-scale beekeepers, so they consider the start of splitting to be the start of the beekeeping season. To reflect this, we defined the seasons for large-scale beekeepers based on this key beekeeping management action. If a large-scale beekeeper made their first split of the calendar year before 1 April (i.e., between 1 January and 31 March), their summer season started on the date of their first split. If no splits were made between 1 January and 31 March, then their start of summer defaulted to 1 April. Regardless of the timing of the start of summer, 1 October was always considered the beginning of the winter season for all beekeepers.

To determine colony losses for each beekeeping operation during the summer, winter, and annual seasons, we asked the same core questions as previous BIP Loss and Management Surveys (e.g., Kulhanek et al., Citation2017; Seitz et al., Citation2015). In brief, respondents reported the number of colonies they owned at the beginning and end of each season, the net number of new colonies made or obtained from within and outside their operation, and the number of colonies they sold or gave away. This allowed us to calculate the number of colonies lost. We also calculated the number of colonies “at risk” for each operation and season (i.e., the maximum number of colonies that could be lost during a season). As in previous studies, the number of colonies at risk was calculated by taking the starting number of colonies, adding any new colonies made or obtained, and subtracting any colonies sold or given away by an operation during a season (Steinhauer et al., Citation2014). We also asked beekeepers to report what they thought were the most “prominent cause(s) of colony death” for the winter and summer seasons. Beekeepers selected up to three options for both “causes of summer loss” and “causes of winter loss.” Finally, we asked beekeepers to report what percentage of loss over the winter they would consider acceptable.

Following previously described methods, valid responses had to come from locations within the United States and have complete data to calculate the respondent’s colony loss rate for either the summer or winter season. Additionally, we used an algorithm to identify and manually review potential duplicate responses (e.g., if someone stopped and later restarted the survey). Since some responses were valid for only one season, we made three datasets—for valid summer, winter, and annual responses. The annual dataset contained responses that were valid for both summer and winter seasons. Data validation, filtering, and per-respondent calculations of colonies lost and at risk were performed in Python 3.8.13 (Python Software Foundation, Citation2022).

Statistics

We performed subsequent data preparation and statistical analyses in R 4.2.1 (R Core Team, Citation2021). We used tools in the tidyverse package for data preparation (Wickham et al., Citation2019), the stats and multcomp package in data analysis (Hothorn et al., Citation2008; R Core Team, Citation2021), the broom (Robinson et al., Citation2022) and gt (Richard Iannone et al., Citation2022) packages to prepare tables, and performed data visualization using the ggplot2 package (Wickham, Citation2016).

To aggregate data from individual respondents into loss estimates for the national honey bee population and loss estimates for subgroups of respondents (e.g., backyard beekeepers, beekeepers in a specific state, or almond pollinating beekeepers), we calculated the sum of the number of colonies lost for all beekeepers in that subgroup, and the sum of the number of colonies at risk for all beekeepers in that subgroup. To assign multi-state operations’ colonies to individual states (for state loss estimates), we used the same method as the USDA National Agricultural Statistical Service Honey Report (USDA-NASS, Citation2022c) and previous survey papers (Steinhauer et al., Citation2014), wherein multi-state operations’ colony numbers are counted in every state that they report keeping bees. However, we counted each colony only once for national estimates. We did not present results for states with fewer than 10 valid responses.

We calculated the loss rate for the national honey bee population and for each subgroup as the total number of colonies lost divided by the total number of colonies at risk. This corresponds to the “total loss” estimate (vanEngelsdorp et al., Citation2013) in previous surveys, in which each colony has the same weight in the estimated loss rate. We repeated this procedure to estimate the loss rate for each season we report on. We used a bootstrap method (n-out-of-n, rep = 1000) to obtain confidence intervals for these loss estimates (at the national level and for all smaller subgroups), which we abbreviate as “b.CI.” The bootstrap is a method that allows us to estimate the uncertainty around statistics inferred from a dataset, without having to make assumptions about their distribution. The standard (non-parametric) bootstrap samples the original n observations from a dataset with replacement, creating b (number of bootstrap resamples, typically over 1000) pseudo-datasets of size n, from which b statistics can be estimated. From the distribution of those estimates, we can calculate the bootstrapped version of the statistic (the mean of the bootstrapped statistics), as well as a nominal 1-α confidence interval for the original statistic (Chernick & LaBudde, Citation2011; Davison & Hinkley, Citation1997). We used the sample function from the base package (R Core Team, Citation2021) to resample our original dataset and create the bootstrapped datasets.

To compare colony loss between different subgroups of beekeepers (e.g., operation types, almond pollinators vs. not), we fitted generalized linear models with a quasibinomial distribution (quasibinomial glm) to account for overdispersion. As our general procedure to determine whether interaction effects or main effects of predictor variables were significant, we fitted a full model and a corresponding reduced model that omitted the term in question. We then performed an analysis of deviance (F test using the anova function) to determine whether the inclusion of the term was warranted, as one would when simplifying models by stepwise deletion (Crawley, Citation2013). For significance tests, we used an alpha level of 0.05.

To determine whether the national loss rate differed between summer and winter, we fitted a separate model for each year, with “season” as a main effect, and compared these models to intercept-only versions. We compared summer vs. winter losses but did not compare summer or winter losses vs. annual losses because an annual loss for a beekeeper is not independent of loss during summer or winter. To compare colony loss rates between operation types (backyard, sideline, commercial) within each season, we fitted a separate model for each combination of survey year and season (summer, winter, annual for each of 2020–21 and 2021–22) with “operation type” as a main effect. For pairwise comparisons between all three operation types within each year and season, we performed post-hoc tests using the glht function in the package multcomp (Hothorn et al., Citation2008) using the single-step method to adjust p-values. To investigate whether summer and winter losses differed within operation type, we fitted a new model for each combination of operation type and survey year (backyard, sideline, commercial for each of 2020–21 and 2021–22), including a “season” main effect. Again, we only compared summer against winter losses. To examine the association of colony migration or almond pollination with winter losses, we fitted a separate model for each year. We performed this analysis only on data for large-scale beekeepers (sideline and commercial) because few backyard beekeepers undertake these activities. Analyzing colony migration separately from almond pollination, we started with models that included main and interaction effects between “operation type” and “management practice” to allow for the possibility that the association between migration or pollination and colony losses could vary between sideline and commercial beekeepers. We did not find an interaction between operation type and migration on winter losses for either year (F = 0.0597, p = 0.808; and F = 0.6455, p = 0.423, respectively for the two years). The interaction between operation type and almond pollination on winter losses was not significant either year (F = 0.2194, p = 0.640; and F = 1.0673, p = 0.304, respectively). Therefore, we analyzed sideline and commercial beekeeper data together, and reported the association of each management practice and winter losses for large-scale beekeepers considered together.

Results

Respondents

Backyard beekeeper respondents were much more numerous in our pool of respondents than sideline and commercial beekeepers. We estimated that the 2020–21 and 2021–22 surveys accounted for 6.7 and 10.9%, respectively, of the honey-producing colonies in the United States on October 1 of that year (USDA-NASS, Citation2021, Citation2022b).

For the 2020–21 survey, 3588 responses were considered valid and included in the analytical dataset (Table S1). This includes responses considered valid for the summer season, the winter season, or both. Of the valid responses, 3444 responses (96.0%) were from backyard beekeepers, 103 (3.0%) from sideline beekeepers, and 41 (1.2%) from commercial beekeepers. On 1 October 2020, these backyard, sideline, and commercial beekeepers managed 23,821 (11.4%), 13,760 (6.6%), and 171,074 (82.0%) colonies, respectively.

The dataset for 2021–22 consisted of 3601 valid responses (Table S1). Of these, 3456 responses (96.0%) were from backyard beekeepers, 93 (2.6%) from sideline beekeepers, and 52 (1.4%) from commercial beekeepers. On 1 October 2021, backyard, sideline, and commercial beekeepers managed 23,319 (7.0%), 10,884 (3.2%), and 301,743 (89.8%) colonies, respectively.

Acceptable losses

In 2020–21, we had 3382 valid responses to the question about acceptable winter loss (3246, 97, and 39 responses from backyard, sideline, and commercial beekeepers, respectively). On average, beekeepers reported that they considered an acceptable level of winter loss to be 23.5% (22.8–24.2; 95% b.CI). The average backyard beekeeper reported 23.7% loss (23.0–24.4; 95% b.CI) to be acceptable. Sideline and commercial beekeepers reported very similar acceptable loss levels as each other—18.8% (16.1–21.7; 95% b.CI) and 18.8% (16.0–21.7; 95% b.CI), respectively.

In 2021–22, we had 3312 valid responses (3178, 90, and 44 responses from backyard, sideline, and commercial beekeepers, respectively). The average level of acceptable winter loss across all operation types was 20.9% (20.3–21.6; 95% b.CI). The average backyard beekeeper reported acceptable winter losses of 21.1% (20.4–21.8; 95% b.CI), while sideline beekeepers reported 18.1% (15.3–21.2; 95% b.CI), and commercial beekeepers reported 16.8% (14.8–19.2; 95% b.CI) as acceptable levels of winter loss.

National losses by season

In 2020–21, we estimated the national annual colony loss to be 50.8% (38.0–63.1; 95% b.CI) (; ). During this year, we recorded a 31.1% (19.2–44.1; 95% b.CI) summer loss and a 33.5% (28.0–40.7; 95% b.CI) winter loss (). Winter losses were significantly greater than summer losses (F = 18.194, p < 0.001), despite the small estimated difference of 2.4 percentage points.

Figure 1. Colony loss rates (%) of Apis mellifera honey bee colonies in the United States. Panes specify losses for (A) all operation types; (B) backyard beekeepers (1–50 colonies); (C) sideline beekeepers (51–500 colonies); and (D) commercial beekeepers (>500 colonies). The points represent the percentage of colonies lost, and the shaded bands and error bars represent the bootstrapped 95% confidence intervals (n-out-of-n, 1000 rep). Yellow, blue, and red represent summer, winter, and annual seasons, respectively. The survey year on the figure denotes the year the survey was conducted. For example, the 2022 survey year includes data for the summer of 2021 and the winter of 2021–22. Data for 2008–2020 are from previous surveys.

Figure 1. Colony loss rates (%) of Apis mellifera honey bee colonies in the United States. Panes specify losses for (A) all operation types; (B) backyard beekeepers (1–50 colonies); (C) sideline beekeepers (51–500 colonies); and (D) commercial beekeepers (>500 colonies). The points represent the percentage of colonies lost, and the shaded bands and error bars represent the bootstrapped 95% confidence intervals (n-out-of-n, 1000 rep). Yellow, blue, and red represent summer, winter, and annual seasons, respectively. The survey year on the figure denotes the year the survey was conducted. For example, the 2022 survey year includes data for the summer of 2021 and the winter of 2021–22. Data for 2008–2020 are from previous surveys.

Table 1. Colony loss rates (%) of Apis mellifera honey bee colonies in the United States for 2020–21 and 2021–22.

In 2021–22, our estimate of annual colony loss was 39.0% (31.5–47.9; 95% b.CI) (; ). The summer loss was 23.8% (16.3–31.9; 95% b.CI), and the winter loss was 24.2% (20.0–30.3; 95% b.CI), which were not significantly different (F = 2.177, p = 0.140).

Losses by operation type within season

During the 2020–21 survey year, backyard beekeepers’ annual losses (56.3%; 54.5–58.1; 95% b.CI) were significantly greater than those of commercial (50.2%; 33.1–63.6; 95% b.CI) or sideline beekeepers (46.8%; 38.3–54.6; 95% b.CI), while commercial and sideline losses were not significantly different from each other (; p-values in Table S2). During summer of 2020–21, commercial beekeepers had the highest estimated colony losses, followed by backyard and sideline (all significantly different). During winter, commercial beekeepers had the lowest losses, sideline losses were intermediate, and backyard beekeepers had the highest colony losses (p < 0.05 for all comparisons except sideline vs. commercial).

In 2021–22, commercial beekeepers had significantly lower annual losses (36.6%; 28.3–47.0; 95% b.CI) than either sideline (53.5%; 44.1–61.8; 95% b.CI) or backyard beekeepers (58.5%; 53.4–64.4; 95% b.CI). Sideline and backyard beekeepers’ losses were not significantly different (; p-values in Table S2). In the summer of 2021–22, commercial beekeepers had significantly lower losses than sideline or backyard beekeepers; sideline and backyard losses were not significantly different. During the winter of 2021–22, commercial beekeepers reported the lowest losses, sideline beekeepers reported intermediate losses and backyard beekeepers reported the highest losses (p < 0.05 for all comparisons).

Losses by season within operation type

In 2020–21, backyard beekeepers experienced 56.3% annual loss (summer 27.1%; winter 42.2%; ). Winter losses were significantly greater than summer losses (F = 411.3, p < 0.001). Sideline beekeepers’ annual losses in 2020–21 were 46.8% (summer 19.7%; winter 35.1%). Winter losses were significantly higher than summer losses (F = 20.048, p < 0.001). Commercial beekeepers in 2020–21 had an annual colony loss rate of 50.2% (summer 32.7%; winter 32.3.%). Summer and winter losses were not significantly different (F = 0.009, p = 0.926).

In 2021–22, backyard beekeepers reported a 58.5% annual loss (summer 35.2%; winter 37.2%). Winter loss was significantly higher than summer loss (F = 6.8367, p = 0.009). Sideline beekeepers reported annual losses of 53.5% in 2021–22 (summer 29.3%; winter 32.0%.). Winter and summer losses were not significantly different (F = 0.478, p = 0.490). Commercial beekeepers’ 2021–22 annual losses were 36.6% (summer 22.7%; winter 22.8%). The difference between summer and winter losses was not significant (F = 0.0197, p = 0.889).

State losses

Annual colony loss estimates for individual states for 2020–21 ranged from 20.7 to 85.8%, whereas annual state losses for 2021–22 ranged from 23.0 to 84.8% (; Table S3). When multi-state operations were considered as a subgroup, these operations had an annual loss of 46.5% in 2020–21 and 38.8% in 2021–22. Annual, summer, and winter losses for individual states and multi-state operations can be found in Table S3.

Figure 2. Colony loss rates (%) of Apis mellifera honey bee colonies in U.S. states and territories in 2020–21 and 2021–22. Orange and green colors represent loss rates above and below 30%, respectively. The grey color indicates that results are redacted for those states that had fewer than 10 respondents during a season. The survey year on the figure denotes the year the survey was conducted. For example, the 2022 survey year includes data for the summer of 2021 and the winter of 2021–22. MSO: multi-state operation; NCU: non-continental U.S. (Guam, Hawaii, Puerto Rico, and the U.S. Virgin Islands). Multi-state operations were counted in each state that they reported owning colonies but were also reported as a separate subgroup.

Figure 2. Colony loss rates (%) of Apis mellifera honey bee colonies in U.S. states and territories in 2020–21 and 2021–22. Orange and green colors represent loss rates above and below 30%, respectively. The grey color indicates that results are redacted for those states that had fewer than 10 respondents during a season. The survey year on the figure denotes the year the survey was conducted. For example, the 2022 survey year includes data for the summer of 2021 and the winter of 2021–22. MSO: multi-state operation; NCU: non-continental U.S. (Guam, Hawaii, Puerto Rico, and the U.S. Virgin Islands). Multi-state operations were counted in each state that they reported owning colonies but were also reported as a separate subgroup.

Self-reported causes of summer and winter loss

When beekeepers listed the “most prominent cause(s) of colony death” in summer of both years, “queen issues” and “Varroa” (“Varroa mites and associated viruses”) were their two most common responses, though the order varied by year and operation type (). In summer of 2020–21, queen issues were the most common response in all operation types. Backyard and sideline beekeepers reported queen issues followed by Varroa in summer of 2021–22, while commercial beekeepers reported Varroa followed by queen issues. Even though it was the third-ranked response among commercial beekeepers, we noticed that “adverse weather” was reported by 41% of commercial beekeepers in summer of 2021—a large proportion compared to the previous 3 years (Bruckner, Wilson, et al., Citation2023).

Figure 3. Causes of Apis mellifera honey bee colony loss as reported by United States beekeepers in 2020–21 and 2021–22. The graph shows the proportion of beekeepers who identified each cause of loss as being in the top three causes of loss for their operation. Operation types are backyard beekeepers (1–50 colonies), sideline beekeepers (51–500 colonies), and commercial beekeepers (>500 colonies). Causes of loss are reported separately for the summer (yellow) and winter (blue) seasons. The points represent the numerical proportion of beekeepers reporting the cause of loss, and the error bars represent bootstrapped 95% confidence intervals (n-out-of-n, 1000 rep). ApiTreat: “Apicultural treatments”; Disaster: “Natural disaster”; EquipFailure: “Equipment failure (e.g., moisture, ventilation); Pollen: “Nutritional stress (pollen deprivation); QueenIssue: “Queen issues”; Scavengers: “Scavenger pests”; Weather: “Adverse weather”; ShedFailure: “Failure of environmental controls in sheds (wintering facilities)”; DK: “Don’t know.”

Figure 3. Causes of Apis mellifera honey bee colony loss as reported by United States beekeepers in 2020–21 and 2021–22. The graph shows the proportion of beekeepers who identified each cause of loss as being in the top three causes of loss for their operation. Operation types are backyard beekeepers (1–50 colonies), sideline beekeepers (51–500 colonies), and commercial beekeepers (>500 colonies). Causes of loss are reported separately for the summer (yellow) and winter (blue) seasons. The points represent the numerical proportion of beekeepers reporting the cause of loss, and the error bars represent bootstrapped 95% confidence intervals (n-out-of-n, 1000 rep). ApiTreat: “Apicultural treatments”; Disaster: “Natural disaster”; EquipFailure: “Equipment failure (e.g., moisture, ventilation); Pollen: “Nutritional stress (pollen deprivation); QueenIssue: “Queen issues”; Scavengers: “Scavenger pests”; Weather: “Adverse weather”; ShedFailure: “Failure of environmental controls in sheds (wintering facilities)”; DK: “Don’t know.”

For winter 2020–21 and winter 2021–22, backyard, sideline, and commercial beekeepers all reported Varroa most commonly. For backyard beekeepers, weather was the second most common response in both years, while for sideline and commercial beekeepers queen issues were the second most common response in both survey years.

Losses by migratory practice and almond pollination

Large-scale beekeepers who moved their colonies across state lines had similar winter losses as those who did not—both in 2020–21 (F = 0.037, p = 0.848) and in 2021–22 (F = 0.0072, p = 0.933). Large-scale beekeepers who participated in California almond pollination had similar winter losses in 2020–21 as those who did not pollinate almonds (F = 0.0874, p = 0.768). In 2021–22, beekeepers who pollinated California almonds lost 23.1% of their colonies, while non-almond pollinators lost 39.3% of their colonies over winter. This estimated 16.2 percentage point difference was statistically significant (F = 6.1455, p = 0.015).

Discussion

Our findings for 2020–22 reaffirm the importance of monitoring colony loss rates in the United States and show that high losses still occur—though not every year. The national annual loss of 50.8% in 2020–21 was the highest estimated loss since 2011 when BIP started to monitor annual U.S. honey bee colony losses. In contrast, the estimated national annual loss rate in 2021–22 (39.0%) was comparable to previously published results (Bruckner, Wilson, et al., Citation2023; Kulhanek et al., Citation2017; Lee et al., Citation2015; Seitz et al., Citation2015; Spleen et al., Citation2013; Steinhauer et al., Citation2014; vanEngelsdorp et al., Citation2007, Citation2008, Citation2010, Citation2011, Citation2012). Backyard beekeepers’ summer losses appear to have risen since they were first surveyed: During the 2020–21 survey year, summer losses for the first time were nearly as high as winter losses. These elevated summer losses contributed to annual losses of above 55% for backyard beekeepers in 2020–21 and 2021–22—the highest annual losses for backyard beekeepers since annual losses were first monitored. We noted a marked increase in the number of commercial beekeepers reporting “weather” as a cause of loss for summer of 2021 and a decrease in the number of backyard beekeepers reporting “starvation” as a cause of winter loss compared to previous years. Nevertheless, the beekeepers overall reported leading causes of loss that align with recent years: Varroa and queen issues (Bruckner, Wilson, et al., Citation2023). While it appears that colony losses in the United States are higher than they have been in past decades, a continued limitation to comparing losses in the U.S. to other parts of the world is that many other survey efforts only report winter losses, whereas annual losses would better allow worldwide comparisons.

Beekeepers of all operation types continue to experience rates of winter colony loss that are above what they consider acceptable, and are above government targets. Winter losses in both years were substantially higher than what beekeepers consider to be acceptable, which aligns with previous survey results (Kulhanek et al., Citation2017). Almost a decade ago, the Pollinator Health Task Force (Citation2015) set a goal to reduce winter losses to no more than 15% by 2025. The winter losses of 34 and 24% estimated by this study are substantially above this target. Although we have not performed a full time-series analysis to analyze trends over time, we do not see evidence that honey bee colony loss rates have decreased since annual losses were first recorded in 2011, which indicates continued poor health of honey bee colonies in the U.S. High rates of loss by individual beekeepers can be costly and discouraging for those who keep bees for pleasure and can have grave economic consequences for professional beekeepers.

Not only was the national annual loss in 2020–21 the highest reported to date, but the estimated colony loss rate during each season during 2020–21 (summer 31.1%, winter 33.5%, annual 50.8%) was higher than the mean of previously reported losses for those seasons (summer 23.5%, winter 29.0%, annual 40.9%) (Bruckner, Wilson, et al., Citation2023; Kulhanek et al., Citation2017; Lee et al., Citation2015; Seitz et al., Citation2015; Spleen et al., Citation2013; Steinhauer et al., Citation2014; vanEngelsdorp et al., Citation2007, Citation2008, Citation2010, Citation2011, Citation2012). It appears that in 2020–21, a summer of high losses was followed by a winter of high losses, resulting in a high loss year. While the annual losses in 2020–21 were high, estimates from USDA-NASS indicate that the number of honey bee colonies in the U.S. decreased by 5% from April 2020 to April 2021 and then returned to previous levels (USDA-NASS, Citation2021, Citation2022b). From these data, there is no indication that the national honey bee population is in decline. It is important to note that our estimates of national colony loss rates and any estimate of the number of colonies in the U.S. mainly reflect the state of the commercial beekeeping sector because commercial beekeepers manage most colonies in the U.S.

Although commercial beekeepers keep most of the national population of colonies, most beekeepers in the U.S. are backyard and sideline beekeepers who report different seasonal patterns of colony loss than their commercial counterparts (e.g., Bruckner, Wilson, et al., Citation2023). Backyard beekeepers tend to have higher annual losses than commercial beekeepers (e.g., Bruckner, Wilson, et al., Citation2023), with sideline beekeepers often but not always having intermediate levels of loss. Furthermore, backyard beekeepers tend to lose substantially more colonies in the winter than they do in the summer, but this pattern changed in the years reported here. While backyard beekeepers’ summer losses in previous years were usually below 20% (Bruckner, Wilson, et al., Citation2023; Kulhanek et al., Citation2017; Lee et al., Citation2015; Seitz et al., Citation2015; Spleen et al., Citation2013; Steinhauer et al., Citation2014), in 2021 they had their highest summer loss so far (35%), and for the first time their summer losses were nearly as high as their winter losses. This apparent increase in summer mortality among backyard beekeepers deserves further investigation to see if it continues—and if so to determine what may be causing an increase in summer loss.

Backyard beekeepers as a group attributed summer and winter colony losses to a broader array of factors than did sideline and commercial beekeepers—and their leading responses included Varroa, queen issues, starvation, weather, and “don’t know.” Commercial beekeepers agreed more with one another in that they most commonly cited two causes of loss: Varroa and queen issues. The broader set of responses from backyard beekeepers may suggest that backyard colonies are subject to a broader set of stressors. Through actions that systematically manage stressors (e.g., sugar feeding) commercial beekeepers’ colonies may experience a reduced set of stressors. Supporting this idea, we know that commercial beekeepers tend to manage their colonies in ways that align more closely with expert recommendations (Steinhauer et al., Citation2021). Additionally, commercial beekeepers often winter their colonies in locations with a mild climate (or indoors) and likely develop their operations in locations that are favorable for colony survival. In contrast, the home location of a backyard beekeeper will strongly influence where they choose to keep their bees, presumably meaning that some of these locations are less suitable than those experienced by commercial colonies. If commercial beekeepers do in fact keep their colonies in more favorable locations, this along with more intensive management could lead to fewer major causes of loss and lead to lower loss overall among commercial beekeepers. Other survey efforts have also found that larger beekeeping operations tend to suffer lower losses (Brodschneider et al., Citation2016; Gray et al., Citation2023) and have also pointed to greater experience among larger operations. Others have also noted the fundamental connection between colony losses and operation size: operations with lower losses would be more able to grow larger. Although commercial beekeepers’ leading self-reported causes of summer loss were queen issues and Varroa, we noted that in summer of 2021, more of them attributed summer losses to “weather” than in the past. In the summer of 2021, 41% of commercial beekeepers cited weather as a cause of summer loss, compared to only about 10% in the preceding three years (Bruckner, Wilson, et al., Citation2023). This increased attribution of summer losses to weather coincides with the widespread drought conditions in the Northern Great Plains and Western U.S. (NOAA-NCEI, Citation2022)—important summertime beekeeping areas for migratory beekeepers (Smart et al., Citation2016). Recent studies further corroborate the potential importance of weather to both summer (Calovi et al., Citation2021) and winter (Insolia et al., Citation2022; Overturf et al., Citation2022) managed honey bee colony losses in the United States.

Our results show a broad consensus across all operation types of U.S. beekeepers that Varroa and associated viruses are the leading cause of colony mortality during winter, and this aligns with our understanding that Varroa infestation has various negative effects on individual bees (e.g., Bruckner, Straub, et al., Citation2023) and is an important driver of colony mortality (e.g., Chauzat et al., Citation2010; Genersch et al., Citation2010; Guzmán-Novoa et al., Citation2010; Steinhauer et al., Citation2018). However, beekeepers’ perception of Varroa as a cause of loss varies with operation type. About 60% of commercial and sideline beekeepers cited Varroa as a prominent cause of winter loss, but only 37–43% of backyard beekeepers did. This difference may indicate that backyard beekeepers are underestimating the effects of Varroa on their colonies. Because large-scale beekeepers more systematically manage non-Varroa stressors, it may be clearer to these respondents that many losses are caused by Varroa and associated viruses, which are difficult to control even with the tools that currently exist (Jack & Ellis, Citation2021). The emergence of amitraz resistance in Varroa (Rinkevich, Citation2020) may be further limiting the usefulness of the existing chemical varroacides, and is one possible contributor to the recent high losses. The report that Canadian beekeepers had an unusually high winter loss of 45.5% in the winter of 2021–22 (CAPA, Citation2022) indicates that the phenomenon of periodic elevated losses is not isolated to the United States and could implicate overlapping causes such as Varroa.

Like in previous BIP surveys, we did not detect an association between colony movement or almond pollination and increased winter losses among large-scale beekeepers. In fact, for the 2021–22 season, almond pollinators lost fewer colonies than non-pollinators. However, we note that we analyzed the association of these activities and losses only to rule out major negative effects of either management practice. A limitation of the single-variable modeling approach is that the actual effects of almond pollination may be confounded by specific characteristics of migratory operations or almond pollinators that affect colony losses but were not included in the analysis. Still, the results presented here and in previous surveys argue against catastrophic direct impacts of either inter-state movement or almond pollination. This suggests that the downsides of migration and almond pollination (e.g., shipping stress, pesticide exposure) (Ahn et al., Citation2012; Simone-Finstrom et al., Citation2016; Wade et al., Citation2019) are balanced by the benefits of migration and almond pollination (e.g., favorable climate, availability of plentiful forage), and/or by the benefits of other management that migratory and almond-pollinator beekeepers perform. One potential hypothesis for the large apparent benefit of almond pollination in 2021–22 is the warm weather favorable for pollination (and bee foraging) in the early part of the California almond bloom in 2022 (USDA-NASS, Citation2022a). If less favorable weather occurred in the locations where non-almond pollinator large-scale operations were based, this could have increased losses among those non-almond pollinators.

It is not straightforward to place U.S. colony losses in the context of present-day losses in other regions (Bruckner, Wilson, et al., Citation2023). To enable comparisons across countries with differing climates, we think it is essential for studies to report annual (year-long) losses, yet many survey efforts report only winter losses (e.g., Brodschneider et al., Citation2010; Gray et al., Citation2023; Zee et al., Citation2014). The focus only on winter losses probably reflects the dominance of winter losses in cooler areas of Europe. However, our results from the U.S. in this and previous BIP surveys show that summer losses can be substantial; this is corroborated by results from Uruguay that showed similar summer and winter losses (Antúnez et al., Citation2017). Additionally, results from a survey in South Africa showed high losses during a 6-month period characterized as the “active season” (Pirk et al., Citation2014), and results from a survey in Saudi Arabia found losses occurring during all seasons of the year (Albarrak & Gray, Citation2023). Even in Europe, we see some evidence that summer losses are currently underappreciated. While few studies report both summer and winter losses, a field study reported both percent summer loss and percent winter loss across 17 European countries ranging from Scandinavia to the Mediterranean (EPILOBEE Consortium et al., Citation2016). From these results, we calculated the ratio of winter losses to summer losses for each country and found a clear effect of latitude on the relative impact of summer and winter losses (Table S4). In the four countries with the highest relative impact of winter losses (Lithuania, Latvia, Poland, and Slovakia), winter losses were 12.9–72.3 times as high as summer losses, while in the four countries with the lowest relative impact of winter losses (France, Italy, Spain, and Greece), winter losses were estimated to be only 1.3–2.3 times as high as summer losses. These calculations may still underestimate the relative impact of summer losses in the Mediterranean region given that: (1) the study tracked already established colonies beginning in autumn (i.e., not nucs/splits), and (2) the colonies were monitored for an abbreviated summer period (i.e., not all the way to the subsequent autumn). Given these facts, we think monitoring colony losses only during the winter season may limit the comparison of colony loss rates even between countries within Europe. In the tropics and subtropics, or in areas with dry summers, phenomena such as precipitation may be more important in determining the active season for honey bee colonies (Requier et al., Citation2018), and in some cases, there is not one season in which losses are concentrated (Antúnez et al., Citation2017).

To bridge the differences in seasonality and in seasons during which the majority of losses occur, we think the best option is to record losses that occur over an entire year—accounting for increases and decreases from causes other than colony loss (Steinhauer et al., Citation2014). Asking beekeepers about their annual losses does not preclude gaining more seasonal resolution by also asking them about specific seasons of interest (e.g., summer and winter seasons in the BIP surveys), and it does not preclude defining the start or end dates in connection to biological events or management actions (e.g., end of winter preparations in the COLOSS surveys; the date that large-scale beekeepers made the first split of a calendar year in this BIP survey). When there may be large differences in the length of a season of interest, such as a beekeeper-defined winter season (van der Zee et al., Citation2012), it would aid interpretation of results if data on the start and end date of the season were collected and descriptive statistics on its length could be presented across different subgroups of respondents.

Conclusions

The 2020–21 and 2021–22 Bee Informed Partnership Loss and Management Survey revealed that high losses of honey bee colonies still occur in the United States. Despite public attention to the importance of pollinators and investment in honey bee research, loss rates do not indicate an improvement in the health of honey bees over the past 10–15 years. Observational studies like this survey cannot on their own conclusively demonstrate causal relationships, but this research points to multiple factors that beekeepers perceive as important causes of colony loss—Varroa, queen issues, and weather. We hope that research efforts to better understand these problems can lead to advances in honey bee breeding, product development, and management practices that can be broadly applied to improve colony health. We should continue to monitor colony losses to gauge progress on colony health.

Supplemental material

Supplemental Material

Download MS Word (59.3 KB)

Acknowledgements

We wish to express our gratitude to the beekeepers who responded to this survey. We are thankful to the following individuals at the Bee Informed Partnership who were essential to making this version of the survey a success: to Esteban Eliza for assisting with custom web development for the online forms at beeinformed.org; to Anne Marie Fauvel and the Tech Transfer Team Field Specialists for gathering responses from commercial beekeepers; and to Annette Meredith and the Board of Directors for supporting the survey as part of BIP’s mission. We also wish to thank two anonymous reviewers for their helpful comments.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

This work was supported by the Bee Informed Partnership, Inc.

References

  • Ahn, K., Xie, X., Riddle, J., Pettis, J., & Huang, Z. Y. (2012). Effects of long distance transportation on honey bee physiology. Psyche: A Journal of Entomology, 2012, 1–9. https://doi.org/10.1155/2012/193029
  • Albarrak, A. B., & Gray, A. (2023). Seasonal colony loss rates and honey bee management in the kingdom of Saudi Arabia: Results of a survey of beekeepers. Insects, 14(6), 513. https://doi.org/10.3390/insects14060513
  • Alger, S. A., Alexander Burnham, P., Lamas, Z. S., Brody, A. K., & Richardson, L. L. (2018). Home sick: Impacts of migratory beekeeping on honey bee (Apis mellifera) pests, pathogens, and colony size. PeerJ, 6, e5812. https://doi.org/10.7717/peerj.5812
  • Antúnez, K., Invernizzi, C., Mendoza, Y., vanEngelsdorp, D., & Zunino, P. (2017). Honey bee colony losses in Uruguay during 2013–2014. Apidologie, 48(3), 364–370. https://doi.org/10.1007/s13592-016-0482-2
  • Brodschneider, R., Gray, A., van der Zee, R., Adjlane, N., Brusbardis, V., Charrière, J.-D., Chlebo, R., Coffey, M. F., Crailsheim, K., Dahle, B., Danihlík, J., Danneels, E., de Graaf, D. C., Dražić, M. M., Fedoriak, M., Forsythe, I., Golubovski, M., Gregorc, A., Grzęda, U., … Woehl, S. (2016). Preliminary analysis of loss rates of honey bee colonies during winter 2015/16 from the COLOSS survey. Journal of Apicultural Research, 55(5), 375–378. https://doi.org/10.1080/00218839.2016.1260240
  • Brodschneider, R., Moosbeckhofer, R., & Crailsheim, K. (2010). Surveys as a tool to record winter losses of honey bee colonies: A two year case study in Austria and South Tyrol. Journal of Apicultural Research, 49(1), 23–30. https://doi.org/10.3896/IBRA.1.49.1.04
  • Bruckner, S., Straub, L., Neumann, P., & Williams, G. R. (2023). Negative but antagonistic effects of neonicotinoid insecticides and ectoparasitic mites Varroa destructor on Apis mellifera honey bee food glands. Chemosphere, 313, 137535. https://doi.org/10.1016/j.chemosphere.2022.137535
  • Bruckner, S., Wilson, M., Aurell, D., Rennich, K., vanEngelsdorp, D., Steinhauer, N., & Williams, G. R. (2023). A national survey of managed honey bee colony losses in the USA: Results from the Bee Informed Partnership for 2017–18, 2018–19, and 2019–20. Journal of Apicultural Research, 62(3), 429–443. https://doi.org/10.1080/00218839.2022.2158586
  • Calovi, M., Grozinger, C. M., Miller, D. A., & Goslee, S. C. (2021). Summer weather conditions influence winter survival of honey bees (Apis mellifera) in the northeastern United States. Scientific Reports, 11(1), 1553. https://doi.org/10.1038/s41598-021-81051-8
  • CAPA (2022). Statement on honey bee wintering losses in Canada (2022).
  • Chauzat, M. P., Martel, A. C., Zeggane, S., Drajnudel, P., Schurr, F., Clément, M. C., Ribière-Chabert, M., Aubert, M., & Faucon, J. P. (2010). A case control study and a survey on mortalities of honey bee colonies (Apis mellifera) in France during the winter of 2005–6. Journal of Apicultural Research, 49(1), 40–51. https://doi.org/10.3896/IBRA.1.49.1.06
  • Chernick, M. R., & LaBudde, R. A. (2011). An introduction to bootstrap methods with applications to R. Wiley.
  • Connor, L. J. (2014). Increase essentials: Nuclei, management, wintering (2nd ed.). Wicwas Press.
  • Crawley, M. J. (2013). The R book (2nd ed.). Wiley.
  • Davison, A. C., & Hinkley, D. V. (1997). Bootstrap methods and their application. Cambridge University Press.
  • EPILOBEE Consortium, Chauzat, M.-P., Jacques, A., Laurent, M., Bougeard, S., Hendrikx, P., & Ribière-Chabert, M. (2016). Risk indicators affecting honey bee colony survival in Europe: One year of surveillance. Apidologie, 47(3), 348–378. https://doi.org/10.1007/s13592-016-0440-z
  • Faurot-Daniels, C., Glenny, W., Daughenbaugh, K. F., McMenamin, A. J., Burkle, L. A., & Flenniken, M. L. (2020). Longitudinal monitoring of honey bee colonies reveals dynamic nature of virus abundance and indicates a negative impact of Lake Sinai virus 2 on colony health. PLOS One, 15(9), e0237544. https://doi.org/10.1371/journal.pone.0237544
  • Frazier, M. T., Mullin, C. A., Frazier, J. L., Ashcraft, S. A., Leslie, T. W., Mussen, E. C., & Drummond, F. A. (2015). Assessing honey bee (Hymenoptera: Apidae) foraging populations and the potential impact of pesticides on eight U.S. crops. Journal of Economic Entomology, 108(5), 2141–2152. https://doi.org/10.1093/jee/tov195
  • Genersch, E., von der Ohe, W., Kaatz, H., Schroeder, A., Otten, C., Büchler, R., Berg, S., Ritter, W., Mühlen, W., Gisder, S., Meixner, M., Liebig, G., & Rosenkranz, P. (2010). The German bee monitoring project: A long term study to understand periodically high winter losses of honey bee colonies. Apidologie, 41(3), 332–352. https://doi.org/10.1051/apido/2010014
  • Goodrich, B. K., Williams, J. C., & Goodhue, R. E. (2019). The great bee migration: Supply analysis of honey bee colony shipments into California for almond pollination services. American Journal of Agricultural Economics, 101(5), 1353–1372. https://doi.org/10.1093/ajae/aaz046
  • Gray, A., Adjlane, N., Arab, A., Ballis, A., Brusbardis, V., Bugeja Douglas, A., Cadahía, L., Charrière, J.-D., Chlebo, R., Coffey, M. F., Cornelissen, B., Costa, C. A. d., Danneels, E., Danihlík, J., Dobrescu, C., Evans, G., Fedoriak, M., Forsythe, I., Gregorc, A., … Brodschneider, R. (2023). Honey bee colony loss rates in 37 countries using the COLOSS survey for winter 2019–2020: The combined effects of operation size, migration and queen replacement. Journal of Apicultural Research, 62(2), 204–210. https://doi.org/10.1080/00218839.2022.2113329
  • Guzmán-Novoa, E., Eccles, L., Calvete, Y., McGowan, J., Kelly, P. G., & Correa-Benítez, A. (2010). Varroa destructor is the main culprit for the death and reduced populations of overwintered honey bee (Apis mellifera) colonies in Ontario, Canada. Apidologie, 41(4), 443–450. https://doi.org/10.1051/apido/2009076
  • Hothorn, T., Bretz, F., & Westfall, P. (2008). Simultaneous inference in general parametric models. Biometrical Journal, 50(3), 346–363. https://doi.org/10.1002/bimj.200810425
  • Iannone, R., Cheng, J., Schloerke, B., Hughes, E., & Seo, J. (2022). gt: Easily create presentation-ready display tables. Retrieved from https://CRAN.R-project.org/package=gt
  • Insolia, L., Molinari, R., Rogers, S. R., Williams, G. R., Chiaromonte, F., & Calovi, M. (2022). Honey bee colony loss linked to parasites, pesticides and extreme weather across the United States. Scientific Reports, 12(1), 20787. https://doi.org/10.1038/s41598-022-24946-4
  • Jack, C. J., & Ellis, J. D. (2021). Integrated Pest Management control of Varroa destructor (Acari: Varroidae), the most damaging pest of (Apis mellifera L. (Hymenoptera: Apidae)) colonies. Journal of Insect Science, 21(5), 1–32. https://doi.org/10.1093/jisesa/ieab058
  • Kulhanek, K., Steinhauer, N., Rennich, K., Caron, D. M., Sagili, R. R., Pettis, J. S., Ellis, J. D., Wilson, M. E., Wilkes, J. T., Tarpy, D. R., Rose, R., Lee, K., Rangel, J., & vanEngelsdorp, D. (2017). A national survey of managed honey bee 2015–2016 annual colony losses in the USA. Journal of Apicultural Research, 56(4), 328–340. https://doi.org/10.1080/00218839.2017.1344496
  • Lee, K. V., Steinhauer, N., Rennich, K., Wilson, M. E., Tarpy, D. R., Caron, D. M., Rose, R., Delaplane, K. S., Baylis, K., Lengerich, E. J., Pettis, J., Skinner, J. A., Wilkes, J. T., Sagili, R., & vanEngelsdorp, D. (2015). A national survey of managed honey bee 2013–2014 annual colony losses in the USA. Apidologie, 46(3), 292–305. https://doi.org/10.1007/s13592-015-0356-z
  • Martin, S. J. (2001). The role of Varroa and viral pathogens in the collapse of honey bee colonies: A modelling approach. Journal of Applied Ecology, 38(5), 1082–1093. https://doi.org/10.1046/j.1365-2664.2001.00662.x
  • Martin, S. J., Highfield, A. C., Brettell, L., Villalobos, E. M., Budge, G. E., Powell, M., Nikaido, S., & Schroeder, D. C. (2012). Global honey bee viral landscape altered by a parasitic mite. Science, 336(6086), 1304–1306. https://doi.org/10.1126/science.1220941
  • NOAA-NCEI (2022). Annual 2021 drought report. Retrieved from https://www.ncei.noaa.gov/access/monitoring/monthly-report/drought/202113
  • Overturf, K. A., Steinhauer, N., Molinari, R., Wilson, M. E., Watt, A. C., Cross, R. M., vanEngelsdorp, D., Williams, G. R., & Rogers, S. R. (2022). Winter weather predicts honey bee colony loss at the national scale. Ecological Indicators, 145, 109709. https://doi.org/10.1016/j.ecolind.2022.109709
  • Pirk, C. W. W., Human, H., Crewe, R. M., & vanEngelsdorp, D. (2014). A survey of managed honey bee colony losses in the Republic of South Africa–2009 to 2011. Journal of Apicultural Research, 53(1), 35–42. https://doi.org/10.3896/IBRA.1.53.1.03
  • Pollinator Health Task Force (2015). National strategy to promote the health of honey bees and other pollinators. The White House. Retrieved from https://obamawhitehouse.archives.gov/sites/default/files/microsites/ostp/Pollinator%20Health%20Strategy%202015.pdf
  • Potts, S. G., Biesmeijer, J. C., Kremen, C., Neumann, P., Schweiger, O., & Kunin, W. E. (2010). Global pollinator declines: Trends, impacts and drivers. Trends in Ecology & Evolution, 25(6), 345–353. https://doi.org/10.1016/j.tree.2010.01.007
  • Python Software Foundation (2022). Python 3.10.5 documentation (3.10.5). Retrieved from https://docs.python.org/3/
  • R Core Team (2021). R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. Retrieved from https://www.R-project.org/.
  • Requier, F., Antúnez, K., Morales, C. L., Aldea Sánchez, P., Castilhos, D., Garrido, P. M., Giacobino, A., Reynaldi, F. J., Rosso Londoño, J. M., Santos, E., & Garibaldi, L. A. (2018). Trends in beekeeping and honey bee colony losses in Latin America. Journal of Apicultural Research, 57(5), 657–662. https://doi.org/10.1080/00218839.2018.1494919
  • Rinkevich, F. D. (2020). Detection of amitraz resistance and reduced treatment efficacy in the Varroa mite, Varroa destructor, within commercial beekeeping operations. PLOS One, 15(1), e0227264. https://doi.org/10.1371/journal.pone.0227264
  • Robinson, D., Hayes, A., & Couch, S. (2022). broom: Convert statistical objects into tidy tibbles. Retrieved from https://CRAN.R-project.org/package=broom
  • Seitz, N., Traynor, K. S., Steinhauer, N., Rennich, K., Wilson, M. E., Ellis, J. D., Rose, R., Tarpy, D. R., Sagili, R. R., Caron, D. M., Delaplane, K. S., Rangel, J., Lee, K., Baylis, K., Wilkes, J. T., Skinner, J. A., Pettis, J. S., & vanEngelsdorp, D. (2015). A national survey of managed honey bee 2014–2015 annual colony losses in the USA. Journal of Apicultural Research, 54(4), 292–304. https://doi.org/10.1080/00218839.2016.1153294
  • Simone-Finstrom, M., Li-Byarlay, H., Huang, M. H., Strand, M. K., Rueppell, O., & Tarpy, D. R. (2016). Migratory management and environmental conditions affect lifespan and oxidative stress in honey bees. Scientific Reports, 6(1), 32023. https://doi.org/10.1038/srep32023
  • Smart, M. D., Pettis, J. S., Euliss, N., & Spivak, M. S. (2016). Land use in the Northern Great Plains region of the U.S. influences the survival and productivity of honey bee colonies. Agriculture, Ecosystems and Environment, 230, 139–149. https://doi.org/10.1016/j.agee.2016.05.030
  • Spleen, A. M., Lengerich, E. J., Rennich, K., Caron, D., Rose, R., Pettis, J. S., Henson, M., Wilkes, J. T., Wilson, M., Stitzinger, J., Lee, K., Andree, M., Snyder, R., & vanEngelsdorp, D. (2013). A national survey of managed honey bee 2011–12 winter colony losses in the United States: Results from the Bee Informed Partnership. Journal of Apicultural Research, 52(2), 44–53. https://doi.org/10.3896/IBRA.1.52.2.07
  • Steinhauer, N., Kulhanek, K., Antúnez, K., Human, H., Chantawannakul, P., Chauzat, M.-P., & vanEngelsdorp, D. (2018). Drivers of colony losses. Current Opinion in Insect Science, 26, 142–148. https://doi.org/10.1016/j.cois.2018.02.004
  • Steinhauer, N., Rennich, K., Wilson, M. E., Caron, D. M., Lengerich, E. J., Pettis, J. S., Rose, R., Skinner, J. A., Tarpy, D. R., Wilkes, J. T., & vanEngelsdorp, D. (2014). A national survey of managed honey bee 2012–2013 annual colony losses in the USA: Results from the Bee Informed Partnership. Journal of Apicultural Research, 53(1), 1–18. https://doi.org/10.3896/IBRA.1.53.1.01
  • Steinhauer, N., vanEngelsdorp, D., & Saegerman, C. (2021). Prioritizing changes in management practices associated with reduced winter honey bee colony losses for US beekeepers. The Science of the Total Environment, 753, 141629. https://doi.org/10.1016/j.scitotenv.2020.141629
  • USDA-NASS (2009). Honey report.
  • USDA-NASS (2020). Honey report.
  • USDA-NASS (2021). Honey bee colonies.
  • USDA-NASS (2022a). 2022 California almond forecast.
  • USDA-NASS (2022b). Honey bee colonies.
  • USDA-NASS (2022c). Honey report.
  • van der Zee, R., Pisa, L., Andonov, S., Brodschneider, R., Charrière, J.-D., Chlebo, R., Coffey, M. F., Crailsheim, K., Dahle, B., Gajda, A., Gray, A., Drazic, M. M., Higes, M., Kauko, L., Kence, A., Kence, M., Kezic, N., Kiprijanovska, H., Kralj, J., … Wilkins, S. (2012). Managed honey bee colony losses in Canada, China, Europe, Israel and Turkey, for the winters of 2008–9 and 2009–10. Journal of Apicultural Research, 51(1), 100–114. https://doi.org/10.3896/IBRA.1.51.1.12
  • van Dooremalen, C., Gerritsen, L., Cornelissen, B., van der Steen, J. J. M., van Langevelde, F., & Blacquière, T. (2012). Winter survival of individual honey bees and honey bee colonies depends on level of Varroa destructor infestation. PLOS One, 7(4), e36285. https://doi.org/10.1371/journal.pone.0036285
  • vanEngelsdorp, D., & Meixner, M. D. (2010). A historical review of managed honey bee populations in Europe and the United States and the factors that may affect them. Journal of Invertebrate Pathology, 103 Suppl 1, S80–S95. https://doi.org/10.1016/j.jip.2009.06.011
  • vanEngelsdorp, D., Caron, D., Hayes, J., Underwood, R., Henson, M., Rennich, K., Spleen, A., Andree, M., Snyder, R., Lee, K., Roccasecca, K., Wilson, M., Wilkes, J., Lengerich, E., & Pettis, J. (2012). A national survey of managed honey bee 2010–11 winter colony losses in the USA: Results from the Bee Informed Partnership. Journal of Apicultural Research, 51(1), 115–124. https://doi.org/10.3896/IBRA.1.51.1.14
  • vanEngelsdorp, D., Hayes, J., Underwood, R. M., & Pettis, J. (2008). A survey of honey bee colony losses in the U.S., fall 2007 to spring 2008. PLOS One, 3(12), e4071. https://doi.org/10.1371/journal.pone.0004071
  • vanEngelsdorp, D., Hayes, J., Underwood, R. M., & Pettis, J. S. (2010). A survey of honey bee colony losses in the United States, fall 2008 to spring 2009. Journal of Apicultural Research, 49(1), 7–14. https://doi.org/10.3896/IBRA.1.49.1.03
  • vanEngelsdorp, D., Hayes, J., Underwood, R. M., Caron, D., & Pettis, J. (2011). A survey of managed honey bee colony losses in the USA, fall 2009 to winter 2010. Journal of Apicultural Research, 50(1), 1–10. https://doi.org/10.3896/IBRA.1.50.1.01
  • vanEngelsdorp, D., Lengerich, E., Spleen, A., Dainat, B., Cresswell, J., Baylis, K., Nguyen, B. K., Soroker, V., Underwood, R., Human, H., Le Conte, Y., & Saegerman, C. (2013). Standard epidemiological methods to understand and improve Apis mellifera health. Journal of Apicultural Research, 52(1), 1–16. https://doi.org/10.3896/IBRA.1.52.1.08
  • vanEngelsdorp, D., Underwood, R., Caron, D., & Hayes, J. Jr. (2007). Estimate of managed colony losses in the winter of 2006–2007: A report commissioned by the Apiary Inspectors of America. American Bee Journal, 147, 599–603.
  • Wade, A., Lin, C.-H., Kurkul, C., Regan, E. R., & Johnson, R. M. (2019). Combined toxicity of insecticides and fungicides applied to California almond orchards to honey bee larvae and adults. Insects, 10(1), 20. https://doi.org/10.3390/insects10010020
  • Wickham, H. (2016). ggplot2: Elegant graphics for data analysis (2nd ed.). Springer-Verlag. https://doi.org/10.1007/978-3-319-24277-4
  • Wickham, H., Averick, M., Bryan, J., Chang, W., McGowan, L., François, R., Grolemund, G., Hayes, A., Henry, L., Hester, J., Kuhn, M., Pedersen, T., Miller, E., Bache, S., Müller, K., Ooms, J., Robinson, D., Seidel, D., Spinu, V., … Yutani, H. (2019). Welcome to the tidyverse. Journal of Open Source Software, 4(43), 1686. https://doi.org/10.21105/joss.01686
  • Zee, R., Brodschneider, R., Brusbardis, V., Charrière, J.-D., Chlebo, R., Coffey, M. F., Dahle, B., Drazic, M. M., Kauko, L., Kretavicius, J., Kristiansen, P., Mutinelli, F., Otten, C., Peterson, M., Raudmets, A., Santrac, V., Seppälä, A., Soroker, V., Topolska, G., Vejsnæs, F., & Gray, A. (2014). Results of international standardised beekeeper surveys of colony losses for winter 2012–2013: Analysis of winter loss rates and mixed effects modelling of risk factors for winter loss. Journal of Apicultural Research, 53(1), 19–34. https://doi.org/10.3896/IBRA.1.53.1.02