Rye-soybean double-crop: planting method and N fertilization effects in the North Central US

Abstract Double-cropping winter rye cover crops (CC) with soybean in the North Central US could help with the global effort to sustainably intensify agriculture. Studies addressing the management of these systems are limited. Therefore, a field study was conducted from 2017 to 2019 in Central Iowa, US to evaluate winter rye CC biomass production, aboveground N accumulation, estimated economics, estimated within-field energy balance and estimated greenhouse gas (GHG) emissions under three N application rates (0, 60, 120 kg N ha−1) and three planting methods (pre- and post-harvest broadcast and post-harvest drilling). Averaged over N rates, all planting methods resulted in >5.0 Mg ha−1 year−1 rye aboveground biomass dry matter. Averaged over the 2-year study and compared with unfertilized treatments, applying 60 kg N ha−1 produced 1.1 Mg ha−1 more aboveground biomass (6.1 vs 5.0 Mg ha−1), accumulated 30 kg ha−1 more N in aboveground biomass (88 vs 58 kg N ha−1), and led to 20 GJ ha−1 more net energy. Biomass production was not significantly higher with 120 kg N ha−1 compared with the 60 kg N ha−1 rate. Even when accounting for an estimated 0.75 Mg ha−1 of above ground rye biomass left in the field after harvesting, more N was removed than applied at the 60 kg N ha−1 rate. The minimum rye prices over the 2-year study needed for double-cropping winter rye CC to be profitable (breakeven prices) averaged $117 and $104 Mg−1 for the 0 and 60 kg N ha−1 rates, which factors in estimated soybean yield reductions in 2019 compared with local averages but not off-site transportation. GHG emissions were estimated to increase approximately threefold between the unfertilized and 60 kg N ha−1 rates without considering bioenergy offsets. While environmental tradeoffs need further study, results suggest harvesting fertilized rye CC biomass before planting soybean is a promising practice for the North Central US to maximize total crop and net energy production.


Introduction
Sustainable intensification of agriculture can be defined as feeding a growing population while minimizing the impact on the environment, which is high on the global policy agenda and one of the grand challenges facing society (Tilman et al., 2011;Garnett et al., 2013;Petersen and Snapp, 2015;NAS, 2021). We must improve our knowledge of management practices that reduce environmental impact while increasing overall production. Harvesting cover crops (CC) has shown promise toward retaining agroecosystem services of CC without harvest, such as erosion control, water quality improvement and weed suppression . Harvesting CC before planting the primary crop may increase overall production by creating a double-cropping system (i.e., two harvestable crops are grown in the same field in succession), which can increase harvestable biomass and provide a revenue stream as livestock forage (Ketterings et al., 2015;Lyons et al., 2019) or bioenergy feedstock (Baker and Griffis, 2009;Feyereisen et al., 2013;Shao et al., 2015). Similar systems, such as relay-cropping, have been discussed as resource-efficient technologies that contribute to improved environmental quality, increased net returns and greater overall crop production per unit area than traditional systems (Tanveer et al., 2017). Optimizing the management of these systems will be critical to maximizing productivity while maintaining agroecosystem services (Tanveer et al., 2017;Blanco-Canqui et al., 2020).
Several studies have shown the potential for fertilizer application to increase rye biomass CC production. Balkcom et al. (2018) reported that in the Southeastern US low rates of commercial N fertilizer resulted in increased rye biomass with similar total costs per Mg of rye compared with unfertilized treatments. Shao et al. (2015) showed that fertilizing and harvesting rye in the Northeastern US can increase biomass and producer revenue compared with unfertilized rye. However, fertilizer applications to CC must be used cautiously since one of the purposes of cover cropping in the drained North Central US is to reduce nutrient loss to the environment. For example, a major benefit of establishing winter rye (Cereal secale L.) CC in corn (Zea mays L.)soybean systems (Glycine max (L.) Merr.) is to scavenge residual soil nitrogen (N) and reduce nitrate levels in drainage water (Kaspar et al., 2007(Kaspar et al., , 2012NAS, 2021). Soil NO3-N can be reduced by harvesting winter rye compared with non-harvested rye (Krueger et al., 2011), and the reduction results in less available N for leaching . A modeling study by Malone et al. (2018) showed that doublecropping fertilized rye with soybean reduced N loss to subsurface drainage compared with unharvested and unfertilized rye while also increasing net revenue and harvested biomass.
Rye CC is typically planted using either a grain drill or surface broadcast seeder, and the planting method can affect biomass production. Drilled rye often shows more rapid and more uniform establishment than broadcast seeding because of improved soil-seed contact and protection from erosion events and extreme fluctuations in temperature and water content (Fisher et al., 2011;Bich et al., 2014). While drilling must occur after harvest of the primary crop, broadcast seeding can be applied into a standing crop to increase the length of the fall growing season (Wilson et al., 2013). Alternatively, broadcast seeding after primary crop harvest can be accompanied by mechanical incorporation of the seed to improve establishment (Fisher et al., 2011;Brennan and Leap, 2014).  reported mixed results when broadcasting rye seeds into standing corn compared with post-harvest drilling. Haramoto (2019) compared broadcasting or drilling rye seeds on the same date and reported that drilling led to higher CC biomass under dry soil conditions. Continued research into planting methods is critical to guide farmers' decision-making across a range of operations and equipment availability to maximize success of systems that include winter rye.
A major barrier to widespread adoption of CC is cost. Singer et al. (2007) reported that 56% of farmers in the Northern US soybean region would plant CC if economic incentives were available. More recent studies have also reported that incentives would likely increase CC use (Roesch-McNally et al., 2018;Plastina et al., 2020). Blanco-Canqui et al. (2020) reported that double-crop CC planting may improve net returns, but more research on management scenarios is needed. While double-cropping soybean systems have been reported to increase production and producer revenue since the 1980s in the Southeastern US (Marra and Carlson, 1986), only recently have such systems been discussed as profitable as far north as Minnesota, US, and Ontario, Canada (Gesch and Archer, 2013;Gesch et al., 2014;Davidson, 2016). The two US states with the largest soybean production, Iowa and Illinois (National Agricultural Statistics Service, 2020;Shahbandeh, 2021), are located in the Northern US soybean region, and each had less than 3% of corn and soybean land in CC in 2015 (Rundquist and Carlson, 2017). Therefore, identifying methods to profitably implement a double-crop CC system within this region could help make the practice more attractive to a large segment of producers.
Finding a combination of management practices that contribute to sustainable intensification must balance productivity, profitability and environmental quality. From a bioenergy perspective, the energy balance of double-cropping systems is an important consideration. For example, Gesch et al. (2014) reported less net energy from double-cropping camelina-soybean systems than soybean-only systems. Similarly, management influences on greenhouse gas (GHG) emissions can be used to evaluate potential tradeoffs between productivity and sustainability goals Malone et al., 2019). Although doublecropping soybean systems show promise in the North Central US, few field studies have addressed agronomic management of these systems. Here, we evaluate how rye planting methods and N fertilizer rates affect rye biomass yields and N content in a winter rye-soybean double-cropping system; and consider potential costs, estimated energy balances and estimated GHG emissions.

Material and methods
We initiated a field study in 2017 on a site approximately 10 km west of Ames, IA in Boone County (42°00 ′ N; 93°47 ′ W), with Spillville loam (Fine-loamy, mixed, superactive, mesic Cumulic Hapludoll) as the predominant soil type. Field plots were laid out in a split-plot randomized complete block design with six replications (blocks). Each year main plots (3.8 × 30.5 m, n = 18) were assigned to three planting methods, and subplots (3.8 × 9.1 m, n = 54) to three N rates (Fig. 1). Two sets of plots (east and west) were managed under a 2-year rotation with one year in corn and the other year in rye-soybean. The rye was planted during corn growth or shortly after harvest (Table 1). In 2018, the west field was in corn and the east field was in rye-soybean. In 2019, the crops were reversed (Fig. 1).

Crop management and analysis
The corn and soybean were managed as a typical Central Iowa no-till system (Table 1). Weeds were controlled through a combination of preemergence and postemergence herbicides for both corn and soybean. Soybean (Pioneer P20T79R, 445,000 seeds ha −1 ) and corn (Viking A81-98R, 86,000 seeds ha −1 ) were planted with a five-row, 0.76 m row width, no-till planter. N fertilizer (32% urea ammonium nitrate, 32-0-0) was applied to corn in each growing season at planting (34 kg ha −1 ) and side-dressed (170 kg ha −1 ) along the row at V4.
Winter rye 'Elbon' was planted during the corn phase of the rotation using three methods: (1) drilled after corn harvest (row width = 19.0 cm; depth = 2.5 cm), (2) overseeded, i.e., broadcast over corn with a spinner spreader mounted on a high-boy tractor at the R6 growth stage, and (3) incorporated, i.e., broadcast and shallow incorporated with a rolling stalk chopper after corn harvest. Rye seeding rates for the three planting methods were 247, 371 and 309 pure live seed (PLS) m −2 ; or 55, 82 and 69 kg total seeds ha −1 using 51,809 total seeds kg −1 and 85% germination. Similarly, seed supplier Green Cover Seed lists a seed weight of 50,700 'Elbon' seeds kg −1 (Green Cover Seed, 2021). Higher rye seeding rates are recommended in Iowa for broadcast compared with drilled seeding. In early April, we fertilized subplots within each planting method with three rates of surface applied granular urea (46-0-0): (1) 0 kg N ha −1 , (2) 60 kg N ha −1 and (3) 120 kg N ha −1 (Fig. 1). Early spring can be the optimum time to apply N to winter crops in the North Central US (Franzen, 2018;Malone et al., 2018).
At the flowering stage, we hand-clipped plants to the soil surface within a 0.38 m 2 rectangular frame. Biomass samples were dried at 60°C until weights were stable, and weights were scaled up based on frame size to Mg ha −1 dry matter (DM). The collected rye biomass was finely ground and analyzed for C and N

446
Robert W. Malone et al.
content using the dry combustion method (FlashSmart C and N analyzer; Thermo Fisher Scientific, Waltham, MA). The remaining rye not collected for analysis was green-chopped and collected as ensilage. Soybean seed yields were determined by harvesting each main plot area, weighing the grain and measuring grain moisture. Although a control treatment without winter rye was not implemented, we compared soybean yields with Boone County, Iowa averages from National Agricultural Statistics Service (2021). We calculated N recovery efficiency (RE N ) by the difference method, where RE N is the difference of N in rye biomass between fertilized and unfertilized treatments divided by the amount of fertilizer N added (Cassman et al., 2002).

Economic analysis
We estimated biomass production costs associated with the different treatments, including establishment, fertilization and harvest costs. We estimated unfertilized rye establishment costs at $92-$110 ha −1 for the different treatments, assuming seed costs of $0.67 kg −1 and $55 ha −1 for custom no-till planting (Plastina and Johanns, 2021). Costs related to rye planting and labor in our analysis differed only according to the seeding rate, although costs for the different methods for on-farm implementation may differ based on field conditions, producer operations and equipment availability. Also, the planting costs used are higher than  those recently reported for drilling or broadcasting CC seeds (Plastina and Johanns, 2021), but are similar to the $99 ha −1 reported by Balkcom et al. (2018) and less than the $151 ha −1 reported by Roley et al. (2016). While Balkcom et al. (2018) and Roley et al. (2016) included rye termination costs, we did not include these costs because harvesting the rye at the flowering stage effectively terminated growth, thereby negating the need for additional management (Clark, 2007). However, CC termination cost using herbicide will increase the production costs, which Roley et al. (2016) estimated at $45 ha −1 . We assume winter rye was harvested as hay, which takes about 3 days under good conditions to cure in the field after mowing (Porter, 2017). In Iowa, the first alfalfa hay cutting can be mid-May (Barnhart, 2010), which is comparable to or earlier than our harvest dates (Table 1). Assuming mowing, raking and tedding at $63.4 ha −1 and baling and moving large round bales to storage at $22.4 Mg −1 DM harvested, total on-farm custom harvest costs not including nutrient replacement were $175.4 ha −1 for 5.0 Mg rye harvested (Plastina and Johanns, 2021). In comparison, if rye harvest costs were similar to harvesting corn stover for biofuel as Baker and Griffis (2009) assumed, estimated rye harvest costs when the producer owns all the necessary equipment (or custom rates) would be $101 ha −1 (or $179 ha −1 ) using the values reported by Edwards (2014) for corn stover harvest in Iowa without including hauling or nutrient replacement costs. Graham et al. (2007) reported $106 ha −1 for 3.5 Mg of shredded/raked/baled corn stover without hauling but including nutrient replacement.
We estimated dry fertilizer N application costs at $15 ha −1 (Plastina and Johanns, 2021), which is similar to Balkcom et al. (2018). Fertilizer prices were estimated at $0.88 kg-N −1 (Mensing, 2019), $0.66 kg-K 2 O −1 and $0.86 kg-P 2 O 5 −1 . We estimated P and K inputs based on 3.5 and 13.0 kg P 2 O 5 and K 2 O removed per Mg rye CC harvested (from FEAT Model described below). We did not include field application costs of P and K, assuming them to be unchanged with or without rye.
Overall, the total production costs for 5.0 Mg of unfertilized harvested rye biomass ranged from $267 to $285 ha −1 . Adding P, K and N fertilizer increased these costs to between $325 and $464 ha −1 depending on seeding and fertilizer rates. We calculated the cost of producing one Mg of biomass DM as the total production cost divided by harvestable biomass. While rye was clipped to the soil surface during collection, we estimated harvestable biomass by subtracting 0.75 Mg from total aboveground biomass. A 10 cm rye height can be assumed equal to 1 Mg rye aboveground biomass, and 7.5 cm height is acceptable to maintain soil ecosystem services when harvesting CC with high soil surface residue after no-till corn (Blanco-Canqui et al., 2020).
The costs above do not include a potential reduction in soybean yield resulting from the double-cropping system, which can occur when harvesting an overwintering crop in late May or early June (Gesch et al., 2014;Nafzinger et al., 2016). Breakeven prices for harvested rye were calculated as the price needed to generate net revenue equal to that without soybean yield reduction (Gesch et al., 2014): where Pr is the breakeven rye price ($ ha −1 ), Ps is the soybean price (assumed $400 Mg −1 DM), Ys0 is the countywide soybean yield (Mg DM ha −1 ) for the year (National Agricultural Statistics Service, 2021), Ysr is the soybean yield (Mg DM ha −1 ) with harvested rye treatments, Yri is the rye yield (Mg DM ha −1 ) for treatment i, and Cri is the rye production costs ($ ha −1 ) for treatment i. For rye prices below this breakeven level, growing soybean alone would be more profitable than the doublecrop system. This approach provides an indication of the minimum rye price that would be needed for a double-crop system to be more profitable than prevailing practices without considering off-site transportation costs.

Energy budgets and GHG with FEAT
We estimated energy balances and GHG emissions for each rye treatment based on inputs and outputs in the field operations using the Farm Energy Analysis Tool (FEAT) version 1.2.7 , unless otherwise noted. Operations associated with corn and soybean management were assumed to remain the same with or without rye. Therefore, included in the energy and GHG balance were only on-farm fuel use for seeding, fertilizing and harvesting the rye, and energy and emissions resulting from the production of seeds, N, P and K fertilizer based on N rate and replacement of P and K exported in crop biomass and transport of inputs. The energy associated with on-farm fuel use (44.8 MJ L −1 ) and transportation associated with inputs (0.64 GJ Mg −1 ) were estimated based on the FEAT database. We did not include energy inputs and emissions resulting from transportation after harvest. Energy inputs and outputs were as follows (GJ Mg −1 ): N 54.8; P 10.3; K 7.0; rye seeds 6.2; rye biomass 21.8; soybeans 23.8. Energy input associated with on-farm fuel use was estimated using 20 L ha −1 diesel for planting and harvesting rye, and 1.7 L ha −1 for N fertilizer application. Fuel use for P and K were not included, and this fuel use was assumed to be the same with or without rye in corn-soybean systems. Net energy of each treatment was calculated by subtracting the energy input from the output of the harvested rye and subtracting the reduced energy output with the estimated soybean yield reduction (Ys0 − Ysr from Equation 1).
IPCC emission factors for wet climates were used to estimate direct emissions of N 2 O from soils following N fertilization (1.6%) and N content in above and below-ground residues (0.6%) resulting from rye management (Hergoualc'h et al., 2019). Calculations in FEAT were followed except where field data for biomass N content and residue remaining after harvest could be input. Estimates of N 2 O emissions were converted to kg CO 2 -e (CO 2 -equivalent) ha −1 using the 100-year estimated warming potential with climate-carbon feedbacks of 298 kg CO 2 -e kg N 2 O −1 (Myhre et al., 2013).

Statistical analysis
Statistical analyses were performed using PROC MIXED in SAS (version 9.4, SAS Institute, Cary, North Carolina, USA) for splitplot designs to determine the impact of N fertilizer and planting method on rye biomass, N content, RE N , and energy and economic parameters. The fixed effects were N rate, planting method and their interactions; the random effects were replication (the blocking factor) and the whole-plot error (the interaction between replication and the whole-plot factorthe planting method). Pairwise comparisons were made using the Tukey-Kramer test when main effects were significant at P < 0.05. We separated data by growing season because of differences in weather

448
Robert W. Malone et al.
conditions and field plot locations. Data were also analyzed using PROC MIXED with year included as a random effect to determine whether treatments were significantly different regardless of year.

Growing conditions
Annual cumulative September through August rainfall during each of the two growing seasons was approximately 1000 mm, above the average of 860 mm (Fig. 2a). The 2017-2018 rainfall was similar to the 30-year average until June 2018, when the monthly rainfall was 296 mm compared with the average of 121 mm. In September 2018, the rainfall was 91 mm greater than the average, and the cumulative rainfall remained above average through August 2019. While rainfall in both years was above average, June through August 2019 received 314 mm less rainfall than 2018 and 110 mm less than the average. Cumulative rye growing degree days (GDD; base 4.4°C) for September through May were 1332°C days in 2017-2018 and 1119°C days in 2018-2019, while the average was 1247°C days (Fig. 2b). The warmer October 2017 and May 2018 compared with 2018-2019 were the largest monthly GDD differences between the two growing seasons.

Soybean yields
Overall soybean yields in 2018 (3.8 Mg ha −1 ) were similar to the county average of 3.7 Mg ha −1 , but the yields in 2019 were 0.6 Mg ha −1 lower than the county average of 3.6 Mg ha −1 (National Agricultural Statistics Service, 2021). Although precipitation and GDD were close to the 30-year averages (Fig. 2), the reduction in overall yield in 2019 may be partly a result of lower 2019 June through August precipitation. Also, in-field variability between the east and west plot areas may be related to the roughly 1 m higher field elevation in 2019 compared with 2018 ( Fig. 1). Kravchenko and Bullock (2000) reported that elevation influenced soybean yield, with higher yields consistently observed at lower landscape positions. Soybean yields were not determined in subplots of N fertilizer rates because of logistical constraints, but adding fertilizer N to rye prior to soybean was only expected to have a small effect, if any, on soybean yields (Salvagiotti et al., 2008;Mourtzinis et al., 2018). While high N fertilizer rates (>200 kg N ha −1 ) to soybean can increase soybean yields in Iowa and Nebraska (Bakhsh et al., 2009;Cafaro et al., 2020), only small N-related effects on soybean yields have generally been detected (Mourtzinis et al., 2018). Soybean yields are typically not reduced following unharvested rye CC (Acharya et al., 2020;, though corn yields are sometimes decreased (Kaspar and Bakker, 2015;Pantoja et al., 2015). However, further investigation into soybean yields in double-cropping systems is warranted, as the later rye termination from harvest and later soybean planting date may reduce yields by shortening the growing season (Egli and Cornelius, 2009;Hu and Wiatrak, 2012). Also, double-or relay-cropping soybean systems can use more water than singlecropping systems in the North Central US (Gesch and Johnson, 2015), and a main limiting factor in double-cropping soybean systems is limited soil water availability (Santos Hansel et al., 2019). Despite some studies in the Northern US soybean region reporting lower soybean yields when double-cropping, these systems have the potential for higher total crop yields (Gesch et al., 2014;Nafzinger et al., 2016;Johnson et al., 2017).

Rye biomass
Growing conditions for both 2018 and 2019 were favorable for rye biomass production, with precipitation higher than the 30-year average and GDD approximately average for both years (Fig. 2). Averaged across planting methods and N rates, Table 2 shows above ground biomass was slightly higher in 2018 (6.1 ± 1.3 Mg ha −1 ) than in 2019 (5.5 ± 1.6 Mg ha −1 ). This may be attributed to (1) more GDD total in September 2017 through May 2018 (Fig. 2), (2) the later rye harvest date in 2018 (Table 1), and (3) rye being N limited in 2019 but not in 2018 as discussed below. For the 60 kg ha −1 N rate, the rye biomass was similar between years and the difference was less compared with the unfertilized treatment (6.3 vs 5.9 and 5.9 vs 4.1 Mg ha −1 , Table 2). Overall, rye biomass production in this study was similar to ranges reported in Alabama (Balkcom et al., 2018) and Maryland (Mirsky et al., 2017), and exceeded production in Nebraska (Koehler-Cole and Elmore, 2020; Koehler-Cole et al., 2020) and Kentucky (Haramoto, 2019).
Above-average precipitation and consistent GDDs in both rye growing seasons likely reduced differences between planting methods. Fall establishment of CC by drilling is typically more uniform and rapid than broadcast (Fisher et al., 2011;Noland et al., 2018), but these advantages are most evident in dry conditions (Haramoto, 2019;. Thus, overseeded and drilled rye produced similar biomass in 2018, but incorporated rye produced significantly higher biomass (Table 2). Seed incorporation likely improved production compared with the overseeded treatment because improved seed-soil contact typically improves germination and emergence, leading to greater biomass production (Fisher et al., 2011;Brennan and Leap, 2014;Peterson et al., 2019). Any advantage of earlier planting and higher seed rate in the overseeded treatment was likely offset by only 36 mm of precipitation in September of 2017 compared with an average of 71 mm (Fig. 2). The higher-than-normal precipitation that facilitated establishment did not occur until October, when the other treatments had also been seeded. Additionally, the broadcast and incorporated rye likely outproduced the drilled treatment because of increased seeding rates, which can lead to more biomass (Wilson et al., 2019; Koehler-Cole and Elmore, 2020), though that is not always a consistent response (Haramoto, 2019). Conversely, September precipitation in the 2019 season led to the early establishment of the overseeded treatment (Wilson et al., 2013); the early establishment allowed a longer growing season that, coupled with the highest seeding rate, contributed to the overseeded treatment producing the most biomass (although not significantly higher than the other two planting methods, Table 2).
Averaged over both years, biomass DM values for the unfertilized and 60 kg N ha −1 plots were 5.0 and 6.1 Mg ha −1 and were significantly different (Table 2). However, in 2018, rye biomass production was not increased by N fertilizer, suggesting that N was not a limiting factor for biomass production (Table 2). Despite rye typically displaying a biomass response to N fertilizer in the Southeast to the Northeast US (Reiter et al., 2008;Mirsky et al., 2017;Balkcom et al., 2018), N demands in 2018 in our North Central US site were likely met by a combination of residual N from the preceding corn crop and soil N mineralization. The east plot area producing rye in 2018 was at the bottom of a slope (Fig. 1). Mamo et al. (2003) reported that field areas non-responsive to N were mostly at lower elevations, likely having increased N availability and mineralization.
Conversely, rye biomass did respond to N fertilization in 2019, suggesting that residual N and soil N mineralization did not meet crop N demand in the higher elevation west plot area (Fig. 1). While adding N fertilizer increased biomass production compared with the unfertilized treatments, the highest rate, 120 kg N ha −1 , did not produce significantly more than the 60 kg N ha −1 treatment (Table 2). This finding roughly agrees with Mirsky et al. (2017), who found that 72.4 kg N ha −1 was needed to reach maximum rye biomass production. Similarly, a study in Canada indicated that N fertilization above 42 kg N ha −1 for rye was not justified (Landry et al., 2019). Blanco-Canqui et al. (2020) reported that while more research is needed, the few studies available indicate that CC harvesting does not generally affect ecosystem services such as soil properties, erosion and weed suppression. In our study, even without fertilization, the rye produced enough biomass to provide agroecosystem services based on yields reported in other studies. Rye reduced nitrate concentrations in subsurface drainage with much lower biomass (<2 Mg ha −1 ) in Iowa corn-soybean systems (Kaspar et al., 2012). Rye biomass production of 3.5-4.0 Mg ha −1 effectively reduced runoff and soil erosion in Wisconsin monocrop corn silage production by more than 50% (Siller et al., 2016), while less than 1 Mg ha −1 of rye biomass reduced erosion rates in 2 out of 3 years in a no-till corn-soybean system (Kaspar et al., 2001). Similarly, CC biomass of about 4 Mg ha −1 suppressed weeds prior to main crop planting in a Pennsylvania corn-oat production system (Finney et al., 2016), although others report that 8 Mg ha −1 biomass or more is needed to reduce weed pressure on subsequent main crops (Ryan et al., 2011;Smith et al., 2011).

Rye N content and potential leaching
The amount of N in rye aboveground biomass differed among planting methods in 2018, but not in 2019 ( Table 2). The incorporated method had the highest amount of biomass N in 2018, reflecting higher total biomass production. In both years, N  (Table 2) and developmental rates (data not shown) among the three N rates, indicating that biomass N concentration was higher in the fertilized treatments. The higher rye biomass N concentration occurred despite little increase in overall biomass production, which is consistent with Lyons et al. (2019) that reported rye biomass N content can increase beyond the N rate that optimized yield. The N content in aboveground biomass and estimated N removed during rye harvest exceeded the additional N applied in the 60 kg ha −1 treatment in both years but not in the 120 kg ha −1 treatment (Table 2). RE N peaked at 63.8% in 2019 and reached only 38.5% in 2018 in the 60 kg N ha −1 treatment. In comparison, Balkcom et al. (2018) reported approximately 50% RE N for rye across three fertilizer rates of spring applications (34, 67 and 101 kg N ha −1 ). These findings suggest some reduction in a CC's ability to limit N loss in the fertilized system, but the fertilized rye double-cropping system may still lead to a reduction in soil N and N losses compared with a no-CC scenario. This finding is critical in Iowa because water quality deterioration caused by excessive nitrate levels is an ongoing concern (Dinnes et al., 2002;Hatfield et al., 2009;Syswerda et al., 2012), so any management practice that reduces the potential for nitrate leaching should be considered. In a modeling study, Malone et al. (2018) reported that fertilized and harvested rye in corn-soybean rotations reduced drainage N loss by 18% compared with rye CC that was neither fertilized nor harvested. In a nearby field with similar soil and other conditions, Kaspar et al. (2012) reported unfertilized and unharvested rye reduced N concentration in drainage water 48-58%. So, even without removing rye and the associated organic N, considerable N reduction in drainage water under similar conditions has been reported when including CC. These findings suggest that harvesting rye at an N application of 60 kg N ha −1 may reduce nitrate leaching and N loss to subsurface drainage and surface water compared with no CC or unharvested and unfertilized CC.
Nonetheless, spring N application must be carefully managed because most tile drainage occurs during spring months (Jaynes et al., 2008;Waring et al., 2020) and some fertilizer N loss would be carried out of tile drainage lines and discharged to surface waters. For example, the lower elevation plots in the current study with rye in 2018 did not require N fertilizer for additional biomass production, and management decisions such as higher or lower fertilizer rates can depend on field variability. Precision N management of crops within a field has been reported to contribute to sustainable food production without degrading the environment (Diacono et al., 2013;NAS, 2021).

Economics
In 2018, planting method and fertilizer rate affected cost required to produce an Mg of harvested biomass (Table 3). The biggest cost difference in 2018 was between N rates of 0 and 120 kg N ha −1 ($26.8 Mg −1 ), where the fertilizer costs were not offset by higher biomass. The overseeded method had the highest cost resulting from the highest seeding rate and low biomass. In 2019, neither fertilizer rate nor planting method significantly affected costs because the higher seed and fertilizer costs were offset by higher Table 3. Data are reported by rye planting method (drilled after corn harvest, broadcast and overseeded into R6 corn, and broadcast and incorporated after corn harvest) and by fertilizer N application rate (0, 60 and 120 kg N ha −1 ). Net energy and breakeven price estimates include soybean yield reduction in 2019. Note that breakeven prices in 2018 are equal to production costs because soybean yield was similar to the county average. Energy output, net energy, costs and prices were based on estimated harvestable biomass (see Table 2).
biomass. Averaged over both years, costs for the 0 and 60 kg N ha −1 rates were $77.4 and $79.7 Mg −1 (and were not significantly different). The costs per Mg of rye biomass are nearly twice the $44 Mg −1 reported by Balkcom et al. (2018) in a study with approximately half the rye biomass (∼2.2 Mg ha −1 ), partly because of including P and K replacement and other harvesting costs in the current study.
With no apparent effect on soybean yield in 2018, the breakeven price would be the rye production costs ($67.4 and $77.5 Mg −1 for 0 and 60 kg N ha −1 ). Although this study did not include a no-rye treatment, we estimated the effects of potential reduction in soybean yields on breakeven prices by comparing with the county average in 2019. The 0.6 Mg ha −1 soybean yield reduction in 2019 contributed to breakeven prices of $166 and $131 Mg −1 for N rates of 0 and 60 kg N ha −1 (Table 3). Thus, the breakeven price was lower for fertilized compared with unfertilized rye in 2019. Averaged over both years, breakeven prices for the unfertilized and 60 kg N ha −1 rates were $117 and $104 Mg −1 (and were not significantly different). Assuming fertilized and unfertilized rye have nearly the same potential price per Mg as reported by Shao et al. (2015), the slightly lower breakeven price and significantly higher yield for the fertilized rye across the current 2-year study indicate that fertilizer input would be neutral to economically viable in these conditions and help sustainably intensify agriculture.
The zero N rate rye production cost of $77.4 Mg −1 averaged over the 2-year study is similar to the estimated $75 Mg −1 for harvested and unfertilized rye biomass reported by Baker and Griffis (2009) and used by Malone et al. (2018). If the fertilizer cost increases from $0.88 to $1.33 kg-N −1 (Malone et al., 2018), the breakeven price for the 60 kg N ha −1 fertilized rye-soybean double-crop system increases from $104 to $110 Mg −1 over the 2-year study, which is still slightly lower than unfertilized rye. While N application in rye production after corn on high organic matter soils is recommended in the Northern US soybean region (Kaiser, 2018), Malone et al. (2019) reported only small reductions in soil organic C and N in Central Iowa long-term cornsoybean rotations with removal of C and N through corn stover harvest without N replacement. Thus, N application may not be required to prevent reduced soil organic matter. Further, economically optimum N fertilizer rates within a field can vary substantially over short distances and year-to-year (Mamo et al., 2003;Scharf et al., 2005), suggesting the lower elevation plots of the current study with rye harvested in 2018 require less N fertilizer on average over several years than the slightly higher elevation plots of rye harvested in 2019.
For comparison with our breakeven prices: (1) prices for good grass hay in Iowa were >$120 Mg −1 in 2021 (AMS, 2021;Hay and Forage Grower, 2021) and (2) value for rye as a bioenergy feedstock plus coproducts were >$150 Mg −1 (Shao et al., 2015) or ∼$400 Mg −1 (Herbstritt et al., 2022). Note, however, that bioenergy feedstock or silage involves different operations (and therefore costs) than hay. If the rye market price were below the breakeven prices calculated here or breakeven prices were higher (e.g., including off-site transportation or higher harvest/operation costs), positive value may still be present factoring in the ecosystem services associated with the harvested rye (Balkcom et al., 2018;Blanco-Canqui et al., 2020).
Winter rye has potential as a cellulosic bioenergy feedstock in the US (Baker and Griffis, 2009;Feyereisen et al., 2013;Shao et al., 2015;Malone et al., 2018), and government initiatives like the US Renewable Fuel Standard (RFS2) call for increasing the production of cellulosic biofuels from crops like winter rye (EPA, 2022). Agricultural biomass similar to winter rye is used for bioenergy production from anaerobic digestion in several countries outside the US (Valli et al., 2017;Theuerl et al., 2019;Kemausuor et al., 2018). Industrial scale renewable natural gas facilities using agricultural residues as their feedstock are operating in the North Central US (e.g., https://www.verbio.us/project/ verbio-nevada-biorefinery/), and interest is growing for these systems to produce biogas from agricultural biomass and convert it into renewable natural gas (Pleima, 2019). Herbstritt et al. (2022) reported the revenue potential of the winter rye collected in this current study as a feedstock for these digestor systems, which included using the digestate as a feed protein supplement. Since the market for winter rye as an energy feedstock is currently limited in the US, our economic analysis focused on breakeven prices rather than revenue and compared these breakeven prices to current market prices for hay and potential prices for the cellulosic bioenergy market.

Energy budgets and greenhouse gas emissions using FEAT
Energy inputs per hectare were the lowest for the unfertilized treatments and the highest for the 120 kg N ha −1 fertilizer treatments in both years (Table 3). In 2018, energy outputs were similar between the fertilizer rates (Table 3), while in 2019 the 60 and 120 kg N ha −1 treatments resulted in higher energy output than the unfertilized treatment. The average annual energy output of approximately 117 GJ ha −1 for the N rate of 60 kg ha −1 is similar to rye silage reported by Camargo et al. (2013) of slightly more than 100 GJ ha −1 . The calculated annual energy input (5.5 GJ ha −1 ) is lower than the Camargo et al. (2013) estimate of approximately 10 GJ ha −1 , partly because we did not include additional lime, herbicide, insecticide or tillage. The current calculated energy inputs were closer to the double-cropping fertilized camelina-soybean system of 7.1 MJ ha −1 reported by Gesch et al. (2014) when subtracting the soybean-related energy inputs.
The net energy associated with the rye management in the double-cropping system was positive in both years for all treatments (greater than 55 GJ ha −1 ) when accounting for the estimated 2019 reduced soybean yield. Averaged over both years, net energy values were 84 and 104 GJ ha −1 for the unfertilized and 60 kg N ha −1 rates (and were significantly different). In contrast, net energy of double-cropping soybean systems with winter camelina were reported 53 GJ ha −1 lower than soybean-only systems (Gesch et al., 2014). The main differences in the net energy of the two studies were the high soybean yield reductions of 2.0 Mg ha −1 and lower camelina yield of approximately 1.2 Mg ha −1 reported by Gesch et al. (2014). They attributed the high soybean yield reduction mostly to the late seeding dates (30 June 2010 and 11 July 2011) relative to the conventional mono-cropping soybean (5 May 2010 and 19 May 2011).
The FEAT model predicted that N fertilization would lead to increased GHG emissions, which ranged from a low of 277 kg CO 2 -e to a high of 1931 kg CO 2 -e ha −1 year −1 in 2019 (Fig. 3). Direct N 2 O emission resulting from fertilizer application was the largest source of GHG, totaling 450 and 899 kg CO 2 -e ha −1 year −1 in the 60 and 120 kg N ha −1 treatments (Fig. 3). This estimate depends heavily on the choice of emission factors and is the largest source of uncertainty in estimating GHG impacts of fertilizer management in the double-cropping system. The geographic location and fertilizer application timing in early spring on wet soils support the higher emission factor of 1.6% for fertilizer

452
Robert W. Malone et al. application in wet climates. However, this effect may be offset by lower soil temperatures in spring, when the growing rye is expected to take up most of the additional N. In balance, the 1.6% emission factor may overestimate emissions as statistical modeling studies of soybean relay-cropping systems suggest a much lower emission factor (Cecchin et al., 2021). Energy required to manufacture the N fertilizer (234 and 469 kg CO 2 -e) was the second greatest source of GHG. Increased biomass production resulting from N fertilization did not offset these emission sources, and GHG intensity (kg CO 2 -e GJ −1 output) increased with each fertilizer rate in both years. In 2018, GHG intensity was 3.1, 9.7 and 18.2 kg CO 2 -e GJ −1 output at 0, 60 and 120 kg N ha −1 , and in 2019 GHG intensity was 4.0, 10.5 and 15.8 kg CO 2 -e GJ −1 output at 0, 60 and 120 kg N ha −1 .
The FEAT model estimates N 2 O emissions resulting from N contained in plant residues returned to the field. In this case, residues included N contained in the rye stubble and the N contained in below-ground residues, which were estimated as a percentage of above-ground biomass based on default values in the FEAT model. Rye residues were predicted to contribute substantially to direct and indirect N 2 O emissions in all systems, averaging 143-206 kg CO2-e depending on biomass production. However, the IPCC emission factors likely overestimate rye residue contribution to N 2 O production. Residues with higher C:N ratio, such as rye, are likely to immobilize inorganic N during decomposition and unlikely to contribute to increased N 2 O emissions (Charles et al., 2017). In a 10-year study on a similar soil type, Parkin et al. (2016) found no increase in N 2 O emission associated with rye that produced less than 2.0 Mg ha −1 in corn-soybean rotations. Rather, rye led to a non-statistically significant decrease in N 2 O emissions. Moreover, rye decreased leaching by 192 kg ha −1 and correspondingly decreased estimated indirect N 2 O emissions by 1.44 kg N 2 O-N over the 10 years. Similarly, Machado et al. (2021) found that residue contribution to N 2 O emissions never exceeded 0.1% in diversified and conventional cropping systems.
Thus, our estimates likely overestimate the burden and underestimate the ecosystem services provided by the double-crop CC in terms of GHG emissions. However, yield benefits from fertilizer application in the rye phase may be evaluated against potential for increased GHG emissions from the system.

Conclusions
Over the 2-year study in the North Central US, we found that double-cropping soybean with winter rye CC showed promise to increase overall production per unit land area. We found that planting method did not consistently affect any metrics included in this study, suggesting that this double-cropping system can be successfully implemented in a variety of ways when conditions are favorable for CC production. Fertilizing the rye double-crop CC at 60 kg N ha −1 was economically competitive with no fertilizer application over the 2-year study, produced more biomass, had higher net energy and had more N removed with harvest than applied. However, increasing the fertilizer rate to 120 kg N ha −1 did not increase biomass production or N accumulation enough to offset increases in N rate. Further research into fertility management of rye within the double-cropping soybean system will improve the efficiency of these systems, as our findings indicated that landscape position and year may have played roles in biomass production response to N fertilizer. Additionally, further research into downstream markets like sustainable bioenergy may offset GHG emissions from N fertilizer and improve these systems overall carbon balance and environmental benefit, as reported by Herbstritt et al. (2022). Sustainable bioenergy systems can remove Fig. 3. Predicted (FEAT model) greenhouse gas emissions as CO 2 -equivalents from management of the rye-soybean double-cropping system under alternate rye planting methods and fertilizer N application rates. For each growing season, 2018 and 2019, data are reported by rye planting method (drilled after corn harvest, broadcast and overseeded into R6 corn, and broadcast and incorporated after corn harvest) and by fertilizer N application rate (0, 60 and 120 kg N ha −1 ). Average rye biomass yield and N content for each treatment from each year were used for this analysis.
carbon from the atmosphere through terrestrial sequestration and bioenergy carbon capture and storage (Sanchez et al., 2015;Field et al., 2020). Overall, including rye in a double-cropping soybean system to intensify agricultural production can be economically and environmentally viable in the North Central US.