Agriculture Canada Logo

1989 - 1994

Predicting Pesticide Migration Through Soils
of the Great Lakes Basin

R. de Jong, W.D. Reynolds, S.R. Vieira and R.S. Clemente
Centre for Land and Biological Resources Research, Research Branch,
Agriculture and Agri-Food Canada, Ottawa, Ontario, Canada


Download Report   (151 KB pdf )



Pesticide contamination of ground water has traditionally been considered to be due primarily to spills, and to improper storage, disposal and application practices. There is increasing evidence, however, that normal agricultural practices can also result in low-level, non-point source contamination of ground water via the downward migration of pesticides through the soil profile (Agriculture Canada, 1990). Although this type of contamination is usually well below Canadian drinking water guidelines, there are growing public concerns over potential health hazards related to long-term exposure to low levels of pesticides (Agriculture Canada, 1990). Consequently, there is a need to determine how important and widespread low-level non-point source pesticide contamination of ground water might be; what the controlling soil, land use and weather factors are; and which agricultural practices are required to mitigate and control this type of pollution at acceptable and sustainable levels.

Essential steps in obtaining this information include identification of the primary mechanisms controlling pesticide movement through the soil profile, and development of the capability to characterize and predict the pesticide movement in space and time with acceptable accuracy. Accordingly, the objectives of this study are: (i) to develop a methodology for predicting, characterizing and quantifying pesticide migration rates through the soil profile, and (ii) to conduct a test and evaluation of the methodology by applying it to atrazine migration through the soil profiles of the Grand River watershed in Southern Ontario, Canada.

The herbicide atrazine was chosen for study because it is the most commonly used pesticide for corn production in Southern Ontario; and because atrazine residues have been detected in the ground waters of many agricultural watersheds in the Great Lakes basin, particularly where there is some combination of high atrazine usage, intensive agriculture, high precipitation, coarse textured and other highly permeable soils, high water tables, and sloping topography (Millette and Torreiter, 1992). The watershed scale was chosen (rather than the field scale, for example) because it is a natural landscape unit; because the most readily available and complete soil and weather databases are applicable on the watershed scale; and because the watershed provides a convenient basis for estimating agrochemical loadings to the Great Lakes. The Grand River watershed (Fig. 1) was chosen as a test case because it is one of the largest in Southern Ontario ( 680,000 ha); it contains a large range of soil textures with a complexity of distribution that is typical for the region; the primary land use is field crop production using standard agricultural practices and "normal" rates of pesticide usage (Shelton et al., 1988); and it empties directly into the Great Lakes (Lake Erie). These features make this watershed ideal for testing and demonstrating the methodology. In addition, the results obtained for the Grand River watershed should be "characteristic" of the entire Southern Ontario region.

Ground water contamination in this study is defined as non-zero values of predicted annual mass loading of pesticide to the 90 cm depth in the soil. Non-zero mass loading was used, rather than pore water concentrations above a specified threshold, because of the need to estimate the quantities and distributions of all pesticide additions to ground water, not just the "high level" additions. The 90 cm depth was selected because it reflects the mean tile drain depth in Southern Ontario, as well as the primary rooting depth for most field crops. It was assumed that if a pesticide reached this depth, it would not be intercepted by tiles and roots, but continue to percolate downward and eventually enter the ground water.



A comprehensive assessment of the accuracy and validity of the predictions is not possible due to a lack of appropriate field measurements. There are, however, sufficient field data available to get a general indication of the plausibility of the predictions; as well as an indication of the sensitivity of the predictions to the quality and quantity of input data.


i) Comparison to Ground Water Survey Data

The Grand River watershed predictions are consistent with the results of a recent survey of ground water quality in Southern Ontario (Agriculture Canada, 1992). The survey, which included over 1000 farm water supply and monitoring wells (approximately 110 of these wells fell within the Grand River watershed), found that non-point atrazine contamination of ground water was infrequent, highly variable spatially, and generally low-level. Atrazine and d-ethyl atrazine were detected in < 7 % of the wells, with maximum concentrations of 18 ppb and mean concentrations 1.1 ppb. The survey also found that the former Canadian drinking water guideline for atrazine (60 ppb) was never exceeded; that the USEPA drinking water standard (3 ppb) was exceeded only occasionally (1 % of the atrazine detections were above the 3 ppb concentration); and that the contamination was not strongly related to soil type or land use, implying control by many interacting soil, weather and land use factors rather than one or two dominant factors (see Reynolds et al., 1994 for further detail). This lends credibility to the LEACHM-Kriging-GIS methodology, notwithstanding that more extensive comparisons with field data are required before definite conclusions can be drawn.

ii) Effect of Scale and Missing Data

The effects of map scale and missing soil hydraulic properties on the predicted atrazine loadings were briefly assessed using a subregion of the Grand River watershed where the soil data were more complete and available at a much more detailed scale. The subregion consists of an approximately 7400 ha section of Haldimand-Norfolk county falling within the watershed. The subregion contained 359 soil polygon centroids at 1:45,000 scale, and 15 different soil types. For 13 of the 15 soils, values of sand, silt and clay content, OC, BD, Ks and 2-3 points on the soil water characteristic were available for 4 soil layers extending to 100 cm depth. Consequently, much less estimation of data via pedotransfer functions was required, and distances between polygon centroids were much smaller, than for the original 1:1 million scale data. The subregion was kriged on a 200 m x 200 m grid to produce 1781 georeferenced grid points of soil, weather and atrazine loading data.

The 1:45,000 scale data yielded a mean atrazine loading of 0.035 mg/m2/yr, a range of 0-0.781 mg/m2/yr and a CV of 323 %. The original 1:1 million scale data within the subregion gave a mean loading of 0.041 mg/m2/yr, a range of 0.037-0.044 mg/m2/yr and a CV of 4.1 %. The closeness of the two mean loadings (within 17 %) suggests that the methodology can still give reasonable overall predictions when map scales are coarse and substantial amounts of data must be estimated. The much smaller range and CV of the 1:1 million scale results indicate, however, that considerable detail is lost when coarse scales are used.


Although the LEACHM - pedotransfer function - geostatistics - GIS methodology is still under development, the results of this study are encouraging. The required input data was extractable, or derivable (via pedotransfer functions) from information archived in the NSDB and AWDB databases. The pedotransfer function, kriging and ILWIS manipulations were effective and sufficiently robust to accommodate small map scales and a high percentage of missing data. Predictions of potential ground water contamination by atrazine for the Grand River watershed are plausible and compare favourably with recent ground water survey results. It is consequently felt that this methodology will prove very useful in the development of agricultural practices and guidelines that maintain agrochemical inputs to the groundwater at acceptable and sustainable levels.


A new methodology, which consists of solute transport modelling in combination with geostatistical analyses and a geographic information system, was developed for predicting, characterizing and quantifying non-point source contamination of ground water due to the migration of agrochemicals through the soil profile. The methodology can be used to develop inventories of "pollution potential" and to develop agricultural land use practices and guidelines that mitigate and control agrochemical inputs to the ground water at acceptable and sustainable levels. It can also be used to predict the impact of changes in land management practices and land use, using so-called "what if" scenarios. The methodology can be applied to a range of agrochemicals, and to virtually any landscape unit (e.g. plot, field, watershed, region), providing the necessary data are also available at that scale.


Predictions were made of atrazine migration through the soil profiles of a "representative" agricultural watershed within the Great Lakes Basin (the Grand River watershed, Southern Ontario). The results indicate that potential non-point atrazine contamination of ground water is highly variable spatially, generally low level, and determined by complex interactions among several soil, weather, crop and solute transport factors. Continuous corn cropping with annual applications of atrazine is predicted to produce a low, but steady, loading of atrazine to the ground water, and hence to the Great Lakes.


Although the solute transport modelling-geostatistics-geographic information system methodology still requires further development and testing, the preliminary results indicate that it should be useful for: pesticide licensing and usage guidelines; land use planning and management; and the development of agricultural practices that control agrochemical pollution of ground water. Consequently, potential users of the methodology include scientists, planners, managers and policy makers.


Several important factors were not considered in the Grand River watershed application, including land use patterns, crop rotations, annual variation in water table depth, topography, and the simultaneous transport of several agrochemicals and metabolites. These factors are likely to be important in the Great Lakes Basin, and should be taken into account. Most watersheds, especially those in the Great Lakes Basin, have substantial non-agricultural areas (e.g. 25 % of the Grand River watershed is used for non-agricultural activities) and this will obviously affect the amounts and distributions of agrochemical inputs to the ground water. Land use and crop rotations not only affect water movement and water content distributions in the soil profile (through crop water use), but also determine the type, amount, timing and frequency of application of fertilizers and pesticides. In humid regions, the depth to the water table can vary from virtually zero at spring thaw to 3 m or more in late summer.

Thus, the distance agrochemicals must travel to enter the ground water varies substantially throughout the year. Run-off and run-on of water, solutes and sediment due to variations in topography have a strong impact on the amount and spatial distribution of water and agrochemical entry into the soil. Any particular agricultural practice (e.g. continuous corn cropping) is likely to contribute a variety of agrochemicals and metabolites to the ground water (e.g. fertilizer nitrate, atrazine plus its main metabolite d-ethyl atrazine, metolachlor, etc.), rather than a single chemical. Except for topography, the methodology in its present form can account for all of the above factors through adjustments and additions to the various input data files. A run-off - run-on based routine that accounts for topography has not been developed.

The representation of "bypass flow" of solute in LEACHP (Eq. 5-7, Appendix A) is simplistic and may be inadequate in soils where extensive bypass flow occurs (e.g. see Table 8, Woodslee field site, Reynolds et al., 1994). Improved, soil property based, representations of bypass flow should be developed and added to LEACHP so that early arrival of agrochemicals to the water table can be detected.

The laboratory column studies suggest that the current form of LEACHP may overestimate the concentration of atrazine in solution, possibly due to an underestimate of effective atrazine dissipation rates (see Reynolds et al., 1994 for details). Further investigations of pesticide - soil interactions should therefore be conducted so that more accurate representations of pesticide transformation and dissipation can be incorporated into LEACHP. Only a very cursory assessment of the accuracy and uncertainty of the predictions has thus far been attempted. Major sources of uncertainty that require further investigation include:

i) NSDB database

As mentioned above, many of the required soil data for the Grand River watershed were missing from the NSDB. In addition, many of the values that are present (e.g. Ks) are estimates made by soil survey personnel, rather than actual measurements. Consequently, the use of the NSDB database may be limited in some watersheds and for certain applications.

ii) Accuracy and precision of the pedotransfer functions.

The accuracy and precision of pedotransfer functions should be clearly established before they are used. For example, the largely texture - OC based relationships used in the Grand River watershed application do not account for soil structure which is known to have a strong impact on soil hydraulic properties. Consequently, the accuracy and/or precision of some of these functions may be rather low.

iii) Use of the dominant soil type in the landscape polygons.

Only the dominant soil in the landscape polygons was used for the LEACHW and LEACHP simulations because the distribution of soil types within the polygons was not available. However, the subdominant soil, which can occupy up to 30 % of the polygon area, may strongly influence, or even control, water and agrochemical movement. Procedures should be developed to account for both the dominant and subdominant soils when determining the solute transport characteristics of a polygon.

iv) Effect of map scale.

The 1:1 million map scale was used for the Grand River watershed predictions because that scale is compatible with the majority of the NSDB data. This scale may not be appropriate for certain applications, because it is too coarse to yield the required detail in soil properties and chemical transport behaviour. Criteria should be developed for matching the appropriate map scale to the intended use of the predictions.

Obtaining the required input data, and validation of the predictions, are important and difficult aspects of applying the LEACHM-Kriging-ILWIS methodology. Appropriate, high quality, field measured input data (e.g. soil hydraulic properties, water table depths, dispersivities, partition coefficients, dissipation rate constants) are very scarce; and appropriate data for ground truthing the predictions (e.g. soil water content profiles, agrochemical concentration profiles) are even more scarce. A catalogue of all such data should be compiled so that this and similar methodologies can be calibrated, tested and assessed before they are used as management, regulatory or policy making tools.



  Table of Contents | GLWQ HOME


Thursday, May 05, 2011 07:19:02 AM