Simulated Biomass, Climate Change Impacts, and Nitrogen Management to Achieve Switchgrass Biofuel Production at Diverse Sites in U.S

Switchgrass (Panicum virgatum L.) is a C4, warm season, perennial native grass that has been strongly recommended as an ideal biofuel feedstock. Accurate forecasting of switchgrass yield across a geographically diverse region and under future climate conditions is essential for determining realistic future ethanol production from switchgrass. This study compiled a switchgrass database through reviewing the existing literature from field trials across the U.S. Using observed switchgrass data, a process-based model (ALMANAC) was developed. The ALMANAC simulation results showed that crop management had more effect on yield than location. The ALMANAC model consists of functional relationships that provide a better understanding of interactions among plant physiological processes and environmental factors (water, soil, climate, and nutrients) giving realistic predictions in different climate conditions. This model was used to quantify the impacts of climate change on switchgrass yields. Simulated lowland switchgrass would have more yield increases between Illinois and Ohio in future (2021–2050) under both Representative Concentration Pathway (RCP) 4.5 and 8.5 pathways with low N fertilizer inputs than high N fertilizer inputs. There was no significant effect of climate variability on upland simulated yields, which means that N fertilization is a key factor in controlling upland switchgrass yields under future climate conditions.


Introduction
In 2015, 80% of the United States' total primary energy consumption was delivered from petroleum, natural gas, and coal, making the U.S. the second largest energy consumer in the world (after China) [1]. Over the last few decades, biofuel has been strongly recommended as an alternative renewable energy source to fossil-derived fuels [2,3]. Biofuels can play a critical role in reducing the emissions of greenhouse gases, especially in the transportation sector [4]. In the U.S., the transportation sector accounted for 27% of carbon dioxide emission in 2015, the second largest emission source after electricity (29%) [5]. The Energy Independence and Security Act (EISA) of 2007 mandated the use of 136 billion yield response to nitrogen (N) rates from 48 different switchgrass yield trials and reported that most findings of the field studies have consistently shown switchgrass yield increases from increasing N fertilizer application when averaged across locations and years. They also found switchgrass responded significantly to precipitation at zero or low N rates. These results revealed that switchgrass yield can vary within the same crop management strategies depending on climate and soil nutrients. Long-term field studies can give more accurate switchgrass yield estimation by properly accounting for the movement of N between soils and plants in different weather conditions. Despite the fact that the long-term field studies are more likely to reflect realistic switchgrass production estimates, such long-term switchgrass yield studies are uncommon because they require long-term commitment, funding, and dedicated study areas. Thus, plant simulation models have been developed to predict the long-term switchgrass yield and biomass production in response to different management systems and environments [28][29][30][31][32].
Simulation models of agricultural system are effective tools that can be used to estimate yield potential as a function of climate and soil conditions as well as crop management. Model simulations can be conducted with sufficient detail in space and time to quantify long-term changes in soil carbon, nutrients (N and P) cycling, water resources, and plant production across diverse environmental conditions. To use models of agriculture systems, calibration and subsequent validations should be performed [33]. To calibrate a model, high quality datasets (e.g., weather, soil, and cropping system) are required. A calibrated model should be validated with independent field data. The ALMANAC (Agricultural Land Management Alternatives with Numerical Assessment Criteria) model is a field-scale agricultural vegetation model [34]. This model simulates plant growth, water balance, erosion, soil organic carbon, and nutrients (N and P) cycling at a daily time step. Simulations are subjected to various management strategies [35]. ALMANAC predicts growth and biomass of switchgrass and has been validated against field data for multiple sites across the southern U.S. [36,37]. Future application of this model with switchgrass will benefit from its improved simulation of switchgrass production for biofuels in diverse environments. Behrman et al. [29] simulated switchgrass yields across the central and eastern U.S. with different climate change scenarios. However, in this study, a measured field data set collected in limited locations were used for model validation, which may increase the uncertainty in yield predictions over large areas. A more extensive study with more cultivars, years and additional location is needed to test how accurately the model simulates long-term yield dynamics of switchgrass across a wider range of soil types, temperature zones, and rainfall zones.
To improve model performance, large amounts of switchgrass productivity data from multiple regions including North Central, South Central, South East, North East regions in U.S. were collected and used to first calibrate and then validate ALMANAC's simulated yields in this study. This will increase confidence of the use of ALMANAC for yield forecasting with long-term changes in environmental factors (e.g., climate, water source, and soil) in various locations. Based on the collected yield data, the major factors influencing yield variability of switchgrass in U.S. will be analyzed using an additive regression model (ARM). The ARM applies functions fitted to data without considering the physiological processes involved in plant growth and morphogenesis and uses key variables that significantly affect yield of switchgrass as inputs. The ARM, based on the data collected over a wide geographic area under diverse management practices, was used to analyze past switchgrass growth data and to describe the yield responses of switchgrass lowland and upland cultivars to environmental covariates and management practices. The significant factor for switchgrass yield will be used for ALMANAC simulation analysis. The objectives of the study are the following: to develop realistic switchgrass yield estimates for diverse regions of the U.S.; to project future scenarios; and to explore the diversity of yield responses of two switchgrass ecotypes to climate change and different fertilizer applications. Each of these objectives was evaluated under three climate scenarios (base climate in 1976-2005 and two future climate scenarios using a representative concentration pathway (RCP 4.5 and RCP 8.5).

Data Collection
In this study, two major switchgrass ecotypes, "Upland" and "Lowland," were parameterized using extensive measured field data. The upland and lowland ecotypes were differentiated by morphological and genetic variations [38][39][40]. The lowland ecotype generally has fewer and larger tillers, erect growth, a compact crown, and blue-green waxy leaves, compared to upland ecotype [38,40]. Lowland ecotype is well adapted in wetter conditions, while upland ecotype is commonly found in drier soils [40]. Upland ecotype has higher cold tolerance than lowland ecotype [41,42]. Moreover, Zhang et al. [39] found that genetic variation in cytoplasm was highly correlated with the divergence of the ecotypes.
The experimental datasets used in this study were compiled from numerous peer-reviewed studies that reported switchgrass productivity in experiments conducted in stands ≥2 years old growing in plots (Table S1, Figure 1). Only biomass yields harvested once per year, either at the end of the growing season, after first killing frost, or in early spring, were used for the analysis. In this study, a total of 66 different locations from 32 peer-reviewed publications were analyzed and included 45 switchgrass cultivars (16 lowland cultivars and 29 upland cultivars). Fertilizer rates ranged from 0 to 896 kg N ha −1 yr −1 . The study locations, switchgrass ecotypes, cultivars, maximum and minimum dry yields of cultivars within each state, and references are listed in Table S1. Based on the database, yield patterns for each cultivar and each ecotype were analyzed. Weather data included total precipitation, minimum temperature, and maximum temperature for the study sites for the years of the study collected from National Oceanic and Atmospheric Administration (NOAA) [43]. Soil information including percentage of clay, silt, and sand in 100 cm soil depth was downloaded from the Natural Resources Conservation Service, United States Department of Agriculture, Soil Survey Geographic Database (USDA NRCS SSURGO database) [44]. Values of measured leaf area index (LAI, total one-sided area of leaf tissue per unit ground surface) for each switchgrass cultivar were collected from seven peer-reviewed studies conducted in latitudes from 31 • N to 45 • N.

Analysis of Factors Determining Switchgrass Yield
Switchgrass growth and development are affected by numerous factors, including temperature [45], rainfall [46], fertilizer inputs [47], and field location (e.g., soil type and latitude) [37,48]. To determine which variable explains most of the variation in switchgrass yields, in first phase, correlations and the importance of 14 variables for yield in switchgrass were analyzed. The

Analysis of Factors Determining Switchgrass Yield
Switchgrass growth and development are affected by numerous factors, including temperature [45], rainfall [46], fertilizer inputs [47], and field location (e.g., soil type and latitude) [37,48]. To determine which variable explains most of the variation in switchgrass yields, in first phase, correlations and the importance of 14 variables for yield in switchgrass were analyzed. The correlation scores for the 14 variables were switchgrass ecotype (−0.09); latitude (−0.40); longitude (0.37); soil type (0.06); components of sand (−0.05), silt (−0.05), and clay (0.18); Growing Degree Days (GDD) (0.29); total precipitation from April to September (0.36); average precipitation (0.36); minimum temperature (0.37); maximum temperature (0.22); average temperature (0.31); and N fertilizer rate (0.35). GDD is the accumulated growing degree days or heat units throughout growing season (April 1st-July 15th) at each study location. Variables with absolute correlation scores less than 0.3 have been removed so that only 8 variables were chosen at the first phase. At the second phase, since a high correlation between predictor variables decreases the precision of the estimated regression coefficient, the highly correlated predictor variables were removed from the model. To be more specific, minimum temperature is correlated with other selected factors, such as average temperature (0.96) and maximum temperature (0.86). Among the three factors, the minimum temperature that has the highest correlation score with the switchgrass crop yield is selected. Similarly, total precipitation is selected and its correlated factor (i.e., average precipitation (0.99)) was removed.
Five predictor variables, including N fertilizer rate (X 1 ), latitude (X 2 ), longitude (X 3 ), minimum temperature (X 4 ), and total precipitation (X 5 ) from April to September, were used to establish the switchgrass yield estimation additive regression model (ARM). The values of these variables were standardized prior to regression (i.e., data point subtracted from the mean and divided by the standard deviation of the distribution).
ARM is expressed as follows [49]: where f j X j = β j1 X j + β j2 X 2 j + · · · + β jL X L j and α = 5 j=1 β j0 . In Equation (1), Y is the response based on the five predictors variables and β j1 is the weight for the jth predictor with lth order, where l = 1, 2, . . . , L. The weight in each term represents significance of each variable on determining a value of Y (switchgrass yield). Error term ε follows the normal distribution with mean zero and variance 5 j=1 σ 2 j . Thus, the ARM can be illustrated as the multivariate normal distribution shown below: In Equation (2), X j = X 1 j , X 2 j , . . . , X L j and β j = β j1 , β j2 , . . . , β jL . Once the maximum order L is given, the maximum likelihood approach can be used to find the weight β jl , where l = 1, . . . , L.
However, since L is unknown, the structure learning approach with the 10-fold cross validation (CV) is utilized. In other words, the CV test will find a structure of f j X j in Equation (1) with the minimum LOOCV error [50].
In Equation (

ALMANAC Model Simulation of Development
The ALMANAC model is a field-based process-level simulator that accurately simulates growth, development, and yield of plants as influenced by environmental conditions and management. In simulation, to determine the effects of environmental condition on yields, switchgrass plants were produced only in rainfed areas. These are specified for the model as input data. Thus, to improve the stability and accuracy of plant growth simulation, high-quality input data for weather, soil, and management are required [51]. For calibrating ALMANAC, three research studies were selected from the total of 45 studies. These studies measured initial soil characteristics (e.g., nutrients, organic matter, N content, etc.) in the fields. The detailed soil information, location, and information on previous land use are listed in Table S2. Management for each study is below. No N, P and K fertilizer were applied in 2007, but 45 kg N ha −1 was applied as urea (N-P-K, 46-0-0) in 2008. Switchgrass was harvested in November in both years. More detailed information is available in Propheter et al. [53].
Study 3: The field study was conducted in 2008-2010 at Ardmore, Oklahoma and Waurika, Oklahoma. A one-year old lowland switchgrass cultivar (Alamo) was transplanted in 2008. Each experiment was a randomized complete block design with a split-plot arrangement and four replicates. N fertilizer rates (0 and 135 kg N ha −1 ) were whole plots, while P fertilizer rates (0, 15, 30, and 45 kg P ha −1 ) were sub-plots. Switchgrass was harvested after first killing frost in 2008-2010 at both locations. More detailed information is available in Kering et al. [54].
Weather data used for each study location were from the nearest and most complete National Oceanic Atmospheric Administration (NOAA) station [43]. Soil data for each study location were obtained from the USDA NRCS SSURGO database [44]. Based on the published values of initial soil data (Table S2), field characteristics were initialized at the start of the run for each site.
In the ALMANAC Plant database, two separate sets of plant parameters named Lowland and Upland were created for switchgrass (Table S3). Since upland cultivars consistently flower earlier than lowland cultivars [55,56], the values of degree days from green up to senescence for upland cultivars are less than for lowland cultivars [57]. The values of degree days varied with latitude [57]. Thus, potential heat units (PHU), degree days from green up to senescence, were adjusted for each of cultivar and location. The degree-days were set to zero when new tiller growth begins to emerge from the base of the plants [29]. The PHU value for each state and ecotype was calculated as accumulated heat units; that involves averaging the maximum and minimum local temperature for each day and subtracting the baseline temperature of 12 • C [57]. Upland cultivar had 1200 PHU for the three locations in OH. Lowland cultivar had 2000 PHU for two locations in KS and 2500 PHU for two locations in OK.
The value of the fraction of seasonal PHU at anthesis (DLAI) was 0.7 for both upland and lowland cultivars. In simulation, after switchgrass approached its maximum yield, the leaf area index started to decrease until harvest. We assumed that the measured yields were reported as maximum yields within the growing season for the study. Thus, in simulation, the potential maximum yields remained unchanged until harvest. The maximum potential leaf area index (DMLA) parameter for lowland and upland was firstly assigned by state based on leaf area index taken from previous studies [58][59][60][61][62][63]. Lowland cultivars have consistently greater values of leaf area index (LAI) than upland cultivars in southern part of U.S. The highest values of LAI were obtained in both southern and northern part of U.S. Unlike lowland cultivars, the lowest LAI for upland was obtained in southern part of U.S. Then, the DMLA parameters for lowland and upland were determined by comparing the ALMANAC simulation results with measured yields. The value for DMLA for lowland was determined to be 12 for all simulated locations. The DMLA for upland cultivars in 26 • -32 • N latitude range was determined to be 2.5, while the LAI in 33 • -48 • N latitude range was determined to be 6.0.
Based on the measured plant height values [55,56], the values of maximum crop height (HMX), for lowland and upland cultivars, were determined. Additionally, the values of HMX for lowland and upland cultivars were 2.5 and 1.3 m, respectively. Other switchgrass parameters were derived from previous studies [36,37].
For model validation, the calibrated ALMANAC model simulated switchgrass biomass across 59 additional sites from previous studies ( Figure 1). For these other sites, lowland and upland cultivars had PHU 3000 and 2500 in 30 • -32 • N latitude range, respectively; PHU 2500 and 2000 in 32 • -37 • N latitude range, respectively; PHU 2000 and 1500 in 37 • -38.5 • N latitude range, respectively; PHU 1800 and 1500 in 38.5 • -39 • N latitude range, respectively; PHU 1600 and 1200 in 39 • -41 • N latitude range, respectively; PHU 1400 and 1100 in 42 • -43.9 • N latitude range, respectively. Upland cultivars had PHU 900 in 44-48 • N latitude range. The same managements (e.g., planting, fertilizer, and harvesting) were used in the model as reported for each field study. Due to lack of information on soil characteristics in the publications, soil initial nutrient contents (N and P) were calibrated for most sites using published biomass data.
Simulated biomass values were compared with measured values that came from field measurements of Studies 1-3. The measured yields that were averaged across cultivars within each ecotype from each study location were also compared with simulated yields of lowland and upland switchgrass cultivars across other 59 study sites. The root-mean square error (RMSE, Equation (4)), RMSE-observations standard deviation ratio (RSR, Equation (5)), Nash-Sutcliffe efficiency (NSE, Equation (6)), and percent bias (PBIAS, Equation (7)) were calculated to quantify the model adequacy.
where i is the ith observation, n is the total number of observations, S i is the ith simulated value, O i is ith observed value, and O mean is the mean of observed data. Additionally, Pearson's correlation coefficient of determination (R 2 ) was estimated using Proc REG in Statistical Analysis Software version 9.3 (SAS 9.3, SAS Institute Inc., Cary, North Carolina, U.S.).

Climate Change Projection
For the future climate projection, APCC Integrated Modeling Solution (AIMS) was used to downscale climate change scenario information using eight global circulation models (GCMs) from CMIP5 (Coupled Model Intercomparison Project Phases) as provided in Table 1. Two core representative concentration pathways (RCP) scenarios including RCP4.5 and RCP8.5 were analyzed. The RCPs represents possible changes in future greenhouse gas (GHG) concentrations. Concentrations in RCP 4.5 peak around 2040 and decline, while concentrations in RCP 8.5 continue to increase throughout the 21st century. As a result, the RCP4.5 scenario projects global temperature increases of 1.1-2.6 • C in the future, while the RCP8.5 scenario leads to a temperature increase of 2-3 • C by the end of the 21st century [30]. The selected GCMs were bias-corrected and downscaled for historical  and the future periods (2021-2050), by using simple quantile mapping methods based on the reproducibility of minimum temperature and precipitation-amount-related extreme climate indices for the past [72]. The observed past weather data containing daily minimum and maximum temperatures and precipitation were obtained from National Oceanic and Atmospheric Administration (NOAA) [43]. Observed data for the weather station closest to each switchgrass field study site were collected as input ( Figure 1).
To quantify the impact of climate change on the switchgrass yield in two different N managements, including N50 and N150, the outputs of climate variables from the selected eight GCMs were fed into the calibrated ALMANAC model, and the switchgrass yield was projected in future period (2021-2050).
Here, the period of 1976-2005 is selected to represent the reference climate conditions. The period of 2021-2050 is selected to represent future climate conditions. According to Li et al. [73], expansion of ethanol production in the U.S. has minimal effect on changes in crop land use. Thus, it was assumed that there will be no changes in land use for the study areas in the future. Climate variables in the future, such as wind speed, relative humidity, and sunshine hours, were also assumed to be the same as the baseline period.

Switchgrass Yield Analysis
Lowland switchgrass cultivars produced consistently higher yields than upland cultivars ( Figure 2). The mean dry yields of lowland and upland switchgrass across locations, management scenarios, cultivars, and years were 14.0 Mg ha −1 and 7.4 Mg ha −1 , respectively. Switchgrass yields varied widely among locations depending on establishment time, management, soil, and climatic factors. Among the states, the highest yields of both groups of switchgrass cultivars were in Alabama, while the lowest yields were in Arkansas (Table S1). Switchgrass biomass had been fertilized with 84 kg ha −1 yr −1 N for 10 years in Alabama [74], while switchgrass had been fertilized with 134 kg ha −1 yr −1 N only for 3 years in Arkansas [41]. Among lowland switchgrass cultivars, Alamo had the highest mean value of 15.9 Mg ha −1 ( Figure S1) and produced as high as 39.3 Mg ha −1 with higher N fertilizer inputs at Texas [75]. When rainfall was low in 1996, among lowland switchgrass cultivars, Kanlow produced the lowest yield of 1.9 Mg ha −1 yr −1 with 134 kg ha −1 yr −1 N fertilizer inputs at Dallas in Texas. NL931 produced the lowest mean yield of 9.15 Mg ha −1 yr −1 among other lowland cultivars ( Figure S1).
Among upland switchgrass cultivars, cave-in-rock can be produced as high as 26.1 Mg ha −1 yr −1 with 200 kg N fertilizer inputs in 2009 at Western in Ohio [52]. However, based on 108 cave-in-rock data values, the mean yield value for this cultivar was 8.97 Mg ha −1 yr −1 . Carthage cultivar produced the lowest mean yield of 2.5 Mg ha −1 yr −1 with no nitrogen input between 2009 and 2011 at Hobet in West Virginia [76]. Among upland cultivars that are commercially available, shelter and Caddo cultivars produced the highest mean yields near 9.7 Mg ha −1 yr −1 ( Figure S1).
while the lowest yields were in Arkansas (Table S1). Switchgrass biomass had been fertilized with 84 kg ha −1 yr −1 N for 10 years in Alabama [74], while switchgrass had been fertilized with 134 kg ha −1 yr −1 N only for 3 years in Arkansas [41]. Among lowland switchgrass cultivars, Alamo had the highest mean value of 15.9 Mg ha −1 ( Figure S1) and produced as high as 39.3 Mg ha −1 with higher N fertilizer inputs at Texas [75]. When rainfall was low in 1996, among lowland switchgrass cultivars, Kanlow produced the lowest yield of 1.9 Mg ha −1 yr −1 with 134 kg ha −1 yr −1 N fertilizer inputs at Dallas in Texas. NL931 produced the lowest mean yield of 9.15 Mg ha −1 yr −1 among other lowland cultivars ( Figure  S1). Among upland switchgrass cultivars, cave-in-rock can be produced as high as 26.1 Mg ha −1 yr −1 with 200 kg N fertilizer inputs in 2009 at Western in Ohio [52]. However, based on 108 cave-in-rock data values, the mean yield value for this cultivar was 8.97 Mg ha −1 yr −1 . Carthage cultivar produced the lowest mean yield of 2.5 Mg ha −1 yr −1 with no nitrogen input between 2009 and 2011 at Hobet in West Virginia [76]. Among upland cultivars that are commercially available, shelter and Caddo cultivars produced the highest mean yields near 9.7 Mg ha −1 yr −1 ( Figure S1).

Analysis of Factors Influencing Yield Variability of Lowland and Upland Ecotypes
Using the latitude, longitude, climatic, and management dataset collected from total 66 study sites, ARM has been developed separately for lowland and upland cultivars. The ARM is: For lowland cultivars,

Analysis of Factors Influencing Yield Variability of Lowland and Upland Ecotypes
Using the latitude, longitude, climatic, and management dataset collected from total 66 study sites, ARM has been developed separately for lowland and upland cultivars. The ARM is: For lowland cultivars, For upland cultivars, In Equations (8) and (9), an absolute coefficient in each term represents the significance of each variable on the potential yield estimation of switchgrass. In other words, a larger value of the absolute coefficient refers to higher correlation between a predictor variable and a response variable. In Equations (8) and (9), because longitude X_ 3 has the highest total coefficients of 1.658 (= 0.827 + 0.424 + 0.406) and 43.614 (= 2.136 + 16.102 + 25.376) for lowland and upland cultivars, respectively, we can intuitively understand that it is the most critical variable to determine the potential yield estimation of switchgrass.
The resulting high correlation between switchgrass biomass and longitude is interesting and perplexing. Casler et al. [77], with latitudes ranging from 94 • W to 101 • W, also found significant differences due to longitude of their evaluation locations. In the study [77], few switchgrass cultivars (e.g., Pathfinder, Blackwell, CIR, WS98-IP) had significant linear relation to longitude of the study locations because their heading dates were changed by longitude. However, most of previous studies have reported a weak relationship among switchgrass yield and longitude [48,[78][79][80]. This may be because most of the previous studies were carried out between 95 • and 100 • W. Longitudes in the present study ranged from 77.97 to 102.8 • W. Looking at mean solar radiation (NOAA, 1977 NOAA, National Oceanic and Atmospheric Administration, 1993. Climatic Atlas of the United States, Asheville, North Carolina), there are no clear trends with longitude in this region in the months of April through August. However, there are definite east-west gradients in relative humidity for the months of June, July, and August. Likewise, there is an obvious east-west gradient in mean annual pan evaporation. It seems that switchgrass as a species favors the lower humidity, higher evaporative demand environments of the more westerly locations.
According to the Equation (8), significant factors on the yields of lowland cultivars are as follows: longitude (1.658), latitude (0.659), N fertilizer rate (0.576), minimum temperature (0.452), and total precipitation (0.006). The total precipitation had the least effect on the yields of lowland cultivars. This result is supported by Abdullahi et al. [81] who reported that the dry yield of upland cultivar was significantly decreased by decreasing irrigation frequency, while the irrigation regimes had a negligible impact effect on dry yield of lowland cultivar. Moreover, Sanderson and Reed [82] reported that the growth of switchgrass lowland cultivar in Central Texas did not respond to additional water. In Equation (9), the significance of predictor variables for the yields of upland cultivars is as follows: longitude (43.614), minimum temperature (0.868), total precipitation (0.685), latitude (0.463), and N fertilizer rate (0.031). Unlike lowland cultivars, N fertilizer rate had the least effect on the yields of upland cultivars, which means that upland yields are more influenced by regional adaptation, rather than N rate. According to Fike et al. [83], upland cultivars had lower yields in poorly drained sites and did not respond to fertilizer rate under low rainfall condition. However, upland cultivars produced higher yields and responded noticeably to N fertilizer rate from moderately well drained soil and under higher rainfall condition.
To evaluate the performance of ARM simulation, the ARM simulated yields were compared with the measured yields. In overall, ARM simulated switchgrass yields agreed moderately well with the measured yields across locations while showing the model failed to identify catchments with higher and lower yields for both ecotypes (Figure 3). Regression analysis for ARM simulated and measured yields of lowland and upland cultivars including all data from total 66 study sites revealed R 2 = 0.30 and 0.37, respectively ( Figure 3). Since ARM was developed based on the few key parameters, R 2 was relatively low. However, the simulated yields and measured yields for both lowland and upland cultivars were not significantly different (Both p ≥ 0.78) (Figure 3).

ALMANAC Evaluation
The ALMANAC model was calibrated against measured data for biomass yields from Study 1-3 ( Table 1). The Studies 1 and 3 were carried out on upland and lowland cultivars, applying various N fertilizer treatments, at multiple locations in Ohio and Oklahoma for 3-5 years ( Table 2). In both studies, measured switchgrass yields increased as N fertilizer input increased.

ALMANAC Evaluation
The ALMANAC model was calibrated against measured data for biomass yields from Study 1-3 ( Table 1). The Studies 1 and 3 were carried out on upland and lowland cultivars, applying various N fertilizer treatments, at multiple locations in Ohio and Oklahoma for 3-5 years ( Table 2). In both studies, measured switchgrass yields increased as N fertilizer input increased.
The ALMANAC simulated yields, which responded to different N fertilizer rate applications, were compared with measured yields using R 2 , RMSE, RSR, NSE, and PBIAS ( Table 3). As shown, the percent bias values for upland and lowland cultivars were under 10%. Based on the general performance rating of Moriasi et al. [84], the simulated yields of upland and lowland cultivars in Ohio, Kansas, and Oklahoma may be evaluated as "very good." However, based on values of RSR (0.97) and NSE (0.05), upland simulated yields can be evaluated as "unsatisfactory." The model underestimated upland switchgrass yields at high soil organic carbon and high total N. For example, high soil organic carbon and high total N were observed in the Northwest site, OH site (Table S2). In this site, the simulated yield differences were not as great as those observed yields from 50 to 200 N kg ha −1 ( Table 2). At high soil organic carbon, in simulation, the optimal plant growing condition (no N stress) was likely to appear with less N application. This may cause failure of model to adequately simulate the yield changes at upper range of N input. Unlike upland simulations, lowland simulations can be estimated as "good," based on values of RSR and NSE ( Table 3). The RMSEs between measured and simulated yields for upland and lowland cultivars were 3.52 and 1.44 Mg ha −1 , respectively.  After initial calibration of the model using detailed soil characteristics from seven different sites, the developed switchgrass model simulated yields in 59 other sites to determine how well it can predict yield under various climatic conditions. Both upland and lowland simulations can be evaluated as "very good" because PBIAS for both ecotypes were under 10% (Table 3). However, based on values of RSR and NSE, both upland and lowland simulations may be evaluated as "satisfactory" according to Moriasi et al. [84]. As the model validation results, regression analysis for simulated and measured yields of lowland and upland cultivars, including all data from total 59 study sites, revealed R 2 = 0.65 ( Figure 4 and Table 3). The upland ecotype yield simulation had a lower RMSE (2.36 Mg ha −1 ) than the lowland simulation (RMSE = 3.74 Mg ha −1 ). After initial calibration of the model using detailed soil characteristics from seven different sites, the developed switchgrass model simulated yields in 59 other sites to determine how well it can predict yield under various climatic conditions. Both upland and lowland simulations can be evaluated as "very good" because PBIAS for both ecotypes were under 10% (Table 3). However, based on values of RSR and NSE, both upland and lowland simulations may be evaluated as "satisfactory" according to Moriasi et al. [84]. As the model validation results, regression analysis for simulated and measured yields of lowland and upland cultivars, including all data from total 59 study sites, revealed R 2 = 0.65 ( Figure 4 and Table 3). The upland ecotype yield simulation had a lower RMSE (2.36 Mg ha −1 ) than the lowland simulation (RMSE = 3.74 Mg ha −1 ).

Switchgrass Yields in the Reference Period
The switchgrass yields in two different nitrogen treatments (N50 and N150) for each study area were estimated for the reference period  using the validated ALMANAC model, as shown in Figure 5. According to ARM regression analysis, the switchgrass yields were significantly related to longitude in study sites (Equations (8) and (9)). The ALMANAC simulated yields were highly related to longitude when lowland and upland cultivars were treated by N50 and N150, respectively ( Figure 5). The yields of lowland and upland cultivars appeared to decline from 80 to 100° W under N50 and N150 treatments, respectively ( Figure 5). Many of the observations were made between 90 and 100 °W, so both high and low yield values were observed in that regions. The simulated yields of lowland cultivars with higher nitrogen treatment were not related to the longitude. The lowland cultivars tended to have average yields of 11.72 Mg ha −1 consistently across locations. The yields of

Switchgrass Yields in the Reference Period
The switchgrass yields in two different nitrogen treatments (N50 and N150) for each study area were estimated for the reference period  using the validated ALMANAC model, as shown in Figure 5. According to ARM regression analysis, the switchgrass yields were significantly related to longitude in study sites (Equations (8) and (9)). The ALMANAC simulated yields were highly related to longitude when lowland and upland cultivars were treated by N50 and N150, respectively ( Figure 5). The yields of lowland and upland cultivars appeared to decline from 80 to 100 • W under N50 and N150 treatments, respectively ( Figure 5). Many of the observations were made between 90 and 100 • W, so both high and low yield values were observed in that regions. The simulated yields of lowland cultivars with higher nitrogen treatment were not related to the longitude. The lowland cultivars tended to have average yields of 11.72 Mg ha −1 consistently across locations. The yields of upland cultivars with low nitrogen treatment were not related to the longitude; they appeared to remain consistent around 5.53 Mg ha −1 across locations. According to ALMANAC simulation results, switchgrass yield's response to longitude varied by crop management, which means that the management may be a more important determinant yield factor than longitude.

Switchgrass Yields in the Future Climate
Comparing historical-observed  and future period (2021-2050) total precipitation (April-September), maximum temperature(Tmax), and minimum temperature (Tmin) values across the study areas, on average, the total precipitation, Tmax, and Tmin values increase from current period to future period in both RCP scenarios ( Figure S2). Maximum temperature values in future period increased by 2.22 and 2.52 • C on average across study locations, years, and GCMs under the RCP 4.5 and 8.5 pathways, respectively. The minimum temperature values in future period increased by 2.26 and 2.67 • C on average across study locations, years, and GCMs under the RCP 4.5 and 8.5 pathways, respectively. Unlike temperature, precipitation changes varied greatly across GCMs within each study area ( Figure S3). In the future period, according to CESM1-CAM5 model, most parts of study areas will experience a significant increase in precipitation under both RCP4.5 and 8.5 pathways, while the HadGEM2-CC model expected that precipitation will mostly decrease across locations ( Figure S3). Agronomy 2020, 10, x FOR PEER REVIEW 10 of 20 upland cultivars with low nitrogen treatment were not related to the longitude; they appeared to remain consistent around 5.53 Mg ha −1 across locations. According to ALMANAC simulation results, switchgrass yield's response to longitude varied by crop management, which means that the management may be a more important determinant yield factor than longitude.

Switchgrass Yields in the Future Climate
Comparing historical-observed  and future period (2021-2050) total precipitation (April-September), maximum temperature(Tmax), and minimum temperature (Tmin) values across the study areas, on average, the total precipitation, Tmax, and Tmin values increase from current period to future period in both RCP scenarios ( Figure S2). Maximum temperature values in future period increased by 2.22 and 2.52 °C on average across study locations, years, and GCMs under the RCP 4.5 and 8.5 pathways, respectively. The minimum temperature values in future period increased by 2.26 and 2.67 °C on average across study locations, years, and GCMs under the RCP 4.5 and 8.5 pathways, respectively. Unlike temperature, precipitation changes varied greatly across GCMs within each study area ( Figure S3). In the future period, according to CESM1-CAM5 model, most parts of study areas will experience a significant increase in precipitation under both RCP4.5 and 8.5 pathways, while the HadGEM2-CC model expected that precipitation will mostly decrease across locations ( Figure S3).
Switchgrass is C4 warm season grass, and therefore an increase in temperature can, in theory, increase its yield in the absence of drought limitations. However, the results show that two switchgrass ecotypes show different yield patterns in the future period ( Figure 6). In the future period, lowland switchgrass yields tended to increase in the areas located between longitude −85 and −80 degrees under high and low nitrogen treatments in both RCP 4.5 and 8.5 pathways. In these study areas, increased temperature may increase the duration of switchgrass growth (according to maximum heat unit approach), which may increase vegetative development (node and leaf appearance rate). However, lowland switchgrass yield variation across GCMs within each study area was becoming greater in RCP 8.5 pathway. As shown in Figure S3, RCP 8.5 pathway showed higher means of the daily maximum and minimum temperatures than RCP 4.5 scenarios. These warmer temperatures may limit switchgrass production increases achieved through the high nitrogen application across all study regions. Under high nitrogen treatment, for example, the lowland Switchgrass is C 4 warm season grass, and therefore an increase in temperature can, in theory, increase its yield in the absence of drought limitations. However, the results show that two switchgrass ecotypes show different yield patterns in the future period ( Figure 6). In the future period, lowland switchgrass yields tended to increase in the areas located between longitude −85 and −80 degrees under high and low nitrogen treatments in both RCP 4.5 and 8.5 pathways. In these study areas, increased temperature may increase the duration of switchgrass growth (according to maximum heat unit approach), which may increase vegetative development (node and leaf appearance rate). However, lowland switchgrass yield variation across GCMs within each study area was becoming greater in RCP 8.5 pathway. As shown in Figure S3, RCP 8.5 pathway showed higher means of the daily maximum and minimum temperatures than RCP 4.5 scenarios. These warmer temperatures may limit switchgrass production increases achieved through the high nitrogen application across all study regions. Under high nitrogen treatment, for example, the lowland switchgrass yield differences between the least predicted yield and the greatest predicted yield among eight GCMs within each study area were between 1.7 and 8.7 Mg ha −1 . Yield simulations based on MIROC5 produced nearly half the yields predicted based on other GCMs in areas located between longitude −101 to −93 degrees. According to MIROC5 model, maximum temperatures will mostly increase about 3.4 • C, which is higher than temperature predicted by other models ( Figure S3). The total precipitation will vary between study areas, but several areas will experience a decrease, or little change in precipitation ( Figure S3). The increase in temperature can cause faster heat unit accumulation; that can lead to shorter duration of growth and increased heat stress.
In comparison with lowland cultivars, yield changes were barely observed for upland cultivars under both high and low nitrogen treatments under RCP 4.5 and 8.5 pathways. Yield variation across GCMs within each study area was much smaller than lowland cultivars. However, yield patterns for upland cultivars vary under different nitrogen treatments. Under high nitrogen treatment, small yield variations for upland cultivars across GCM were observed between longitude −102 and −97 degrees, while yield variations were observed between −89 and −83 degrees under low nitrogen treatments. This may be related to the potential heat units. As temperature increases under RCP 4.5 and 8.5 pathways, in simulation, the maximum potential heat units of upland ecotypes may be reached earlier, which prevents their growth. Due to this early maturity, upland ecotypes may not respond to the fertilizer application in most study locations.

Conclusions
This study analyzed switchgrass productivity data from multiple sites across the U.S. Each observation has been identified by cultivar, ecotype, stand age, harvest system, N application rate, planting location, soil, and year of harvest. The additive regression model, ARM, was developed based on few key parameters (e.g., climate, fertilizer, and location) and used to evaluate relationship among the key parameters and yield through analyzing past growth data. According to ARM

Conclusions
This study analyzed switchgrass productivity data from multiple sites across the U.S. Each observation has been identified by cultivar, ecotype, stand age, harvest system, N application rate, planting location, soil, and year of harvest. The additive regression model, ARM, was developed based on few key parameters (e.g., climate, fertilizer, and location) and used to evaluate relationship among the key parameters and yield through analyzing past growth data. According to ARM simulation results, the yields of switchgrass lowland upland cultivars were highly related to longitude of the study location. This result was used for analysis of ALMANAC climate change simulation results. The collected yield data was used to develop a process-based model and to improve switchgrass yield simulation accuracy. ALMANAC model simulates key physiological processes linking weather to plant yield production. Many of the essential parameters in the model have been established through experiments or literature reviews. As a result, the ALMANAC model can predict switchgrass yield relatively well for both switchgrass cultivars. For climate change simulations, the validated ALMANAC was used to quantify the climate change impacts on switchgrass yields. Based on the simulation results in reference period (1976−2005), overall, switchgrass yield response to longitude varied according to nitrogen application rates. For example, yields of switchgrass lowland cultivar responded well to variation in longitude with low nutrient application, while upland cultivars appeared to suffer from a nutrient deficiency across locations. Under high nutrient application rates, upland cultivars appeared to relate to longitude and had greater dry yield at the more eastern locations. The validated ALMANAC model was used to quantify the impact of climate change on switchgrass yield with two different nitrogen fertilizer applications. In the future period (2021−2050), both RCP4.5 and 8.5 of eight climate models showed an increase of temperature across study locations. Changes in precipitation differed across GCMs and locations. According to CESM1-CAM5 model, most areas will experience an increase in precipitation, while the HadGEM2-CC model tends to have a decrease in precipitation. ALMANAC simulated yields of lowland cultivars are projected to increase in areas located in longitude between −88 and −80 degrees under high and low nitrogen treatments in both RCP 4.5 and 8.5 pathways. Under RCP 8.5 pathway, the simulation results for lowland cultivars show substantial variation across GCMs. One of the GCMs, MIROC5, predicted nearly half size of yields that were simulated based on other GCMs across each study area located between longitude −101 to −93 degrees. Unlike lowland cultivars, climate variability barely affected production of upland cultivars. On the basis of the overall simulation results, yields of lowland and upland switchgrass increased with increasing nitrogen application. However, lowland and upland ecotypes responded differently to climate variability. While upland cultivars appeared to be showing no responses from climate change, the warmer temperatures may limit productions of lowland cultivars even though high nitrogen was applied to the plots. This paper shows that the fertilizer management may be the most important factor for upland cultivars, which will control switchgrass yield under increasing temperature. This finding will be useful for overcoming the impacts of climate change on switchgrass yield and make switchgrass more adapted and perform better in a wide range of locations in the U.S.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2073-4395/10/4/503/s1, Table S1: Study locations (states), switchgrass ecotype (L: Lowland, U: Upland), cultivars, minimum and maximum biomass dry yield (Mg ha −1 ), and references. Table S2. Summary of soil characteristics (0−30 cm depth) measured for switchgrass production studies conducted in Ohio, Kansas, and Oklahoma, USA. SOC indicates soil organic carbon in soil; TN indicates total nitrogen in soil. Table S3. Key ALMANAC plant parameters for switchgrass lowland and upland ecotypes. Figure S1: Biomass yields for lowland and upland cultivars of switchgrass. Grey bar indicates the mean yield value for each cultivar. Upper bar indicates standard deviation for each cultivar. Number above cap of bar indicates the number of observations for each cultivar. Figure  Funding: This work was also supported in part by an appointment to the Agricultural Research Service administered by the Oak Ridge Institute for Science and Education through interagency agreement between the U.S. Department of Energy (DOE) and the U.S. Department of Agriculture (USDA), Agricultural Research Service Agreement #60-3098-5-002.