- Open Access
A GIS-based method for household recruitment in a prospective pesticide exposure study
© Allpress et al; licensee BioMed Central Ltd. 2008
- Received: 20 September 2007
- Accepted: 30 April 2008
- Published: 30 April 2008
Recent advances in GIS technology and remote sensing have provided new opportunities to collect ecologic data on agricultural pesticide exposure. Many pesticide studies have used historical or records-based data on crops and their associated pesticide applications to estimate exposure by measuring residential proximity to agricultural fields. Very few of these studies collected environmental and biological samples from study participants. One of the reasons for this is the cost of identifying participants who reside near study fields and analyzing samples obtained from them. In this paper, we present a cost-effective, GIS-based method for crop field selection and household recruitment in a prospective pesticide exposure study in a remote location. For the most part, our multi-phased approach was carried out in a research facility, but involved two brief episodes of fieldwork for ground truthing purposes. This method was developed for a larger study designed to examine the validity of indirect pesticide exposure estimates by comparing measured exposures in household dust, water and urine with records-based estimates that use crop location, residential proximity and pesticide application data. The study focused on the pesticide atrazine, a broadleaf herbicide used in corn production and one of the most widely-used pesticides in the U.S.
We successfully used a combination of remotely-sensed data, GIS-based methods and fieldwork to select study fields and recruit participants in Illinois, a state with high corn production and heavy atrazine use. Our several-step process consisted of the identification of potential study fields and residential areas using aerial photography; verification of crop patterns and land use via site visits; development of a GIS-based algorithm to define recruitment areas around crop fields; acquisition of geocoded household-level data within each recruitment area from a commercial vendor; and confirmation of final participant household locations via ground truthing. The use of these procedures resulted in a sufficient sample of participants from 14 recruitment areas in seven Illinois counties.
One of the challenges in pesticide research is the identification and recruitment of study participants, which is time consuming and costly, especially when the study site is in a remote location. We have demonstrated how GIS-based processes can be used to recruit participants, increase efficiency and enhance accuracy. The method that we used ultimately made it possible to collect biological samples from a specific demographic group within strictly defined exposure areas, with little advance knowledge of the location or population.
- Aerial Photography
- Pesticide Exposure
- Ground Truthing
- Residential Proximity
- Recruitment Area
Pesticide exposures and health outcomes
The U.S. Environmental Protection Agency has defined pesticides as substances that are intended to prevent, destroy, repel or mitigate pests. Pests include insects, animals, weeds, fungi and microorganisms such as bacteria and viruses . Although pesticides are associated primarily with the agricultural industry, many household products are also pesticides. These include flea and tick treatments, cockroach sprays and baits, disinfectants and sanitizers, and lawn and garden products.
Animal bioassays have shown many pesticides to be carcinogenic and some are known or suspected human carcinogens. Agricultural Health Study results have indicated that prostate and ovarian cancer incidence rates are elevated for male and female pesticide applicators, respectively . Pesticides have been linked with cancers in children including leukemia, neuroblastoma, non-Hodgkins lymphoma, Wilm's tumor, Ewing's sarcoma, soft-tissue sarcoma, and cancers of the brain, colorectum and testes . The link to cancer in children is of particular concern as young children are at greater risk from pesticide exposure than adults are, due to their higher skin to body mass ratio, hand-to-mouth activities, crawling and floor play that bring them into increased contact with soil and dust [4–6]. A high tendency to spend time indoors also increases possible exposure as residue is likely to persist longer indoors because sun, rain and soil microbial activity are not present to aid degradation . With non-point-source pesticide contamination from agricultural applications occurring in water supplies, soil and airborne particulate matter [8–12] and only 10–15% of pesticides applied reaching target crops  it is important to develop accurate methods of measuring pesticide exposure. Recent advances in GIS technology and remote sensing (e.g. satellite images and aerial photography) have provided new opportunities to collect ecologic data on pesticide exposure.
Geographic information systems, pesticide exposure studies and household recruitment
A number of studies have used GIS and remote sensing technologies to aid in agricultural pesticide exposure research. Many of these incorporated spatial functions such as distance measurement, buffering, and overlay analysis. Most used historical data on crops and their associated pesticide applications and measured residential proximity to agricultural fields.
Schreinemachers  used county-level data on wheat acreage to examine birth defects in the upper Midwest. Other studies used residential proximity to agricultural fields as a surrogate for exposure [14–16]. Area-based measures have been used to estimate crop acreages within specified buffer zones and, in some cases, model probabilities of pesticide use [17–20]. Data from California's pesticide use reporting (PUR) database have been used to examine a number of health outcomes, using a combination of distance- and area-based measures [21–26]. In some studies, pesticide use data has been linked to available land use maps [24, 25]. Brody et al. have used GIS-based exposure modelling to examine historical patterns of pesticide use on Cape Cod [27, 28].
Due to time and financial constraints, or an historical timeframe, very few of these studies collected environmental and/or biological samples from household members. The challenges faced by researchers in the analysis of household samples include: 1) identification and recruitment of study participants, 2) estimation of exposure and/or verification of pesticide use, and 3) accurate measurements of distance between study homes and pesticide application sites. Many pesticide study populations are drawn from disease registries or health outcomes databases, such as: state cancer registries [21, 23, 28], birth records or birth defects registries [14, 17, 20, 25], Medicare databases , and population-based disease studies [18, 19, 24].
There are some exceptions, however. For their study of organophosphate metabolite concentration, Royster et al.  recruited children at local clinics where measles, mumps, rubella (MMR) shots were being administered. Koch et al.  contacted potential participants for their pesticide exposure study at Women, Infants and Children (WIC) clinics in central Washington State.
The identification of study fields and recruitment of study participants is a complex process that is often glossed over in reports of pesticide research. In this paper, we present a cost-effective, GIS-based method for crop field selection and household recruitment in a prospective pesticide exposure study. For the most part, our multi-phased approach was carried out in a research facility, but involved two brief episodes of fieldwork for ground truthing purposes. The method presented here is part of a larger study designed to examine the validity of indirect pesticide exposure estimates by comparing measured exposures in household dust, water and urine with records-based estimates that use crop location, residential proximity and pesticide application data. The timing of field identification, household recruitment and collection of biological samples was critical to our study, due to the half-life of atrazine, which is about 60 days in topsoil.
The use of GIS was central to the process of household recruitment. In the initial stages, we used it for the identification of potential study fields, the spatial definition of household recruitment areas, and mapping of study households. In later stages (not reported here), we used it for distance and directional measurements. The methods that we used to recruit household members ultimately made it possible to collect biological samples from a specific demographic group within strictly defined exposure areas, with little advance knowledge of the location or population.
Study area selection
Atrazine Use in Illinois, 1999–2003
% of Corn Acreage with Atrazine
Identification of potential study (candidate) fields
For inclusion in the study, a field had to be: 1) planted in corn during the spring of 2003, as biological and environmental samples would be collected after atrazine applications; 2) at least 1000 meters away from a "conflict field," i.e. any other field planted in corn that could also be treated with atrazine, thus potentially making determination of exposure source problematic, and 3) within 200 meters of a moderately populated residential area. Our study team used aerial photographs from the Illinois Natural Resources Geospatial Data Clearinghouse  to visually identify potential field locations. We downloaded and examined 284 Digital Orthophoto Quadrangles (DOQs), taken in 1998 and 1999, from 29 counties. Because aerial photographs are taken during the winter months, when leaves are off deciduous trees, it was not possible to identify specific crops. Fields were identified as large clearings with little or no vegetation.
Our preliminary review indicated that the best field candidates would be on the edge of medium to large towns, thus allowing for an appropriate housing density moving away from the field's edge. From these photographs, we identified 95 potential fields from 17 counties. Given the large number of cornfields in Illinois, one might ask why our examination included so many counties. The reason is that cornfields are so ubiquitous that it was very difficult to find "isolated" fields, i.e. fields far enough away from conflicting fields.
Without ground truthing, there was no way to determine which crops were being grown in each field. Field selection was reliant on aerial photography taken in 1998 and 1999 because processed, up-to-date satellite imagery was not available for our large study area and our project budget did not include funds for image processing. Additionally, many of the areas identified for possible study were on the edges of towns experiencing population growth and new housing development. Because such construction could completely change the study area as it had been defined, or remove the selection of a field as a candidate altogether, it was impossible for the study team to predict which fields would actually be eligible for the study without traveling to Illinois to confirm crop patterns and current land use.
We conducted fieldwork in fall 2002 and spring 2003. The purpose of the 2002 visits was to verify crop and land use patterns and eliminate fields that were not candidates. This would provide us with adequate planning time to recruit study participants prior to the spring 2003 planting. Due to the fact that corn and soybeans are typically rotated annually as a "best management practice" for better crop yields [33, 34], we sought to identify 2002 fields with soybean residue, assuming that they would be planted in corn the following spring. Hence, we used crop rotation practices as a predictor of field utilization for 2003.
At the end of the two initial field trips, 34 potential fields in 10 counties remained in the study. Those removed from consideration either had corn residue on the fields at the time of observation, suggesting that soybeans would be the crop for the following year, or the location was no longer a crop field.
Spatial definition of recruitment areas
The next step was to identify households eligible for participant recruitment. Based on information from the fall 2002 site visits, boundaries for some fields were redrawn and new buffers created. Our goal was to develop an appropriate sample size of children living within each of the buffer zones (200, 500, 800, 1000 meters). This required the following information: 1) household composition – occupants' ages, type of dwelling (e.g. single family, apartment); 2) household location, i.e. latitude and longitude; and 3) contact information – address and phone number. After contacting marketing list vendors, we chose a vendor that provided information about all households fitting the study requirements that were within a specified radius from any known address that we provided. This method was determined to be the most cost-effective, as we could locate an address at the center of the study area and then create a radius to incorporate all the homes within the study buffers. Because fees were based on the number of addresses provided by the vendor, one of our challenges was to develop a method to capture all of the homes within the study areas while eliminating areas that did not meet our requirements, thus keeping costs down.
Next, we had to find a real address as close as possible to this center point. This was a multi-step process. First we overlaid TIGER 2000 streets on the study area and identified the road closest to the center point. Using the address range along the road segment as a guide, we selected a fictitious address on the street and found a corresponding ZIP+4 for the address on the U.S. Postal Service web site . When this was completed for all fields, we sent the ZIP+4 data to the vendor, who provided geocoded addresses for all homes in each ZIP+4. These addresses were mapped, enabling us to identify the address closest to the center point.
Finally, we determined the size of the radius needed to capture all potential candidate households in the study area by buffering the central address in increments of one-tenth of a mile (a requirement set by the vendor) until all of the homes in the study area were contained within the radius. On completion, we provided the central address and radius length to the data vendor.
Mapping and evaluation of potential participants
The vendor provided geographic coordinates and other data about all homes in the specified study area. Using these data we identified all single-family homes with at least one child aged 4 or younger. For each study field, we created a GIS layer consisting of homes that met the study criteria and overlaid this with the aerial photography, buffer zones, and conflict fields. From this procedure we were able to determine the approximate number of potential candidate households in the area, as well as how these households were dispersed across the four buffer zones. In accordance with our study design, we determined that the optimum number of participating households in a study area was six, with two in each of the 0 to 200 m and 800 to 1000 m buffers, and one household in each of the other two buffers. From this initial overlay, it was possible to determine whether any of the fields were unacceptable due to a lack of potential household participants or because of inconsistently-spaced households throughout the study area. Some homes were also excluded because of their proximity to conflict fields. At the end of this process, 24 candidate fields remained, located between Rockford and East St. Louis, Illinois.
Verification of house location and crop
In late spring 2003, we made two field trips to Illinois with the primary goal of verifying that corn was being grown in the candidate fields (and not in any conflict fields). We also needed to ascertain that the geocoded homes of potential recruitment households were located accurately. Mapping the residential location of study participants is an important component of pesticide exposure studies, as it allows researchers to measure distance and direction from agricultural fields. In most cases, address geocoding, which interpolates where an address falls along a street centerline, is used to assign geographic coordinates to residences. Much has been written about the accuracy of geocoding methods for epidemiologic studies [36–38]. When tax parcel databases are available, they can be used with a higher level of accuracy.
Using address lists of all housesholds that met our study criteria, we navigated to each household and compared the geocoded point to the actual house location. Many of the geocoded points were inaccurate. To resolve this, we edited the shapefile using ArcPad and moved each point to its corresponding rooftop on the aerial photograph. Making edits in the field reduced the chance of inaccuracies during data transfer from paper to digital formats. In several cases, we discovered that the correct location of the house put it into a different buffer zone or even outside the study area. The inaccuracy of the original geocoded houses confirmed that the fieldwork verification process was integral to the accuracy of the study's final results. During this process we visited and verified the locations of 304 houses. At the conclusion of this stage of the investigation 21 fields remained viable for the study.
Study Fields and Household Distribution Across Buffer Areas
Field Size in Acres
In recent years, there has been a proliferation of epidemiologic research that utilizes GIS and spatial analysis. One of the greatest challenges of these studies is the accurate assessment of exposure. Too often, proximity is used as a surrogate for exposure, without reliable information to back up its use. The two main objectives of our larger research have been to 1) determine whether there is a clear exposure gradient with distance from corn fields and 2) evaluate the correspondence between predicted vs. actual exposures by comparing record-based methods with biological and environmental samples. The results of our larger study are being reported elsewhere, but the initial first steps of identifying study fields and participant recruitment areas have been reported here.
We successfully implemented several procedures designed to meet the requirements of locating candidate fields and homes. Through aerial photograph interpretation, we were able to identify potential fields and recruitment areas; however, ground truthing was necessary to confirm these selections. Without fieldwork, this portion of the investigation would have been impossible. Once in the field, the utilization of digital technologies with ArcPad aided in data collection and increased accuracy of field delineation and household location.
Using geocoded points provided by marketing list vendors was integral to the success of our task. However, geocoding as a means for precise address location is not infallible. Without ground truthing the geocoded results, we would have been unable to accurately locate homes that possessed the desired demographic characteristics in relation to candidate fields, which would have significantly compromised the project.
Fieldwork always has its challenges and, in this case, timing was critical since it relied on the planting and harvest cycles of crops. We did not anticipate the loss of so many potential candidate fields to surburban housing development. Despite these challenges, we were able to accurately and efficiently place homes within the study area and locate those that met the criteria imposed by the project. None of this would have been possible without the use of GIS and digital spatial technologies.
We have presented a novel approach to the complex process of selecting study fields and identifying households for the recruitment of study participants. However, in the past two-three years, new technologies and data products have become available that could have enhanced our multi-step process. In October 2004, Google Earth was launched. With up-to-date aerial photography, the use of Google Earth may have precluded the need to ground truth geocoded addresses.
In 1996, the U.S. Department of Agriculture (USDA) implemented the Cropland Data Layer (CDL) program . Each year, this program produces raster, georeferenced, categorized land cover data layers for states that are major crop producers. The focus is on cropland and the purported accuracy rate for major crops ranges from 80% to the high 90s. These data layers were not available at the time of our study, but are useful for historical or landuse change analyses.
In order to evaluate the accuracy of Cropland Data Layer corn and soybean classifications for our study, we conducted a retrospective analysis using DOQs of 46 fields, on which we had identified and marked crop types during our 2002 and 2003 field visits. Of these 26 (57%) were classified correctly and 4 fields (9%) were covered by clouds. In general, accuracy was greater for larger fields. The two most common misclassifications were: 1) soy or corn was classified as pasture/grassland and 2) soy and corn pixels were interspersed. Given the rate of accuracy and the small size of some of our study fields, ground truthing would have been necessary to our study even if the satellite imagery data had been available in a more timely manner.
Corn/Soy Rotation, 1999–2003, Illinois
SOY TO CORN
CORN TO SOY
While the Cropland Data Layer focuses on cropland, it also contains information about other landuses, including urbanization. Between the years 1999–2003, Illinois lost an average of 4–5% of its corn/soy acreage annually to urbanization, much of it in the form of new housing. We certainly noted this trend during fieldwork and lost many potential study fields to it.
We have demonstrated how GIS-based processes can be used to recruit participants, increase efficiency and enhance accuracy, resulting in the collection of biological samples from a specific demographic group within strictly defined exposure areas, with little advance knowledge of the location or population. In addition to being cost-effective, our iterative, multi-pronged method provided us with the quality control needed for our study. If we had not used this methodology, we would have comprised quality and accuracy. Inaccurate geocoding would have placed some potential study houses in the wrong buffer and invalidated our sampling strategy. Without field checks, time would have been wasted in geocoding homes in invalid study fields, such as those that were not planted in corn or had disappeared due to suburban development. While our methodology may require refinement, in light of new data and technologies, it provides a valuable framework for study field identification and household recruitment.
This project was made possible by funding from the National Institutes of Health, Grant #5R21CA094717-02.
- U.S. Environmental Protection Agency (EPA): What Is a Pesticide?. http://www.epa.gov/pesticides/about/index.htm#what_pesticide
- Alavanja M, Sandler D, Lynch C, Knott C, Lubin J, Tarone R, Thomas K, Dosemeci M, Barker J, Hoppin J, Blair A: Cancer incidence in the Agricultural Health Study. Scand J Work Env Health. 2005, 31 (Suppl 1): 39-45.Google Scholar
- Zahm S, Ward M: Pesticides and childhood cancer. Environ Health Perspect. 1998, 106 (Suppl 3): 893-908. 10.2307/3434207.PubMedPubMed CentralView ArticleGoogle Scholar
- Lewis R: Human exposure to pesticides used in and around the household. The Effect of Pesticides on Human Health. Edited by: Baker S, Wilkinson C. 1989, Princeton, NJ, Princeton Scientific PublishingGoogle Scholar
- National Research Council (NRC): Pesticides in the Diets of Infants and Children. 1993, Washington DC, National Academy PressGoogle Scholar
- Lewis R, Fortmann R, Camann D: Evaluation of methods for monitoring potential exposure of small children to pesticides in the residential environement. Arch Environ Contam Toxicol. 1994, 26 (): 1-10. 10.1007/BF00212792.View ArticleGoogle Scholar
- Camann D: Investigating the effects of pesticide exposure in the home. Technology Today. 7-1994, June 1Google Scholar
- Chester G, Ward RJ: Occupational exposure and drift hazard during aerial application of paraquat to cotton. Arch Environ Contam Toxicol. 1984, 13: 551-563. 10.1007/BF01056333.PubMedView ArticleGoogle Scholar
- Gaynor J, MacTavish D, Findlay W: Atrazine and metolachlor loss in surface and subsurface runoff from three tillage treatments in corn. J Environ Qual. 1995, 24: 246-256.View ArticleGoogle Scholar
- Simcox N, Fenske R, Wolz S, Lee I-C, Kalman D: Pesticides in household dust and soil exposure pathways for children of agricultural families. Environ Health Perspect. 1995, 103: 1344-1353. 10.2307/3432609.View ArticleGoogle Scholar
- Cohen S: Pesticides in ground water in the United States: Monitoring, modeling, and risks from the U.S. perspective. J Environ Sci Health B. 1996, 31: 345-352.PubMedView ArticleGoogle Scholar
- Wauchope R: Pesticides in runoff: Measurement, modeling, and mitigation. J Environ Sci Health B. 1996, 31: 337-344.View ArticleGoogle Scholar
- Schreinemachers DM: Birth malformations and other adverse perinatal outcomes in four U.S. wheat-producing states. Environ Health Perspect. 2003, 111 (9): 1259-1264.PubMedPubMed CentralView ArticleGoogle Scholar
- Bell EM, Hertz-Picciotto I, Beaumont JJ: Case-cohort analysis of agricultural pesticide applications near maternal residence and selected cases of fetal death. Am J Epidemiol. 2001, 154: 702-710. 10.1093/aje/154.8.702.PubMedView ArticleGoogle Scholar
- Lu C, Fenske RA, Simcox NJ, Kalman D: Pesticide exposure of children in an agricultural community: Evidence of household proximity to farmland and take home exposure pathways. Environ Res. 2000, 84: 290-302. 10.1006/enrs.2000.4076.PubMedView ArticleGoogle Scholar
- Royster M, Hilborn E, Barr D, Carty C, Rhoney S, Walsh D: A pilot study of global positioning system/geographical information system measurement of residential proximity to agricultural fields and urinary organophosphate metabolite concentrations in toddlers. J Expo Anal Environ Epidemiol. 2002, 12: 433-440. 10.1038/sj.jea.7500247.PubMedView ArticleGoogle Scholar
- Xiang H, Nuckols R, Stallones L: A geographic information assessment of birth weight and crop production patterns around mother's residence. Environ Res (Section A). 2000, 82 (2): 160-167. 10.1006/enrs.1999.4009.View ArticleGoogle Scholar
- Ward MH, Nuckols JR, Weigel SJ, Maxwell SK, Cantor KP, Miller RS: Identifying populations potentially exposed to agricultural pesticides using remote sensing and a geographic information system. Environ Health Perspect. 2000, 108: 5-12. 10.2307/3454289.PubMedPubMed CentralView ArticleGoogle Scholar
- Ward MH, Lubin J, Giglierano J, Colt JS, Wolter C, Bekiroglu N, Camaan D, Hartge P, Nuckols JR: Proximity to crops and residential exposure to agricultural herbicides in Iowa. Environ Health Perspect. 2006, 114: 893-897.PubMedPubMed CentralView ArticleGoogle Scholar
- Meyer KJ, Reif JS, Veeramachaneni DNR, Luben TJ, Mosley BS, Nuckols JR: Agricultural pesticide use and hypospadias in eastern Arkansas. Environ Health Perspect. 2006, 114 (10): 1589-1595.PubMedPubMed CentralView ArticleGoogle Scholar
- Reynolds P, Von Behron J, Gunier RB, Goldberg DE, Hertz A, Harnly ME: Childhood cancer and agricultural pesticide use: An ecologic study in California. Environ Health Perspect. 2002, 110 (3): 319-324.PubMedPubMed CentralView ArticleGoogle Scholar
- Reynolds P, Hurley SE, Goldberg DE, Yerabati S, Gunier RB, Hertz A, Anton-Culver H, Bernstein L, Deapen D, Horn-Ross PL, Peel D, Pinder R, Ross RK, West D, Wright WE, Ziogas A: Residential proximity to agricultural pesticide use and incidence of breast cancer in the California Teachers Study cohort. Environ Res. 2004, 96: 206-218. 10.1016/j.envres.2004.03.001.PubMedView ArticleGoogle Scholar
- Reynolds P, Hurley SE, Gunier RB, Yerabati S, Quach T, Hertz A: Residential proximity to agricultural pesticide use and incidence of breast cancer in California, 1988–1997. Environ Health Perspect. 2005, 113 (8): 993206-1000.View ArticleGoogle Scholar
- Ritz B, Costello S: Geographic model and biomarker-derived measures of pesticide exposure and Parkinson's disease. Ann NY Acad Sci. 2006, 1076: 378-387. 10.1196/annals.1371.074.PubMedPubMed CentralView ArticleGoogle Scholar
- Rull RP, Ritz B, Shaw GM: Neural tube defects and maternal residential proximity to agricultural pesticide applications. Am J Epidemiol. 2006, 163: 743-753. 10.1093/aje/kwj101.PubMedView ArticleGoogle Scholar
- Marusek JC, Cockburn MG, Mills PK, Ritz BR: Control selection and pesticide exposure assessment via GIS in prostate cancer studies. Am J Prev Med. 2006, 30 (2S): S109-S116. 10.1016/j.amepre.2005.09.002.PubMedView ArticleGoogle Scholar
- Brody JG, Vorhees DJ, Melly SJ, Swedis SR, Drivas PJ, Rudel RA: Using GIS and historical records to reconstruct residential exposure to large-scale pesticide application. J Expo Anal Environ Epidemiol. 2002, 12: 64-80. 10.1038/sj.jea.7500205.PubMedView ArticleGoogle Scholar
- Brody J, Aschengrau A, McKelvey W, Rudel R, Swartz C, Kennedy T: Breast cancer risk and historical exposure to pesticides from wide-area applications assessed with GIS. Environ Health Perspect. 2004, 112: 889-897.PubMedPubMed CentralView ArticleGoogle Scholar
- Koch D, Lu C, Fisker-Anderson J, Jolley L, Fenske R: Temporal association of children's pesticide exposure and agricultural spraying: report of a longitudinal biological monitoring study. Environ Health Perspect. 2002, 110: 829-833.PubMedPubMed CentralView ArticleGoogle Scholar
- Illinois Agricultural Statistics Service.http://www.agstats.state.il.us/website/reports.htm
- Illinois Natural Resources Geospatial Data Clearinghouse.http://www.isgs.uiuc.edu/nsdihome/
- Woods N, Craig IP, Dorr G, Young B: Spray drift of pesticides arising from aerial application in cotton. J Environ Qual. 2001, 30: 697-701.PubMedView ArticleGoogle Scholar
- David Pike, University of Illinois Pesticide Impact Assessment Program, personal communication.Google Scholar
- University of Minnesota Extension: Tillage Best Management Practices for Corn-Soybean Rotations in the Minnesota River Basin.http://www.extension.umn.edu/distribution/naturalresources/DD6676.html
- United States Postal Service.http://www.usps.gov
- McElroy J, Remington P, Trentham-Dietz A, Robert S, Newcomb P: Geocoding addresses from a large population-based study: lessons learned. Epidemiology. 2003, 14: 399-407.PubMedGoogle Scholar
- Bonner M, Han D, Nie J, Rogerson P, Vena J, Freudenheim J: Positional accuracy of geocoded addresses in epidemiologic research. Epidemiology. 2003, 14: 408-412.PubMedGoogle Scholar
- Ward MH, Nuckols JR, Giglierano J, Bonner MR, Wolter C, Airola M, Mix W, Colt JS, Hartge P: Positional accuracy of two methods of geocoding. Epidemiology. 2005, 16: 542-547. 10.1097/01.ede.0000165364.54925.f3.PubMedView ArticleGoogle Scholar
- U.S. Department of Agriculture, National Agricultural Statistics Service: Cropland Data Layer.http://www.nass.usda.gov/research/Cropland/SARS1a.htm
- Babcock BA, Hennessy DA: Getting more corn acres from the Corn Belt. Iowa Ag Review Online. 2006, 12 (4):http://www.card.iastate.edu/iowa_ag_review/fall_06/article3.aspx
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.