Introduction
The United States is ranked fifth among countries producing pumpkin, squash, and gourds (FAO 2022). In 2020, pumpkin production in the United States totaled $194 million; Indiana ranked fifth among the top pumpkin-producing states, with 2,428 ha valued at approximately $16 million (USDA NASS 2021). Although production practices can vary widely, pumpkins are usually direct-seeded into bare-ground rows placed 1.2 to 1.8 m apart. Pumpkins can be bushy or vining. In-row seed spacing is determined based on this distinction and ranges from 46 to 240 cm (Phillips Reference Phillips2021). The wide row and plant spacing required for this crop’s growth allows weeds to establish easily.
Weeds can reduce pumpkin yield by as much as 67% (Walters and Young Reference Walters and Young2010). Common, difficult-to-control weeds in Midwestern cucurbit production are Eastern black nightshade (Solanum ptychanthum Dunal), horseweed (Erigeron canadensis L.), morningglory spp. (Ipomoea spp. L.), pigweeds (Amaranthus spp. L.), giant ragweed (Ambrosia trifida L.), wild buckwheat [Fallopia convolvulus (L.) Á. Löve], Canada thistle [Cirsium arvense (L.) Scop.], dandelion (Taraxacum officinale F.H. Wigg.), field bindweed (Convolvulus arvensis L.), Johnsongrass [Sorghum halepense (L.) Pers.], and yellow nutsedge (Cyperus esculentus L.) (IPMdata 2005). Weed management in conventional pumpkin production generally includes chemical control. Preemergence herbicides combined with shielded postemergence row-middle application of nonselective herbicides are those most often used by pumpkin producers (Phillips Reference Phillips2021).
Herbicides can significantly reduce production costs by helping farmers overcome labor scarcity and elevated costs associated with other weed management practices. It has been estimated that US crop production would decrease by 20% without herbicides (Gianessi and Reigner Reference Gianessi and Reigner2007). Farmers can only use state-registered herbicides for tolerant crops. In Indiana, few herbicides have been registered for use in pumpkin production (Phillips Reference Phillips2021), so farmers have to rely on the same herbicides year after year, a practice that can contribute to herbicide resistance (Evans et al. Reference Evans, Tranel, Hager, Schutte, Wu, Chatham and Davis2016; Gressel Reference Gressel, Casley, Cussans and Atkin1991).
Fomesafen is a Group 14 herbicide that inhibits protoporphyrinogen oxidase. It is registered with 24C Special Local Need labels for use preemergence after pumpkin seeding but before crop emergence in Illinois, Kansas, Michigan, Minnesota, and Ohio, where it successfully controls several of the problematic weeds found in pumpkin production. However, it is not registered in Indiana (Phillips Reference Phillips2021). Farmers have noted this inconsistency in extension meetings and want fomesafen to be registered in Indiana as well. In-state tolerance data are desirable for an herbicide to be registered with a 24C label. A crop is considered tolerant when the applied herbicide does not cause any toxicity (Pitty Reference Pitty1995) or when it shows some injury but completely recovers by the end of its growing cycle (Seefeldt et al. Reference Seefeldt, Jensen and Fuerst1995).
Dose-response studies can be used to derive a model from the biological effect of an herbicide, or multiple herbicides on a crop, or multiple crops (Streibig Reference Streibig1980). Dose-response curves are often sigmoidal and constrained by an upper and a lower limit. The upper and lower limits are defined by the response from nontreated plants (control) and the highest dose applied (Knezevic et al. Reference Knezevic, Streibig and Ritz2007). Our objective was to fit fomesafen dose-response curves to evaluate the biological response of two pumpkin cultivars. With this, we can determine possible outcomes regarding crop tolerance at other fomesafen rates within the range of rates used in our study.
Materials and Methods
Fomesafen dose-response field trials were conducted in 2020 and 2021 at the Southwest Purdue Agricultural Center (SWPAC), Vincennes, IN (38.73°N, 87.48°W) and in 2021 at the Pinney Purdue Agricultural Center (PPAC), Wanatah, IN (41.44°N, 86.93°W). At SWPAC, soil types were a Conotton gravelly loam (loamy-skeletal, mixed, active, mesic Typic Hapludalfs, 65% sand) with 0.8% organic matter (OM) and pH 6.6 in 2020, and a mixture of Lomax loam (coarse-loamy, mixed, superactive, mesic Cumulic Hapludalfs, 32% sand) and Lyles sandy loam (coarse-loamy, mixed superactive, mesic Typic Endoaquolls, 67% sand) with 0.9% OM and pH 6.4 in 2021. At PPAC, the soil type was a mixture of Tracy sandy loam (coarse-loamy, mixed, active, mesic Ultic Hapludalfs, 59% sand) and Bourbon sandy loam (coarse-loamy, mixed, active, mesic Aquultic Hapludalfs, 56% sand) with 1.7% OM and pH 6.8.
Fields were prepared with tillage prior to the formation of raised beds. Raised beds with subsurface drip tape were prepared on June 17, 2020 and June 15, 2021 at SWPAC and June 2, 2021, at PPAC. The experimental design was a randomized complete block design with a split-plot treatment arrangement and four replications. The main plots consisted of the fomesafen rate and the subplots of the pumpkin cultivar randomly placed within each main plot. Subplots were 27 m2 and contained three 4.9-m-long rows, 1.8 m apart. Fomesafen (Reflex®; Syngenta Crop Protection, LLC, Greensboro, NC) rates were 0, 280, 560, 840, and 1,120 g ai ha–1, where 0 g ai ha–1 was the nontreated control. Pumpkin cultivars were ‘Bayhorse Gold’ and ‘Carbonado Gold’ (Rupp Seeds, Inc., Wauseon, OH). Crop fertilization, irrigation, and diseases and insect management followed recommendations by Phillips (Reference Phillips2021).
In each 27-m2 subplot, two pumpkin seeds were hand-planted into the same hole 1.2 m apart in-row, totaling 24 seeds per subplot, on June 18, 2020 and June 16, 2021 at SWPAC, and on June 2, 2021 at PPAC. Two pumpkin seeds were planted per planting hole in the event that one seed failed to germinate. Subplots were thinned to one plant per hole 2 to 4 wk after planting. Fomesafen was broadcast-applied on top of the bed and respective row middles (Figure 1A) within 1 d of planting. To help manage weeds, S-metolachlor (Dual Magnum®; Syngenta Crop Protection, LLC, Greensboro, NC) at 1,070 g ai ha–1 was broadcast-applied in a separate application across all plots within 1 d after applying fomesafen. At SWPAC, both herbicides were applied using a tractor-mounted PTO-driven Hypro 7560 C roller pump with four TeeJet XR 8003 VS nozzles (Spraying Systems Co., Wheaton, IL) calibrated to deliver 187 L ha–1 at 207 kPa. At PPAC, fomesafen was applied using a CO2-pressurized backpack sprayer equipped with four TeeJet XR 11004 VS nozzles calibrated to deliver 187 L ha–1 at 165 kPa and S-metolachlor was applied using PTO-driven Hypro model 6500 C roller pump with four TeeJet XR 8003 VS nozzles calibrated to deliver 187 L ha–1 at 138 pKa.
Data collection included counting the number of emerged pumpkin plants out of the 24 seeds that were planted 2 wk after treatment (WAT). Visible crop injury was rated using a scale of 0 (no injury) to 100% (crop death) at 2, 4, 6, and 8 WAT. Weed control was rated 4 WAT on a scale of 0 (no control) to 100% (complete control) relative to the 0 g ha–1 fomesafen treatment. After the 4-WAT weed control rating, weeds were removed either by hand or with hoes or cultivators to maintain plots weed-free and avoid yield loss due to weed interference. Pumpkin harvest was performed on September 12, 2020 [86 d after planting (DAP)] and September 17, 2021 (93 DAP) at SWPAC, and on September 1, 2021 (91 DAP) at PPAC. All fruits were harvested from each plot, individually weighed, and the color of each fruit recorded. A fruit was classified as marketable if it weighed ≥1.5 kg. Marketable fruits were categorized as orange (≥50% of the surface area was orange), green (<50% of the surface area was orange), and immature (green, tender rind). Individual marketable fruit weight average was calculated by dividing the marketable yield by the marketable fruit number of each category.
Emergence and marketable orange pumpkin yield and fruit number data were converted to a percent of the nontreated control using Equation 1:
where M was the average of the nontreated control variable value pooled across the four repetitions within a location-year for each pumpkin cultivar and B was the variable value of each data point for each location-year.
Data were subjected to statistical analysis using R software (RStudio®; PBC, Boston, MA). Data were first analyzed for each location-year with a linear model and subjected to ANOVA to determine if the models were statistically significant for each trial. If models were significant, data were combined across all three or only two location-years to check if the normality of the data was affected and determine if statistically interactions (P ≤ 0.05) existed between fomesafen rate, pumpkin cultivar, and location-year for each response variable. If data normality was affected or interactions between the explanatory variables existed, data are presented separately. Response variables were emergence as a percent of the nontreated control, visible pumpkin injury at 2, 4, and 6 WAT, weed control 4 WAT, marketable orange pumpkin yield and fruit number as a percent of the nontreated control, and marketable pumpkin yield (kg 27 m–2), fruit number, and average individual fruit weight (kg fruit–1) for the green and immature fruits. Visible pumpkin injury and weed control data were arcsin-squareroot transformed for analysis and are presented as back-transformed data. Data from the nontreated check were excluded from the visible pumpkin injury and weed control data analysis due to zero variance.
Significant response variables’ models were then subjected to nonlinear regression analyses using the package drc in R software and fit to either a three-parameter log-logistic model using Equation 2:
where d is the upper limit, b is the growth rate, e is the inflection point, and x is the fomesafen rate in g ai ha–1, or a three-parameter logistic model using Equation 3:
where d is the upper limit, b is the relative slope, e is the inflection point, and x is the fomesafen rate in g ai ha–1. Nonlinear models fit were analyzed with a lack-of-fit test, where a P > 0.05 indicates that the nonlinear model provides adequate description of the data. If data did not fit a model, a Tukey’s HSD means separation test was performed at a P ≤ 0.05 significance level.
Results and Discussion
Pumpkin Emergence
The SWPAC-2020 emergence as a percent of the nontreated control data was separated from the 2021 data because there was significant fomesafen rate-by-location-year interaction (F8, 89 = 7.32, P = 1.98 × 10–7) when pooled across all three location-years. However, emergence data from both locations in 2021 were pooled. Data were pooled across cultivars because there were no significant fomesafen rate-by-cultivar interactions. A three-parameter log-logistic model (Equation 2) was fit to the SWPAC-2020 and the pooled 2021 data. At SWPAC-2020, as fomesafen rate increased from 280 to 1,120 g ha–1, predicted emergence decreased from 85% to 25% of the nontreated control at SWPAC-2020, but only from 99% to 74% in 2021 at both locations (Figure 2). After thinning the subplots to one plant per planting hole, the nontreated subplots had 11 (SWPAC-2020) or 12 (SWPAC and PPAC-2021) plants per subplot. At SWPAC-2020, as the fomesafen rate increased from 280 to 1,120 g ha–1, the pumpkin density decreased to 8 (280 and 560 g ha–1), 4 (840 g ha–1), and 3 (1,120 g ha–1) plants per subplot. At the other two location-years, subplots treated with fomesafen averaged 11 plants per subplot.
We attributed the reduction in emergence to excessive rainfall. Cumulative rainfall within 2 WAT in 2020 was 120 mm, but in 2021 it rained only 44 mm at SWPAC and 17 mm at PPAC (Table 1). Soil-applied herbicide uptake happens mainly in the root for dicotyledonous plants via diffusion, interception, or mass flow. Herbicide uptake via mass flow, the process where the herbicide moves due to the hydrostatic gradient, accounts for the majority of the herbicide uptake (Menendez et al. Reference Menendez, Rojano-Delgado and De Prado2014). Rain is necessary to incorporate preemergence herbicides into the soil profile (Figure 1B). However, excessive rain moves the herbicide deeper in the soil profile into the crop’s root zone (Figure 1C), enhancing the hydrostatic gradient, thus increasing herbicide absorption. Fomesafen is highly mobile under water-saturated soil conditions, especially in soils with low OM content, high pH, and a high proportion of sand content (Guo et al. Reference Guo, Zhu, Shi and Sun2003; Li et al. Reference Li, Grey, Price, Vencill and Webster2019; Weber et al.Reference Weber, Strek and Sartori1993, Reference Weber, Wilkerson and Reinhardt2004). Low soil OM content and the excessive rain through the first 2 wk at SWPAC-2020 increased fomesafen available for uptake in the crop root zone in an early, vulnerable stage, thus reducing emergence.
a Data from the Midwest Regional Climate Center, West Lafayette, IN.
Peachey et al. (Reference Peachey, Doohan and Koch2012) reported that fomesafen at 560 g ha–1 reduced the emergence of ‘Eureka’ cucumber (Cucumis sativus L.), ‘Golden Delicious’ Hubbard squash (Cucurbita maxima Duchesne), ‘Dickinson’ pumpkin, and ‘Ultra’ butternut winter squash (Cucurbita moschata Duchesne ex. Poir.), and ‘Elite’ zucchini, ‘Yellow Crookneck’ summer squash, and ‘Small Sugar’ pumpkin (C. pepo L.) on average from 2.8 to 2.1 plants m–1 (25% reduction). However, the pumpkin cultivar ‘Small Sugar’ emergence was not affected by the fomesafen rate of 280 g ha–1 and was reduced only by 8% at the fomesafen rate of 560 g ha–1. Our results differ from their result in that we found a somewhat wide range of emergence reduction (5% to 21%) even at the lowest fomesafen rate of 280 g ha–1, presumably because of the soils’ OM content. Peachey et al. (Reference Peachey, Doohan and Koch2012) reported OM content ≥2.1% for all soil types, which could have increased fomesafen sorption to the soil. Also, other environmental conditions like rainfall must be taken into consideration. Similar to our results, Ferebee (Reference Ferebee2018) reported that fomesafen at 280 g ha–1 reduced the plant stand of ‘Kratos’ (C. moschata) and ‘Cougar’ (C. pepo L.) pumpkin by 20% and 63% to 75%, respectively. These trials were also conducted in soils with low OM content (1%). Likewise, they attributed the reduction in plant stand to rainfall (7 to 26 mm) shortly after planting (2 to 3 DAP).
Pumpkin Injury
We observed necrosis (Figure 3A), small white and brown spots (Figure 3B, C), chlorosis (Figure 3D, E), and stunting injury (Figure 4). Injury data were analyzed separately by location-year because of a significant fomesafen rate-by-location-year interaction. Injury was pooled across both cultivars in each location-year because of a nonsignificant fomesafen rate-by-cultivar interaction. Injury data 4 WAT at SWPAC-2020 and PPAC-2021 were fit a three-parameter logistic model (Equation 3; Figure 5). All other injury data were subjected to a Tukey’s HSD mean comparison test (Table 2). At 2 WAT, as the fomesafen rate increased from 280 to 1,120 g ha–1, injury increased from 6% to 28% at SWPAC-2021 and 5% to 36% at PPAC-2021 (Table 2). At 4 WAT, predicted injury increased from 7% to 26% (SWPAC-2020) and 4% to 50% (PPAC-2021) (Figure 5), and observed injury data at SWPAC-2021 increased from 0 to 13% (Table 2). At 6 WAT, injury ranged from 1% to 11% at SWPAC-2020 and 2021 and from 1% to 21% at PPAC-2021 (Table 2). Injury 8 WAT increased from 0 to 4% at SWPAC-2021 and 5% to 33% at PPAC-2021 (Table 2). With the exception of PPAC-2021, injury decreased from 2 to 8 WAT.
a Injury was arcsin transformed for analysis and back-transformed for the table. Scale: 0% = no injury, 100% = crop death.
b Data for SWPAC-2020 and PPAC-2021 at 4 WAT were fit a three-parameter logistic model (Figure 5).
c Means separation using Tukey’s HSD test P ≤ 0.05. Means followed by the same letter are not significantly different. Lack of letters indicates that the F statistic was not significant at a α = 0.05.
Heavy rainfall events increase the chance of injury due to the splashing of fomesafen from the soil onto the leaves (Peachey et al. Reference Peachey, Doohan and Koch2012). This could explain the necrosis and chlorosis scattered patterns on cotyledons and leaves lying close to the ground. Injury 8 WAT was mainly stunting. Lingenfelter and VanGessel (Reference Lingenfelter and VanGessel2016) also mentioned stunting up to 8 WAT with fomesafen applied at 175 and 350 g ha–1 to five pumpkin cultivars (C. pepo and C. maxima). As mentioned before, fomesafen persistence in the soil varies with OM, sand content, and pH. For this reason, fomesafen half-life values in diverse soil types ranged variably from 4 to 66 d (Li et al. Reference Li, Grey, Price, Vencill and Webster2019; Mueller et al. Reference Mueller, Boswell, Mueller and Steckel2014). Pumpkin injury inconsistency across trials was possibly due to its variable persistence depending on soil characteristics and other environmental factors such as microbial degradation (Feng et al. Reference Feng, Li, Zhang, Zhang, Huang, Lu and Li2012; Mielke et al. Reference Mielke, Mendes, de Sousa and Medeiros2022) and rainfall pattern. PPAC-2021 had the most prolonged injury (Table 2), probably because there could have been more available herbicide as a result of less leaching. Herbicide was probably less likely to leach at PPAC-2021 because of a higher OM content (1.7%) than the other two location-years and less rainfall within the first 8 WAT (Table 1).
Weed Control
Weed control data were analyzed separately by location and year. Increasing fomesafen rates did not affect weed control at SWPAC-2021 (F7,21 = 2.01, P = 0.102). Relative to the 0 g ha–1 fomesafen rate treatment that only received S-metolachlor, weed control was above 90% for all fomesafen rates in all three location-years 4 WAT (Table 3). Fomesafen controlled carpetweed (Mollugo verticillata L.), common ragweed (Ambrosia artemisiifolia L.), morningglory spp., pigweeds, and prickly sida (Sida spinosa L.), and grass species at SWPAC-2020; carpetweed, common purslane (Portulaca oleracea L.), and grass species at SWPAC-2021; and carpetweed, common lambsquarters (Chenopodium album L.), giant ragweed, morningglory spp., velvetleaf (Abutilon theophrasti Medik.), volunteer soybean [Glycine max (L.) Merr.], and grass species at PPAC-2021. Weed control during the first 4 wk after emergence is ideal in pumpkin production because of its critical weed-free period of 4 to 6 wk (Dittmar and Boyd Reference Dittmar and Boyd2019; Schonbeck Reference Schonbeck2015). Because plots were maintained weed-free after 4 WAT, we cannot determine from this study how fomesafen rate-related weed control would have affected crop growth and yield.
a Weed control was arcsin transformed for analysis and back-transformed for the table. Scale: 0% = no weed control, 100% = complete weed control.
b Means separation using Tukey’s HSD test P ≤ 0.05. Means followed by the same letter are not significantly different. Lack of letters indicates that the F statistic was not significant at a α = 0.05.
Pumpkin Yield
Because there was a significant fomesafen rate-by-location-year interaction, marketable orange pumpkin yield (F8, 88= 4.78, P = 6.67 × 10–5) and fruit number (F8, 89 = 5.32, P = 1.81 × 10–5) as a percent of the nontreated control data were analyzed separately by location-year. There were no differences in yield, nor fruit number among treatments at SWPAC-2021, where the average marketable orange pumpkin yield was 109 kg 27 m–2 and fruit number was 16 fruits 27 m–2 pooled across all treatments (data not shown). Marketable orange pumpkin yield data as a percent of the nontreated control at SWPAC-2020 and PPAC-2021 fit a three-parameter log-logistic model (Equation 2, Figure 6A). As the fomesafen rate increased from 280 to 1,120 g ha–1, marketable orange pumpkin yield decreased from 95% to 24% of the nontreated control (102 kg 27 m–2) at SWPAC-2020 and 99% to 66% of the nontreated control (119 kg 27 m–2) at PPAC-2021 (Figure 6A). Marketable orange pumpkin fruit number as a percent of the nontreated control fit a three-parameter log-logistic model at SWPAC-2020 (Equation 2) and a three-parameter logistic model at PPAC-2021 (Equation 3) (Figure 6B). As the fomesafen rate increased from 280 to 1,120 g ha–1, the marketable orange pumpkin fruit number decreased from 94% to 21% of the nontreated control (15 fruits 27 m–2) at SWPAC-2020 and 98% to 74% of the nontreated control (17 fruits 27 m–2) at PPAC-2021 (Figure 6B).
Fomesafen rate did not significantly influence the individual marketable orange pumpkin fruit weight nor the marketable green and immature pumpkin yield, fruit number, and individual fruit weight (data not shown).
Although predicted pumpkin marketable orange pumpkin yield and fruit number decreased as the fomesafen rate increased from 280 to 1,120 g ha–1, the values for the lowest fomesafen rate used were not statistically different from the nontreated control. These results confirm the results of Lingenfelter and VanGessel (Reference Lingenfelter and VanGessel2016) and Peachey et al. (Reference Peachey, Doohan and Koch2012), who reported that pumpkin yield was not affected by fomesafen rates of 175 and 350 g ha–1, and 280 and 560 g ha–1, respectively. Lingenfelter and VanGessel (Reference Lingenfelter and VanGessel2016) noted no effect on individual fruit weight as well. Because the individual marketable orange pumpkin fruit weight average was not affected by any fomesafen rate, marketable yield loss at high fomesafen rates was attributed only to the reduced plant stand.
Overall, the recommended, labeled fomesafen rate for use preemergence in other Midwestern states of 280 g ha–1 was safe for use preemergence in Jack O’Lantern pumpkin cultivars ‘Bayhorse Gold’ and ‘Carbonado Gold’ at SWPAC and PPAC. Despite the impact on emergence at SWPAC-2020, the pumpkins recovered, and predicted yield loss was only 5%, and visible injury was less than 7% in all the ratings at all locations. Also, adding this fomesafen rate to a blanket application of S-metolachlor improved weed control (>90% compared to the nontreated control). However, in soils with a low OM content and a high portion of sand, heavy rainfall events shortly after planting are expected to move the herbicide to the crop root zone and affect emergence. Consequently, it is necessary to plan its application carefully. If emergence reduction happens, significant yield loss due to reduced plant stand is expected only at fomesafen rates higher than 280 g ha–1. When considering its use, Indiana pumpkin growers must also consider that fomesafen can only be applied in alternate years and has a maximum total use rate of 343 g ha–1 year–1 in all areas north of Interstate 70 and at 410 g ha–1 year–1 in all areas south of of Interstate 70 in Indiana.
Acknowledgments
We thank the following individuals for assistance with this research: Dennis Nowaskie, Bill Davis, Ange Thompson, Tristand Tucker, Paul Howard, Jacob Imbus, Jessica Conner, Brandi Woolam, Nathaly Vargas, Laura Rodríguez, and Luz Cardona. We thank Syngenta Crop Protection, LLC for supplying herbicides used in this study, and Rupp Seeds, Inc. for supplying pumpkin seeds used in the study. This project was funded by the Indiana State Department of Agriculture through grant A337-19-SCBG-19-003. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the ISDA. This work is supported by the US Department of Agriculture–National Institute of Food and Agriculture Hatch project 7000862. No conflicts of interest have been declared.