Mapping Intra-Field Yield Variation Using High Resolution Satellite Imagery to Integrate Bioenergy and Environmental Stewardship in an Agricultural Watershed

Biofuels are important alternatives for meeting our future energy needs. Successful bioenergy crop production requires maintaining environmental sustainability and minimum impacts on current net annual food, feed, and fiber production. The objectives of this study were to: (1) determine under-productive areas within an agricultural field in a watershed using a single date; high resolution remote sensing and (2) examine impacts of growing bioenergy crops in the under-productive areas using hydrologic modeling in order to facilitate sustainable landscape design. Normalized difference indices (NDIs) were computed based on the ratio of all possible two-band combinations using the RapidEye and the National Agricultural Imagery Program images collected in summer 2011. A multiple regression analysis was performed using 10 NDIs and five RapidEye spectral bands. The regression analysis suggested that the red and near infrared bands and NDI using red-edge and near infrared that is known as the red-edge normalized difference vegetation index (RENDVI) had the highest correlation (R = 0.524) with the reference yield. Although predictive yield map showed striking similarity to the reference yield map, the model had modest correlation; thus, further research is needed to improve predictive capability for absolute yields. Forecasted impact using the Soil and Water Assessment Tool model of growing switchgrass (Panicum virgatum) on under-productive areas based on corn yield thresholds of 3.1, 4.7, and 6.3 Mg·ha showed reduction of tile NO3-N and sediment OPEN ACCESS Remote Sens. 2015, 7 9754 exports by 15.9%–25.9% and 25%–39%, respectively. Corresponding reductions in water yields ranged from 0.9% to 2.5%. While further research is warranted, the study demonstrated the integration of remote sensing and hydrologic modeling to quantify the multifunctional value of projected future landscape patterns in a context of sustainable bioenergy crop production.


Introduction
Higher corn grain yields due to the expansion of production acreage and to better production efficiencies have led to an increase in availability of corn stover as a readily available feedstock for cellulosic biofuels.Corn stover consists of leaves and stalks of a corn plant (Zea mays L.) left in the field after grain harvesting.However, there are constraints on how much corn stover can be removed without negatively impacting subsequent soil productivity [1][2][3][4][5].For example, Thompson and Tyner [6] used a linear programming model to project that at a price of about $66 per metric ton of corn stover, land use in the U.S. Midwest could shift from about 12% to 35% continuous corn.Such cost and farmer supply response projections would lead to more fertilizer use because of the nutrient requirements of a continuous corn agricultural system.The majority of current nutrient contributions to the Gulf of Mexico hypoxic zones is indeed fertilizer application in Midwestern agricultural systems (≈1500-3100 kg• N• km −2 • year −1 ) [7][8][9].The U.S. Department of Agriculture (USDA), the U.S. Department of Energy (DOE), and the U.S. Environmental Protection Agency (EPA) also project cellulosic bioenergy crops such as perennial grasses as important components of the U.S. biomass feedstock portfolio, even in the U.S. Corn Belt [10].
The interest and focus on perennial herbaceous grasses for cellulosic bioenergy feedstock is due to the accrued beneficial impacts on the environment and ecosystem biodiversity as integral components of sustainable biomass.However, the integration of cellulosic bioenergy crops into existing agricultural systems requires a design approach that enhances environmental sustainability and minimally affects current net annual food, feed, and fiber production to trigger indirect land use change (iLUC).One option explored at plot and field scales is the use of marginal agricultural lands [11][12][13][14].Definitions of marginal agricultural lands include sub-optimal production areas for commodity crops [14][15][16][17][18], areas susceptible to environmental degradation (e.g., erosion and nitrate leaching) when put into agricultural production [15,16], and small or irregular fields that are not suited for optimal farm machinery operations (based on discussions with farmers).Intra-field sub-optimal production areas are economically burdensome for farmers and may coincide with nutrient sink hotspots [19] and are therefore ideal for the integration of bioenergy crops into current agricultural systems to support sustainable landscapes and the provisioning of ecosystem services.Other potential uses of marginal lands not explored in this manuscript include sites for solar panels or wind turbines.However, patches of marginal lands are likely too small for those types of energy development.
Remote sensing is a reliable and cost-effective method to forecast crop yields over large areas based on the optical properties of crops such as chlorophyll absorption and mesophyll structure of foliage [20][21][22].Spectral vegetation indices (SVIs) derived from satellite imagery collected from various platforms have been used to predict crop yields across a range of spatial scales.The platforms have included Advanced Very High Resolution Radiometer (AVHRR), the Moderate Resolution Imaging Spectroradiometer (MODIS) [23][24][25], Landsat Thematic Mapper (TM) [23,26,27], Satellite Pour l'observation de la Terre (SPOT) [28,29], IKONOS [30,31], WorldView 2, and RapidEye imagery.For an exhaustive list of satellites and other remote sensing platforms, refer to Konecny [32].
Shanahan, et al. [33] compared the strength of association between SVIs and corn yields in Nebraska.The crops in the study area were irrigated during the growing season.Three SVIs-the normalized difference vegetation index (NDVI, Tucker [34]), the transformed soil adjusted vegetation index (TSAVI, Baret, et al. [35]), and the green normalized difference vegetation index (GNDVI, [36])-were computed using 50-cm resolution multispectral images collected from six different growth-stages: the V6 stage (6-leaves; mid-June), VT stage (Tasseling; early July), R2 stage (blister stage; mid-July), R3 stage (milk stage; early August), mid-August, and in early September.Of all SVIs tested, the GNDVI exhibited the highest correlation (R 2 ; coefficient of determination) with grain yield for the 1997 and 1998 growing seasons.The highest correlation (R 2 = 0.70 in 1997 and R 2 = 0.92 in 1998) were observed during the mid to late August period (R3 stage; during the mid-grain filling stage).Cicek, et al. [37] used NDVI and GNDVI derived from Landsat TM and SPOT at four corn and soybean growth-stages to predict respective yields.GNDVI exhibited more accurate predictions of field aggregated yield (0.81 ≤ R 2 ≤ 0.91) than that of NDVI (0.47 ≤ R 2 ≤ 0.85).These studies demonstrated the effectiveness of high resolution multispectral remote sensing for predicting crop and the possibility of identifying under-productive croplands under corn grain by estimating and mapping yields at a sub-field scale.
The goal of this study was to demonstrate the utility of high resolution remote sensing to inform farmers' decision making on agricultural land use planning at the sub-field scale, in order to facilitate environmental sustainability production of biofuel feedstock.The objectives were to: (1) formulate a regression model for predicting corn yields in a U.S. Midwest watershed using SVIs, more specifically normalized difference band ratio indices, and spectral bands; (2) demonstrate an example of identifying under-productive areas within a plot; and (3) forecast biomass yield and impact on water quality if the low corn yield areas were converted to an alternative bioenergy crop (e.g., switchgrass for this study) using hydrologic model simulations.The study was motivated by the need to develop techniques to obtain data predictive of yields at a spatial scale fine enough to resolve variability of yields and its spatial patterns at the sub-field level.The techniques also require cost-effective and repeatable for predicting crop yield over a long period of time across a watershed.This information would allow farmers to plan at high resolution where they should concentrate their highest production and where they could optimally manage their land for conservation and ecosystem services.For this study, biomass yield refers to cellulosic bioenergy crop yield (e.g., switchgrass yield) and crop yield refers to corn or soybean yield.

Study Area
The Indian Creek watershed comprises three 12-digit hydrologic unit code (HUC) subwatersheds (HUCs: 071300020203, 071300020204 and 071300020205) which drain part of the Vermillion River watershed (8-digit HUC; USGS 07130002) that flows into the Illinois River.The watershed covers parts of three Illinois counties of Livingston, McLean, and Ford (Figure 1) and drains approximately 207.4 km 2 of mostly corn and soybean cropland (>87%) on 104 farms [38].The major environmental concern in the watershed is excessive nitrate from upland cropland and its impact on public water supply use and aquatic life use (>10 mg• L −1 as monitored by USGS gage 05554300).The mean annual rainfall and mean annual total evapotranspiration are 887 mm and 661 mm, respectively [19].The three dominant soils in the watershed included (1) Drummer silty clay loam (16.3% of the watershed); (2) Reddick clay loam (9.5%); and (3) Saybrook silt loam (8.8%).Some of the soils in the watershed are susceptible to nitrate leaching.These constituted 6% of the watershed and included the Symerton silt loam and the Varna silty clay loam.Major crop rotations included continuous corn (10.7%), corn-soybean (31.3%), soybean-corn (31.2%), and corn-corn-soybean (11.0%).Crop rotations were generated by overlaying the 2010 to 2012 National Agricultural Statistics Service (NASS) data [39].The watershed is predominantly under tile drainage due to its low gradient and poorly drained soils; however, there is no readily available data on spatial coverage of the tile-drainage network.

Data
Two remotely sensed data types, RapidEye and National Agriculture Imagery Program (NAIP) imagery, were used for this study (Table 1).Both images were collected over the Vermillion River watershed.The RapidEye image was collected on 26 August 2011, and the NAIP image consists of image frames collected on 26 August 2011 for our study watershed.The RapidEye satellite imagery consists of blue, green, red, red-edge, and near infrared (NIR) spectral bands at 5 m ground sampling distance at nadir.The RapidEye image was orthorectified by the vendor.The NAIP aerial imagery consists of three visible bands and has 1 m spatial resolution.Both image sets contained less than 10% cloud cover within the scenes.Reference yield data were obtained from corn yield map a 6.5 ha site near the City of Fairbury located in the northern part of the watershed (Fairbury field).The site is under mostly Comfrey loam soils with 0-2% slope.The 2011 corn yields at this site varied between 0.6 and 14 Mg• ha −1 with an average of 9 Mg• ha −1 .The georeferenced corn yield map at 15 cm resolution was generated by the John Deere Apex TM Farm Management Software based on the volumetric grain flow data from the GreenStar TM yield monitoring system (John Deere Company; Moline, IL, USA) mounted on a combined harvester.Average yields and intra-field variability at this site were consistent with observed watershed values.Observed average corn productivity in the watershed was 9.1 Mg• ha −1 with a 95% confidence interval (95% CI) of 6.2 to 12 Mg• ha −1 .These were derived from 2010 to 2013 National Agricultural Statistics Service data (USDA-NASS, 2012) county data, weighted by spatial extent of the three counties in the watershed (Livingston, Ford, and McLean).

Image Processing and Model Development
The RapidEye image was spatially registered to the NAIP image using 1000+ ground control points to achieve sub-pixel positional accuracy (RMSE < 0.8 pixel).The image-to-image registration was performed using affine transformation with a nearest neighbor resampling method in the ERDAS IMAGINE 2013 AutoSync workstation.The three visible bands of the NAIP image were aggregated to simulate a 1 m resolution panchromatic image.A pansharpening technique was applied to the 5-m RapidEye multispectral image and the 1-m panchromatic image derived from the NAIP image in order to generate a 1-m resolution multispectral image.Using the 1-m pansharpened multispectral image, pixel values of the five spectral bands-blue, green, red, red-edge, and NIR-were extracted from 1500 randomly selected locations within the Fairbury field.A total of 1060 locations that had reference yield values were used as a training dataset for subsequent analyses.Using all possible two-band combinations in the training data, normalized difference band ratios (or normalized difference indices; NDIs) were computed [40] as: where Ra and Rb are pixel values of spectral band a and b, respectively.NDIs were selected for this study because the indices such as NDVI are widely known to correlate with crop yield [24,41,42].
To develop a predictive model for crop yield, a multiple linear regression analysis was performed on a total of 15 predictor variables-10 NDIs and five spectral bands of the pansharpened RapidEye image.Predictor variables that were statistically significant (p < 0.01) were included in the model.It should be noted that the spectral bands of remotely sensed imagery are expected to be highly correlated.In our dataset, a principal component analysis suggested that the first two bands explain 99.4% of the variance.Spectral transformations such as a normalized difference band ratio provide information that is less correlated than the original bands via nonlinear transformation of two-band combination.Upon developing a regression model, correlation between predictor variables was examined to account for multicolinearity.When variables are highly correlated (r > 0.90), the one that has higher correlation to yield was kept in the model and the rest of variables were discarded.
The predictive model was applied to the pansharpened RapidEye image, and a predictive yield map was generated.The crop yield map was examined for its reliability based on visual assessment for overall agreement in spatial patterns between the reference and predicted yields.Quantitative reliability in prediction was also examined using F-statistics with p-value, coefficient of determination (R 2 ), and standard error (SE) based on 300 samples selected from random locations.To visualize abundance and spatial configuration of productivity within a field, the predictive yield map was classified into productive and under-productive areas by applying arbitrary percentage of total area of the field (e.g., lowest 10%, 20%, and 30%) for a sample study plot.

Forecasting of Bioenergy Crop Impact
A calibrated Soil and Water Assessment Tool (SWAT) model for the Indian Creek watershed, IL was used to forecast the impact of growing switchgrass on sub-optimal corn yield areas in the watershed.Watershed corn yields at a subfield scale were estimated using the above developed predictive corn yield model.For details on the calibrated SWAT model refer to Ssegane and Negri [43].Three assessed scenarios of sub-optimal corn productivity included (1) yields below or equal to 3.1 Mg• ha −1 ; (2) yields below or equal to 4.7 Mg• ha −1 ; and (3) yields below or equal to 6.3 Mg• ha −1 .The three thresholds are equivalent to approximately 30%, 50%, and 70% of observed average corn productivity in the watershed of 9.1 Mg• ha −1 .Next, SWAT was modified to model an alternative landscape where switchgrass replaced current land use and land cover at a subfield level on areas meeting each of the three scenarios.The forecast evaluated the impact of conversion to switchgrass on annual corn and soybean yields, NO3-N and sediment exports, and water yield.Statistical assessments of pairwise differences in nutrient exports under different scenarios were achieved using a non-parametric Mann-Whitney test [44].

Regression Model for Predictive Yield
The regression analysis of 10 NDIs and five spectral bands indicated that four variables-red and NIR bands, a NDI using red-NIR combination (known as NDVI; [34]), and a NDI using red-edge-NIR combination (known as red-edge normalized difference vegetation index or RENDVI; [45])-were significant predictors for crop yield (p < 0.01).Among the four variables, correlation between NDVI and RENDVI were 0.94; thus, RENDVI that had higher correlation with yield (R 2 = 0.733) was selected for the model instead of NDVI (R 2 = 0.711).The final multiple linear regression model for predicting crop yield was as follows: Yield (Mg • ha −1 ) = 0.0009107 × 3 + 0.0001729 × 5 + 4.531 × RENDVI − 1.301 (2) where B3 and B5 are red and NIR spectral bands, respectively.The model showed a large F-value (F = 418) with statistical significance (p < 0.001) and the coefficient of determination (R 2 = 0.543) with standard error (SE) of 1.635.The regression model was applied to the pansharpened RapidEye image to map predictive crop yield.
The predictive yield map of the study plot derived from the model was shown in Figure 2. The predicted yield showed striking similarity to the reference yield estimated at harvest.High yield appeared to be concentrated in the central portion of the field with a couple of small isolated distributions in the southern part of the field.Considerably lower yield was observed in the southern portion section of the field particularly near the mid-section.Overall the model is statistically significant (p < 0.001) with a large F-value (F = 328).Correlation (R 2 ) between the predicted and reference yield was 0.524, and estimated error (SE) in prediction was 1.61.The predictive yield map was classified into two areas-productive and under-productive.Sample classification of the study field shown in Figure 3 were based on the three areal proportions-10%, 20% and 30% of total area.Areas of low productivity are primarily concentrated in the southern portion of the field in a spatially contiguous manner across the three classification maps.The 10% of total area exhibiting low productivity coincides with part of the eroded Symerton silt loam soil series with slope gradients of 5%-10%.With the 30% threshold, under-productive areas coincide with most of the eroded Symerton silt loam soil series with slopes of 5%-10% and parts of the same soil series with slopes of 2%-5%.These soils have less organic matter and are more susceptible to NO3-N leaching.The high yield areas are on the Comfrey loam soils with slopes of 0%-2% and are not susceptible to NO3-N leaching.In addition to the distribution, the maps also indicate size and connectivity of under-productive areas within the field.

Simulated Grain and Switchgrass Yields
The total area in the watershed meeting the three corn yield thresholds of 3.1 Mg• ha −1 , 4.7 Mg• ha −1 , and 6.3 Mg• ha −1 , were 646 ha (3.2% of the watershed), 1487 ha (7.4%), and 2892 ha (14.3%), respectively.When areas whose corn yield productivity were less or equal to 3.1 Mg• ha −1 were converted to switchgrass (scenario 1), annual total grain yields were projected to reduce by 6.9% (Figure 4).However, this land use conversion projected 10,664 metric tons of switchgrass biomass.Conversion of areas at or below corn productivity thresholds of 4.7 Mg• ha −1 and 6.3 Mg• ha −1 were projected to decrease annual total grain yields by 9.4% and 16.6%, respectively.For example, under the second scenario, 19,642 metric tons of grain (Corn: 15,876 and Soybean: 3766) were lost from the baseline 209,810 metric tons.Projected switchgrass biomass was 16,797 metric tons.Yield projections under the third scenario (6.3 Mg• ha −1 threshold) predicted a reduction of 34,747 metric tons (Corn: 24,919 and Soybean: 9828) from baseline annual total grain yields.Projected switchgrass biomass yield under the third scenario was 31,522 metric tons.The figure also shows substantial decrease in sediment yield under all three scenarios (25%-39% reduction).

Forecasted Impact on Water Quality and Quantity
Conversions of watershed areas whose corn productivity was less or equal to 3.1 Mg• ha −1 , 4.7 Mg• ha −1 , and 6.3 Mg• ha −1 , on average were projected to reduce tile NO3-N leachate by 15.9%, 18.4% and 25.9%, respectively (Figure 5a).For sediment, the corresponding respective reductions were 25%, 35.5% and 39% (Figure 5b), while, for water yield, the reductions were 0.9%, 1.1% and 2.5% (Figure 5c). Figure 5 shows that the reductions in tile NO3-N are not significantly different (p > 0.05) for Cases 1 and 2 (3.2% and 7.4% of the watershed).However, the corresponding reductions in sediment export are.This was attributed to a relatively higher reduction in annual surface runoff (92.8 mm vs. 90.7 mm) compared to the increase in annual tile-flow (118.4 mm vs. 120.1 mm).Average reductions in annual NO3-N export and sediment yield were projected to be over 10% and therefore considered hydrologically significant for all three land cover change scenarios because these changes are higher than typical probable errors associated with data measurement and sample collection [46].Projected reductions in water yield under all three scenarios were below 10% and therefore not hydrologically significant.3).The corresponding areal coverages are 3.2% of the watershed, 7.4%, and 14.3%, respectively.The relative reductions are comparisons between current land use (baseline condition) and conversion of areas below the above yield thresholds to switchgrass.

Discussion
The multiple regression model developed using a high-resolution multispectral image showed modest performance for predicting crop yield, yet the predictive yield map showed considerable resemblance with the reference yield map.The coefficient of determination (R 2 = 0.52) and SE = 0.16 suggests that approximately 52% of variability in yield is explained by the three terms in the model (i.e., reflectance values of the red and NIR bands and RENDVI), and a typical deviation of predicated yield from the reference yield is 1.61 Mg• ha −1 .Timing of image collection was considered to have substantial impact on the correlation between SVIs and crop yield.Based on the analysis of SVIs and crop phenology metrics derived from the 250-m MODIS products, Bolton and Friedl [24] determined that approximately 65-75 days and 80 days after greenup were optimal timing of image collection for predicting yield for corn and soybean, respectively.Typical greenup of corn in our study watershed is in late May, and the image data collected on 26 August, was approximately a few weeks later than the optimal date.While high resolution satellite imagery has promise for providing sub-field level information, it is challenging to acquire imagery having desired spectral bands at the optimal time of the year with good sky condition.The use of an aerial platform would permit control over image collection dates and sky conditions, but mission planning and standardization of data over time would be challenging and costly.
The difference between predicted and reference yield most likely stemmed from the indirect nature of yield prediction using remote sensing.Spectral reflectance signatures of crops are primarily influenced by plant's physical properties (e.g., chlorophyll content, mesophyll structure, canopy architecture, and biomass).We predicted yield indirectly by determining a relationship between the spectral reflectance signatures and crop yield [47,48]. Positional misalignment and the difference in spatial resolution between the image and the reference yield data were also sources of uncertainty.The reference yield data was reprojected to improve spatial alignment with the RapidEye image.Some discrepancy could have been introduced during the transformation of the reference data, which could have influenced the predictive power of the model.The reference yield data also exhibited a blocky pattern that was created by multiple spatially contiguous pixels having a common yield value although the yield was recorded at 15 cm resolution.The block size varies across the field ranging from approximately 4 m 2 to 12 m 2 .This pattern appeared as multiple index values falling in a single yield value in Figure 2. The use of appropriate spatial scales in conjunction with careful spatial registration between field measurement and input image is expected to improve accuracy in predictive yield.
Normalized difference indices based on ratios of two multispectral bands were tested for their effectiveness for predicting crop yield in this study.More exhaustive analysis of spectral indices (such as simple ratios and differences) and angle indices derived from the angle formed by three spectral bands in a multispectral signature plot [49], in conjunction with partial least square regression [50] or random forest regression [40] may improve accuracy of yield estimates.Recent examples for yield estimates using remote sensing are provided in Rembold et al. [48].
Theoretically, computing predictive yields across a small watershed using remote sensing is feasible.Because of heterogeneity of soil conditions within the watershed, however, further investigation is warranted to determine how well we can predict yields using a single-date image.Additional information such as crop types, planting date, and fertilizer application for each field [51,52] in conjunction with careful selection of the image collection timing would improve the accuracy of yield prediction.If the use of a single image is found insufficient, it is important to learn about the minimum number of images [24,53] needed for reliable yield prediction in order to keep manageable operational cost.
Spatial comprehensive information of predictive yield, particularly with sub-field yield, provides valuable information for farmers to plan their land use at a very fine spatial scale that meets their objectives and for land and watershed planners to design landscapes to deliver specific conservation targets.The demonstration presented in this study is one example based on the proportional area characterized by low predictive yield.Farmers could use the information to determine a certain proportion of the field for alternative or experimental crops without sacrificing yield.Spatial representation of predictive yield also provides valuable information such as size, proximity, spatial arrangement, and connectivity of under-productive areas (Figure 3), which simple statistics or numerical information does not provide.The predictive yield maps at a resolution that can resolve variability of yields and its spatial patterns would be valuable for farmers in planning their sub-field scale land use in order to achieve their production goals while mitigating impacts on ecosystem services and environment.The information can be used to examine accessibility of areas by large equipment, efficiency of management for multiple areas, and potential impacts on primary crop adjacent to areas that experience alternative treatment/management.There are multiple ways to define under-productive areas such as based on predictive yield or multi-year yield record. Criteria for under-productive areas should be determined based on the objectives.
For this study, we do not make the assumption that the financial loss due to the acreage put out of grain production is directly compensated by the switchgrass yield.Because current market valuation frameworks do not reflect the full benefits of the ecosystem services associated with such land cover conversions.However, current and future work on valuation frameworks of the accrued ecosystem services may reflect the true cost and benefit ratio of these integrated systems and thus make them more acceptable to farmers.

Conclusions
Biofuels derived from renewable biological materials are important alternatives for meeting our future energy needs.A sustainable integration of new bioenergy crops within existing cropping systems requires a landscape design that identifies sub-field scale details to correctly allocate land based on its fertility and productivity to the best matching crops.This design allows the bioenergy crops to perform important ecosystem services in balance with the production of bioenergy feedstock and commodity crops.To facilitate sustainable bioenergy crop production, this study demonstrated the utility of single date, high resolution remote sensing and hydrologic modeling for determining under-productive areas within a plot and examining impacts of using those areas to grow bioenergy crops.The findings are summarized as follows:  The regression model using the red and near infrared spectral bands and the red-edge normalized difference vegetation index showed the best performance for predicting crop yield (R 2 = 0.524; p < 0.001; standard error = 1.61) of all variables tested. The predictive yield map showed that under-productivity areas largely coincide with the eroded Symerton silt loam soil series with slope gradient of 5%-10%. When testing three scenarios of corn-to-switchgrass conversion for areas with predicted corn yield less than or equal to 3.1, 4.7 and 6.3 Mg• ha −1 , total grain yield was projected to decrease from baseline by 6.9%, 9.4% and 16.6%, respectively, and the corresponding switchgrass biomass production was projected to be 10,664, 16,797, and 31,522 metric tons, respectively. Between scenarios 1 and 2, the reduction in tile NO3-N transport was not significantly different (15.9% vs. 18.4% reduction, respectively) while the reduction in sediment yield was (25% vs. 35.5%,respectively).
While there is more work to be conducted, this study is a first attempt to demonstrate the promise of remote sensing technology to deliver actionable data for landscape design, by resolving at an unprecedented resolution, issues related to land marginality, spatially constrained productivity and optimal land allocation to meet specific functional objectives at the decision-making scale.It has been shown that land use decisions regarding privately owned land are made at the field, farm, and watershed scales.Such a methodology would be an invaluable tool in allowing this decision making to prioritize crop selection, specific yield and conservation objectives that are the basis for sustainable agricultural practices.In addition to determining the robustness of the model using multiple plots across the watershed, future studies should include determining the limitation of a single date, high resolution multispectral image for accurately predicting crop yield at a subplot scale, and testing exhaustive spectral indices and other techniques (e.g., partial least square regression and random forest regression) to improve prediction accuracy.

Figure 1 .
Figure 1.Location of the Indian Creek watershed in the State of Illinois showing the dominant land use and land cover, and three year crop rotations.The watershed drains parts of three Illinois counties of Livingston, Mclean, and Ford into the Vermilion River.

Figure 2 .
Figure 2. Predicted corn yield maps generated based on the regression model of the watershed (left) and the study field (center); and reference yield map at harvest (right).

Figure 3 .
Figure 3. Examples of stratification of under-productive land by percent area in Fairbury Field.

Figure 4 .
Figure 4. Grain and biomass yields (corn, soybean, and switchgrass biomass) for the baseline (current land use and land cover) and three suboptimal corn yield productivity scenarios of less or equal to 3.1 Mg• ha −1 , 4.7 Mg• ha −1 , and 6.3 Mg• ha −1 .The error bars are standard errors.

Figure 5 .
Figure 5. Relative reductions in: (a) NO3-N export; (b) sediment yield; and (c) water yield under different areal extents of land use change derived from applying three corn productivity thresholds of 3.1 Mg• ha −1 (Case 1), 4.7 Mg• ha −1 (Case 2), and 6.3 Mg• ha −1 (Case 3).The corresponding areal coverages are 3.2% of the watershed, 7.4%, and 14.3%, respectively.The relative reductions are comparisons between current land use (baseline condition) and conversion of areas below the above yield thresholds to switchgrass.

Table 1 .
Specification of Remotely Sensed Images.