Characterizing Dominant Field-Scale Cropping Sequences for a Potato and Vegetable Growing Region in Central Wisconsin

: Crop rotations are known to improve soil health by replenishing lost nutrients, increasing organic matter, improving microbial activity, and reducing disease risk and weed pressure. We characterized the spatial distribution of crops and dominant ﬁeld-scale cropping sequences from 2008 to 2019 for the Wisconsin Central Sands (WCS) region, a major producer of potato and vegetables in the U.S. The dominant two- and three-year rotations were determined, with an additional focus on assessing regional potato rotation management. Our results suggest corn and soybean are the two most widely planted crops, occurring on 67% and 36% of all agricultural land at least once during the study period. The most frequent two- and three-year crop rotations include corn, soybean, alfalfa, sweet corn, potato, and beans, with continuous corn being the most dominant two- and three-year rotations (13.2% and 8.5% of agricultural land, respectively). While four- and ﬁve-year rotations for potato are recommended to combat pest and disease pressure, 23.2% and 65.9% of potato ﬁelds returned to that crop in rotation after two and three years, respectively. Furthermore, 5.6% of potato ﬁelds were planted continuously with that crop. Given potato’s high nitrogen (N) fertilizer requirements, the prevalence of sandy soils, and ongoing water quality issues, adopting more widespread use of four- or ﬁve-year rotations of potato with crops that require zero or less N fertilizer could reduce groundwater nitrate concentrations and improve water quality.


Introduction
Coarse textured or sandy soils can present significant challenges to agricultural producers due to their low nutrient and water holding capacities [1]. However, with precision nitrogen fertilizer applications [2], irrigation [3], and other management options such as crop rotations to combat pest, weed, and disease pressure [4], growers can overcome key limitations that can limit agricultural production on sandy soils. Irrigation is often necessary on sandy soils in temperate crop growing regions of the United States where the climate provides an ideal growing season length and growing degree days (GDD) and approximately 700-900 mm of precipitation, with about 50% or more falling during the growing season.
In the U.S., one of the leading potato and vegetable (e.g., green peas, sweet corn, snap beans, carrots) production regions is found in the Midwest across central Wisconsin on sandy soils [5]. The Wisconsin Central Sands (WCS) is a region located east of the Wisconsin River and is characterized by its deposits of sand and gravel, often >30 m thick above sandstone and crystalline bedrock [6,7] (Figure 1). The highly permeable sand and gravel layer contains a regional aquifer that sits at a typical depth of 3 to 20 m below the surface [6,8]. Typical recommended management of potato and vegetable crops in this region includes significant nitrogen (N) inputs from inorganic fertilizer (e.g., 145-250 kg N/ha; [1]), and approximately 100-380 mm of supplemental irrigation water (average of 274 mm from N/ha; [1]), and approximately 100-380 mm of supplemental irrigation water (ave 274 mm from 2008 to 2012 for potato [9]) from late May through early September the growing season, which is dependent on the seasonal weather [9]. Addit because of the high value of potato and vegetable crops and significant weed a pressure, a diverse crop rotation is also recommended for the WCS region [4,10]. . This boundary is defined as a contiguous are the Wisconsin River with sand and gravel deposits greater than 50 feet deep (Wisconsin DN Implementing a diverse crop rotation can help build soil organic matter and e soil health by replenishing lost nutrients, while reducing the risk of significant we pest pressure [4]. Crop rotation is important to weed management because it disco the dominance of individual weeds [10] and is the oldest and most effective nonch method in limiting populations of plant diseases caused by soilborne pat [11,12,13,14]. However, rotation effectiveness is dependent on crop sequence a length of time between host crops [15,16].
For high-value crops such as potato, the number of years between plantings same field is extremely important in maintaining a high level of plant health by re pest and disease pressure [4]. Potatoes are vulnerable to many different soilborne d that can cause significant economic losses each year by infecting the roots, tube foliage and causing common diseases such as potato early dying, early blight, co scab, black scurf, and Fusarium wilt [10,17,18]. Previous research has provid necessary information for determining which crops should be planted in rotatio potato, as well as how many years should occur between successive potato crop pl on the same field to maximize suppression of natural diseases and [4,10,12,14,15,16,17,19,20,21]. Implementing a diverse crop rotation can help build soil organic matter and enhance soil health by replenishing lost nutrients, while reducing the risk of significant weed and pest pressure [4]. Crop rotation is important to weed management because it discourages the dominance of individual weeds [10] and is the oldest and most effective nonchemical method in limiting populations of plant diseases caused by soilborne pathogens [11][12][13][14]. However, rotation effectiveness is dependent on crop sequence and the length of time between host crops [15,16].
For high-value crops such as potato, the number of years between plantings on the same field is extremely important in maintaining a high level of plant health by reducing pest and disease pressure [4]. Potatoes are vulnerable to many different soilborne diseases that can cause significant economic losses each year by infecting the roots, tubers, and foliage and causing common diseases such as potato early dying, early blight, common scab, black scurf, and Fusarium wilt [10,17,18]. Previous research has provided the necessary information for determining which crops should be planted in rotation with potato, as well as how many years should occur between successive potato crop plantings on the same field to maximize suppression of natural diseases and insects [4,10,12,[14][15][16][17][19][20][21]. When determining crops to rotate with potato, other root or tuber crops that are susceptible to the same plant pathogens should be avoided as they will continue to harbor disease causing pathogens in the soil [16,22,23]. Incorporating crops that are non-hosts to the same plant pathogens as potato will break the host-pathogen cycle and cause a natural decline in the root zone [4]. The absence of a host for three or more years will allow pathogen populations to naturally decrease, improving soil and microbiome conditions for a host to be grown again [4,19,23]. Recommendations for non-host crops in the WCS include snap beans, peas, sweet corn, field corn, small grains, and alfalfa [10].
Across large regions, documentation of specific agricultural management techniques such as tillage/residue management [24] and crop rotations at the individual farm field level or better has been difficult [25]. Improving our understanding of cropping systems management and land use change can help further research and the use of numerical models to understand the impacts of these land use/land management choices on ecosystems services such as food production, carbon sequestration, or water quality. Most recently, the annual U.S. Department of Agriculture (USDA) Cropland Data Layer (CDL) has afforded scientists a new source of land use data and been used in multiple studies to track largescale agricultural land use and land cover change in great detail [25][26][27][28][29]. The CDL was developed in 1997 by the USDA's National Agricultural Statistical Service (NASS) with the primary focus to use satellite imagery to provide crop acreage estimates for select states [30]. Since 2008, the CDL has been updated annually for the continuous United States at 30-and 56-m resolutions and includes over 100 different crop types.
In this study, we use the USDA CDL dataset to document agricultural land use patterns and dominant crop rotations in the WCS region from 2008 to 2019. The objectives of this work were to: (1) document the spatial distribution (and total hectares) of potato and other vegetable crops each year and their frequency of occurrence over the study period; (2) document the dominant crop rotations of the region, the amount of land area in each, and if any trends occurred; and (3) assess potato management preferences by quantifying the number of hectares of potato that returned after 0 (i.e., continuous potato), 1, 2, and 3 years between successive plantings.

Methods
National USDA CDL files for years 2008-2019 were obtained from the USDA CDL website [31] at 30-m resolution. A shapefile of the WCS, defined as a contiguous area east of the Wisconsin River with sand and gravel surficial deposits greater than 50 feet deep [32], was used to define the entire area of interest for this study and applied to each year using ESRI ArcGIS 10.7.2 [33] during image processing. Through continuous product improvement, the accuracy of the CDL has gradually increased over time and others have suggested that a significant greater increase of cropland expansion has happened recently compared to other data sources [29]. When using the CDL, certain precautions and adjustments should be taken to account for product biases when calculating total acres and cropland changes over time. A complete list of CDL limitations and recommended methods of analysis and use can be found in Lark et al. [29].
One important limitation that is pertinent to this study given crop types grown in the WCS, is that during 2017 and 2018 the CDL's crop identification model incorrectly classified all sweet corn as field corn (Patrick Willis, NASS, personal communication, 4 September 2019). This error resulted in sweet corn not being recorded for those two years and its hectares were combined with field corn. To account for this error, rotations that included sweet corn were averaged over the remaining years and that average value was used to replace zero-hectare values in 2017 and 2018. This average was then subtracted from the equivalent corn rotation. These adjusted values are used in the current study.

Classification of Dominant Crop Types in Individual Fields
Additional common misclassification errors in the CDL can be caused by field boundaries and 30m pixels often not aligning, creating a sub-pixel area bias or a false edge. Low spots in fields, contour strips, and field corners outside of irrigation tracks-due to differences in reflective properties-can also create non-homogeneous classification for entire fields that can contribute to over-or under-estimation of crop acreage. To address this issue, we integrated additional datasets in combination with the baseline CDL [29]. We applied the 2008 Common Land Unit (CLU) data layer that outlines field boundaries and assumes that the majority of pixels within the boundary is the land cover for the entire field. Due to restrictions in the 2008 Farm Bill, 2008 is the most recent CLU version available for scientists to use. Changes to field boundaries, as well as use of buffer and contour strips since 2008 would not be accounted for from a single snapshot. However, despite possible changes, the CLU still adds valuable boundary information that can help minimize uncertainty surrounding field edges as well as improve the categorization of crop type (and pixel counts) within similarly managed fields [25,28,29]. By combining the CLU with the CDL, changes in crop type can be evaluated on a field-by-field basis each year with higher confidence [25]. The majority crop type was calculated using zonal statistics from the spatial analysis toolbox within ArcMap.
Additional challenges exist in the CDL's inability to distinguish among crop types that exhibit spectral signatures that are similar to grassland vegetation. By consolidating raster classifications into categories, this can help eliminate misclassification errors [25,29,30]. For this study, we combined within ArcMap the classes of barley, spring wheat, winter wheat, rye, oats, millet, other hay/non-alfalfa, and camelina and labeled these as cereal (Ce) crops. The CDL also does not categorize snap bean and other edible beans. Snap beans are commonly grown in the WCS region, but for this study, the category of beans (B) will be used to represent snap, kidney, and all other edible beans found in the region. The other crops we were most interested in for this study region include corn (C), soybean (S), sweet corn (Sc), potato (Po), alfalfa (A), and peas (Pe). All other crops identified in the region were consolidated and labeled other (O). These eight main crops were chosen because of their dominance within the WCS region.

Classification of Multiyear Crop Rotations
Within any given year, a field is simultaneously embedded in several multiyear rotation sequences across years (depending on starting and ending points) that can be broken down into shorter rotations that may repeat over a long time period as mathematically demonstrated by Hennessy [34]. For this study, we have 12 consecutive years of crop classification data from 2008 to 2019. We adapted a method by Hennessy [34] that uses a moving n-year rotation window where we break down the continuous 12-year crop sequence into multiple year rotations including two-and three-year rotations. For example, the 12-year rotation sequence, PeScBPoScBScPoScPePoSc, can be broken down into eleven two-year rotations (PeSc, ScB, BPo, PoSc, ScB, BSc, ScPo, PoSc, ScPe, PePo, and PoSc) and ten three-year rotations (PeScB, ScBPo, BPoSc, PoScB, ScBSc, BScPo, ScPoSc, PoScPe, ScPePo, and PePoSc). These simpler rotations can show which crops are most likely to be grown in combination together and how common these are adapted across the region. The process to calculate when potato returns in a rotation used a similar moving window approach but also includes rotations of four and five years in the analysis. Calculations of when potato returns in rotation only considered rotations that began and ended with potato. Analysis of four-and five-year rotations that began and ended with potato were a focus of this study given there is preference and recommendation to have at least three years of other crop types planted before returning to potato to break the plant pathogen cycle [4]. All possible rotation combinations with potato that could be created using the remaining crop classes, including the other crop category, were calculated, and totaled for each of the two-, three-, four-, and five-year combined datasets. All rotations in the multiyear crop rotation analysis were identified using the Combine tool within the Spatial Analyst toolbox in ArcMap. The tool combines the years together and creates a table of crop combinations, or rotation sequence, and the corresponding pixel counts where unique rotations occurred on the landscape. With the improvement of the CDL each year and varied identification accuracies across crop types, converting pixels to hectares requires the use of an independent dataset to adjust the value of individual pixels so that when totaled, they agree with county-or statelevel statistics [29,30,35]. Published NASS statistics were used to determine an adjustment multiplier unique to each crop and year in this study. The first crop in the desired rotation determined the multiplier which is then applied to the pixel count to calculate the ending total of hectares for that rotation.

Ecological Landscapes
Wisconsin is divided into 16 ecological landscapes, with 5 occurring within the study area. These landscapes were created by the Ecosystem Management Planning Team (EMPT) based on the National Hierarchical Framework of Ecological Units (NHFEU) to create ecological sections that are similar in climate, geology, topography, soils, and vegetation [36]. These ecological boundaries available for public use by the Wisconsin DNR were used in ArcMap to determine the spatial distribution of each crop.

Crop Type Spatial Distributions, Area, and Frequency
In the WCS, 40.2% of the total land area (285,300 ha of 708,200 ha) was in agricultural use. On average, from 2008 to 2019 the eight crops reported in Table 1 account for 64% of the agricultural land in the region. Distribution is greatly influenced by the five ecological landscapes that make up the region: Central Sand Hills, Central Sand Plains, Forest Transition, Central Lake Michigan Coastal, and Southeast Glacial Plains. The spatial distribution of several of the most dominant individual crop types (corn, potato, sweet corn, and alfalfa) in the WCS and their planting frequency during the 2008-2019 period are shown in Figure 2. Corn, soybean, cereal, and alfalfa were the only crops that are found throughout all five ecological landscapes. Potato, sweet corn, beans, and peas are found concentrated in Central Sand Plains, Central Sand Hills, and Forest Transition ( Figure 2). Table 1     Corn was planted at least one time during the study period on 190,189 hectares, or 66.7% of agricultural land in the WCS, and had an annual average of 70,616 ha (24.8% of total agricultural land) ( Table 2). One percent of agricultural land (2025 ha) was planted with continuous corn. Corn was found in all five ecological landscapes, with 75% occurring between Central Sand Plains and Central Sand Hills regions. Common frequencies range from three to eight years in all regions except for Central Sand Plains which has the highest amount of land with corn only grown once during the study period (16,336 ha or 24% of corn grown in the region). Sweet corn was planted at least one time on 78,482 hectares or 27.7% of agricultural land, with an average of 20,915 ha planted annually or 7.3% of all agricultural land ( Table 2). Sweet corn is found in three ecological landscapes, with 70.9% (40,536 ha) of all sweet corn grown in the Central Sand Plains. Individual land unit frequencies are as high as nine total years (frequency totals were not adjusted for missing acres) but are commonly between one through four years. Potato was planted at least once on 66,257 hectares, or 23.4% of agricultural land, with 72.4% (47,988 ha) occurring in the Central Sand Plains. Potato was also grown in the Central Sand Hills (22.2% of total potato plantings, 14,679 ha) and Forest Transition (5.4%, 3556 ha). Potato is annually planted on only 6.5% (18,417 ha) of agricultural land in the WCS ( Table 2). Fields that are in potato commonly return between one through four years during the study period, with less than 1% returning seven out of the twelve years. Alfalfa was planted at least once on 62,772 hectares, or 22.2% of agricultural land, and has an annual average of 23,457 ha (or 8.2% of agricultural land) ( Table 2). Alfalfa is found throughout all five ecological landscapes and has the most hectares planted in the Central Sand Hills (27,660 ha). Frequency values most commonly range from one to six years with 1.4% in continuous alfalfa. While soybean is not shown in Figure 2, it is most often grown in rotation with corn and was planted at least once on 101,717 ha (35.7% of agricultural land), with an annual average of 21,935 ha (or 7.7% of agricultural land) ( Table 2). Soybean frequencies are most common between one and four years (out of twelve years), with 34.5% grown only once during the study period. Table 2. The dominant crops found in the Central Sands and their respective mean hectares and percent coverage of agricultural land. Standard deviation indicates the extent of deviation from the mean each rotation has; the CV is the coefficient of variation of annual average variability across all years; the slope, R 2 , and p-value are from linear regression and indicate the change in hectares planted during the study period, how well this change represents the data, and whether this change is significant. A concentrated area of agricultural land is found within 10 km to either side of the Central Sand Plains and the Central Sand Hills ecological boundaries. Ninety-four percent of potato plantings are found within Central Sand Plains and Central Sand Hills, with 73.1% of that occurring within this boundary. Sweet corn, beans, and peas also have a strong presence within this boundary, accounting for 75.6%, 75.3%, and 73.3% of plantings, respectively. Corn, soybean, alfalfa, and cereal are also planted but account for 47.1%, 40.2%, 36.3%, and 33.3% of plantings, respectively.

Common Two-and Three-Year Rotations
Two-and three-year rotations were determined for the eight crops and the other crops category and were compared using their annual average land area totals over the study period. Fourteen dominant two-year rotations were identified to occur on 130,801 hectares each year, or 45.9% of WCS agricultural land (Table 3). Continuous corn was the most common two-year rotation, which occurred on 13.2% of agricultural land (Table 3). Other dominant two-year rotations included alfalfa-alfalfa (5.6%), soy-corn (5.3%), corn-soy (5.0%), and sweet corn-potato (2.8%) ( Table 3). Continuous corn reached a max imum extent of 51,514 ha in 2012 and then steadily decreased thereafter (Figure 3). The two-year rotation of continuous alfalfa declined at the same time continuous corn and corn-soy were increasing, suggesting that corn and soy were replacing alfalfa in some of these rotations (Figure 3). At the time that continuous corn was declining after 2012, corn-soy and soycorn exhibited a steady increase in planted area (Figure 3). The increases in corn-soy and soy-corn were significant (p < 0.05) across the entire study period (Table 3). There was also a significant increase in potato-corn rotations (p = 0.019) and a significant decline in potato-sweet corn (p = 0.034) during the study period. The remaining two-year rotations remained relatively consistent during the study period with no significant trends. Table 3. Most common two-year rotations and their respective average hectare totals and percent coverage of agricultural land in the WCS. Standard deviation indicates the extent of deviation from the mean each rotation has; the CV is the coefficient of variation of annual average variability across all years; the slope, R 2 , and p-value are from linear regression and indicate the change in hectares planted during the study period, how well this change represents the data, and whether this change is significant. The following abbreviations correspond to the crop grown; C = corn, A = alfalfa, S = soybean, Sc = sweet corn, Po = potato, B = beans.

Rotations
Mean We identified sixteen three-year rotations that were found on 98,481 ha, or 31.4% of agricultural land in the WCS (Table 4). Continuous corn was the most common three-year rotation (8.5%), followed by continuous alfalfa (4.0%) and various combinations of corn and soy including corn-soy-corn (4.0%), corn-corn-soy (2.3%), soy-corn-soy (2.3%), and soycorn-corn (1.8%). Continuous corn reached a peak in 2011 of 31,298 ha and then exhibited a long-term decreasing trend through 2017 (Figure 4). Corn-soy-corn and soy-corn-soy rotations exhibited a significant increasing trend through the period (p < 0.05; Table 4). Two of the three-year rotations including potato-beans-corn and potato-corn-beans exhibited significant increases in area during the study period (p < 0.05; Table 4). None of the other dominant three-year rotations showed significant trends from 2008 to 2019 (Table 4).  (Table 4).

Return Time to Potato and Most Common Rotations
The analysis of multiyear rotations that began and ended with potato show th average 5.6% (1025 ha) of potato returned to potato immediately, 16.6% (3050 ha) retu the second year after initial planting, 43.7% (8072 ha) returned in the third year, and (2109 ha) returned after four years ( Figure 5). This suggests that the most popular y a potato rotation to come back to potato is year three (two years after the initial plan A total of 64.8% and 77.4% of fields planted with potato have returned to that crop b fourth and fifth years in the rotation, respectively (three-and four-year gaps, respect ( Figure 5). The most common crops that are rotated with potato include sweet corn, and beans. Sweet corn most commonly is grown between potato plantings, accountin 47.5% of all two-year rotations, followed by corn (17.5%), and beans (15.8%) (Tab Rotations with two years of sweet corn grown in between potato crops account for of all three-year rotations, and combinations of beans and sweet corn comprise 21.5% ble 5). The cropping sequences that are most frequently found in four-year potato tions include sweet corn-beans-sweet corn (16.2%), beans-sweet corn-beans (7.0% corn-corn-beans (6.0%) ( Table 5). Only one potato rotation (potato-corn-beans-p

Return Time to Potato and Most Common Rotations
The analysis of multiyear rotations that began and ended with potato show that on average 5.6% (1025 ha) of potato returned to potato immediately, 16.6% (3050 ha) returned the second year after initial planting, 43.7% (8072 ha) returned in the third year, and 11.6% (2109 ha) returned after four years ( Figure 5). This suggests that the most popular year in a potato rotation to come back to potato is year three (two years after the initial planting). A total of 64.8% and 77.4% of fields planted with potato have returned to that crop by the fourth and fifth years in the rotation, respectively (three-and four-year gaps, respectively) ( Figure 5). The most common crops that are rotated with potato include sweet corn, corn, and beans. Sweet corn most commonly is grown between potato plantings, accounting for 47.5% of all two-year rotations, followed by corn (17.5%), and beans (15.8%) ( Table 5). Rotations with two years of sweet corn grown in between potato crops account for 12.7% of all three-year rotations, and combinations of beans and sweet corn comprise 21.5% (Table 5). The cropping sequences that are most frequently found in four-year potato rotations include sweet corn-beans-sweet corn (16.2%), beans-sweet corn-beans (7.0%), and corn-corn-beans (6.0%) ( Table 5). Only one potato rotation (potato-corn-beans-potato) exhibited a significant trend (increasing area; p = 0.047) over the study period (Table 5).

Figure 5.
Bar graph represents the average hectares of potato that return after one year (continuous potato), two years (one crop between potato plantings), three years (two crops between), and four years (three crops between), and their respective standard deviations marked by error bars. The line graph represents the cumulative percent of the initial potato crop that has returned to potato. Table 5. Common rotations when rotating with potato after two-(one year off potato), three-(two years off potato), and four-years (three years off potato), mean number of hectares that rotation would occur, and the percent weight compared to all other two-, three-, and four-year potato rotations. Standard deviation indicates the extent of deviation from the mean each rotation has; the CV is the coefficient of variation of annual average variability across all years; the slope, r2 and p-value are from linear regression and indicates the change in hectares planted during the study period, how well this change represents the data, and whether this change is significant. The following abbreviations correspond to the crop grown; C = corn, Sc = sweet corn, Po = potato, B = beans.

Rotation
Mean  Figure 5. Bar graph represents the average hectares of potato that return after one year (continuous potato), two years (one crop between potato plantings), three years (two crops between), and four years (three crops between), and their respective standard deviations marked by error bars. The line graph represents the cumulative percent of the initial potato crop that has returned to potato. Table 5. Common rotations when rotating with potato after two-(one year off potato), three-(two years off potato), and four-years (three years off potato), mean number of hectares that rotation would occur, and the percent weight compared to all other two-, three-, and four-year potato rotations. Standard deviation indicates the extent of deviation from the mean each rotation has; the CV is the coefficient of variation of annual average variability across all years; the slope, r2 and p-value are from linear regression and indicates the change in hectares planted during the study period, how well this change represents the data, and whether this change is significant. The following abbreviations correspond to the crop grown; C = corn, Sc = sweet corn, Po = potato, B = beans.

Spatial Distribution and Frequency of Common Crops Grown in the WCS Is Dictated by Soil Formation Processes
The five ecological landscapes that are located within the WCS (Central Sand Hills, Central Sand Plains, Forest Transition, Central Lake Michigan Coastal, and Southeast Glacial Plains) have landforms left behind by glacial activity that influence spatial distribution of agricultural production. Over 24,000 years ago, the Green Bay Lobe ice sheet entered the WCS from the east and reached its maximum extent between the Central Sand Hills and Central Sand Plains ecological landscapes, creating series of moraines and leaving glacial outwash behind as it advanced and retreated before fully retreating from the area over 12,000 years ago [7]. The outermost moraine is known as the Johnstown Moraine and is exemplified by low hills that run north and south along the border of these two ecological landscapes and into Forest Transition [7]. About 19,000 years ago, the Green Bay Lobe advanced into and blocked the ancient river now known as the Wisconsin River, creating Glacial Lake Wisconsin within the Central Sand Plains [37]. Layers of silt were deposited on the lakebed and then covered with sand, leading to poor drainage and wetlands throughout this region when the glacial lake was drained. The eastern part of the Wisconsin Sand Plains near the Johnstown Moraine has outwash sand deposits that are thicker and higher above the water table, allowing for soils to be drained enough to grow vegetable crops with the use of irrigation [37]. The spatial distribution east of the Johnstown Moraine is influenced by the location of moraines and glacial outwash left between them after the glacier fully retreated. Deposits range from poorly to excessively well drained, making agricultural production dependent on irrigation and drainage [7,[37][38][39][40].

A High Number of Growers Use Multiple Year Rotations with Potato, but a Significant Fraction Return to Potato Immediately or by the Third Year
Potato is a high-value crop and a major component to the WCS's economy, contributing USD 522 million to Wisconsin's economy annually [41], thereby creating incentives for farmers to keep it in rotation. However, potatoes require large quantities of water and nutrients in a region that is susceptible to leaching or fertilizer, which has contributed to the region's groundwater quality issues [42]. Finding a balance between economic return and environmental sustainability is needed to keep the region profitable and reduce impacts on land and water resources. Implementing a diverse crop rotation for highvalue potato can be beneficial in replenishing soil nutrients and decrease the likelihood of disease and weed persistence [4], but also can be economically beneficial. According to Halloran et al. [43], there are three sources of potential economic benefits to crop rotations: the rotation either increases yields of the target crop (potato) at current or reduced input levels, the rotation crop itself is profitable, and the economic risk is reduced by income variability. In Halloran et al. [43], barley, sweet corn, green bean, soybean, canola, and continuous potato were simulated as two-year rotation sequences to estimate the expected mean net income and income variability. Potato-potato, green bean-potato, and sweet corn-potato rotations were found to be most profitable due to high-valued commodities. However, economic risk associated with potato-potato is high due to a loss in diversity, whereas the economic risk with green bean-potato is low, although a decline in potato yield was found the year following green bean from an increase in powdery mildew. In another study by Watkins et al. [44], they used the EPIC model to study two-, three-, and four-year rotations and their economic impact and environmental outcomes for seed potatoes in Idaho. A three-year rotation of potatoes followed by two grain crops of either spring wheat, barley, oats, and canola (potato-grain crop-grain crop) was found to have reduced nitrogen loss and soil erosion and obtained the most economic return.
These studies generally reflect recommended best practices for potato to spend two or more years out of rotation before returning to potato and to include snap beans, peas, sweet corn, field corn, small grains, and alfalfa [10]. In the Wisconsin Central Sands, these rotation recommendations with potato were generally observed, with 44% of fields returning to potato after three years off, and that sweet corn, corn, and beans are the most frequently rotated crops. Beans, potato, and sweet corn come from different plant families and differ in term of disease and insect host, making this combination an environmentally sustainable rotation [10]. Although most fields are implementing rotations, there were still a significant fraction of fields that returned to potato immediately in the following year or after only one year of an alternate crop planting. Rowberry et al. [45] compared profitability between a seven-year rotation of barley, corn, soybean, and potato to a rotation of continuous potato and found the continuous potatoes to be more profitable with a slight decline in yield. Some producers in the WCS region might be sacrificing yield, reduced soil health, and increased input costs to continually produce a higher-valued crop each year. Furthermore, given the high amounts of nitrogen fertilizer required for potato in the WCS, this may be exacerbating the ongoing water quality issues. However, the most common rotation of potato in the WCS is where two years of other crops are planted in between potato crops.

Changes in Continuous Corn Area and Increase in Soybean Reflect Changes in Energy Policy and Market Prices
The overall trend of declining continuous corn area after 2012 coupled with increases in soybean in rotation are reflective of changes in markets and prices and a push towards more sustainable management practices in the Midwest U.S. that lead to improved soil and water health. The desire for cleaner fuels and reduced greenhouse gas emissions have led to changes in policy and standards for biofuels [46], which has impacted land use change in agricultural landscapes such as the Wisconsin Central Sands. In the Energy Policy Act of 2005, the Renewable Fuel Standard (RFS) was established, requiring fuel to contain specific amounts of ethanol and the 2008 Farm Bill promoted the processing and production of this ethanol through subsidies [47,48]. In the years following, corn and soybean saw an intensification in acres planted in the United States coming from a reduction in pasture, CRP enrollment, and cropland expansion on land typically deemed not suitable for row crop production [35]. From 2009 to 2013 in the WCS, continuous corn area sharply increased until it reached a peak in 2012 before decreasing thereafter. Countering this rise, hectares planted in cereals and beans declined during these years. The rise and fall of the amount of continuous corn grown also reflects a pattern of changing corn prices. In August 2012, a drought impacted much of the corn belt and caused corn prices to reach a maximum of USD 8.24/bushel [49]. At the beginning of 2013, corn prices were forecasted to be within USD 6.95-8.25/bushel, likely causing the spike in corn hectares planted within the WCS [50]. During the 2013 growing season, the price of corn started to drop and was likely a major factor in the decline of corn planted in years following.
Common rotations with corn in the WCS region include soybean and alfalfa, with significant increases in those rotations that include corn and soybean during the study period. Soybean is a desired crop due to its high economic value and its ability to fix nitrogen from the atmosphere through determinate nodules in its root system which can be used by the plant as well as replenishing the soil nitrogen content after the plant dies and decomposes [51]. Adding soybean to a rotation can increase crop yield and enhance soil fertility, while reducing the need for additional sources of nitrogen [52][53][54]. These are particularly important for a region such as the WCS with its sandy soils and generally lower organic matter content. Alfalfa is commonly left in rotation for several years, allowing roots to penetrate deeper, improving soil structure and increasing organic matter. Rotating alfalfa with corn or other small grains has been found to be more profitable for farmers. Younger alfalfa stands (three years and younger) provide a higher yield than older stands and provide a boost in corn yield the year following alfalfa by 10% compared to corn following corn [55].

Conclusions
The WCS is a very diverse and important agricultural region that is faced with challenges of growing high-value crops on soils with low nutrient and water holding capacities. By implementing best management practices through more precise fertilizer applications, deficit irrigation, and crop rotations, the region can address some of the current environmental challenges around water resources facing the region. For this study, we looked at the spatial distribution of the dominant crops, common rotation combinations, and the amount of land area in each of these, the typical time between consecutive potato plantings, and agricultural land use trends during the 2008-2019 time period.
The spatial distribution and frequency of crops grown within the region is partially influenced by the ecological landscapes and soil formations within them left from the last glaciation from 12,000 years ago. The Johnstown Moraine, which roughly divides the WCS north to south, marks the furthest extent of the Green Bay Lobe which carved moraines and deposited glacial outwash as it advanced and retreated in the area [7]. Outwash deposits left between moraines range from poorly to excessively well drained and depth to water table influences the spatial distribution of agricultural land in the WCS. High frequencies of plantings are observed out to twelve years for potato, corn, alfalfa, and beans, but due to the diverseness of the region, continuous cropping is not commonly observed.
Corn is the most planted crop in the region and covers on average 24.8% of agricultural land each year. We defined thirteen two-year and sixteen three-year rotation combinations that cover 45.9% and 31.4% of agricultural land each year. These rotations include corn, soybean, sweet corn, potato, alfalfa, and beans. Continuous corn was the most dominant rotation followed by continuous alfalfa and corn-soy combinations for both two-and three-year rotation combinations. The amount of agricultural land in the WCS planted in continuous corn has declined since 2012, whereas there has been a significant increase in soybean planted in rotation with corn since 2008.
A majority of producers are following recommended management practices to reduce disease and weed pressure by introducing two or more annual crops into potato rotations; 16.6% of fields growing potato return to that crop after two years, with an additional 43.7% of fields returning to potato after three years off. On average each year, 5.6% of fields with potato are planted with that crop continuously. We suspect a key reason for growers to not exercise longer 4-and 5-year rotations for potato is the high value of the crop; this likely contributes to a significant fraction of producers to keep potato in continuous rotation (or return after only one year off), but this comes with the increased risk of disease and lower yields on average.
Given the N fertilizer requirements of potato and ongoing water quality issues in the WCS, where a significant fraction of private wells are greater than the EPA safe drinking water standard threshold of 10ppm [56], rotating crops with potato that require zero or less inorganic N fertilizer, and extending all potato rotations to at least 4 or 5 years could be a management option that helps to reduce nitrate leaching to groundwater in the long term.