SciELO - Scientific Electronic Library Online

vol.48 número3Interferência da variação vertical de fósforo no solo e do estresse hídrico no crescimento de milho, soja e girassolQualidade do corte basal de cana-de-açúcar efetuado por facas de diferentes angulações e revestimentos índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados




Links relacionados


Revista Ciência Agronômica

versão impressa ISSN 0045-6888versão On-line ISSN 1806-6690

Rev. Ciênc. Agron. vol.48 no.3 Fortaleza jul./set. 2017 

Agricultural Engineering

Correlation between 2,4-D herbicide residues and soil attributes in southern of Brazil1

Correlação espacial do residual do herbicida 2,4-D com atributos do solo

Denilson Baumgartner2 

Eduardo Godoy de Souza3  * 

Silvia Renata Machado Coelho3 

Marcio Furlan Maggi3 

2Departamento de Engenharia Agrícola, Universidade Estadual do Oeste do Paraná/UNIOESTE, Campus Cascavel, Paraná-PR, Brasil,

3Centro de Ciências Exatas e Tecnológicas, Universidade Estadual do Oeste do Paraná/UNIOESTE/CCET/PGEAGRI, Campus Cascavel, Rua, Universitária 2069, Jardim Universitário, Cascavel-PR, Brasil, 85.819-110,,,


The 2,4-D herbicide is widely used in Brazilian and global agriculture to control weeds. It is potentially toxic to mankind, and one of the herbicides final destinations is the soil. The aim of the study was to determine the spatial correlation of 2,4-D herbicide in soil attributes after 12 months of application. The experiment was conducted in an agricultural area of 5.04 ha in the municipality of Nova Santa Rosa/PR, in 2011 and 2012. Two soil samples collections were performed: collection 1 (09/03/2011) and collection 2 (08/25/2012), both of which were made 12 months after 2,4-D application. The herbicide was applied in a single dose of 1.5 L ha-1 in 09/2010 and in two doses in 09/2011 (1.0 L ha-1or 1.5 L ha-1 2,4-D). The correlation between soil attributes and 2,4-D was determined with the use of a spatial correlation matrix. The 2,4-D soil residues varied widely from one location to another in the same area, with points of no detection of the active ingredient. It was not possible to detect differences between 1.0 L ha-1 and 1.5 L ha-1 2,4-D treatments. Several factors may have interfered in the herbicide residues in the soil, including: climate, product characteristics, physical and chemical soil characteristics, soil management, among others. The spatial correlation of the soil attributes and 2,4-D residues in the soil proved to be an important tool to verify one parameter interference over another.

Key words: Agrochemicals; Weeds; Spatial distribution


O herbicida 2,4-D é amplamente utilizado na agricultura brasileira e mundial para o controle de plantas espontâneas, sendo potencialmente tóxico ao homem e um dos destinos finais dos herbicidas é o solo. O objetivo do trabalho foi determinar a correlação espacial do herbicida 2,4-D com os atributos do solo após 12 meses de aplicação. O experimento foi desenvolvido em área agrícola de 5,04 ha, do município de Nova Santa Rosa/PR, nos anos de 2011 e 2012. Duas coletas de solo foram realizadas: coleta 1 (03/09/2011) e coleta 2 (25/08/2012), sendo que ambas foram feitas após 12 meses da aplicação de 2,4-D. O herbicida foi aplicado em dose única de 1,5 L ha-1 em 09/2010 e em duas dosagens em 09/2011 (1,0 L ha-1 e 1,5 L ha-1 de 2,4-D). A correlação entre os atributos do solo e o 2,4-D foi determinada com a utilização de matriz de correlação espacial. O residual de 2,4-D no solo variou muito de um local para outro na mesma área, apresentando pontos sem detecção do ingrediente ativo. Não foi possível detectar diferenças entre os tratamentos 1,0 L ha-1 e 1,5 L ha-1 de 2,4-D. Diversos fatores podem ter interferido no residual do herbicida no solo, entre eles: o clima, características do produto, características físico-químicas do solo, manejo do solo, entre outros. A correlação espacial dos atributos do solo e do residual de 2,4-D no solo mostrou-se ser uma importante ferramenta para verificar a interferência de um parâmetro sobre o outro.

Palavras-chave: Agrotóxico; Plantas espontâneas; Distribuição espacial


Concern about effects caused by agrochemicals on the environment has increased due to its intensive use in agricultural areas (RIBEIRO et al., 2013). In addition, the method of soil tillage can have considerable influence on soil erosion, rain infiltration, runoff and leaching Associated with this movement of soil and water are agrochemicals, either bound to soil particles or in a soluble form. The contamination of surface waters with silt, pesticides and nutrients have been frequently found to damage these ecosystems (HOLLAND, 2004; URI; ATWOOD; SANABRIA, 1998). Among widely used agrochemicals, herbicides can be mentioned as being one of the most used methods to control weeds.

The 2,4-D (2,4-dichlorophenoxyacetic acid) is one of the widely used herbicide to control broadleaf weed sin the Brazilian and world agriculture, which is potentially toxic to humans (BOIVIN et al., 2005), thus it is important to know your behavior in the environment because it has the ground as one of the final destinations.

"The distribution and degradation of herbicides in the soil are dynamic and unique processes for each herbicide-soil relation" (ROCHA et al., 2013). In studies performed by several authors, different 2,4-D herbicide half-life periods were verified in soils, which can remain for days up to months (SILVA; VIVIAN; D'ANTONINO, 2009). For example, 2,4-D + picloram residues can persist in the soil up to 2.5 years (NASCIMENTO; YAMASHITA, 2009), according to soil and climate conditions.

Agrochemicals residues can be found in different environment places, such as in soil (VIEIRA et al., 2009), and its environmental behavior can be influenced by sorption, drift, biological and chemical decomposition, volatilization, leaching and runoff processes (SPADOTTO; GOMES, 2013; VIEIRA et al., 1999), besides these processes interaction (SPADOTTO; GOMES, 2013; VIVIAN et al., 2006). Its behavior in the soil is very complex and results of several factors, such as interaction with soil constituents (SILVA et al, 2013; VIEIRA et al., 1999; VIVIAN et al., 2006).

Among attributes that typically have soil-agrochemicals interaction is soil texture, thus it is very important to have its knowledge before agrochemical application. Clay soils have a greater herbicide molecules sorption potential (SILVA; VIVIAN; D'ANTONINO, 2009). In addition to texture, Vieira et al. (1999) argued that organic matter presence also helps in increasing 2,4-D adsorption.

Understand the interaction between herbicides and soil attributes is of fundamental importance in the knowledge of herbicides persistence in the environment. It is important to know the residual potential of the herbicide into the soil, or soil properties which have correlation with the herbicide 2,4-D. Since it is potentially toxic to humans and can cause negative effects on plants, such as interruption of growth (SILVA et al., 2011). Among other problems, Nascimento and Yamashita (2009) reported that soil contaminated with 2,4-D + picloram caused reduction in emergency, height and toxicity in vegetable crops. To understand which attributes may have more interaction, the study of the spatial correlation between herbicide residues and soil attributes is shown as an important tool for its knowledge.

Usually, georeferenced collected samples are used to know soil attributes spatial distribution. For Cavalcante et al. (2007a) and Silva et al. (2003), this sampling allows to spatialize chemical attributes with geostatistics use, which determine the correlation between soil attributes and 2,4-D herbicide residues. In assessing the spatial dependence of agronomic attributes, Deutsch and Journel (1998) emphasized that can make use of geostatistics to describe phenomena that fluctuate over time and/or space. Thematic maps are used for displaying the attributes evaluated, however, the interpolation of values is required for those places where the attribute has not been sampled (VILELA; CAETANO, 2009).

The half-life of the herbicide in the soil can also be considered as an interference factor, because herbicide permanence time in the soil is influenced by many factors, such as temperature, precipitation, soil moisture, among others. Authors mention that the half-life of 2,4-D in the soil may vary from seven days up to several years (MCCALL; VRONA; KELLEY, 1981;. NASCIMENTO; YAMASHITA, 2009; SILVA; VIVIAN; D'ANTONINO, 2009) and this time is dependent on climate factors, soil management conditions, chemical product characteristics, physico-chemical soil characteristics, among others.

Cavalcante et al. (2007a) consider the study of the spatial variability of soil chemical attributes important when the soil is subjected to different managements. The study area has soybeans and corn in rotation systems its main crops, with different NPK formulations for each culture. The system of no-tillage presented significant accumulation of organic matter, phosphorus, potassium and elevation of CEC in relation to the other studied systems, besides the improvement in the chemical conditions of the soil, the organic matter was larger in relation to the natural system.

As the soil is considered a final herbicides reservoir, it is important to determine herbicide residues in it. In this context, the aim of the study was to determine the correlation between 2,4-D herbicide and soil attributes in an agricultural area located in Nova Santa Rosa/PR municipality, 12 months after herbicide application.


The experiment was conducted at Nova Santa Rosa, western part of Paraná state, in an area of 5.04 ha (800 x 63 m) (Figure 1a) that use no-tillage system with corn and soybeans crop rotation for more than 15 years, with central geographic coordinates of 24°22'54" S and 53°56'14" O (datum WGS-84), and an average altitude of 363 m. The lowest part of the area is located in the northeast region and the highest part in the southwest region (Figure 1b).

Figure 1 Experimental area image (Google Earth, 2014) (a) and altimetry map (m) (b) 

The municipality predominant soil is the typical eutrophic RED LATOSOL (EMBRAPA, 2013) and the climate, according to Köppen classification, is Cfa, humid temperate climate with hot summers and infrequent frosts. The monthly precipitation data from 2011 and 2012 were provided by SIMEPAR Technology Institute (Figure 2a).

Figure 2 Monthly precipitation chart of 2011 and 2012 in Nova Santa Rosa/PR region, in mm (a), and spatial distribution of the sampling points (b) 

In the experimental area, 40 sampling points (Figure 2b) were georeferenced using the Trimble Geo Explorer XT 2005 Topographic GPS, with post-processed differential correction. The sampling grid was of 40 x 33 m. Forty soil composed samples were collected for each year (2011 and 2012, Tables 1 and 2) using a Dutch auger at depths of 0 to 20 cm. 8 subsamples (Figure 3a) were collected in each of the 40 points and homogenized, removing approximately 300 g of soil from each sample and accommodating it in an identified plastic bag.

Table 1 Descriptive statistics of granulometry and soil chemical attributes 

Attribute Min Mean Median Maximun SD CV(%) Kur Ass Normal
% Clay 59.0 68.7 68.5 78.0 4.7 6.9 -0.4 (A) -0.1 (a) Yes
% Silt 14.0 21.0 21.0 29.0 3.6 17.3 -0.2 (A) -0.1 (a) Yes
% Sand 7.0 10.3 10.0 15.0 1.8 17.9 0.1 (A) 0.2 (a) Yes
pH in CaCl2 4.4 4.8 4.8 5.3 0.2 4.8 -0.6 (A) 0.1 (a) Yes
Cu (mg dm-3) 6.4 12.4 11.7 20.2 3.6 28.7 -0.6 (A) 0.4 (a) Yes
Zn (mg dm-3) 1.0 4.3 3.5 12.9 3.0 69. 6 0.5 (A) 1.0 (b) No
Fe (mg dm-3) 14.0 22.1 21.0 33.0 4.6 20.7 -0.3 (A) 0.2 (a) Yes
Mn (mg dm-3) 142.0 259.6 261.5 361.0 42.8 16.5 1.3 (A) -0.4 (a) Yes
P (mg dm-3) 4.0 15.4 13.1 60.0 10.0 65 9.3 (C) 2.6 (b) No
H+Al(cmolc dm-3) 4.6 6.2 6.2 7.8 0.8 12.6 -0.6 (A) 0.2 (a) Yes
Ca (cmolc dm-3) 4.4 6.1 5.9 8.5 1.0 15.8 0.8 (A) 0.8 (b) Yes
Mg (cmolc dm-3) 0.8 1.2 1.2 1.8 0.2 17.6 0.5 (A) 0.7 (a) Yes
Al3+ (cmolc dm-3) 0 0.1 0.1 0.4 0.1 118.9 0.6 (A) 1.1 (b) No
K (cmolc dm-3) 0.2 0.5 0.4 1.1 0.2 45.8 0.2 (A) 0.9 (b) Yes
CTCpH7 (cmolc dm-3) 11.7 14.0 13.8 17.2 1.1 7.7 1.2 (A) 0.5 (a) Yes
V% 43.5 55.6 55.4 67.7 6.0 10.8 -0.8 (A) -0.1 (a) Yes
MO (g dm-3) 20.8 34.5 34.9 53.1 5.6 16.1 2.4 (C) 0.5 (a) Yes
pH em CaCl2 4.2 4.6 4.6 5.0 0.2 4.4 -0.7 (A) 0.2 (a) Yes
Cu (mg dm-3) 7.4 16.6 16.7 26.7 5.1 30.6 -0.9 (A) 0.1 (a) Yes
Zn (mg dm-3) 2.2 7.6 6.0 31.5 5.7 74.7 8.7 (C) 2.7 (b) No
Fe (mg dm-3) 22.0 33.7 33.0 51.0 6.3 18.7 0.3 (A) 0.4 (a) Yes
Mn (mg dm-3) 186.0 297.8 301.0 394.0 44.4 14.9 0.1 (A) -0.3 (a) Yes
P (mg dm-3) 9.4 21.9 20.0 51.2 10.5 48.0 0.6 (A) 1.1 (b) No
H+Al(cmolc dm-3) 5.8 7.2 7.2 9.0 0.8 10.9 -0.7 (A) 0.3 (a) Yes
Ca (cmolc dm-3) 3.5 4.9 4.6 8.5 1.2 24.1 2.2 (C) 1.5 (b) No
Mg (cmolc dm-3) 0.8 1.2 1.2 1.6 0.2 14.3 -0.2 (A) 0.4 (a) Yes
Al3+ (cmolcdm-3) 0.0 0.1 0.1 0.5 0.1 81.8 0.1 (A) 0.9 (b) No
K (cmolc dm-3) 0.3 0.5 0.5 1.0 0.2 35.1 -0.1 (A) 0.8 (b) Yes
CTCpH7 (cmolc dm-3) 11.8 13.8 13.5 17.1 1.2 8.5 1.3 (A) 1.2 (b) Yes
V% 35.9 47.7 47.7 63.1 6.5 13.7 -0.2 (A) 0.50 (a) Yes
MO (g dm-3) 22.8 31.7 30.9 39.7 3.9 12.3 0.1 (A) 0.2 (a) Yes

Min: minimum; Max: maximum; SD - Standard Deviation; CV - Coefficient of Variation; ASS - Symmetry: Symmetric (a); positive asymmetry (b); negative asymetry (c); KUR - Kurtosis: Mesokurtic (A); Platikurtic (B); Leptokurtic (C). Normal: Normality

Table 2 Detection frequency and descriptive statistics of 2,4-D residues (µg kg-1) for the 40 points 

Year Min Mean Median Maximum SD CV(%) Kur Ass Normal Detection Frequency
Application Dose... 1.0 L ha-1 1.5 L ha-1
2011 0 83 59 609 128 155 7 (C) 2 (b) No - 53%
2012 0 92 0 652 191 207 3 (C) 2 (b) No 40% 10%

Min: minimum; Max: maximum; SD - Standard Deviation; CV - Coefficient of Variation; ASS - Symmetry: Symmetric (a); positive asymmetry (b); negative asymetry (c); KUR - Kurtosis: Mesokurtic (A); Platikurtic (B); Leptokurtic (C). Normal: Normality

Figure 3 Soil sampling for 2,4-D residual analysis in the soil (a) and herbicide 2,4-D application interval in 2011 and 2012 (b) 

Samples 1 (2011) and 2 (2012) were made 12 months after 2,4-D herbicide application. In Sept/2010 a single 1.5 L ha-1 dose of 2,4-D herbicide (DMA*806 BR) was applied throughout the area, and 12 months after (09/03/2011), collection 1 was held in the 40 sampling points. In Sep/2011, after collection 1, two herbicide doses (Figure 3b) were applied: Interval 1, in sampling points 1-20, with a 1.5 L ha-1(1.005 g e.a ha-1) herbicide dose application; and Interval 2, in sampling points 21-40, with a 1.0 L ha-1(670 g e.a. ha-1) 2,4-D dose application, and 12 months after herbicide application, collection 2 was performed, in 2012 (08/25/2012) (Figure 3a).

Soil samples for 2,4-D residual analysis were dug, air dried and protected from light for a period of 5 to 7 days, being then sieved with a 1.2 mm mesh and frozen until 2,4-D liquid chromatography (HPLC) determination. For 2,4-D extraction, the following steps proposed by Amarante Junior (2003) and Kashyap et al. (2005) were adapted: 1) sieve the soil to 1.2 mm and weigh 5 g; 2) add acidified water with H2SO4(1 mL, pH=1.0; 9 M) + dichloromethane (15 mL); 3) stay in the ultrasound system for 1 hour; 4) Collect and reserve the extract; 5) add dichloromethane back to the ground (15 mL) for a new extraction; 6) collect the extract along with the previously reserved; 7) Filter the extract by suction with anhydrous sodium sulfate; 8) dry evaporate the solvent on a rotary evaporator up to 1 mL; 9) dry with N2; 10) Dissolve the residue in 1 mL methanol, accommodate in Eppendorf and refrigerate until injection on HPLC.

Chromatographic conditions were adapted from Amarante Junior (2003) and Kashyap et al. (2005), according to the following conditions: oven temperature: 40 °C; mobile phase flow: 1.0 ml min-1; run time: 23.01 min; wavelength used: 225 nm; mobile phase: 60% MilliQ H2O and 40% TFA acetonitrile; Sample injection volume 20:µL.

After HPLC 2,4-D herbicide determination in the soil, herbicide frequency detection was verified, and with Minitab version 16 software, residues classical descriptive analysis was performed. The coefficient of variation (CV) was considered low when CV ≤ 10% (homogeneity); Average when CV was between 10 and 20%; high when CV was between 20 to 30%; and very high when CV> 30% (heteroskedasticity) (PIMENTEL-GOMES, 2009). Asymmetry and kurtosis coefficients were compared with the confidence interval provided by Jones (1969) at 5% probability level. The normal data probability distribution was verified using Kolmogorov-Smirnov and Anderson-Darling tests at a significance level of 5%. Data that presented normality in at least one of the tests were considered normal.

To evaluate herbicide 2,4-D residues spatial distribution, residual thematic maps were generated with Surfer 10 software, using the inverse interpolator of the square of the distance.

The spatial correlation of 2,4-D herbicide residues with soil attributes (Table 1) was calculated using Moran's bivariate spatial autocorrelation statistics Equation 1, (BONHAM; REICH; LEADER, 1995) using the Software for Definition of Management Zones - SDMZ (BAZZI et al., 2013.). The spatial correlation matrix checks at 1% and 5% level and identifies the attributes that influence the yield positively or negatively, and checks whether a sample is correlated spatially (spatial autocorrelation).

IYZ=i=1nj=1nWij*Yi*ZjWmY2*mZ2 (1)

where, IYZ - Degree of association between Y and Z variables, ranging from -1 to1, as it is followed: positive correlation IYZ > 0 and negative correlation IYZ < 0; Wij - correspond to the ij element of spatial association matrix, calculated by Wij = (1/(1 + Dij)), so that Dij is the distance between i and j points; Yi - transformed Y value at point i; W - the sum of degrees of spatial association, obtained from Matrix Wij, for ≠ ji;mY2 - sample variance of variable Y; mZ2 - sample variance of variable Z.


The 2,4-D residues in soil found in 2012 (collection 2) presented a smaller number of points with 2,4-D detection than residues found in soil in 2011, possibly due to the soil being drier in function of low precipitation in the collection 2 period (Figure 2a). This difference in detection points can be better evaluated by the median (Table 2). In 2012, the median was 0 µg kg-1, i.e., herbicide residues were not detected in most sampling sites. Besides the period of hydric stress, other factors that may have interfered with the detection of 2,4-D residues in the soil can be surface runoff, leaching, and the degradation of the active ingredient. For Schaffrath et al. (2008), runoffs do not occur evenly, as it runs across preferential paths.

Another factor to consider is soil collection samples depth (0-20 cm), because the herbicide could be displaced by mobility to deeper soil regions. Pinheiro, Moraes and Silva (2011) detected 2,4-D herbicide with the highest frequency (56%) at a depth of 80 to 100 cm, showing mobility according to application time.

According to Malone, Shipitalo and Meek (2004), macropores produce preferential paths. Hence, water runs through empty spaces, carrying pesticides to soil deeper regions. For Pinheiro, Moraes and Silva (2011) this factor may explain concentration differences found for the herbicide.

Vieira et al. (1999) consider the soil as one of herbicides final reservoirs. In a study of cucumber plants poisoning, Santos et al. (2006) found that the higher the dose of mixed 2,4-D + picloram and the shorter the time after herbicide application, the greater the residual effect in the soil. The average 2,4-D residues in the soil was of 83 µg kg-1in collection 1 (2011) and 92 µg kg-1in collection 2 (2012). In 2012, higher point values than in 2011 were found, despite detection frequency being higher in 2011. In both years, the minimum residual values found were 0 0 µg kg-1(nd) and the maximum values were of 609µg kg-1in 2011 (collection 1) and 652 µg kg-1in 2012 (collection 2).

Residues CV found in soil is considered very high (155% and 207%) (Table 2). One of the factors that contributed in having a very high CV is the large number of points without residues detection. The maximum residues amount of 2,4-D in the soil encountered by Pinheiro, Moraes and Silva (2011) was of 845.5 µg kg-1in the depth of 0 to 20 cm. In Brazil, the CONAMA Resolution No. 357 of March the 17 th 2005 and the CONAMA Resolution No. 396 of April the 3 rd 2008 determine that the acceptable limits of the final 2,4-D residues in water for human consumption are 30 µg L-1.

2,4-D residues in soil in collection 1 presented detection frequency (Table 2) of 40% in points 21 to 40, and of 65% in points 1 to 20. In collection 2 (two 2,4-D doses), herbicide frequency detection (Table 2) was of 40% in points 21 to 40 and of 10% in points 1 to 20. Pinheiro, Moraes and Silva (2011), while evaluating 2,4-D residues in Haplumbrept soil (clayey and medium texture) had also found concentration values variable at the same depth (0-20 cm), with active ingredient detection frequency of 33%. In the same soil class and at the same depth, Pinheiro et al. (2009) found a 2,4-D detection frequency of 21% in the soil, i.e., herbicides residues values in soil may vary depending on several factors, including soil characteristics and climatic conditions. An ANOVA could not been made in order to test if there were statistical differences between 1.0 L ha-1 and 1.5 L ha-1 2,4-D treatments because both of samplings presented spatial dependence (Table 3).

Table 3 Correlation matrix of the chemical attributes, texture and 2,4-D herbicide residues in soil for 2011 and 2012 

Al3+ Sand Clay CTC Cu H+Al K Mg Mn MO P pH RI-2,4-D Sil Zn

Legend: RI-2,4-D = = 2,4-D residuesin soil

Negative correlation in 2011 and 2012

Positive correlation in 2011 anda 2012

By means of the average and the maximum values of residual 2,4-D was not possible to detect differences between 1.0 L ha-1 and 1.5 L ha-1 2,4-D treatments.

By mapping the spatial distribution of 2,4-D herbicide residues in soil for 2011 (collection 1) (Figure 4a), it was found that most of the experimental area had residual amounts of 2,4-D ranging from 0 (nd) to 122 µgkg-1. However, there was a residual concentration of 2,4-D in the northeastern part of the experimental area, with intervals ranging from 244 to 609 µgkg-1, in the property lower part. According to Spadotto and Gomes (2013), after agrochemicals application, its environment behavior may be influenced by several processes, including transport (drift, volatilization, leaching and runoff), i.e., the 2,4-D applied at the property higher part may have shifted to the lower part by runoff. Schaffrath et al. (2008) consider soil as heterogeneous, thus, presenting high soil physical properties spatial variability.

Figure 4 Spatial distribution of 2,4-D residues (µg kg-1) in 2011 (collection 1) (a), 2012 (collection) (b) and residues difference between 2012 and 2011 (c). Map coordinates are in meters (Universal Transverse Mercatur - UTM) 

In 2012 (collection 2) (Figure 4b), there were residues concentration on the extreme and central parts of the experimental area, with maximum intervals of residual herbicide ranging 391 to 652 µg kg-1 and minimum interval ranging from 0 (na) and 261µg kg-1, or concentration limits are found in the search. During this period there was a lower precipitation rate (Figure 2a). In the 2,4-D residues difference map (Figure 4c) it is observed that there was an increase in herbicide residues at the extremes and central parts of the experimental area in 2012, although, 12 months before collection 2, a different 2,4-D herbicide application was performed. This spatial variability was commented by Pinheiro, Moraes and Silva (2011), saying that environmental conditions in agricultural areas are dynamic and heterogeneous, influenced by the soil spatial and temporal distribution. Hence, agrochemicals may be lost by degradation, evaporation or organism absorption.

Analyzing the spatial correlation between soil attributes, it was found that many attributes have different correlation coefficients from one year to another, and many attributes have not presented spatial correlation. Table 3 presented the attributes which have in both years, 2011 and 2012, the same behavior. According to Deon (2007) attributes may have different correlation levels from one year to another, probably due to nutritional disorders that may have been induced by improper and excessive fertilization management, among which toxicity and induced deficiency of one nutrient by another.

According to Cavalcante et al. (2007b), in the no-tillage system, the V% in depths of 0-10 cm and 10-20 cm, and the pH and Ca in 10-20 cm depths showed a strong spatial dependence. This detected spatial dependence may be derived from any of the soil formation factors, especially source material and relief, while moderate spatial dependence would be derived from soil homogenization that fertilization and liming provided to it. In this study the spatial dependence was found for Cu, Al H +, Zn and 2,4-D soil residues.

Binding of agrochemicals, and other attributes such as Al3+ + with clay is related to the particle size, then the smaller the particle of the soil is, the higher the reactivity, and clay is the soil reactive fraction, as it is the smallest granulometric fraction of the soil. According to Procopio et al. (2009), clay favors 2,4-D adsorption. In accordance with Silva et al. (2011) Al3+ showed positive spatial correlation with the clay content (2011 and 2012). In addition to the particle size, electric charges may also influence the binding, i.e., clays have negative charges and Al3+ positive charges, thus, it tends to bind. In 2011 and 2012, Al3+correlation was negative with silt granulometric content.

Cu spatial correlation was positive with sand, Mn, P, and Zn. Mn showed a positive spatial correlation with sand and Zn. As Mn is an element that was found above the soil necessary limit to plant development (30.0 mg dm-3) (COSTA; OLIVEIRA, 2001), the influence of it on the organic MO was higher. It was also found that the chemical attribute K showed a negative spatial correlation with Cu, Mg, and P, and Mg and had a positive spatial correlation with Zn.

Evaluating 2,4-D residues spatial correlation in the soil, there was a negative correlation with sand and positive with clay, in agreement with Silva, Vivian and D'Antonino (2009), who stated that, with the increase of clay, 2,4-D residues increase in the soil may occur, as 2,4-D sorption in clay particles is higher, since clay soils normally have a higher herbicides molecules sorption potential, possibly because it is the reactive granulometric fraction of soil. No spatial correlation of 2,4-D herbicide with pH was found, although Johnson, Lavy and Gbur (1995) state that 2,4-D adsorption decreases with pH.

Although the soil MO of the experimental area has not presented spatial correlation with 2,4-D herbicide residues in the soil, this attribute in conjunction with soil texture is one of the main attributes to consider in herbicides choosing and recommended doses, as it can favor an increase in adsorption of 2,4-D to the soil (VIEIRA et al., 1999). In both years, 2,4-D residues present spatial correlation (positive or negative) with at least one granulometric attribute evaluated, demonstrating the importance of knowing soil texture in the selection of the herbicide to be applied in the area. Agreeing with Silva et al. (2011), who evaluated 2,4-D residual activity in soybean emergence in soils of different textures and verified that there were residual activity variations due to dosage and time period between soybean planting and application, presenting a more pronounced effect on plants grown on medium textured soils (20% clay, 4% silt and 76% sand).

With 2,4-D herbicide determination results in the soil in 2011 and 2012 in depths of 0 to 20 cm, it was found that with the time this is a herbicide that may remain in the soil although residues were not detected at certain points of the area. Depending on each soil characteristics, permanence time will be higher or lower, i.e., detection may be related to chemical, physical, and management soil factors.


  1. In both years it was found only cross-correlation between 2,4-D herbicide and clay and sand;

  2. 2,4-D detection in the soil varied greatly from one location to another in the same area, presenting points without the active ingredient detection;

  3. It was not possible to detect differences between 1.0 L ha-1 and 1.5 L ha-1 2,4-D treatments;

  4. The cross-correlation matrix using the Moran's bivariate spatial autocorrelation statistics proved to be an important tool to verify the interference of one attribute over another. To do so, a chemical and physical soil properties continuous verification must be performed, besides monitoring agrochemicals residual activity, as these can be toxic to plants.

1Trabalho extraído da Dissertação de Mestrado do primeiro autor


The authors acknowledge the support received from the State University of West Paraná (Unioeste) and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).


AMARANTE JUNIOR, O. P. et al. Determination of 2,4-dichlorophenoxyacetic acid and its major transformation product in soil samples by liquid chromatographic analysis. Talanta, v. 60, p. 115-121, 2003. [ Links ]

BAZZI, C. L. et al. Management zones definition using soil chemical and physical attributes in a soybean area. Engenharia Agrícola, v. 33, n. 5, p. 952-964, 2013. [ Links ]

BOIVIN, A. et al. 2,4-Dichlorophenoxyacetic acid (2,4-D) sorption and degradation dynamics in three agricultural soils. Environmental Pollution, n. 138, p. 92-99, 2005. [ Links ]

BONHAM, C. D.; REICH; R. M.; LEADER, K. K. Spatial cross-correlation of Boutelouagracilis with site factors. Grassland Science, v. 41, n. 1, p. 196-201, 1995. [ Links ]

CAVALCANTE, E. G. S. et al. Variabilidade espacial de atributos químicos do solo sob diferentes usos e manejos. Revista Brasileira de Ciência do Solo, n. 31, p. 1329-1339, 2007b. [ Links ]

CAVALCANTE, E. G. S. et al. Variabilidade espacial de MO, P, K e CTC do solo sob diferentes usos e manejo. Ciência Rural, v. 37, n. 2, p. 394-400, 2007a. [ Links ]

COSTA, J. M.; OLIVEIRA, E. F. de. Fertilidade do solo e nutrição de plantas: culturas: soja-milho-trigo-algodão-feijão. 2 ed. Campo Mourão: COAMO; Cascavel: COODETEC, 2001. 93 p. [ Links ]

DEON, M. D. Crescimento e nutrição mineral da soja submetida a excesso de P, S, K, Ca e Mg em solução nutritiva. 2007. 71 f. Dissertação (Mestrado em Agronomia)- Universidade de São Paulo, Escola Superior de Agricultura "Luiz de Queiroz", Piracicaba, 2007. [ Links ]

DEUTSCH, C. V.; JOURNEL, A. G. GSLIB: geostatistical software library and user's guide. New York: Oxford University Press, 1998. 369 p. [ Links ]

EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA. Serviço Nacional de Levantamento e Conservação de Solos. Sistema brasileiro de classificação de solos. 3. ed. Brasília: Embrapa Produção de Informação, 2013. 353 p. [ Links ]

HOLLAND, J. M. The environmental consequences of adopting conservation tillage in Europe: reviewing the evidence. Agriculture, Ecosystems & Environment, v. 103, p.1-25, 2004. [ Links ]

JOHNSON, W. G.; LAVY, T. L.; GBUR, E. E. Sorption, mobility, and degradation of triclopyr and 2,4-D and four soils. Weed Science, v. 43, p. 678-684, 1995. [ Links ]

JONES, T. A. Skewness and kurtosis as criteria of normality in observed frequency distributions. Journal Sedimentary Petrology, v. 39, p. 1622-1627, 1969. [ Links ]

KASHYAP, S. M. et al. Rapid Analysis of 2,4-D in soil sample by modified Soxhlet apparatus using HLPC with UV detection. Journal of Chromatographic Science, v. 43, p. 81-86, 2005. [ Links ]

MALONE, R. W.; SHIPITALO, M. J.; MEEK, D. W. Relationship between herbicide concentration in percolate, percolate breakthrough time, and number of active macropores. Transactions of the ASAE, v. 47, p. 1453/1456, 2004. [ Links ]

MCCALL, P. J.; VRONA, S. A.; KELLEY, S. S. Fate of uniformly carbon-14 ring labeled 2,4,5-Trichlorophenoxyacetic acid and 2,4-Dichlorophenoxyacetic acid. Journal of Agricultural Food Chemistry, v. 29, p. 100-107, 1981. [ Links ]

NASCIMENTO, E. R.; YAMASHITA, O. M. Desenvolvimento inicial de olerícolas cultivadas em solos contaminados com resíduos de 2,4-D + picloram. Semina: Ciências Agrarias, v. 30, n. 1, p. 47-54, 2009. [ Links ]

PIMENTEL-GOMES, F. Curso de estatística experimental. 15. ed. Piracicaba: FEALQ, 2009. 451 p. [ Links ]

PINHEIRO, A. et al. Estudo da presença de pesticidas no perfil do solo, sob diferentes tipos de culturas. Revista Brasileira de Recursos Hídricos, v. 14, n. 3, p. 51-59, 2009. [ Links ]

PINHEIRO, A.; MORAES, J. C. S.; SILVA, M. R. Pesticidas no perfil de solos em áreas de plantação de cebolas em Ituporanga, SC. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 15, n. 5, p. 533-538, 2011. [ Links ]

PROCÓPIO, S. O. et al. Utilização do herbicida 2,4-D na dessecação de manejo em lavoura de soja no sistema de plantio direto. Magistra, v. 21, n. 3, p. 187-193, 2009. [ Links ]

RIBEIRO, A. C. A. et al. Resíduos de pesticidas em águas superficiais de área de nascente do Rio São Lourenço-Mt: validação de método por extração em fase sólida e cromatografia líquida. Química Nova, v. 36, n. 2, p. 284-290, 2013. [ Links ]

ROCHA, P. R. R. et al. Meia-vida do diuron em solos com diferentes atributos físicos e químicos. Ciência Rural, v. 43, n. 11, p. 1961-1966, 2013. [ Links ]

SANTOS, M. V. et al. Eficácia e persistência no solo de herbicidas utilizados em pastagem. Planta Daninha, v. 24, n. 2, p. 391-398, 2006. [ Links ]

SCHAFFRATH, V. R. et al. Variability and correlation among physical properties of an Oxisol under no-till and conventional tillage. Revista Brasileira de Ciência do Solo, v. 32, p. 1369-1377, 2008. [ Links ]

SILVA, A. A.; VIVIAN, R.; D'ANTONINO, L. Dinâmica de herbicidas no solo. 2009. Disponível em: <>. Acesso em: 31 maio 2012. [ Links ]

SILVA, F. M. L. et al. Atividade residual de 2,4-D sobre a emergência de soja em solos com texturas distintas. Revista Brasileira de Herbicidas, v. 10, n. 1, p. 29-36, 2011. [ Links ]

SILVA, G. R. et al. Sorption of fomesafen in brazilian soils. Planta Daninha, v. 31, n. 4, p. 971-977, 2013. [ Links ]

SILVA, V. R. et al. Variabilidade espacial das características químicas do solo e produtividade de milho em um Argissolo Vermelho-Amarelo distrófico arênico. Revista Brasileira de Ciência do Solo, v. 27, n. 6, p. 1013-1020, 2003. [ Links ]

SPADOTTO, C. A.; GOMES, M. A. F. Perdas de agrotóxicos. Disponível em: <>. Acesso em: 27 jul. 2013. [ Links ]

URI, N. D.; ATWOOD, J. D.; SANABRIA, J. The environmental benefits and costs of conservation tillage. Science of the Total Environment, v. 216, n. 1/2, p. 13-32, 1998. [ Links ]

VIEIRA, E. M. et al. Estudo da adsorção/dessorção do ácido 2,4 diclorofenoxiacético (2,4d) em solo na ausência e presença de matéria orgânica. Química Nova, v. 22, n. 3, p. 305-308, 1999. [ Links ]

VIEIRA, E. O. et al. Avaliação da contaminação do carbofuran nos solos do Distrito de Irrigação do Gorutuba. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 13, p. 250-256, 2009. [ Links ]

VILELA, L. C.; CATANEO, A. Análise do método de interpolação krigagem ordinária aplicado à pesquisa agronômica. Revista Energia na Agricultura, v. 24, n.1, p.85-104, 2009. [ Links ]

VIVIAN, R. et al. Persistência de sulfentrazone em Argissolo Vermelho-Amarelo cultivado com cana-de-açúcar. Planta Daninha, v. 24, n. 4, p. 741-750, 2006. [ Links ]

Received: December 02, 2014; Accepted: October 03, 2016

*Autor para correspondência

Creative Commons License This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.