Effects of Field and Landscape Scale Habitat on Insect and Bird Damage to Sunflowers

Agriculture-dominated landscapes harbor significantly diminished biodiversity. Woody vegetation along field margins can provide farmers with ecosystem services and benefit biodiversity. However, when crops are damaged by the biodiversity harbored in such vegetation, farmers are reluctant to incorporate field margin habitat onto their land and may even actively remove such habitats. We investigated how damage by both insect pests (sunflower moth, Homoeosoma electellum) and avian pests to sunflower (Helianthus annuus) seed crops varied as a function of field-margin and landscape-scale habitat, as well as by bird abundance and diversity. Surveys for insect damage, avian abundance, and bird damage were carried out over 2 years in 30 different fields. The mean percentage of moth-damaged sunflowers sampled was nearly four times higher in fields that had bare or weedy margins (23.5%; $877/ha) compared to fields with woody vegetation (5.9%; $220/ha) and was positively associated with landscape-scale habitat complexity. Birds damaged significantly fewer sunflower seeds (2.7%) than insects, and bird damage was not affected by field margin habitat type, landscape-scale habitat variables, or avian abundance, but was significantly higher along field edges compared to ≥ 50 m from the field edge. Avian species richness nearly doubled in fields with woody margin habitat compared to fields with bare/weedy margins in both the breeding season and in fall. These results indicate that the benefits of planting or retaining woody vegetation along sunflower field margins could outweigh the ecosystem disservices related to bird damage, while simultaneously increasing the biodiversity value of intensively farmed agricultural landscapes.


INTRODUCTION
In the face of significant losses of both diversity and abundance of avian species (Rosenberg et al., 2019), farming agroecosystems represent a critical frontline for improving vast tracts of land for the conservation of biodiversity beyond the reserve system (Kremen and Merenlender, 2018;Grass et al., 2019). Establishing and protecting agroecosystems that harness functional diversity to provide ecosystem services at the farm and landscape level may simultaneously decrease chemical inputs and increase biodiversity (Daily et al., 2000;Bommarco et al., 2013;Weier et al., 2018;Kleijn et al., 2019). For example, establishing or maintaining strips of woody vegetation along field margins can increase the diversity, abundance, and corresponding ecosystem services of pollinators (Garibaldi et al., 2011;Sardiñas et al., 2016;M'Gonigle et al., 2015), arthropod predators (Eilers and Klein, 2009;Gareau et al., 2013), and birds (Heath et al., 2017;Gonthier et al., 2019). Likewise, higher amounts of natural habitat within agricultural landscapes (landscape-level complexity) can also increase biodiversity and associated ecosystem services (Chaplin-Kramer et al., 2011;Klein et al., 2012;Karp et al., 2018;Gonthier et al., 2019).
Increasing structural complexity within the working agricultural landscape can enhance the conservation value of larger tracts of protected land situated on the periphery of agricultural lands (Heath et al., 2017;Kremen and Merenlender, 2018). By integrating habitat complexity into agricultural landscapes, dispersal corridors, access to food sources, and greater genetic connectivity may occur across a wide diversity of taxa (Isbell et al., 2017). Such changes to the agricultural landscape could therefore have dual benefits for both landscape-scale biodiversity (DeClerck et al., 2010) and on-farm sustainability in the form of ecosystem services including biological pest control (Mitchell et al., 2013). However, natural habitat in agroecosystems can also harbor pests, which is true for avian pests of many seed and fruit crops worldwide (De Grazio, 1978;Gebhardt et al., 2011;Kross et al., 2012;Schäckermann et al., 2014). Farmers that perceive birds as detrimental to their crops will take action to deter birds (Kross et al., 2018), often by removing field margin habitat (Gennet et al., 2013) or utilizing commercially available bird deterrents such as gas guns, reflective tape, or netting (Baldwin et al., 2013), all of which can be costly for both farmers and non-target wildlife. Bird depredation of crops therefore not only has direct economic implications for growers, but can lead farmers to oppose conservation programs within agricultural communities and on their own properties (Kross et al., 2018).
Farmers are the primary decision makers for land management choices within agricultural regions, and their decisions, including those related to implementation of integrated pest management (IPM), are mostly based on direct economic returns (Kleijn et al., 2019). Ecosystem services, when enhanced by integrating habitat complexity into the working agricultural landscape, can serve as a part of IPM and contribute to a more sustainable and holistic preventative pest management paradigm (Stenberg, 2017). However, the effects of natural vegetation on biological control are complex and can vary with crop type, seasonality, farm management, and the demographic effects of interactions between natural enemies and pests (Karp et al., 2018;Settele and Settle, 2018). There is a strong need to provide clear, balanced empirical information to better inform habitat-modification based IPM strategies on the working, farm-scale level. This goal can be confused at times when, for example, studies into the detrimental behaviors of birds rarely focus on potentially beneficial impacts, and similarly, studies into beneficial pest-control services from birds rarely focus on damage that the same species may cause to crops (Pejchar et al., 2018; but see a few recent exceptions: Peisley et al., 2016;Gonthier et al., 2019). Therefore, disentangling the relationships between landscape-and field-level habitat complexity and crop damage from insect and avian pests has critical implications for both habitat management and preventative pest management in agroecosystems.
In California, one of the world's most productive and intensive farming regions, <4% of potential field margins have been planted with woody vegetation such as hedgerows (Brodt et al., 2009); field margins therefore have significant potential for increasing the biodiversity conservation value of California's working agricultural landscape. However, farmer surveys in California showed that major obstacles to hedgerow implementation included uncertainty around the potential ecosystem service benefits of hedgerows, along with concern that these hedgerows could harbor plant, insect, and vertebrate pests (Brodt et al., 2009). Research to provide empirical documentation of the costs and benefits of planting (and in some cases, retaining) such habitats is therefore critical to inform land management decisions.
Avian species and their abundance and diversity relationships to landscape structure can, in particular, create a still-unresolved conundrum in the analysis of costs and benefits to farmers (Pejchar et al., 2018). Local-and landscape-scale habitat influences both pests and their potential predators at varying scales and effect sizes (Karp et al., 2018). At the same time, some species of birds can provide ecosystem services, in the form of insect pest control, during the breeding season, and then can become major pests of the same crops when they switch to a more omnivorous diet in the fall (Figure 1).
In the present study, we analyzed the potentially conflicting roles of avian species within the unique sunflower (Helianthus annuus) seed growing region of California's Central Valley. We investigated the effects of both field-margin and landscape-scale habitat complexity on the occurrence of (1) potential benefits to farmers in the form of (A) avian insectivory leading to reduced occurrence of major invertebrate pests of sunflower seeds; and also (2) of potential costs to farmers in the form of (B) insect damage to sunflower seeds, and (C) bird damage to sunflower seeds.

Study Area and Crop
California's Central Valley is one of the most productive agricultural landscapes in the world, producing over 25% of the fresh produce consumed in the United States (USDA 2015), and valued at over $45 billion (USD) per year. Over 95% of the Central Valley's riparian and wetland ecosystems have been replaced by highly intensive agriculture and urban development (Katibah, 1984;Frayer et al., 1989), with remnant native habitat existing only in fragmented and isolated patches. Nevertheless, some native biodiversity in this region persists despite the highly human-modified landscape (Heath et al., 2017). FIGURE 1 | A conceptual framework for the potential effects of landscape-and local-scale habitat complexity on avian diversity (represented by bird silhouettes), and potential benefits (ecosystem services in the form of insect pest control) or costs (ecosystem disservices in the form of bird damage to crops), which can also vary with seasonality and diet, that affect yields for farmers.
Each year, sunflower is grown for hybrid seed production on an average of 20,234 ha (50,000 acres) across California's Sacramento Valley, producing over 31,750 tons, valued at ∼$70 million/year (Long et al., 2019). California's Central Valley produces over 95% of the United States' hybrid sunflower seeds, and over 25% of global sunflower seeds (Long et al., 2019). Sunflowers grown for seed are valued 5-10 times higher than the commercial oil crops for which they are used (Long et al., 2019), and growers therefore have a low threshold for damage. All sunflower fields in our study were grown for the same seed company and therefore were grown using the same standard sunflower production field-management practices (Long et al., 2019). This study was conducted within conventional fields (i.e., non-organic fields), but no growers reported utilizing insecticides on their fields over the duration of this study.
The sunflower moth (Homoeosoma electellum) is a major pest of sunflowers in North America and is the predominant insect pest of the crop in California (Long et al., 2019). Female sunflower moths lay eggs among the florets of sunflowers in early bloom, and eggs take 2-5 days to hatch. After hatching, larvae remain on the face of flowers for 8 days before boring into the developing seeds where they can cause losses of 30-60% of a crop (Long et al., 2019). Birds are the predominant vertebrate pest of sunflower crops around the world (De Grazio, 1978;Schäckermann et al., 2014;Ernst et al., 2019;Long et al., 2019). In North America, and in our study area, Icterid birds and the non-native European starling (Sturnus vulgaris) cause the most damage to sunflowers because of their omnivorous diets and behavior of foraging in large flocks. In Israel, bird damage to sunflowers is often concentrated to the edges nearest to habitat (specifically, large trees) that can act as shelter for birds, and can also be affected by landscape-scale habitat (Schäckermann et al., 2014), however flocking birds often cause damage within the interior of fields as well (e.g., European starlings in New Zealand vineyards; Kross et al., 2012).

Local-Scale Habitat Complexity
The presence or absence of woody habitat along field margins has been demonstrated to affect avian diversity and abundance in our study area (Heath et al., 2017), and at the time of our study the effects of habitat on sunflower moth damage was not described. To quantify the influence of habitat complexity on the bird community and sunflower damage, we included fields with woody margin habitat (n = 6 in 2014 and n = 12 in 2015) and fields with bare or weedy field margins in (n = 7 2014 and n = 5 in 2015), for a total of 30 fields (Figure 2). To quantify local (field) habitat complexity, at 5 evenly spaced locations along each 200-m long field margin transect (used for bird counts, see 2.6 below), we collected data on the maximum height of vegetation within 10 m of the of the transect, we estimated the number of canopy layers (out of a possible 7 predetermined canopy layers), and used satellite imagery to measure the width of field margin vegetation perpendicular to the focal field (see Heath et al., 2017 for details).
Because the variables describing field margin habitat (height, width, and number of vegetation layers) were highly correlated, we used a Principle Components Analysis (PCA) to reduce these into two orthogonal axes that explained over 95.5% of the variance among them. The two axes, PC1 and PC2, were included as predictor variables in our candidate models for sunflower damage and for bird abundance and richness. PC1 explained 86.2% of the variability among habitat variables and was negatively associated with all three variables, whereas PC2 was positively associated with habitat width and height, and negatively associated with habitat layers. Therefore, if PC1 is a positive predictor of damage, we would expect less damage at sites with habitat that is taller, wider and has more layers (because of the inverse relationship). If PC2 is a positive predictor of damage, we would expect less damage at sites with more habitat layers and more damage at sites with taller/wider habitat.

Landscape-Scale Habitat Complexity
To quantify and incorporate landscape habitat complexity into our study design, we selected fields at varying distances from natural habitat, which in our study area consists mainly of remnant and restored riparian areas (Figure 2). We used preexisting habitat data for our study area (CA DWR, 2008;Geographic Information Center, 2009), and added polygons for any trees within 800 m of each transect that were not included in the existing dataset (e.g., trees lining driveways, trees around homesteads). To calculate the distance to riparian areas, we first created a distance raster that encompassed the entire study area by using the Euclidean distance algorithm in ArcGIS 10.1 (ESRI, 2010). We used the riparian vegetation GIS dataset (habitats classified as native riparian, blue oak woodland, valley foothill riparian, fresh emergent wetland, saline emergent wetland, and valley foothill riparian) as the "source" input for the algorithm and set the output grid cell size to 10 m. Each field's transect center point was then buffered by 50 m, and we calculated the distance from each grid cell within the buffer to the nearest riparian vegetation polygon. The mean distance for all cells within each buffer was calculated as the distance value for each field. We also calculated the mean proportion area consisting of natural habitat at concentric buffer distances of 100, 200, 400, and 800 m, which have been shown to be relevant scales for riparian bird species in the Central Valley (Seavy et al., 2009).
To account for collinearity among landscape-scale habitat complexity variables, we ran separate models using the predictor variables for landscape-scale habitat complexity (either distance to the nearest natural habitat or percent natural cover at varying distances) for both the insect damage (Table S1) and the bird damage (Table S2) models.

Vertebrate Exclosures
In 2015, we created exclosures to prevent vertebrates (birds and bats) from accessing sunflowers (see Maas et al., 2019 for a review of exclosure methods). Exclosures consisted of nylon bird netting (No-Knot Bird Netting ¾ polypropylene mesh, Bird B Gone Inc R , Irvine, CA) draped over an area 4 sunflowers in width and ∼20 flowers in length (for a total of ∼80 plants) and secured to cover the tops of the flowers to a height of ∼2-4 feet above the ground. Birds were never observed inside our exclosures. Exclosures were installed in late spring, prior to the onset of bloom (which is when sunflower moth typically lay eggs on the flowers), and were checked and maintained over the entire growing season until final damage estimates were made in August-September. We set up four exclosures in each field, with the closest end of each exclosure located 5, 10, 50, and 100 m from the edge of the field. Due to last minute changes in the harvest schedule at some fields, we were able to collect damage data from the exclosures at nine different fields.

Sunflower Damage
We sampled from 10 sunflowers at distances from 0 to 200 m from the field edge. In 2014, we collected observations of both insect and bird damage from each site at 0, 10, 20, 30, 40, 50, 75, 100, 150, and 200 m from the field edge. In 2015, we collected observations from each site at 5, 10, 50, and 100 m from the field edge because we found in 2014 that bird damage dropped to close to 0 at distances beyond 50 m (see Figure 3C), and that insect damage was largely unchanged by distance from the field edge. Estimates for insect and bird damage in 2015 were taken from sunflowers within exclosures and from sunflowers that were ∼10 m from the exclosures (parallel to the field margin), but only data from non-enclosed sunflowers was used in our comparative analysis of insect damage.
We quantified both bird and insect damage by visually inspecting 10 sunflowers within each sampling area. Sunflowers were chosen by reaching out to select a plant stalk, so the seedbearing area of each plant was not seen until after the plant was selected. Observers moved a few steps along and between rows to select each new flower. Bird damage was characterized by missing seeds. We were careful to avoid classifying winddamaged seeds as bird damage. These seeds were generally removed from larger continuous areas of the sunflower head, whereas seeds removed by birds were in patchy sections or removed singularly. Insect damage was characterized by an area of visible frass (insect excrement and webbing) on the surface of multiple sunflower seeds. Seeds under the frass were often shrunken or visibly damaged. All areas that were under frass were classified as insect-damaged.
To estimate the percent of seeds on each sunflower that were damaged, we used a pre-cut circular piece of galvanized steel chicken-wire that was marked to allow for easy measurement of the flowers. Sunflower heads were classified into different size classes based on the diameter (to the nearest 1.3 cm, or 0.5 inches) of the seed-bearing area on each plant. We then estimated the number of hexagons on the wire (to the nearest ¼ hexagon) that was damaged by birds or damaged by insects on each sunflower head. Using the flower circumference and the known area within each hexagon of our grid, we were then able to calculate the percent of each sunflower head that was damaged by birds, and the total that was damaged by insects. For each sampling location, we aggregated the data from the 10 flowers for a single mean for percent damage. To estimate yield, damage from insects and damage from birds were summed for a total percent damage to each sunflower, since both types of damage result in a direct loss of yield for growers.

Economic Estimates
We used published data on mean sunflower yields and economic value for the Sacramento Valley from 2015 to 2018  to calculate the reduction in gross earnings for farmers as a result of insect and bird damage in response to significant predictor variables. Mean sunflower yields were 1,260 lbs/acre (1,412 kg/ha; range 1,076-1,748 kg/ha) after seed companies cleaned and removed non-viable seeds and non-seed material from field harvests . Seeds were valued at a FIGURE 3 | Model estimates (and 95% confidence intervals) of percent of sunflower seeds damaged (primary axis) as a function of the presence (darker colored lines) or absence (lighter colored lines) of woody vegetation along field edges and, (A) increasing field margin height and width (PC2), (B) the distance to the nearest natural habitats; and (C) percent seeds damaged by birds as a function of the distance of sampling points within each field from the nearest field margin. Raw data points are shown, with darker points indicating multiple overlapping points. The secondary axis for economic yield (right side of figure) applies to all 3 panels and shows the mean value (from 2015 to 2018) of sunflower seed crop in the region ($3,736 USD) equivalent to 0% damage, and decreasing as percent damage increases on the primary axis. mean value of $1.2/lb [$0.54/kg; range of $0.41-0.68/kg ]. We calculated the economic effect size of insect or bird damage by multiplying the scaled effect sizes from our model estimates with these mean values, assuming that yields were valued at $3,736/ha if they sustained no bird or insect damage.

Bird Counts
We conducted four bird surveys at each site, two in summer (June 9-July 2) and two in fall (August 5-September 16). All bird surveys were conducted by trained observers and timed to coincide with sunflower bloom in the summer (when sunflower moths typically lay eggs on the flowers), and immediately prior to the seed harvest in the fall. All counts were conducted between dawn and 10 a.m. and were not conducted in very cold (<3 C) or very hot weather (>24 C), in high winds, or in heavy precipitation. Counts were also re-scheduled if there were any farm workers or machinery in our focal field. The order in which fields were visited for counts was randomized within geographical groupings of fields. We conducted two counts per visit at each field: one to quantify the birds utilizing the field margin habitat, and another to quantify the birds utilizing the field interior. To count birds utilizing field margin habitat, observers walked a 200 m transect over 10 min, counting all birds detected by sight or sound within 20 m of the field margin, but not within the field itself. To count birds utilizing the field interior, observers returned to the mid-point of the transect, allowed 5 min for birds to settle, and then conducted a 10-min point count focused only on birds that were observed within the field (a half circle with an ∼200 m radius from our center point). We truncated the detection distance at field margins because margin habitat varied across sites. We assumed that intra-species detectability was the same within all sunflower fields, since sunflowers were at similar levels of maturation and height at the time of our surveys. Sunflowers are a tall (2-3 m) and densely planted crop, so most birds detected within the crop were of individuals flying into or out of fields, or singing/calling within a field.

Statistical Analyses
For all statistical models, we included as predictor variables in our maximal models the continuous variables for the distance from the nearest natural habitat (or proportion natural cover within concentric distance bands), PC1, PC2, as well as the categorical variable for whether the field had a weedy or bare edge (simple edge habitat) or had woody field margin habitat (complex edge habitat). For insect and bird damage to sunflowers, we also included the distance into the field the sample was collected (as both a linear and quadratic predictor, to account for potentially non-linear effects of distance on insect and bird damage). We simplified the maximal models by removing interactions, then main effects, until no further reduction in residual deviance (measured using Akaike's Information Criterion) was obtained. For all regression analyses, we considered candidate models with AIC ≤ 2 and chose the most parsimonious model. All data was analyzed in R (R Core Team, 2015).

Avian Consumption of Sunflower Moth
We used a Wilcoxon rank-sum test to compare the total insect damage observed inside exclosures and in adjacent nonexclosure locations.

Bird and Insect Damage to Sunflowers
Only the data from the non-exclosure sampling locations were used for investigating the effects of habitat variables on sunflower damage. Sunflower moth damage and bird damage were analyzed in separate models. We assessed our data for appropriate model distributions using quantile comparison plots ("qqplot" in Fox and Weisberg, 2019), then the fitdistr function (Venables and Ripley, 2002). We also used likelihood ratio tests to compare the final models to alternative likely distributions. For both insect and bird damage to sunflowers, we used generalized linear models with a negative binomial family of errors ("glm.nb" in MASS package; Venables and Ripley, 2002).

Seasonal Avian Species Richness and Abundance
We ran eight separate generalized linear regressions with a poisson family of errors for avian species richness and abundance along the field edge and within the field interior for data collected in summer and in fall ( Table S3).

Sunflower Damage
Sunflower moth damage was almost four-times higher at sites with bare or weedy field margin habitat (23.46 ± 1.41%) compared to sites with woody vegetation (5.89 ± 1.16%; z = 7.12, p < 0.001). There was a slight decrease in sunflower moth damage as habitat height and width (PC2) increased (z = −2.75, p = 0.005; Figure 3A). Model selection revealed that the variable for mean distance to natural habitat was the most parsimonious landscape-scale habitat variable in our insect damage models ( Table 1) and had a significant reduction in damage as distance from natural habitat increased (z = −2.25, p = 0.02; Figure 3B). The most parsimonious model for bird damage included the proportion of natural habitat within 800 m of the field as the landscape-scale habitat variable (Table 1) and had a marginally significant increase in bird damage (z = 1.96, p = 0.05). Bird damage was highest at the edge of fields, regardless of the presence of field margin habitat, and dropped quickly to near 0% within 50 m of the field edge ( Figure 3C). This effect was driven primarily by distance from field edge, with the linear (z = −4.38, p < 0.001) and quadratic values (z = 2.93, p = 0.003) for distance from field edge retained in the final model.

Economic Estimates
Our damage models estimate that the presence of woody field margin habitat results in significant changes to yield TABLE 1 | Model selection for candidate models explaining sunflower moth damage to sunflower seeds using the distance to nearest natural habitat as a measure of landscape-scale habitat complexity, and for candidate models explaining bird damage to sunflower seeds using the proportion natural habitat within 800m of each site. We used a principle components analysis to consolidate field margin habitat complexity, with PC1 negatively associated with field margin height, width, and number of canopy layers; and with PC2 positively associated with field margin height and width, and negatively associated with number of canopy layers. Field margins for each site were categorically defined based on the presence or absence of woody vegetation along the field margin. The "Distance into Field" measure is the number of meters within the field for each sampling location from the nearest field edge, and Distance into field 2 is the quadratic term for this variable.

Model parameters
and therefore economic value of sunflower crops. Multiplying these results by the mean value of sunflower seeds in the region ($3,736/ha) allows for a coarse estimate of the economic implications of damage to sunflowers. For example, at sites adjacent to natural vegetation, farmers would expect to lose $877/ha in yields due to sunflower moth damage at sites with bare/weedy vegetation along the field margin, compared to $220/ha in lost yields due to sunflower moth damage at sites with woody vegetation, but this difference would lessen with increasing distance from natural areas as overall insect damage also declined ( Figure 3B, secondary axis). To put this into perspective, the mean cost of applying insecticides to treat for sunflower moth is $292/ha (the equivalent of losses of ∼7.8% damage), so our results suggest that fields with more complex, woody margins would be likely to remain under an economic threshold that would trigger growers to apply insecticides, whereas sites with bare/weedy margins mostly incur damage above that threshold (Figures 3A,B, secondary axis). Bird damage at the field edge would result in $100/ha in lost yields but that would decline to negligible damage within 50 m of the field edge ( Figure 3C, secondary axis).

Avian Species Richness and Abundance
We observed 70 different avian species during our summer counts, and 74 species during our fall counts. These included California 'Bird Species of Special Concern ′ (Shuford and Gardali, 2008) like northern harrier (Circus hudsonius), yellow warbler (Setophaga petechia), and California "Threatened" species like Swainson's hawk (Buteo swainsoni), and tri-colored blackbird (Agelaius tricolor, 13 individuals observed at one site). During our summer counts, 64 different bird species utilized sunflower field edges and 49 species utilized field interiors. During our fall counts, we observed 69 species utilizing sunflower field edges and 46 species utilizing field interiors. The most abundant birds observed in fields interiors during the fall were Icterid species including Brewer's (Euphagus cyanocephalus) and red-winged blackbirds (Agelaius phoeniceus), European starling (Sturnus vulgaris), American goldfinch (Spinus tristis), house finch (Haemorhous mexicanus), and mourning dove (Zenaida macroura), all of which are species that are likely to damage sunflower seeds-although mourning doves most likely feed on seeds that have been knocked to the ground.

Summer Avian Species Richness and Abundance
For our summer counts, avian species richness ( Figure 4A) and abundance ( Figure 4B) along field edges were positively associated with increasing field margin habitat complexity ( Table 2, Table S4). Increasing field margin habitat was associated with lower avian species richness within fields in summer (PC1, Figure 4C, Table 2). For summer field interiors, avian species richness was negatively associated with increasing field margin habitat complexity (PC1, Figure 4C, Table 2). Avian abundance within field interiors in the summer was negatively associated with increasing height and width of field margin habitat (PC2, Figure 4D, Table 2). Model results for landscape scale habitat complexity, measured as the distance to nearest natural habitat, indicate that fields located further from natural habitat had lower avian species richness and abundance along field edges in the summer (but not field interiors, Table 2).

Fall Avian Species Richness and Abundance
In the fall, avian species richness ( Figure 4E) and abundance ( Figure 4F) along field edges, and avian abundance within field interiors ( Figure 4H) were positively associated with increasing field margin (local) habitat complexity (PC1; Table 2, Table S4). Avian species richness within field interiors in the fall was not   Table S5 for model selection), and blank spaces indicate parameters that were removed in model simplification. We report the intercept for each model, and the predicted difference per unit increase for each predictor variable (slope). Predictor variables include a categorical measure of "Field Margin" habitat (intercept level = woody vegetation), the results of a PCA for margin habitat complexity (PC1 negatively associated with field margin height, width, and number of canopy layers; and PC2 positively associated with field margin height and width, and negatively associated with number of canopy layers), and a measure of landscape-scale habitat complexity (distance to natural habitat).
significantly correlated with PC1 ( Figure 4G). However, avian abundance along field edges and within field interiors was negatively associated with PC2, which is associated positively with margin habitat height and width ( Table 2). Fall avian abundance was significantly lower along bare/weedy field edges (estimated mean reduction of 21.23 birds/transect) compared to along edges with woody vegetation (estimated mean of 44.39 birds/transect), but the opposite was true within field interiors, where there was an estimated increase of 100.59 more birds per transect at sites with bare/weedy edges compared to fields with woody vegetation along their edges (intercept = 56.63, Table 2). Only the abundance of birds within field interiors was significantly driven by landscape-scale habitat complexity, with fewer birds in fields further from natural habitat patches (at a rate of ∼2 birds per 100 m, Table 2). Avian species richness in field interiors was not correlated with any of our predictor variables for local or landscape scale habitat complexity.

DISCUSSION
Our results suggest that sunflower growers would benefit from planting or maintaining woody vegetation alongside their fields since insect damage was significantly higher at sites without field margin vegetation, while bird damage was not driven by field margin habitat. In fact, bird damage at our 30 fields was similar across sites with and without field margin habitat. Furthermore, within sunflower fields across all distances from the field margin, sunflower moth damage was significantly higher than bird damage, and was the main source of yield loss for sunflower growers in our area. The pest control service benefits that farmers receive from field margin vegetation therefore outweigh the potential ecological disservices associated with bird damage to sunflowers. Our results also indicate a clear benefit for biodiversity, with significantly higher species richness and avian abundance along field margins that had woody habitat. Combined, these results support the assertion that diversified farming systems can provide both farmers and broader society with multiple additive ecosystem services (Kremen and Miles, 2012). Our exclosures did not reveal an effect of bird foraging on sunflower moth damage. This could be the result of small sample size (n = 36 exclosures in 9 fields), or these results could indicate that foliage-gleaning birds and bats are not a major predator of sunflower moth. We hypothesize that the patterns of sunflower moth damage we observed were driven by either increased predation pressure from invertebrates, or from aerially-hunting bats and birds (which would not have been affected by the exclosures). Because of their nocturnal nature, adult sunflower moths are likely to be targeted more by nocturnal arthropod predators and/or bats (both of which would not be affected by the presence of exclosures) than by the predominantly diurnal avian predators. Studies in California have shown that the presence of habitat along field margins is associated with increased diversity and abundance of beneficial insects including natural enemies (Eilers and Klein, 2009;Gareau et al., 2013;Morandin et al., 2014), and with increased bat activity (Kelly et al., 2016), suggesting that our results could be driven by either or both of these groups of predators. Alternative hypotheses to explain our findings could be that woody vegetation along field margins present physical barriers to sunflower moths, or that increasing sunflower resources further from natural habitat dilute concentrations (and therefore damage) of sunflower moths, especially if natural habitats are the source of sunflower moths (e.g., see Tscharntke et al., 2016;O'Rourke and Petersen, 2017). Further research is clearly needed in this system.
Our results demonstrate that both insect and bird damage increased with landscape-scale habitat complexity. California's Central Valley is largely dominated by agriculture, with few corridors of remnant natural habitat along riparian areas (Figure 2). In this landscape, such corridors of natural habitat may be a source of migrating sunflower moths, and may also provide roosting habitat for large flocks of icterid birds in the autumn. Natural areas may therefore be a greater source of pests than they are a source of natural enemies (e.g., Hypothesis 2 in Tscharntke et al., 2016). However, our results show that regardless of how complex the landscape a farm is embedded within is, retaining or planting woody vegetation along field margins leads to a decrease in insect damage, and has no effect on total bird damage compared to sites with low levels of local habitat complexity. Importantly, this is also the scale at which farmers make decisions about planting, and therefore has the largest implications for rapid changes on private lands. Native hedgerows, the primary method for farmers to plant new woody vegetation along field margins, are also an important sources of pollination services for sunflowers (Sardiñas et al., 2016) and support pollinator metacommunity dynamics (Ponisio et al., 2019).
Habitat loss, largely a result of agriculture, is a primary driver of alarming trends of decreasing avian abundance in North America (Rosenberg et al., 2019). Both local (e.g., Hinsley and Bellamy, 2000;Batary et al., 2012;Heath et al., 2017;Gonthier et al., 2019) and landscape (e.g., Railsback and Johnson, 2014;Heath et al., 2017;Karp et al., 2018;Gonthier et al., 2019) habitat complexity have been linked to increased diversity of avian communities in farmlands, and may boost the conservation value of intensive agricultural landscapes. We found that avian species richness was positively associated with local-scale and landscapescale habitat complexity. The presence of woody vegetation also led to higher avian abundance along field edges, but was associated with decreased avian abundance within sunflower fields in the fall, although this did not drive a reduction in bird damage in our study.
Birds on farms can provide multiple, overlapping, and complex benefits and costs for farmers at multiple scales (Pejchar et al., 2018). Individual species can be beneficial to a crop in some seasons and detrimental in others, or may benefit one crop and cause damage to another. Birds may also disrupt other natural trophic cascades that benefit farmers (Grass et al., 2017). Importantly, while our results indicate a net benefit of woody vegetation along field margins for both sunflower yields and avian diversity in California, sunflowers in other regions (Peer et al., 2003;Schäckermann et al., 2014;Ernst et al., 2019) suffer from economically significant bird-damage to the same crop. Therefore, we caution that land managers and scientists should consider local climate, habitat availability, agricultural practices, and avian communities before translating our findings into management changes in other regions.
Increasing natural habitat in intensive agricultural landscapes can provide numerous ecosystem services and support biodiversity. However the risk of also attracting pests is a major cause of farmer reluctance to plant or retain such habitat (Brodt et al., 2009). Finding solutions that lead to landscapes that benefit both wildlife and farming is essential to ensuring food security and a thriving biodiversity in the future (Bommarco et al., 2013;Kremen and Merenlender, 2018). Our study demonstrates that while landscape habitat complexity leads to slight increases in insect and bird damage, but also increases bird diversity and abundance. Regardless of landscape-scale habitat, fields with local habitat complexity have higher yields compared to fields with bare or weedy edges, and fields with woody vegetation along field margins also harbor a greater diversity and abundance of birds. To move into a more preventative pest management paradigm, encouraging farmers in California to plant or retain woody vegetation along field margins will simultaneously increase sunflower seed yields and increase the diversity and abundance of birds.

DATA AVAILABILITY STATEMENT
Data will be made available on the Columbia University Library Digital Repository, https://doi.org/10.7916/d8-k1sp-wd06.

ETHICS STATEMENT
The animal study was reviewed and approved by Institutional Animal Care and Use Committee (IACUC #18033) at the University of California, Davis.

AUTHOR CONTRIBUTIONS
SK, TK, and JT conceived the ideas and designed methodology. SK, BM, and RB field assistants collected the data. CM performed the landscape analysis. SK analyzed the data and led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.

FUNDING
Fieldwork was funded by the David H. Smith Conservation Research Fellowship to SK, who was hosted by J. Eadie at UC Davis.