Hostname: page-component-586b7cd67f-2plfb Total loading time: 0 Render date: 2024-11-24T00:24:04.991Z Has data issue: false hasContentIssue false

Impacts of repeated glyphosate use on growth of orchard crops

Published online by Cambridge University Press:  18 August 2020

O. Adewale Osipitan
Affiliation:
Postdoctoral Researcher, Department of Plant Sciences, University of California, Davis, CA, USA
Bahar Yildiz-Kutman
Affiliation:
Assistant Professor, Institute of Biotechnology, Gebze Technical University, Kocaeli, Turkey
Seth Watkins
Affiliation:
Staff Research Associate, Department of Plant Sciences, University of California, Davis, CA, USA
Patrick H. Brown
Affiliation:
Professor, Department of Plant Sciences, University of California, Davis, CA, USA
Bradley D. Hanson*
Affiliation:
Cooperative Extension Specialist, Department of Plant Sciences, University of California, Davis, CA, USA
*
Author for correspondence: Bradley D. Hanson, Department of Plant Sciences, One Shields Avenue, University of California, Davis, CA95616 (Email: [email protected])
Rights & Permissions [Opens in a new window]

Abstract

Glyphosate is an important component of herbicide programs in orchard crops in California. It can be applied alone or in tank-mix combinations under the crop rows or to the entire field and often is used multiple times each year. There has been speculation about the potential impacts of repeated use of glyphosate in perennial crop systems, because of uptake from shallow root systems or indirectly because of effects on nutrient availability in soil. To address these concerns, research was conducted from 2013 to 2020 on key orchard crops to evaluate tree response to glyphosate regimens. Almond, cherry, and prune were evaluated in separate experiments. In each crop, the experimental design was a factorial arrangement of two soil types, four glyphosate rates (0, 1.1, 2.2, and 4.4 kg ae ha−1, applied three times annually), and two post-glyphosate application irrigation treatments. In the first 2 yr of the study, there was no clear impact of the glyphosate regimens on shikimate accumulation or leaf chlorophyll content, which suggested no direct effect on the crop. In the seventh year of the study, after six consecutive years of glyphosate application to the orchard floors, there were no negative impacts of glyphosate application on leaf nutrient concentration or on cumulative trunk growth in any of the three orchard crops. Lack of a negative growth impact even at the highest treatment rate, which included 18 applications of glyphosate totaling nearly 80 kg ae ha−1 glyphosate over the course of the experiment suggest there is not likely a significant risk to tree health of judicious use of the herbicide in these production systems. Given the economic importance of orchard crops in California, and grower and industry concerns about pesticides generally and specifically about glyphosate, these findings are timely contributions to weed management concerns in perennial specialty crops.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2020. Published by Cambridge University Press on behalf of Weed Science Society of America

Introduction

Tree fruit, tree nut, and vineyard production systems are economically important in California. In 2017, orchards and vineyards accounted for more than 1.6 million irrigated ha in the state and had an aggregate farm-gate value of more than $21 billion (CDFA 2018). Of the top 10 commodity groupings in California in 2017, grape (Vitis vinifera L.), almond, pistachio (Pistacia vera L.), oranges [Citrus x sinensis (L.) Osbeck] alone had a collective value of $16 billion (CDFA 2018). Although production practices in orchard systems vary among crops and growing regions, weed management is an important component in these intensely farmed high-value crops. As in most crops, weeds can directly compete with trees for limited resources, especially during the establishment period. In addition to direct competition, weeds can interfere with cultural operations such as irrigation, pruning, harvesting, and application of fertilizers and pesticides. In some crops, understory vegetation is managed to reduce the risk of frost during critical periods in the spring. Because of planting arrangements and irrigation infrastructure, weed management in these crops often uses different approaches within the tree row compared with the area between crop rows (Hanson et al. Reference Hanson, Roncoroni, Hembree, Molinar, Elmore, Fennimore and Bell2014). Commonly, California orchards are managed using relatively intense herbicide programs to manage weeds in the “strips” within the tree row, whereas less-intense chemical and/or physical methods are used in the “middles.” The width of the herbicide-treated strip varies by crop and among growers but can range from as narrow as 0.5 m in grapevines and young trees to 5.5 m in large-tree crops such as walnut. Although less common, some growers use full orchard-floor herbicide programs instead of mowing or tillage to reduce management time and dust generated by these mechanical weed-control operations.

In most orchard crops in California, herbicide programs are highly dependent on glyphosate as part of a tank mix applied with PRE herbicides and for POST weed control at multiple times during the year. In all tree crops for which data are available, glyphosate is by far the most widely used herbicide (CDPR 2020). Based on treated acres and gross crop acres, the average orchard is treated two to three times each year with glyphosate and potentially more often if PRE herbicides are not used as part of the year-round weed control program. For several years in the late 2000s, several publications suggested there are nontarget effects of glyphosate, including interactions with crop nutrient status, plant disease interactions, and soil microbial-community effects (reviewed by Duke et al. Reference Duke, Lydon, Koskinen, Moorman, Chaney and Hammerschmidt2012). Although these reports were largely in the context of glyphosate-tolerant soybean systems, California orchard-crop growers and industries expressed concern about micronutrient availability issues or direct glyphosate effects in orchard systems in which glyphosate is regularly used. In particular, one trade publication (Huber Reference Huber2007) was widely distributed among the crop-input supply chain and used to support sales of manganese and other micronutrient products in tree nut crops.

In their review, which focused on annual cropping systems, Duke et al. (Reference Duke, Lydon, Koskinen, Moorman, Chaney and Hammerschmidt2012) suggested that “significant effects of glyphosate on soil mineral content or availability to plants are highly unlikely.” Although relatively less well explored in perennial cropping systems, available crop production statistics in California do not suggest there are negative indirect effects of glyphosate on orchard crops despite decades of use (CDFA 2018). However, given the significant investment involved in establishing and maintaining these long-lived orchard crops, their high value, and the conflicting information available from trade channels, many growers remain concerned about the potential for subtle cumulative direct or indirect effects of glyphosate use in orchard crops. To address these concerns, a 7-yr study was conducted on three orchard crops to evaluate growth and various plant-health metrics related to glyphosate treatments. To create a worst-case scenario, the experimental design included rates up to 4.4 kg ae ha−1, multiple applications per year, a coarse soil in some planting sites, and, in the first 2 yr of the study, herbicide treatments were immediately followed with a simulated flood irrigation event to facilitate downward movement of the herbicide into the root zone.

Materials and Methods

Experimental Description and Design

This research was conducted for the 7 yr during 2013 to 2020. Experimental orchards were planted at the University of California, Davis Plant Sciences Field Facility (38.5382°N, 121.7617°W) for almond, cherry, and prune in spring 2013 (Table 1). Half of the trees in each orchard were planted in the soil native to the field in Yolo County (Rincon silty clay loam), with 1.74% organic matter (OM), pH of 7.7, and 28.1 mEq 100 g−1 cation exchange capacity (CEC) (Andrews Reference Andrews1972); and half were planted in soil imported from Merced County (Delhi sandy loam) with 0.87% organic matter (OM), pH of 6.9, and 3.8 mEq 100 g−1 CEC (Qin et al. Reference Qin, Ga and Ajwa2013). Before planting, each tree site was prepared using a tractor-mounted auger to create a 90-cm diam by 60-cm deep hole; every second hole was refilled with either the native silty clay loam soil or the imported sandy loam soil; these were considered paired plots. The orchard was planted in March 2013 with dormant, bareroot nursery stock in a 3 by 6 m spacing irrigated with microsprinklers at each tree.

Table 1. Planting, glyphosate application, and trunk diameter measurement dates in three orchard experiments conducted to evaluate the cumulative effects of glyphosate over six growing seasons in California.

Each experiment was managed according to local practices (e.g., Micke et al. Reference Micke1996) with microsprinkler irrigation, fertilization, insect and disease management, and pruning practices. Weeds were managed with glufosinate several times each season in addition to the experimental treatment regimens. In 2017, the entire almond experiment was inadvertently treated with one application of glyphosate in addition to the glyphosate and no-glyphosate treatments in the experimental design.

Beginning in spring 2014, glyphosate (Roundup PowerMAX®; Monsanto, St. Louis, MO) was applied at 0, 1.1, 2.2, and 4.4 kg ae ha−1 (equivalent to 0X, 1X, 2X, and 4X, respectively, of a common use rate in orchard crops) three times each growing season (Table 1). Ammonium sulfate (1% vol/vol, Bronc® Max; Wilbur-Ellis Agribusiness, San Francisco, CA) was included with each glyphosate treatment. Glyphosate was applied to a 2 by 2 m area around the base of each tree using two passes (one on either side of the tree to simulate a grower “strip” treatment) with a CO2-pressurized backpack sprayer at 186 kPa, equipped with two flat-fan XR11002 nozzles (TeeJet® Technologies; Spraying Systems Co., Wheaton, IL), calibrated to deliver 187 L ha−1 of spray solution. In the first season, the trunks were protected with wax-paper cartons, but during the remainder of the experiment, the lower 20 cm of trunk was exposed to spray solution when the spray patterns overlapped from the two passes, which is consistent with production practices in the region. Herbicide-application equipment setup, lower-limb pruning, and sucker-removal practices minimized the risk of foliar exposure.

In the first 2 yr of glyphosate treatment, a post-glyphosate irrigation treatment (drench vs. none) was included as a split-plot factor. In the drench treatments, a shallow earthen berm was built around the tree and, immediately after each glyphosate application, 20 L of water was poured into the basin—a volume that approximated a 2.5 ha-cm irrigation within the confined area. The intent of this drench treatment was to increase the potential for leaching of the just-applied herbicide into the root zone of the young tree.

The experimental design was a factorial arrangement of two planting site soils, four glyphosate rates, and two post-treatment irrigation regimens (Table 2). Each treatment was replicated four times in single-tree plots, making a total of 64 experimental units. The almond, cherry, and prune crops were managed separately and considered independent experiments.

Table 2. Treatment structure in three California orchard experiments conducted to evaluate the cumulative effects of three annual glyphosate applications over six consecutive years.

a Experiments were conducted in a site with a silty clay loam soil, but a split-plot factor included half of the trees planted in a 90-cm diam planting site filled with imported sandy loam soil.

b Each treatment was applied three times each year during 2014–2019. Ammonium sulfate was used as adjuvant.

c In the first 2 yr of the experiments, an additional split-plot factor included a 20-L flush of water immediately after each glyphosate application to facilitate movement of the herbicide into the tree root zone.

Data Collection and Analysis

In 2014, the first year of glyphosate application, shikimate assays were conducted 14 d after each glyphosate application to evaluate direct herbicidal effects of glyphosate. Five young leaves were collected from different parts of each tree and this composite leaf sample was prepared for shikimate accumulation assays using the procedure described by Ozturk et al. (Reference Ozturk, Yazici, Eker, Gokmen, Römheld and Cakmak2008) and Hanson et al. (Reference Hanson, Shrestha and Shaner2009). Additionally, in the first 2 yr of glyphosate application (2014 and 2015), the relative chlorophyll content in six randomly selected, youngest, fully expanded leaves from each tree was measured at 30 d after each glyphosate application using a SPAD-502 meter (Spectrum Technologies Inc., Plainfield, IL). These values provided a general approximation of the health of the photosynthetic apparatus and would be affected by either direct herbicidal effect of glyphosate or by indirect effects of micronutrient limitations. On the basis of the initial results, which indicated no treatment-related effects, the shikimate and chlorophyll analyses were discontinued after the first and second growing seasons, respectively. In October 2019, after 6 yr of glyphosate treatments, five young, fully expanded leaves were collected from different parts of each tree to measure the following nutrients: nitrogen (N); phosphorus (P); potassium (K); calcium (Ca); magnesium (Mg); sulfur (S); boron (B); iron (Fe); zinc (Zn); manganese (Mn); and copper (Cu). Leaf samples were collected from all trees but composited over soil type and soil drench split plots to represent only the glyphosate rate main effect in the four experimental replicates in each orchard crop. The nutrients were quantitatively determined by the University of California Davis Analytical Laboratory using combustion (AOAC, 2005) or nitric acid digestion (Sah and Miller, Reference Sah and Miller1992) methods as appropriate for each element.

Trunk diameter was used as a measure of growth of each orchard crop. Trunk diameter 45 cm above the soil surface was measured before the first glyphosate application in 2014 and then in each subsequent year during the dormant season (between January and March of each year), from 2015 to 2020, accounting for 6 yr of post-glyphosate observation. A three-parameter sigmoidal regression model was used to characterize the growth of the orchard crops over 6 yr for the different glyphosate treatments; Equation 1 was used:

([1]) $$T = {\rm{ }}{a \over {\left\{ {1 + {\rm{exp}}\left[ { - \left( {x - h} \right)/b} \right]} \right\}}}$$

where T was the trunk diameter; a was the trunk diameter at the final observation (i.e., sixth year after glyphosate application); x was the year of observation; h was the year at which half of the observed final trunk diameter was achieved, as a measure of growth speed at the log stage; and b was the slope around h. ANOVA was used to determine the effects of soil type, glyphosate rate, post-glyphosate application irrigation, and their interactions on the orchard crop response. A treatment was considered to be significant if P ≤ 0.050, and this was followed by the Tukey honestly significant difference test for mean comparison. All statistical analyses and graphs were performed with SigmaPlot® 14 software (Systat Software Inc., San Jose, CA).

Results and Discussion

Shikimate, Chlorophyll, and Plant Nutrients

There were no differences in shikimate accumulation among the glyphosate-treated and nontreated plants in the almond experiment irrespective of the soil type or postapplication irrigation (Tables 3 and 4). For example, the shikimate accumulation in untreated almond trees (0.20 µmol g−1) 14 d after the third glyphosate application was not significantly different from the 0.28, 0.14, and 0.15 µmol g−1 accumulated shikimate in trees grown in a coarse soil (sandy loam) treated with glyphosate rates of 1.1, 2.2, and 4.4 kg ae ha−1, respectively, with postapplication drench (Tables 3 and 4). Similar shikimate accumulation results were observed in the cherry and prune experiments, except for a significant interaction between soil type (P = 0.035) or postapplication irrigation (P = 0.047) and glyphosate rate, 14 d after the first glyphosate application in cherry (Tables 3 and 4). However, because there was no clear association with the higher-risk treatments (i.e., high application rates, coarse soil, post-treatment drench), these few significant interactions in cherry may be spurious rather than an indication of cherry sensitivity.

Table 3. Mean shikimate concentrations by soil type, postapplication irrigation (drenched vs. none), and glyphosate rate in 2014.

a Experiments were conducted in a site with a silty clay loam soil, but a split-plot factor included half of the trees planted in a 90-cm diam planting site filled with imported sandy loam soil.

b Each treatment was applied three times each year during 2014–2019.

c Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment; FW, fresh weight.

Table 4. P values based on ANOVA of soil type, post-application irrigation (drenched vs. none), glyphosate rate, and their interactions on orchard crops treated with glyphosate three times in 2014.

a Sandy or clay loam soil in original 2013 planting site.

b Post-glyphosate irrigation treatment conducted in the first 2 yr of the experiments.

c Glyphosate applied three times y−1 at 0, 1.1, 2.2, or 4.4 kg ae ha−1.

d Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment.

* P ≤ 0.050.

When absorbed and translocated in a sensitive plant, glyphosate binds with the 5-enolpyruvylshikimate 3-phosphate (EPSP) enzyme in the shikimate pathway, thereby inhibiting the conversion of shikimate-3-phosphate (and phosphoenolpyruvate) to EPSP and inhibiting the biosynthesis of important aromatic acids (Duke and Powles Reference Duke and Powles2008; Wiersma et al. Reference Wiersma, Gaines, Preston, Hamilton, Giacomini, Buell, Leach and Westra2015). A relative increase in shikimate accumulation has been widely used as a measure of the presence and phytotoxicity of glyphosate in a plant (Gaines et al. Reference Gaines, Patterson and Neve2019; Hanson et al. Reference Hanson, Shrestha and Shaner2009; Hernandez et al. Reference Hernandez, Garcia-Plazaola and Becerril1999; Osipitan and Dille Reference Osipitan and Dille2017). At sublethal levels, shikimate accumulation would be expected to increase with increase in glyphosate dose (Shaner et al. Reference Shaner, Nadler-Hassar, Henry and Koger2005; Wilson et al. Reference Wilson, Takano, Van Horn, Yerka, Westra and Stoltenberg2020). In this study, similar shikimate levels in trees grown in glyphosate-treated and untreated plots and no dose-response trend suggest little or no direct herbicidal impact of glyphosate due to root uptake and movement into the aboveground portions of the tree. This is consistent with previous reports of extremely low soil activity of glyphosate due to its strong binding to soil particles (Duke et al. Reference Duke, Lydon, Koskinen, Moorman, Chaney and Hammerschmidt2012). Even in the worst-case scenario in these three orchard experiments, which had young trees growing in sandy soil, treated three times with 4.4 kg ae ha−1 glyphosate, and immediately followed with a flood event (postapplication irrigation), no detectable differences in shikimate accumulation were observed (Tables 3 and 4).

Once taken up by a sensitive plant, glyphosate reduces chlorophyll content, with a consequential effect on photosynthesis and growth in plants (Gomes et al. Reference Gomes, Le Manac’h, Maccario, Labrecque, Lucotte and Juneau2016; Reddy et al. Reference Reddy, Hoagland and Zablotowicz2001; Ye et al. Reference Ye, Huang, Qiu, Wu and Xu2019). Chlorophyll content of the orchard trees was evaluated 30 d after each of the three glyphosate applications in 2014 and 2015, accounting for six observation times for each orchard experiment. Only one of these six observations showed a treatment-related influence on chlorophyll content in each of the three orchard experiments; however, timing and form (increase or decrease) of the treatment effect were inconsistent across experiments (Tables 510). In almond, there was a significant interaction (P = 0.043) between postapplication irrigation and glyphosate rates, whereas in cherry, glyphosate rate was significant as a main effect (P = 0.006) 30 d after third glyphosate application in 2015 (Tables 58). Meanwhile, in prune, the significant influence (P = 0.031) of glyphosate rate was at 30 d after third glyphosate application in 2014 (Tables 9 and 10). If there was a biologically meaningful effect, a greater response would be expected from treatments likely to result in the greatest amount of herbicide in the tree root zone. The inconsistent pattern with regard to the highest glyphosate rate, coarse soil, and postapplication irrigation suggest the relatively few statistically significant results may be due to random variation or experimental artifacts.

Table 5. Mean chlorophyll content (SPAD values) based on soil type, postapplication irrigation, and glyphosate rates in 2014 and 2015.

a Experiments were conducted in a site with a silty clay loam soil, but a split-plot factor included half of the trees planted in a 90-cm diam planting site filled with imported sandy loam soil.

b Each treatment was applied three times each year during 2014–2019.

c Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment; SPAD, Soil Plant Analysis Development.

Table 6. P values based on ANOVA of soil type, postapplication irrigation, glyphosate rate, and their interactions on almond treated with glyphosate three times in 2014 and 2015.

a Sandy or clay loam soil in original 2013 planting site.

b Postglyphosate irrigation treatment conducted in the first 2 yr of the experiments.

c Glyphosate applied three times yr−1 at 0, 1.1, 2.2, or 4.4 kg ae ha−1.

d Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment.

* P ≤ 0.050.

Table 7. Mean chlorophyll content (SPAD values) based on soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on cherry treated with glyphosate three times in 2014 and 2015.

a Experiments were conducted in a site with a silty clay loam soil, but a split-plot factor included half of the trees planted in a 90-cm diam planting site filled with imported sandy loam soil.

b Each treatment was applied three times each year during 2014–2019.

c Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment; SPAD, Soil Plant Analysis Development.

Table 8. P values based on ANOVA of soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on cherry treated with glyphosate three times in 2014 and 2015

a Sandy or clay loam soil in original 2013 planting site.

b Postglyphosate irrigation treatment conducted in the first 2 yr of the experiments.

c Glyphosate applied three times yr−1 at 0, 1.1, 2.2, or 4.4 kg ae ha−1.

d Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment.

* P ≤ 0.050.

Table 9. Mean chlorophyll content (SPAD values) based on soil type, postapplication irrigation (drenched vs. none), glyphosate rate and their interactions on prune treated with glyphosate three times in 2014 and 2015.

a Experiments were conducted in a site with a silty clay loam soil, but a split-plot factor included half of the trees planted in a 90-cm diam planting site filled with imported sandy loam soil.

b Each treatment was applied three times each year during 2014–2019.

c Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment; SPAD, Soil Plant Analysis Development.

Table 10. P values based on ANOVA of soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on prune treated with glyphosate three times in 2014 and 2015.

a Sandy or clay loam soil in original 2013 planting site.

b Postglyphosate irrigation treatment conducted in the first 2 yr of the experiments.

c Glyphosate applied three times yr−1 at 0, 1.1, 2.2, or 4.4 kg ae ha−1.

d Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment.

* P ≤ 0.050.

It has been hypothesized that cumulative effects of glyphosate use in orchard production systems could affect plant growth indirectly by limiting micronutrient availability in the soil and, ultimately, the mineral nutrition of plants (Duke et al. Reference Duke, Lydon, Koskinen, Moorman, Chaney and Hammerschmidt2012). However, evaluation of leaf nutrients (N, P, K, Ca, Mg, S, B, Fe, Zn, Mn, and Cu) of the orchard crops after 6 yr of repeated application of glyphosate at extreme rates did not provide evidence of negative impacts of the treatments on crop nutrient status in three orchard crop experiments (Table 11). This may be explained by the relatively weak chelation of these nutrients by glyphosate (Duke et al. Reference Duke, Lydon, Koskinen, Moorman, Chaney and Hammerschmidt2012; Mertens et al. Reference Mertens, Höss, Neumann, Afzal and Reichenbecher2018) by relatively rapid degradation in the soil environment (Zablotowicz et al. Reference Zablotowicz, Accinelli, Krutz and Reddy2009) and by the relatively small amount of glyphosate applied compared with the concentration of the nutrients in soil (Duke et al. Reference Duke, Lydon, Koskinen, Moorman, Chaney and Hammerschmidt2012).

Table 11. Leaf nutrient analysis in three California orchard crops in October of 2019 after 18 glyphosate applications made during 2014–2019. a

a Data are reported as mean ± SE; n = 4 for each rate and crop. Within-crop and block composite leaf samples included two soil types in the original planting site and two postapplication irrigation treatments for each glyphosate application rate.

b Treatments were applied three times each year during 2014–2019.

c Abbreviations: B, boron; Ca, calcium; Cu, copper; Fe, iron; K, potassium; Mg, magnesium; Mn, manganese; N, nitrogen; P, phosphorus; S, sulfur; Zn, zinc.

* P ≤ 0.050.

Trunk Diameter Growth

The use of trunk diameter has been widely used as a robust measure of orchard crop growth (Hernandez-Santana et al. Reference Hernandez-Santana, Fernández, Cuevas, Perez-Martin and Diaz-Espejo2017; Martín-Palomo et al. Reference Martín-Palomo, Corell, Girón, Andreu, Trigo, López-Moreno, Torrecillas, Centeno, Pérez-López and Moriana2019; Moriana et al. Reference Moriana, Orgaz, Pastor and Fereres2003). From 2014 to 2020, average trunk diameter (across treatments) increased from 28 to 169 mm in almond, 28 to 217 mm in cherry, and 32 to 127 mm in prune (Figure 1). There were no differences in final trunk diameter among treatments in the almond or prune experiments; in cherry, there was a slight but significant increase in trunk diameter with the highest glyphosate rate. ANOVA revealed no significant impact of soil type in the planting hole, post-glyphosate application irrigation, or their interactions on cumulative trunk-diameter change in three orchard crops over 7 yr of observation that included 6 yr of glyphosate treatment (Table 12). Similarly, the ANOVA results suggested that glyphosate rates applied on the soil had no negative impact on growth of the orchard crops and, in the case of cherry, a slight increase in trunk diameter in the highest glyphosate rate (Table 13). The regression model adequately fitted the cumulative trunk growth data over the 6-yr period of observation (Figure 1) with small root mean square error (≤6%) in all three experiments. In terms of growth rate, it took 2 to 3 yr for the orchard crops to attain half of their final trunk diameter, but these growth-rate data were not influenced by glyphosate regimens (Table 8).

Figure 1. Increase in almond, cherry, and prune trunk diameter over time with or without glyphosate applied around the base of the tree three times per year at up to 4.4 kg ha−1 over 6 yr. Data were averaged over the soil type in the original planting hole and a postapplication irrigation conducted after each glyphosate treatment in the first two growing seasons.

Table 12. P values based on ANOVA of soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on trunk diameter of orchard crops in 2014 to 2020.

a Sandy or clay loam soil in original 2013 planting site.

b Postglyphosate irrigation treatment conducted in the first 2 yr of the experiments.

c Glyphosate applied three times yr−1 at 0, 1.1, 2.2, or 4.4 kg ae ha−1.

d Abbreviations: DAT1, days after first glyphosate treatment; DAT2, days after second glyphosate treatment; DAT3, days after third glyphosate treatment.

* P ≤ 0.050.

Table 13. Trunk diameter increase in California orchard crops after 18 glyphosate applications made during 2014–2019.

a Each treatment was applied three times each year during 2014–2019.

b Regression parameters are final trunk diameter (a), slope (b), and number of years (h) to achieve half of the final trunk diameter for different rates of soil-applied glyphosate in three orchard crops.

c Data are reported as mean ± SE.

In apple (Malus spp.), glyphosate has been linked to trunk injury and scaffold death (e.g., Rosenberger et al. Reference Rosenberger, Watkins, Miranda-Sazo, Kahlke, Fargione, Nock and Rugh2013). However, in these California experiments, no trunk cankers or other trunk and limb malformations were observed in almond, cherry, or prune during the 7-yr evaluation period (data not shown). The lack of negative impact of glyphosate on the growth rate or total growth of the orchard trees over time supports the finding of no measurable effects on shikimate accumulation in the first year, leaf chlorophyll content in the first 2 yr, or leaf nutrient levels after 6 yr of treatment. The highest glyphosate rate (4.4 kg ae ha−1) evaluated in this study resulted in a total annual soil glyphosate load of 13.3 kg ae ha−1, or nearly 80 kg ae ha−1 over the life of the experiment, which is well beyond what would typically occur in a commercial orchard in California.

Several studies that simulated drift cases in sensitive plants or at above-label doses in glyphosate-tolerant plants have shown that foliar-applied glyphosate can have negative impacts on plant nutrient uptake, photosynthetic apparatus, and plant productivity (Al-Khatib et al. Reference Al-Khatib, Parker and Fuerst1992; Cakmak et al. Reference Cakmak, Yazici, Tutus and Ozturk2009; Foshee et al. Reference Foshee, Blythe, Goff, Faircloth and Petterson2008; Gomes et al. Reference Gomes, Le Manac’h, Maccario, Labrecque, Lucotte and Juneau2016; Huang et al. Reference Huang, Silva, Shen, Jiang and Lu2012; Su et al. Reference Su, Ozturk, Cakmak and Budak2009). In a glyphosate-resistant soybean, foliar application of glyphosate at 2.4 kg ae ha−1 substantially reduced chlorophyll content (Soil Plant Analysis Development (value), photosynthetic rate, plant nutrients and growth of the plant (Zobiole et al. Reference Zobiole, Kremer, de Oliveira and Constantin2012). Root uptake of glyphosate, with resulting impacts on plant nutrient status, is possible, such as reported by Ozturk et al. (Reference Ozturk, Yazici, Eker, Gokmen, Römheld and Cakmak2008), who reported reduced ferric reductase activity in sunflower exposed to glyphosate in soil-free hydroponic conditions. However, this has not been observed from glyphosate residues in soil under field conditions (Bromilow et al. Reference Bromilow, Evans, Nicholls, Todd and Briggs1996; Duke et al. Reference Duke, Lydon, Koskinen, Moorman, Chaney and Hammerschmidt2012; Liphadzi et al. Reference Liphadzi, Al-Khatib, Bensch, Stahlman, Dille, Todd, Rice, Horak and Head2005) and does not appear to be the case in California orchard crops in the current study.

California orchard-production systems are unique compared with many other U.S. crops in terms of potential for direct or indirect impacts of glyphosate use on crop safety, due to the long lifespan of the crop, potential for repeated and high-rate applications of the herbicide, minimal soil disturbance, and the Mediterranean climate conditions. Three experiments with treatments intended to create a worst-case scenario for cumulative glyphosate soil residues to directly or indirectly affect the growth of California orchard crops over 7 yr did not suggest any negative impact. These findings are timely and add value to current conversations about glyphosate use in perennial specialty crops.

Acknowledgements

This research was partially supported with funding from the Almond Board of California and the California Dried Plum Board, and a donation of orchard nursery stock from Sierra Gold Nurseries. No conflicts of interest have been declared

Footnotes

Associate Editor: Darren Robinson, University of Guelph

References

Al-Khatib, K, Parker, R, Fuerst, EP (1992) Sweet cherry (Prunus avium) response to simulated drift from selected herbicides. Weed Technol 6:975979 CrossRefGoogle Scholar
Andrews, WF (1972) Soil Survey of Yolo County, California. U.S. Department of Agriculture, Soil Conservation Service. 102 pGoogle Scholar
[AOAC] AOAC International (2005) Protein (crude) in animal feed, combustion method. Chapter 4, pages 30-31, in Official Methods of Analysis of AOAC International . 18th ed. Gaithersburg, MD: AOAC International Google Scholar
Bromilow, RH, Evans, AA, Nicholls, PH, Todd, AD, Briggs, GG (1996) The effect on soil fertility of repeated applications of pesticides over 20 years. Pestic Sci 48:6372 3.0.CO;2-I>CrossRefGoogle Scholar
Cakmak, I, Yazici, A, Tutus, Y, Ozturk, L (2009) Glyphosate reduced seed and leaf concentrations of calcium, manganese, magnesium, and iron in non-glyphosate resistant soybean. Eur J Agron 31:114119 CrossRefGoogle Scholar
[CDFA] California Department of Agriculture (2018) California agricultural production statistics. https://www.cdfa.ca.gov/statistics/. Accessed: April 2, 2020.Google Scholar
[CDPR] California Department of Pesticide Regulation (2020) California Pesticide Information Portal application. https://calpip.cdpr.ca.gov/main.cfm. Accessed: April 2, 2020.Google Scholar
Duke, SO, Powles, SB (2008) Mini-review. Glyphosate: a once-in-a-century herbicide. Pest Manag Sci 64:319325 CrossRefGoogle Scholar
Duke, SO, Lydon, J, Koskinen, WC, Moorman, TB, Chaney, RL, Hammerschmidt, R (2012) Glyphosate effects on plant mineral nutrition, crop rhizosphere microbiota, and plant disease in glyphosate-resistant crops. J Agric Food Chem 60:1037510397 CrossRefGoogle ScholarPubMed
Foshee, WG, Blythe, EK, Goff, WD, Faircloth, WH, Petterson, MG (2008) Response of young pecan trees to trunk and foliar applications of glyphosate. HortScience 43:399402 CrossRefGoogle Scholar
Gaines, TA, Patterson, EL, Neve, P (2019) Molecular mechanisms of adaptive evolution revealed by global selection for glyphosate resistance. New Phytol 223:17701775 CrossRefGoogle ScholarPubMed
Gomes, MP, Le Manac’h, SG, Maccario, S, Labrecque, M, Lucotte, M, Juneau, P (2016) Differential effects of glyphosate and aminomethylphosphonic acid (AMPA) on photosynthesis and chlorophyll metabolism in willow plants. Pestic Biochem Phys 130:6570.CrossRefGoogle ScholarPubMed
Hanson, BD, Shrestha, A, Shaner, DL (2009) Distribution of glyphosate-resistant horseweed (Conyza canadensis) and relationship to cropping systems in the Central Valley of California. Weed Sci 57:4853 CrossRefGoogle Scholar
Hanson, BD, Roncoroni, J, Hembree, KJ, Molinar, R, Elmore, CL (2014) Tree, vine, and soft-fruit crops. In Fennimore, SA and Bell, C (eds). Principles of Weed Control . Salinas CA: California Weed Science Society Google Scholar
Hernandez, A, Garcia-Plazaola, JI, Becerril, JM (1999) Glyphosate effects on phenolic metabolism of nodulated soybean (Glycine max L. Merr.) J Agric Food Chem 47:29202925 CrossRefGoogle ScholarPubMed
Hernandez-Santana, V, Fernández, JE, Cuevas, MV, Perez-Martin, A, Diaz-Espejo, A (2017) Photosynthetic limitations by water deficit: effect on fruit and olive oil yield, leaf area and trunk diameter and its potential use to control vegetative growth of super-high density olive orchards. Agric Water Manag 184:918 CrossRefGoogle Scholar
Huang, J, Silva, EN, Shen, Z, Jiang, B, Lu, H (2012) Effects of glyphosate on photosynthesis, chlorophyll fluorescence and physicochemical properties of cogongrass (Imperata cylindrical L.). Plant Omics 5:177 Google Scholar
Huber, DM (2007) What about glyphosate-induced manganese deficiency? Fluid J 15:2022 Google Scholar
Liphadzi, KB, Al-Khatib, K, Bensch, CN, Stahlman, PW, Dille, JA, Todd, T, Rice, CW, Horak, MJ, Head, G (2005) Soil microbial and nematode communities as affected by glyphosate and tillage practices in a glyphosate-resistant cropping system. Weed Sci 53:536545 CrossRefGoogle Scholar
Martín-Palomo, MJ, Corell, M, Girón, I, Andreu, L, Trigo, E, López-Moreno, YE, Torrecillas, A, Centeno, A, Pérez-López, D, Moriana, A (2019) Pattern of trunk diameter fluctuations of almond trees in deficit irrigation scheduling during the first seasons. Agric Water Manag 218:115123 CrossRefGoogle Scholar
Mertens, M, Höss, S, Neumann, G, Afzal, J, Reichenbecher, W (2018). Glyphosate, a chelating agent—relevant for ecological risk assessment? Environ Sci Pollut Res Int 25:52985317 CrossRefGoogle ScholarPubMed
Micke, WC (1996) Almond Production Manual . Davis, CA: University of California Division of Agriculture and Natural Resources. Publication 3364. 289 p Google Scholar
Moriana, A, Orgaz, F, Pastor, M, Fereres, E (2003) Yield responses of a mature olive orchard to water deficits. J Am Soc Hortic Sci 128:425431 CrossRefGoogle Scholar
Osipitan, OA, Dille, JA (2017) Fitness outcomes related to glyphosate resistance in kochia (Kochia scoparia): what life history stage to examine? Front Plant Sci 8:1090 CrossRefGoogle ScholarPubMed
Ozturk, L, Yazici, A, Eker, S, Gokmen, O, Römheld, V, Cakmak, I (2008) Glyphosate inhibition of ferric reductase activity in iron deficient sunflower roots. New Phytol 177:899906 CrossRefGoogle ScholarPubMed
Qin, R, Ga, S, Ajwa, H (2013) Emission and distribution of fumigants as affected by soil moistures in three different textured soils. Chemosphere 90:866872 CrossRefGoogle ScholarPubMed
Reddy, KN, Hoagland, RE, Zablotowicz, RM (2001) Effect of glyphosate on growth, chlorophyll, and nodulation in glyphosate-resistant and susceptible soybean (Glycine max) varieties. J Seed Sci 2(3):3752 Google Scholar
Rosenberger, D., Watkins, C, Miranda-Sazo, M, Kahlke, C, Fargione, M, Nock, J, Rugh, A (2013) Effects of glyphosate on apple tree health. New York Fruit Quarterly 21(4):2327.Google Scholar
Sah, RN, Miller, RO (1992) Spontaneous reaction for acid dissolution of biological tissues in closed vessels. Anal Chem 64:230233 CrossRefGoogle ScholarPubMed
Shaner, DL, Nadler-Hassar, T, Henry, WB, Koger, CH (2005) A rapid in vivo shikimate accumulation assay with excised leaf discs. Weed Sci 53:769774 CrossRefGoogle Scholar
Su, YS, Ozturk, L, Cakmak, I, Budak, H (2009) Turfgrass species response exposed to increasing rates of glyphosate application. Eur J Agron 31:120125 Google Scholar
Wiersma, AT, Gaines, TA, Preston, C, Hamilton, JP, Giacomini, D, Buell, CR, Leach, JE, Westra, P (2015) Gene amplification of 5-enol-pyruvylshikimate-3-phosphate synthase in glyphosate-resistant Kochia scoparia . Planta 241:463474 CrossRefGoogle ScholarPubMed
Wilson, CE, Takano, HK, Van Horn, CR, Yerka, MK, Westra, P, Stoltenberg, DE (2020) Physiological and molecular analysis of glyphosate resistance in non-rapid response Ambrosia trifida from Wisconsin. Pest Manag Sci 76:150160 CrossRefGoogle ScholarPubMed
Ye, J, Huang, C, Qiu, Z, Wu, L, Xu, C (2019) The growth, apoptosis and oxidative stress in Microcystis viridis exposed to glyphosate. Bull Environ Contam Toxicol 103:585589 CrossRefGoogle ScholarPubMed
Zablotowicz, RM, Accinelli, C, Krutz, LJ, Reddy, KN (2009) Soil depth and tillage effects on glyphosate degradation. J Agric Food Chem 57:48674871 CrossRefGoogle ScholarPubMed
Zobiole, LHS, Kremer, RJ, de Oliveira, RS Jr, Constantin, J (2012) Glyphosate effects on photosynthesis, nutrient accumulation, and nodulation in glyphosate-resistant soybean. J Plant Nutr Soil Sci 175:319330 CrossRefGoogle Scholar
Figure 0

Table 1. Planting, glyphosate application, and trunk diameter measurement dates in three orchard experiments conducted to evaluate the cumulative effects of glyphosate over six growing seasons in California.

Figure 1

Table 2. Treatment structure in three California orchard experiments conducted to evaluate the cumulative effects of three annual glyphosate applications over six consecutive years.

Figure 2

Table 3. Mean shikimate concentrations by soil type, postapplication irrigation (drenched vs. none), and glyphosate rate in 2014.

Figure 3

Table 4. P values based on ANOVA of soil type, post-application irrigation (drenched vs. none), glyphosate rate, and their interactions on orchard crops treated with glyphosate three times in 2014.

Figure 4

Table 5. Mean chlorophyll content (SPAD values) based on soil type, postapplication irrigation, and glyphosate rates in 2014 and 2015.

Figure 5

Table 6. P values based on ANOVA of soil type, postapplication irrigation, glyphosate rate, and their interactions on almond treated with glyphosate three times in 2014 and 2015.

Figure 6

Table 7. Mean chlorophyll content (SPAD values) based on soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on cherry treated with glyphosate three times in 2014 and 2015.

Figure 7

Table 8. P values based on ANOVA of soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on cherry treated with glyphosate three times in 2014 and 2015

Figure 8

Table 9. Mean chlorophyll content (SPAD values) based on soil type, postapplication irrigation (drenched vs. none), glyphosate rate and their interactions on prune treated with glyphosate three times in 2014 and 2015.

Figure 9

Table 10. P values based on ANOVA of soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on prune treated with glyphosate three times in 2014 and 2015.

Figure 10

Table 11. Leaf nutrient analysis in three California orchard crops in October of 2019 after 18 glyphosate applications made during 2014–2019.a

Figure 11

Figure 1. Increase in almond, cherry, and prune trunk diameter over time with or without glyphosate applied around the base of the tree three times per year at up to 4.4 kg ha−1 over 6 yr. Data were averaged over the soil type in the original planting hole and a postapplication irrigation conducted after each glyphosate treatment in the first two growing seasons.

Figure 12

Table 12. P values based on ANOVA of soil type, postapplication irrigation (drenched vs. none), glyphosate rate, and their interactions on trunk diameter of orchard crops in 2014 to 2020.

Figure 13

Table 13. Trunk diameter increase in California orchard crops after 18 glyphosate applications made during 2014–2019.