1. Introduction
Conventional farming and integrated pest management (IPM) regimes relate to the utilization of pesticides for the protection of agricultural crops against diseases and pests. However, pesticides are potentially toxic to humans and can have both acute and chronic health effects, depending on the quantity and ways in which a person is exposed [
1]. Hence, controlled use of agricultural pesticides is important, modelling pesticide distribution in crops as a tool in limiting the excessive use of pesticides [
2]. Maximum residue limits (MRLs) for pesticide/commodity combinations have been established by the European Union in the Regulation of European Commission [
3] as the highest level of pesticide residues that are legally tolerated in food or feed. MRLs are based on Good Agricultural Practice (GAP) data and must meet requirements for pesticide registration. A uniform MRL as low as 0.01 mg·kg
−1 has been established by Commission Directives [
4,
5] for any pesticide residue in baby foods and processed cereal-based foods [
6]. To meet the legal limits for pesticide residues in various food crops at a given harvest time, a pre-harvest interval (PHI) is officially established for particular pesticides defined as the time between the last pesticide application and the harvest of the treated crop. This value is based on FAO (Food and Agriculture Organization of the United Nations) recommendation [
7] and is usually fixed at a country level [
8]. The aim of “zero-residue” vegetable production (“pesticide residue free production”) is achieving pesticide residues in respective crops as low as 0.01 mg·kg
−1 (or even lower). This limit corresponds to the limit used for baby food production [
4]. Adopting MRLs and PHIs helps to keep residue-related health risks below an “acceptable risk level”.
The most recent EU report summarizes the study of pesticide residue contamination of food carried out in the European Union (EU) Member States, Iceland, and Norway [
9]. It summarizes the results of both the 2017 EU-coordinated control program (EUCP) and the national control programs (NP). To allow the assessment of representative consumer exposure to pesticide residues by food commodity, the same pattern of commodities is monitored for the presence of pesticide residues in 3-year cycles in the same countries. Regarding the 2017 EUCP, results were compared with the ones of 2014 for those commodities sampled in both years. The results of this showed some improvement in the overall situation in comparison with the results of 2014. No quantifiable residues were reported in 54.1% and 53.6% of the samples in 2017 and 2014, respectively. Quantified residues at or below the MRLs were reported in 41.8% and 43.4% of the samples in 2017 and 2014, respectively [
9]. Overall, 95.9% of the samples analyzed in 2017 fell within the legal limits; nevertheless, 41.8% of the samples contained quantified residues at or below the MRLs. Nevertheless, the dietary risk assessment indicated that, for the samples analyzed, the probability of European citizens being exposed to pesticide residue levels that could lead to negative health outcomes is relatively low [
9].
However, supermarkets, and thus also food crop producers, are continuously under pressure due to consumer concern about the health risks associated with the potential occurrence of pesticide residues in their diet. When seeking to rationalize pesticide use, both government [
10] and supermarkets have tended to make the implicit assumption that any rationalization in this field is primarily an issue of decreasing the quantity of pesticides used, coupled with the banning of certain persistent or highly toxic substances [
11,
12]. Some retail chains apply for the assessment of contamination of products, mainly fruits and vegetables, originating from low-residue production. So called “action thresholds” declare the highest acceptable concentration of pesticide residues in the product corresponding to a percentage of the respective MRL, e.g., 25%. The required action thresholds commonly range from 25% to 75% MRL depending on the retail chain. In all cases, the action thresholds are fairly below the MRL. Achieving action thresholds for low-residue production or zero-residue production is based on a setting of an action pre-harvest interval (APHI), which is the minimal time between the last pesticide application and crop harvest within which the required drop in pesticide residue levels occurs. The knowledge of APHIs, which is necessary for setting relevant pesticide treatment regimes in low-residue and/or zero-residue production, needs the development of dissipation models of the active substances in pesticide preparations, based on which calculations are performed. In this context, information on the dissipation kinetics of pesticide residues in food crops and other plants is an additional key aspect of current risk and impact assessment practices. Considering more than 740 published studies, Fantke and Juraske [
13] identified nine different calculation models, corresponding residual pesticide concentration curves
Ct and related half-lives
t1/2. The authors concluded that pesticide dissipation in plants generally follows first-order kinetics, although some exceptions exist. Currently, experimental data are not available for all relevant pesticide−plant combinations; nevertheless, Fantke et al. [
14] created a list of corrected geometric means of dissipation half-lives for 333 pesticides applied to an average plant under reference air temperature (20 °C). Half-lives
t1/2 ref,i were calculated using model II based on dissipation data obtained for 346 pesticide residues in 183 treated plant species [
13,
14]. Model II also provides estimates for the influence of specific plants, which can be used to correct
t1/2 ref,i for pesticides applied to a particular crop under given study conditions and temperature. Model III was designed to estimate dissipation half-lives from physico-chemical properties of individual pesticides of 14 substance classes for which no reference half-lives are available, i.e., they are not listed in
t1/2 ref,i [
14]. The model III is more uncertain as it tries to predict half-lives for chemicals not included in the fitting dataset.
Pesticide residues were analyzed in Brassica vegetables with application until harvest in order to establish mathematical models of residue dissipation and forecast the action pre-harvest interval for a given action threshold as a percentage of MRL for low-residue production in Brassica vegetables [
15]. This study describes the dissipation of 22 active substances of pesticide formulations applied on Chinese cabbage, head cabbage, and cauliflower.
Experimental studies usually present aggregated dissipation data estimated from measuring the changes of the overall residual pesticide concentration in the respective plant over time. However, the overall dissipation of residues involves several processes. Jacobsen et al. [
2] showed that degradation of a pesticide molecule is, in many cases, the most relevant dissipation component, followed by growth dilution and volatilization. Dissipation rate and half-lives of pesticides vary according to vegetable species and depend on plant characteristics, environmental conditions, and substance physico-chemical properties [
7]. Hence, the models of pesticide dissipation are applicable in regions with similar climatic conditions and under comparable conditions of vegetable growing [
2].
The main objective of the study was to evaluate dynamics of pesticide residues in iceberg lettuce, onion, leek, carrot, and parsley, and, based on generated data, to suggest a procedure enabling targeted (low) levels of pesticide residues in respective vegetable at harvest. The specific objectives of this work were (1) to quantify the dissipation rates of pesticide residues in tested vegetables, (2) to determine the APHI25 and APHI0.01 that would guarantee achieving “action thresholds” of 25% MRL and 0.01 mg residue kg−1, respectively, for a low-residue and zero-residue production of tested vegetables, and (3) to critically assess differences in the dissipation rates and dissipation half-lives of particular pesticide residues in tested vegetable species and compare the results with available literature data.
2. Materials and Methods
Pesticide residues were analyzed in samples collected from semi-field experiments on iceberg lettuce (
Lactuca sativa L.), leek (
Allium ampeloprasum L.), carrot (
Daucus carota L. subsp. sativus (Hoffm.) Arcang), and parsley (
Petroselinum crispum (Mill.) Fuss) at the Crop Research Institute in Prague (GPS: 50.0864428N, 14.2985553E, soil type: illimerized luvisol, soil pH: 6.8) and on onion (
Allium cepa L.) at the Czech University of Life Sciences Prague (GPS: 50.1267258N, 14.3770625E, soil type: black earth, pH: 7.3). The specific information of the geographic coordinates, soil type and soil pH are presented according to current reporting recommendations [
16]. The studied vegetables were usually treated with pesticides registered for the control of vegetable diseases and pests in the Czech Republic. However, some additional pesticides permitted in neighboring countries were also selected for the case of their future use in the Czech Republic. Iceberg lettuce and both root vegetables were planted in a single season, while onion and leek were cultivated in two consecutive years (
Table 1). Cultivation technologies (row and in-row plant spacing, fertilization) respected general field practices. The same pesticide application and sampling methodology were used for all tested vegetable crops. The average air temperature (min/max range) from the first chemical treatment until the last harvest was 20 °C (13.8/26.8 °C) in lettuce 2011; 18.3 °C (13.7/29.7 °C) and 17.3 °C (12.8/24.7 °C) in onion in 2012 and 2013, respectively; 9 °C (1.7/13.9 °C) in leek in 2008; 10.5 °C (1.7/19.3 °C) in leek as well as in carrot in 2009; and 12 °C (6.4/15.2 °C) in parsley in 2010. The crops were cultivated in three plots (A–C); three replications per plot were set up in a randomized complete block design. Each plot was treated with three different pesticide mixtures (
Supplementary Table S1 online) on three different dates (
Table 1). Thus, pesticide mixture 1 was sprayed in plot A on the first date, in plot B on the second date and in plot C on the third date of experiment. Pesticides in each mixture are marked with a superscript in
Supplementary Table S1 online. The pesticide mixtures were different for each crop. The mixture composition depended on the efficacy of the pesticides in each crop. Pesticides were applied from the dates selected based on plant phenophase and the estimated harvest date. Plant stages at the time of application of pesticides were BBCH 42, 45, 45, 45, and 45 for iceberg lettuce, leek, carrot, parsley, and onion, respectively [
17]. All the pesticides were administered at maximum label application rates. In all variants, 500 L·ha
−1 spraying with 0.1 L·ha
−1 Silwet was used. Sample collection started three days after the last spraying. In each plot, the samples were collected on four different dates; thus, twelve samples of each vegetable were analyzed per year except for onion, where only three sample collections were performed in 2012. Crops, planting dates, varieties, pesticide application dates, and harvest dates are listed in
Table 1.
Analyses of pesticide residues were performed by the testing laboratory of the University of Chemistry and Technology accredited according to the EN ISO/IEC 17025 standard. The standardized analytical method used in this study is based on the EN 12393-2 standard (Foods of plant origin-Multiresidue methods for the determination of pesticide residues by GC or LC-MS/MS-Part 2: Methods for extraction and clean-up) and EN 12393-2 (Foods of plant origin-Multiresidue methods for the determination of pesticide residues by GC or LC-MS/MS-Part 3: Determination and confirmatory tests).
Approximately 1 kg of representative vegetable samples from field experiments were delivered to the laboratory at the University of Chemistry and Technology in Prague. The edible part of the crop and the top leaves were separated and homogenized using a Retsch GM200 blender (Retsch GmbH, Haan, Germany) and stored in a freezer (−18 °C) until analysis. An ISO 17025 [
18] accredited multiresidue method based on QuEChERS extraction followed by high-performance liquid chromatography coupled with tandem mass spectrometry (HPLC-MS/MS) was employed, encompassing most of the examined pesticides (abamectin, acetamiprid, azoxystrobin, boscalid, chlorantraniliprole, chlorpyrifos, cymoxanil, cyprodinil, cypermethrin, deltamethrin, difenoconazole, dimethoate, dimethomorph, fludioxonil, fluoxastrobin, indoxacarb, mandipropamid, metalaxyl-M, methoxyfenozide, pirimicarb, propamocarb-hydrochloride, prothioconazole, pyraclostrobin, pyridaben, spinosad, tebuconazole, thiacloprid, thiamethoxam); for pymetrozine, a single residue method with pH adjustment of the sample was performed. Sample preparation and HPLC-MS/MS analysis were described in detail in our previous study [
19], and the multiple reaction monitoring (MRM) conditions optimized for the pesticides evaluated in this study are summarized in
Supplementary Table S2 online. These analytical procedures were introduced earlier by Ticha et al. [
20] and subsequent validation protocol as well as the Internal Quality Control (IQC) measures followed document SANTE/11945/2015 (latest consolidated version SANTE/12682/2019) “Method Validation & Quality Control Procedures for Pesticide Residues Analysis in Food & Feed” [
21], obligatorily used by the EU official control laboratories. MS measurement parameters of the optimized method are given in
Table S2 (LC-based) and Table S3 (GC-based, confirmatory for GC amenable pesticides) in Supplementary Materials. The performance characteristics for all pesticides involved in this study are summarized in
Supplementary Table S4 online. The mean recovery, repeatability standard deviation, and limit of quantitation (LOQ) were evaluated for each pesticide. This means that the entire analytical procedure (started by the handling of the analytical portion) is represented by the metrological characteristics given in
Table S4. The performance of the accredited laboratory was verified (external quality control, QC) through participation of the laboratory in regular proficiency tests: (i) Food Analysis Performance Assessment Scheme (FAPAS
®) and (ii) European Commission’s Proficiency Testing Program (EU–PT).
The following first-order kinetic equation was used to characterize the dissipation rate of active substances of the pesticides in crop products:
where
Ct is the residual pesticide concentration (mg·kg
−1) at time
t (days) after application;
C0 is the initial pesticide concentration; and
kdiss is the pesticide dissipation rate constant (day
−1). ANCOVA was used to analyze the effect of the year and time after application on the residual pesticide concentration. If there was no significant difference between the years, the data from two years were combined. The corresponding dissipation half-life equation is:
where
t1/2 is the pesticide dissipation half-life (days), which was used for all pesticides with significant models (
R2 > 0.5). The model parameters of
C0 and
kdiss were calculated from all experimentally determined residue values of each active pesticide substance in the products (
Ct) and the corresponding number of days after pesticide application (
t). Calculations were performed by Kruskal–Wallis test in the XLSTAT 2009 program (Addinsoft, New York, NY, USA). Our
t1/2 values were compared with corrected geometric means of estimated half-lives applied to specific plants in each set of production conditions and temperature using model II [
14]. Two active substances (fluoxastrobin [
22] and prothioconazole [
23]) were not involved in the pesticide list needed for the
t1/2 calculation according to model II, so we used model III, which Fantke et al. [
14] designed to predict the half-lives of individual pesticides for which no reference half-lives are available.
Maximum residue limits were obtained in EU Pesticides database [
24]. Because the MRL of active substances used for carrot and parsley roots was the same, the missing MRL for carrot leaves was replaced by the MRL of parsley leaves. MRLs for carrot leaves are not established because carrot is grown almost exclusively for the root. However, raw young carrot leaves have recently been used not only as food for animals but also to enrich dishes such as salads or soups. “Action thresholds” for the 25% MRL (APHI
25) and 0.01 mg kg
−1 limit (APHI
0.01) were established for 31 active substances in the pesticide formulations. To calculate the APHI, the following equation can be used:
where
Ct is the residual pesticide concentration in mg·kg
−1 at time
t (days) after application;
C0 is the initial pesticide concentration; and
kdiss is the pesticide dissipation rate constant (day
−1). For zero-residue production,
Ct = 0.01 mg·kg
−1; for low-residue production,
Ct corresponds to a given percentage of the MRL (e.g., 25% MRL). The APHIs calculated from the pesticide dissipation models were extended by one-third based on a confidence interval of the model for active substances with the highest variability of dissipation rate to increase the reliability of APHIs, i.e.,
The parameters
C0 and
kdiss differ according to the active pesticide substance and the agricultural commodity. When the calculated APHI prolonged by one-third was shorter than the obligatory PHI according to a list of registered products [
25], the PHI was used.
4. Discussion
In the tested vegetables, the pesticide dissipation rate
kdiss and the pesticide dissipation half-life
t1/2 varied depending on the active substances and vegetable species. The fastest pesticide dissipation occurred overall in lettuce (
Table S5A;
Figure 1a). Our results corresponded with the findings of Song et al. [
27], who compared dissipation curves and half-lives of six pesticides on six leafy vegetables planted in field trails. The dissipation of dimethoate, chlorpyrifos, beta-cypermethrin and deltamethrin in leaf lettuce expressed in
kdiss and
t1/2 was in the range of −0.793 to −0.176 and 0.87 to 3.94 days, respectively.
The dissipation half-lives of the active substances were compared with the predicted geometric mean of the dissipation half-lives at 20 °C (
t1/2 ref,i) and the corrected
t1/2 ref,i of pesticides applied to a specific vegetable under a given average temperature (
t1/2 plant, active subst.) calculated according to model II [
14]. The results are given in
Table S7. In many cases, the corrected
t1/2 plant, active subst of the active substances correspond more to our results; e.g., the corrected
t1/2 in lettuce and onion were shorter than
t1/2 ref,i, whereas in leek, the corrected half-lives were prolonged for several days (
Table S7). Two active substances, fluoxastrobin and prothioconazole, sprayed on onion did not appear in the pesticide list needed for
t1/2 calculation according to model II. We compared our results with
t1/2 ref,i and
t1/2 plant, active subst. computed according to model III [
14]. Although the calculation according to model III has its limitations, such as higher uncertainty compared to half-lives based on model II [
14], the results correlated with our values (
Table S7). Computed predicted half-lives (
t1/2 ref,i;
t1/2 plant, active subst.) of fluoxastrobin 2012–13 (6.14; 6.73 days) and prothioconazole 2013 (7.19; 8.14 days) corresponded to our results of 7.20 and 6.78 days, respectively.
Eight of the thirty-two active substances were sprayed on all vegetable species. The results showed differences in concentrations of active substances in underground and aboveground parts of the plants after foliar application of the pesticides. The initial deposition of the active substance in aboveground part could be influenced by the foliar surface area, which can be defined as vegetation cover (VC), e.g., the proportion of soil area covered by leaves [
8]. Leaves of carrot and parsley are overlapping, while the compact shape of iceberg lettuce head and leek could allow less deposition of the active substance. The role of foliar area in the initial deposition of active substances has been mentioned, e.g., by Lu et al. [
28] and Song et al. [
27]. Two fungicides (azoxystrobin, tebuconazole) and three insecticides (cypermethrin, spinosad, thiacloprid) were selected to demonstrate differences in pesticide dissipation.
The systemic fungicide azoxystrobin had robust models for all vegetables and its parts. The dissipation rate of azoxystrobin was as follows (ranked from the fastest dissipating): lettuce, parsley leaves, carrot root, onion, carrot leaves, leek, and parsley root. The coefficient
kdiss ranged from −0.725 to −0.093, and
t1/2 ranged from 0.96 to 7.42 days. Azoxystrobin belonged to eight active substances most often detected in vegetables and mushrooms obtained during official inspection of pesticide residues in Czech supermarkets in 2018 [
29]. Similarly, azoxystrobin was detected in eight vegetable species from southeastern Poland but never exceeded the MRL [
30]. According to our results, azoxystrobin never exceeded the MRL and in all tested vegetables produced APHI
25 values that were equal to or lower than the PHI, but it was not suitable for zero-residue production, except for lettuce.
The second fungicide, tebuconazole, had significant effects on all vegetables except carrot root and parsley leaves. The dissipation rate of this systemic triazole fungicide was as follows (ranked from the fastest): lettuce, leek, onion, carrot leaves, and parsley root. The dissipation constant kdiss and t1/2 ranged from −0.523 to −0.064 and 1.32 to 10.88 days, respectively. Tebuconazole is not allowed on any of the tested vegetables in the Czech Republic, and our results confirm that it is not suitable for low-residue or zero-residue production.
Spinosad, a fermentation product of the actinomycete bacterium
Saccharopolyspora spinosa, dissipated in descending order: lettuce, carrot leaves, carrot roots, onion, parsley leaves, and leek. In parsley roots, the dissipation model was not significant. Based on the
kdiss (−0.614 to −0.132) and
t1/2 (1.13 to 5.23 days) values, spinosad is a rapidly dissipating substance. In the Czech Republic, spinosad is permitted in lettuce (PHI 14 days), onion, and leek (PHI 7 days). In carrot and parsley, we used a longer PHI of 14 days. In all tested vegetables, spinosad was suitable for low-residue production with APHI
25. In lettuce, carrot root and onion, production with less than 0.01 mg kg
−1 spinosad was feasible. Sikorska-Zimny et al. [
31] calculated the half-life of spinosad in onion, carrot and cabbage and found the values to be 5.2, 3.6 and 2.9 days, respectively. Our calculated half-lives of spinosad in onion (3.58 days), carrot root (2.55 days), and cabbage (3.85 days; data shown in Kocourek et al. [
15]) confirm the fast dissipation process of this substance.
Similarly, thiacloprid, a neonicotinoid insecticide, dissipated fast in almost all tested vegetables (ranked from the fastest): lettuce, onion, carrot root, leek, carrot leaves, and parsley roots, with kdiss values of −0.660 to −0.102 and t1/2 values of 1.05 to 6.81. The parsley leaf data did not enable the generation of any significant model. In the Czech Republic, thiacloprid is permitted in lettuce, carrot, and parsley with a PHI of seven days. In our study, we used the longest PHI of 21 days for leek and onion, as recommended for other vegetable species. In onion, the 21-day PHI was sufficient for 25% MRL as well as for zero-residue production. In the case of leek, a 21-day PHI was suitable for 25% MRL as well as for zero-residue production in one year only. In 2008, the calculated APHI25 and APHI0.01 prolonged the PHI for 6 and 18 days, respectively.
Cypermethrin, a synthetic pyrethroid, dissipated in descending order: lettuce, carrot root, onion, parsley leaves, leek, and carrot leaves, with
kdiss values of −0.507 to −0.051 and
t1/2 values of 1.37 to 11.53. In parsley roots, cypermethrin was below the detection limit. Its dissipation strongly depended on the vegetable species. For example, no cypermethrin residues were detected in carrot roots after foliar application or later [
32]. On the other hand, Yuan et al. [
33] reported a high concentration of cypermethrin in some vegetable samples from supermarkets, where residue levels varied in different vegetables, high levels being found in radish and cauliflower. Cypermethrin was the insecticide observed frequently at concentrations above the MRL in leafy vegetables [
34]. Cypermethrin is the most used pesticide in vegetables, but the residue levels were always below the European MRL [
35].
Two others tested pyrethroids (deltamethrin and lambda-cyhalothrin) dissipated according to vegetable species and their parts. Both were decomposed very fast in onion, carrot, and parsley roots without the ability to create a dissipation curve. Similarly, Riplay et al. [
32] did not detect deltamethrin residues in onion from the first day of foliar application. In lettuce, the
kdiss coefficients −0.539 and −0.404 were obtained for deltamethrin and lambda-cyhalothrin, respectively. On the other hand, in leek, carrot, and parsley leaves, these two active substances dissipated very slowly and are thus not suitable for low-residue and zero-residue production.
Pesticides sprayed at the highest recommended dosages for vegetables in the Czech Republic dissipated under the MRL except for iprodione and methoxyfenozide (onion), acetamiprid, chlorpyrifos and thiamethoxam (leek) and tebuconazole (parsley leaves).
The MRLs of active substances in carrot and parsley roots can even be one thousand times lower than in leaves, e.g., the MRL of spinosad in parsley roots is 0.02 mg kg
−1 and in parsley leaves is 60 mg kg
−1 [
24]. Our model for parsley roots was not robust enough, but the highest detected amount of spinosad was 0.03 mg kg
−1 and in parsley leaves was 0.45 mg kg
−1 three days after spraying. Similarly, the maximum level of spinosad detected in carrot roots and carrot leaves was 0.03 mg kg
−1 and 1.18 mg kg
−1 three days after spraying, respectively. The limit of 60 mg kg
−1 appears to be very mild compared to the lettuce (10 mg kg
−1) and leek (0.2 mg kg
−1) MRLs. Leaves of parsley are often used raw in many dishes. High initial deposits of active substances in combination with mild MRLs allow PHI compliance.
Gonzáles-Rodríguez et al. [
34] showed this problem in lettuce compared to other leafy vegetables. The MRLs for Swiss chards and spinaches can even be one hundred times lower than for lettuce because lettuce is highly sensitive to pests and needs successive applications of pesticides, consequently leaving a higher level of residues that are tolerated. Recently, there was a high decrease in the MRL of dimethoate in lettuce from 25 mg kg
−1 (Reg. (EU) 2015/400) to 0.01 mg kg
−1 (applicable from 31.7.2019; Reg. (EU) 2019/38) [
24]. However, the MRLs of many active ingredients evaluated in lettuce were still much higher than those in leek (
Table 2 and
Table 4).
Models of pesticide dissipation in vegetables are suitable for the regulation of pesticide residues in harvested crops. In this paper, significant models were established for 114 pesticide/crop combinations. According to these models, it is possible to predict the APHI for the requested “action threshold” of residues in crops at the time of harvest. For pesticides with a very low incidence of residues in the short term after treatment, it is not necessary to prolong the PHI. Fast dissipation of residues was found for deltamethrin and lambda-cyhalothrin in onion, carrot roots, and parsley roots and for cypermethrin, metalaxyl-M and pirimicarb in parsley roots. Similarly, a very fast drop in residue levels was found for abamectin in onion and leek. For these pesticides, it is not necessary to establish APHI
25 or APHI
0.01 because the pre-harvest interval is achieved in any case. A strict limit for zero-residue production (0.01 mg kg
−1) is achieved when pyrethroids are applied in root vegetables. In contrast, high deposits of pyrethroids in the leaves of carrot and parsley after treatment do not enable us to achieve this limit. Very low incidence of residues of several pesticides in the short term after treatment was also detected in cauliflower (beta-cyfluthrin) and head cabbage (acetamiprid, beta-cyfluthrin, cypermethrin, deltamethrin, lambda-cyhalothrin, pymetrozine, thiamethoxam) [
15].
The difficulty of vegetable growing for low-residue or zero-residue production depends on the crop and pesticides used. The order of crops for which the established APHI
25 or APHI
0.01 are longer than the PHI corresponds to the order of crops descending from the lowest to the fastest dissipation of pesticides. For low-residue production, the order of crops is based on the number of pesticides with APHI
25 values exceeding the PHIs from the total number of evaluated pesticides for which a significant model was established: leek (7/13), carrot leaves (2/9), iceberg lettuce (3/21), onion (4/29), parsley leaves (0/7), carrot roots (0/7), and parsley roots (0/8). In comparison, the number of pesticides with APHI
25 values exceeding the PHIs was 7 of 20, 8 of 17, and 1 of 18 evaluated pesticides in Chinese cabbage, cauliflower, and head cabbage, respectively [
15]. The growth of Chinese cabbage, leek and cauliflower for low-residue production is difficult, while the growth of iceberg lettuce, onion, head cabbage, parsley and carrot for roots seems to be easier. Growing of parsley and carrot for leaves in low-residue production is also achievable owing to the very mild MRLs discussed above. For zero-residue production, the order of crops is based on the number of pesticides with APHI
0.01 values exceeding the PHIs from the total number of evaluated pesticides for which a significant model was established: leek (13/13), carrot leaves (9/9), parsley leaves (6/7), iceberg lettuce (9/21), onion (11/29), parsley roots (2/8), and carrot roots (1/8). In comparison, the number of pesticides with APHI
0.01 values exceeding the PHIs were 19 of 20, 13 of 17 and 1 of 18 evaluated pesticides in Chinese cabbage, cauliflower, and head cabbage [
15]. The growth of Chinese cabbage, leek, cauliflower, parsley, and carrot for leaves for zero-residue production is difficult, while the growth of head cabbage, iceberg lettuce, onion, parsley, and carrot roots is feasible.
According to the exceeding of the APHI25, it is possible to select the pesticides tebuconazole (onion, iceberg lettuce, carrot leaves, leek), chlorpyrifos (leek, iceberg lettuce), iprodione (onion, iceberg lettuce), cypermethrin (carrot leaves), methoxyfenozide (onion), azoxystrobin, pyridaben, thiacloprid, and thiamethoxam (leek) as unsuitable for low-residue production.
According to the exceeding of the APHI
0.01, the number of “risky” pesticides is higher. Chlorpyrifos, acetamiprid, thiacloprid, and thiamethoxam were the insecticides that most often exceeded the APHI
0.01 in vegetables. Except for acetamiprid, these pesticides were recently restricted by EU regulation [
36,
37,
38]. Azoxystrobin, boscalid, iprodione, spinosad, pyridaben, propamocarb-hydrochloride, and tebuconazole were the fungicides that most often exceeded the APHI
0.01 in vegetables. However, the limited number of pesticides available for the control of vegetables such as leek, Chinese cabbage or cauliflower did not enable the exclusion of these unsuitable pesticides from pest management.
The half-lives of pesticides set employed in this study for treatment of lettuce were lower than half-lives of those when used in onion, leek, and parsley and carrot roots, including the carrot leaves (
Figure 1a). This means that dissipation rates of pesticides in lettuce were higher compared to other vegetable species (
Supplementary Table S5 online). In this context it is rather surprising that relatively high levels of some fungicides (e.g., azoxystrobin) were found in samples of lettuce collected in a Spanish market, and in some of them the MRL was exceeded [
34]. It is worthwhile to point out that MRLs for the same pesticides in other crops are in some cases lower than in lettuce [
24]. For instance, the MRL for azoxystrobin in lettuce is 1.5-times higher than in onion and 15-times higher than in carrot root. The MRL for difenoconazole in lettuce is 10-times higher than in carrot and MRL for mandipropamid in lettuce is even 250-times higher than in onion [
24]. MRLs for leafy vegetable as Swiss chard or spinach can be one hundred times lower than for lettuce [
34]. The reason of such apparent discrepancies in MRLs setting is a high sensitivity of lettuce to pests and thus the need to apply more treatments, which, consequently, leave higher residues. However, in an earlier Spanish study, the residues of azoxystrobin in greenhouse lettuce harvested at PHI day were in all cases below the officially set MRL mentioned in the study. The results indicate that the time between application and harvest is at least as important as the application dose [
8].
It is necessary to mention that high half-life values of some pesticides do not necessarily lead to occurrence of residues at harvest. For instance, thiamethoxam had in lettuce a relatively high level of half-live in comparison to onion (
Table S5A,B). However, in lettuce the thiamethoxam level lowered at harvest below a limit of detection when PHI was kept [
39].
The occurrence of pesticide residues in lettuce in a Spanish market was reported more often than in other vegetables [
34] and may potentially pose a threat to consumer health. Interestingly, fungicides were most detected, especially in lettuce. Although concentrations of insecticides used to be lower, residues of cypermethrin, chlorpyrifos, difenoconazole, or lambda-cyhalothrin were found in 50% of the lettuce samples, while no detectable pesticide residues were present in potatoes and onions. In 20% of the lettuce samples, residues exceeded MRL; moreover, in some samples, this was the case for two or three pesticides [
40].
Pesticide dissipation half-lives in plants are largely dependent not only on their characteristics (e.g., surface morphology) but also on environmental conditions [
13]. Moreover, dissipation rate and half-lives of pesticides are variable for various vegetable species. Therefore, the mathematic models are suitable only for the vegetable species for which they were developed. One of the key factors affecting pesticide half-lives is temperature in the period between application and harvest, with higher temperatures resulting in higher dissipation rates and thus shorter half-lives. Generally, the PHI forecast according to developed models will be more reliable under conditions with higher temperatures. The application of the predictive models on vegetables grown in glasshouses is rather complicated as the temperature conditions are much higher than those for which they were developed; half-lives are fairly lower, what can be also documented by our earlier data comparing the fate of azoxystrobin under field and glasshouse conditions [
8].
It is noteworthy that the developed models of pesticide dissipation for a forecast of PHI and incidence of residues in plant products before harvest are applicable in regions with similar climatic conditions for the growing of that vegetable. For expanded use of dissipation models developed in this study, parameters such as plant characteristics, substance properties and environmental conditions including temperature are summarized in regression models developed by Fantke et al. [
14]. As an example, corrected values of half-live
t1/2 ref,i(day) and
t1/2 plant, active subst. (day) according to these regression models of Fantke et al. [
14] are presented in
Table S7.
Degradation is for many pesticides the most relevant dissipation process, followed by growth dilution [
2]. Growth of plants, depending on water accessibility, is another limiting factor of this model’s reliability. A better forecast power can be expected when water access for vegetable is sufficient, regardless of whether it is due to irrigation or precipitation in the growing locality. On the other hand, limited plant growth increases the share of growth dilution on the overall dissipation.
The major outcome of this study is 114 regression model equations of tested pesticide dissipation of in five vegetable species. These models can be used in practice for monitoring of pesticide residues in products at harvest in dependence on term of application; the knowledge of dissipation rate enables establishing APHI for targeted limits of pesticide residues in product at harvest, e.g. for a limit of 25% MRL or 0.01 mg·kg−1. In addition, based on the models, it is possible to differentiate between the pesticides for which prolonging of the PHI is not necessary and those for which some adjustment of PHI is needed. Models of pesticide dissipation and the procedure of APHI calculation can also serve as a basis for expert systems aimed at the regulation of pesticide residues in vegetables. Under these conditions, the key users might be vegetable growers, extension services or retail chains that require vegetable products with pesticide residues lower than official MRLs. Overall, protection of consumer health will occur due to a reduction of pesticide vegetable contamination by residues. The other challenging application of the developed models is their use for registration purposes and updating of respective MRLs and PHIs. However, it remains unclear if reduced pesticide application amounts will remain effective against the relevant target pests. Hence, understanding the efficacy of reduced pesticide application for controlling relevant pest organisms requires further research.
The outcomes of this study might nevertheless help to identify a selection of candidate pesticides for stricter pesticide legislation, specifically in the case of pesticides with a higher incidence of residues in food crops, longer half-lives and longer APHI0.01. Considering the focus of future research in this field, the monitoring of levels of rapid dissipation of active substances from the first day after application would be useful as it provides important data for comparison of initial deposits in particular crops with the literature data (18). Evaluation of relationships between MRL, PHI, and acute reference doses (ARfDs) is a challenge for further research of pesticide residues in products. Finally, the models characterizing pesticide dissipation can also be used for screening residue-related consumer exposure, thus informing related risk assessments.