Assessing Remote Sensing Vegetation Index Sensitivities for Tall Fescue ( Schedonorus arundinaceus ) Plant Health with Varying Endophyte and Fertilizer Types: A Case for Improving Poultry Manuresheds

: Tall fescue ( Schedonorus arundinaceus ) is a common perennial forage in cattle pastures of the southeastern United States. A mutualistic fungal endophyte normally infects the grass and produces ergot alkaloids toxic to livestock, but fungal biotypes that have no ergot alkaloid production have been developed. Here remote sensing methods were used to assess plant health in 1 ha grazed paddocks with application amongst different combinations of fertilizer sources (inorganic and broiler litter) and endophyte associations (wild, novel–tall fescue MaxQ type with novel endophyte, and free). Broiler litter fertilization is common in the region due to the presence of many chicken farms. Moreover, broiler litter costs are comparable to inorganic fertilizer depending on distance from source to application. Incorporating remote sensing, we tested the sensitivity of three indices: normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and land surface water index (LSWI) to assess fescue plant health. Indices were obtained from satellite imagery provided by Landsat 7 ETM+ between the years 2005 and 2009. Sensitivity analytics suggested that LSWI was the optimum index to determine fescue plant health. The year experiencing drought (determined by annual precipitation) showed signiﬁcant difference between fertilizer types ( p = 0.05) and a nearly signiﬁcant difference between endophyte associations ( p = 0.08). There was no signiﬁcant difference in years with normal or wet precipitation rates due to tall fescue endophyte association or type of fertilization. Limited availability of satellite imagery during parts of the ﬁve years of study might have inﬂuenced outcomes of statistical analyses. Nevertheless, the data and ﬁndings point to the potential use of satellite imagery in assessing grazingland tall fescue health and advancing the concept of poultry manureshed in the region or elsewhere where poultry manure production is extensive. Y.Z.; visualization, K.C.F., S.M.; supervision, K.C.F., D.E., A.F., Y.Z.; project K.C.F., T.L., D.E., A.F., S.M., Y.Z.;


Introduction
Agricultural systems in the United States have become more specialized in recent decades, but few systems have been developed to promote the recycling of surplus manure nutrients to nearby croplands [1]. There are many consequences when large amounts of nutrients accumulate in and around feeding operations, including detrimental impacts to soil, water, and air quality, as well as quality of life for livestock [2][3][4]. Historically, agricultural producers in the United States have used livestock manure nutrients as a source of crop fertilizer [5][6][7][8]. Regions that support a particular livestock type do not always overlap with regions that produce feed for that livestock, which adds to the difficulty in transporting excess nutrients in the form of fertilizers to agricultural producers [9]. The emerging concept of "manuresheds" expands the understanding and meaning of returning livestock manure nutrients to agricultural lands where they are needed, whether it be across fences, counties, or regions [1]. Utilizing our growing understanding of manuresheds will aid in preventing "sources" of excess manure that can critically harm the environment.
Those regions identified as manureshed sources can provide nutrients for regions identified as manureshed "sinks" (those lacking natural fertilizers to apply to agricultural fields due to a lack of surplus manure). Opportunities to help create a balance between these sources and sinks lie within geographical distribution of grazing lands, but successful application will depend on proper management to minimize hazardous environmental impact [10]. For instance, unlike typical current practice of using poultry litter as nutrient source for crops or forages within close proximities of production facilities, partly as a disposal mechanism. Instead it can be used to meet phosphorus (P) and nitrogen (N) demands of tall fescue (Schedonorus arundinaceus) paddocks located within the poultry manureshed region including source and sink areas in a manner that takes optimal nutrient need of the grass and environmental protection into account. However, the biomass production of the tall fescue under varying climatic conditions may be affected by whether the paddock was fertilized with poultry litter or synthetic nutrients [11]. Other studies with napiergrass have shown that there is not a significant difference in crop yield between inorganic fertilizer and poultry litter [12]. However, biomass N concentration and total N removal were greater with inorganically fertilized fields, while biomass P concentration and P removal were greater in fields fertilized by poultry litter [12]. Manure can be lightly applied to enhance grazing land soil properties, but heavy rates of disposal are considered unsafe [13].
Determining biomass production can be completed by in situ measurements. However, remote sensing techniques present opportunities to conduct this through analysis of satellite imagery. Satellites provide a multitude of imagery in different spectral bands providing opportunity to analyze vegetation from an aerial view. The spectral signature of tall fescue in combination with vegetation indices can supply spatial information on biomass, drought stress, and the nutritional status of crops [14]. Multiple studies have reported the use of the normalized difference vegetation index (NDVI) and other multi-spectral indices for assessing the performance of turfgrass similar to tall fescue [15,16]. NDVI has been shown effective in measuring leaf chlorophyll content, which can provide information concerning the physiological state of a plant [17] and can serve as a proxy for biomass measures [18]. NDVI presents a negative correlation with soil brightness and is influenced by the atmosphere, but the Enhanced Vegetation Index (EVI) simultaneously corrects for both [19], potentially resulting in greater sensitivities to grazingland physiological states. Moreover, light is greatly absorbed by water in the shortwave infrared (SWIR) band, which is used by the Land Surface Water Index (LSWI) to detect the amount of liquid water in vegetation and its soil background [20][21][22], again potentially serving as a proxy to grazingland physiological states.
Thus, the main objective of this paper is to compare sensitivities/performance of three common vegetation indices (NDVI, EVI, and LSWI), derived from the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) as a proxy for plant productivity and health of grazed tall fescue with different endophyte associations (endophyte free, novel endophyte, and wild endophyte) managed under differing fertilizer types (poultry litter vs. inorganic) and under various climatic conditions (dry, normal, and wet years) in northeast Georgia, USA. The second objective is to highlight how manuresheds can benefit from remote sensing technology in grazinglands research within reasonable proximity of poultry manure production facilities. Franzluebbers et al. [23] and Endale et al. [24,25] give details of the parent study that looked at pasture and cattle performance as well as runoff amount and quality in response to fertilization and tall fescue endophyte associations mentioned above. Our paper is focused on evaluating satellite derived indices to assess tall fescue plant health during this parent study.

Study Area and Experimental Design
Georgia is an exemplary location of poultry manureshed potential as it lies within a region (southeastern US) that contributes 60% of the United States annual 8.6 billion broiler production, resulting in nearly 10 million Mg y −1 of byproduct broiler litter [26]. The experimental site is a~20 ha parcel of land near Watkinsville in northeast Georgia (33 • 53 N, 83 • 26 W). At the time of the study, the site was owned and operated by the United States Department of Agriculture-Agricultural Research Service (USDA-ARS) J. Phil Campbell Sr. Natural Resource Conservation Center. According to data from the Soil Climate Analysis Network (SCAN) operated by the USDA Natural Resources Conservation Service (NRCS), the long-term (30-year) mean annual rainfall and temperature are 1261 mm (49.6 in) 120 and 16.5 • C (61.7 • F), respectively [27].
The site is comprised of 14 tall fescue paddocks, 0.9 to 1.1 ha each. All but two were grazed throughout the year; the remaining two were hayed. Each year in September, paddocks were stocked with weaned Angus (Bos taurus) heifers (2002)(2003)(2004)(2005)(2006)(2007) and steers (2008)(2009)) at rates of~3 to 6 head per paddock adjusted every 28 days to maintain 1 to 3 Mg ha −1 of available forage (Franzluebbers et al. 2009). Weather conditions influencing soil water and forage production dictated when cattle were stocked and removed. During periods when estimated forage biomass fell <1 Mg ha −1 for most paddocks, cattle were moved to a nearby holding area of endophyte-free fescue and supplemented with hay until sufficient forage regrowth occurred (1 to 3 Mg ha −1 ) to allow restocking. The site is comprised of 14 tall fescue paddocks, 0.9 to 1.1 ha each. Cattle were allowed to graze in each paddock for variable time periods in both summer and winter. Stocking was in September and used weaned Angus (Bos taurus) heifers (2002)(2003)(2004)(2005)(2006)(2007) and steers (2008)(2009) at rates of~3 to 6 head paddock-1 adjusted every 28 days to maintain 1 to 3 Mg ha −1 of available forage (Franzluebbers et al. 2009). For this study, only 12 of the 14 paddocks were used because two were subject to hay harvest instead of grazing. Each paddock was named based on type of fertilizer (broiler litter (B) or inorganic (I)) and tall fescue endophyte association (free of endophytes (F), novel-tall fescue MaxQ type with novel endophytes (N), or wild type endophytes (W)), with two paddocks of each combination type (i.e., BF, IF, BN, IN, etc.; Figure 1). Stocking rates (Equation (1)) were calculated (Table 1)

Climate Data
Precipitation data were provided by a NRCS SCAN weather station located 739 m (2425 ft) west of the study area. Daily data were aggregated into monthly and yearly data and compared with the long-term (30-

Stocking Rate
A one-way analysis of variance (ANOVA) test was conducted on the stocking rates among each of the paddock types. The analysis was performed within RStudio using the

Stocking Rate
A one-way analysis of variance (ANOVA) test was conducted on the stocking rates among each of the paddock types. The analysis was performed within RStudio using the stats package. One-way ANOVA tests with significant (p < 0.05) p-values indicate that some of the group means are different.

Vegetation Indices Derived from Landsat 7 ETM+
Employing the use of Landsat 7 ETM+ satellite imagery, we evaluated the plant health proxies of different combinations of fertilization regime (broiler litter or inorganic) and tall fescue endophyte association (free, novel, or wild). This satellite was chosen due to the limited availability of satellites at the timeframe of the investigation (2005-2009). The use of Landsat 5 TM was considered, however, due to the differences among the sensors (e.g., band widths etc.) and potential effects on the sensitivity analysis, we ultimately determined it best to utilize only Landsat 7 ETM+ imagery. We cover available satellite options further in the discussion section. Nevertheless, Landsat 7 ETM+ has a spatial resolution of 30 m and a temporal resolution of 16 days. The spectral bands used by Landsat 7 ETM+ for this study are blue (0.45-0.52 µm), red (0.63-0.69 µm), near infrared (0.77-0.90 µm), and shortwave infrared 1 (1.55-1.75 µm). Moreover, the quality assessment (QA) bands were incorporated, represented as bands 1, 3, 4, 5, and QA, respectively. A well-documented issue with Landsat 7 ETM+ is the failing of the scan line corrector (SLC) mechanism on May 31, 2003 [28]. All imagery collected since then exhibit scan-to-scan gaps, resulting in missing data for some paddocks during some data collection.
Google Earth Engine Code Editor was used to access the collection of satellite imagery, which was then filtered to only show images within the study area during the desired timeframe. The imagery utilized was reflectance data from tier 1 bottom-of-atmosphere (BOA) data. These data have been atmospherically corrected using the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS), and includes a cloud, shadow, water and snow mask produced using C Function of Mask (CFMASK), as well as a per-pixel saturation mask. The QA band was used to identify and mask pixels that had a high confidence of cloud coverage to minimize the risk of interference with vegetation index values. The vegetation indices calculated for this study include NDVI (Equation (2)), EVI (Equation (3)), and LSWI (Equation (4)): where NIR, Red, Blue, and SWIR represent the reflectance values of the spectral bands of Landsat 7 ETM+. The vegetation index image collections were then exported to run zonal statistics within ArcGIS v.10.3 (ESRI, Redlands, CA, USA). The zonal statistics were used to spatially average the pixels within a given paddock. The spatial averages for each paddock over time provide the basis for sensitivity analyses and the subsequent spatiotemporal analytics.

Relative Sensitivity Analysis
The three spectral vegetation indices used in this paper have been shown to accurately estimate biophysical parameters of vegetation, but alone they are not capable of quantifying the detailed relationship between the vegetation indices and the biophysical parameters they are measuring [29]. The sensitivity of vegetation indices cannot be assumed to be a constant value, but instead a function of the biophysical parameter. Thus, a sensitivity function is required to clarify the change in sensitivity over the range of the biophysical pa-rameter. Similar to Ji and Peters [29], here we employ a "relative sensitivity function" used to compare the sensitivities between two vegetation indices where in situ measurements of biophysical parameters are unavailable. To begin, we utilize y and x as the dependent and independent variables, respectively, to assume a y on x regression (Equation (5)).
whereŷ' represents the first derivative of the y on x regression function which provides the rate of change. Both the function and the error term can be linear or nonlinear; thus, there are two potential equations for the error term σ that include the linear (Equation (6)) and nonlinear (Equation (7)) error calculations.
where σŷ represents the mean square error of a linear function and X represents a matrix of independent variables with X i representing the ith row of X.
where σŷ represents the mean square error of a nonlinear function and F represents a matrix of derivatives for approximation of a least squares estimation with F i representing the ith row of the F matrix (see [30] for more details). We then inverted y on x regression to an x on y regression (Equation (8)) where we again calculated the first derivative of the x variable while employing the use of a mean square error utilizing similar linear and nonlinear methodologies as seen in the equations above (Equations (6) and (7)).
Following the calculations exhibited above, plotting of both S x|y and S y|x provides a demonstration of the relative sensitivities of the two indices being compared. It is important to note that we are comparing three indices, thus each unique pairing (e.g., NDVI vs. EVI, EVI vs. LSWI, and NDVI vs. LSWI) will be explored. The unique pairings will be analyzed across all years (normal, wet, and dry) to identify the most sensitive index to be utilized in determining tall fescue plant health among varying climate conditions.

Time Series Analyses of Tall Fescue Plant Health
Following the most sensitive index identification, we then used that index throughout for the remainder of the analytics. The remaining methods included ANOVA tests for the stocking rate and the highest sensitivity index values among entire years and grazing periods for normal, wet, and dry climates to identify significant (p < 0.05) differences.
Using RStudio, a script was written to import the values of each zonal statistics file along with its respective paddock number and index type, and then combine all of them into a single data frame usable by the ggplot2 package for plotting and analysis. The index values were plotted for each separate paddock with the precipitation data and grazing period overlaid. A point was plotted for each available image to show the temporal distribution.

Boxplot Comparison of Fertilizer Type with Endophyte Association
Ggplot2 was also used to create a box plot of the most sensitive vegetation index values, comparing fertilizer types across each endophyte-type. Box plots require only a minimum sample size of five, compared to the 30 with histograms, but box plots are also more effectively compared across three or more samples [31]. With the variation in sample size between years, this made box plots a viable choice amongst graphing tools to illustrate the spread and difference between samples. A separate box plot was created for the normal, wet, and dry years. In addition to the standard box plot format, the mean is included in each and is represented by a hollow circle.

Analysis of Variance (ANOVA)
Two-way and interaction ANOVA tests were ran with the most sensitive vegetation index, with fertilizer type and endophyte association as the independent variables. The data were sorted into three categories based on precipitation: normal, wet, and dry. Findings that can be considered significant were subject to a Tukey post hoc test to determine which groups the differences occurred between.

Stocking Rate
Each paddock had a different number of grazing cattle. Stocking rates were calculated for each paddock in each year ( Figure 3). ANOVA tests showed statistical significance (p < 0.00) among the difference in stocking rates between paddock types. Wild-type paddocks (4) had greater stocking rates than other paddock types.

Landsat Vegetation Index Collection
Due to the revisit period (~16 days) of Landsat and the removal of refle using the cloud mask, the number of images collected varied annually dur period (Table 2). Moreover, the presence of two paddocks replicating the e fertilizer type provides twice the sampling for each satellite revisit, except

Landsat Vegetation Index Collection
Due to the revisit period (~16 days) of Landsat and the removal of reflectance values using the cloud mask, the number of images collected varied annually during the study period (Table 2). Moreover, the presence of two paddocks replicating the endophyte and fertilizer type provides twice the sampling for each satellite revisit, except in the occurrence of one paddock being excluded because of clouds. The values for each image were sorted by paddock type and averaged to create the mean annual values for each of the indices ( Table 3). The mean annual NDVI values range from 0.3879 to 0.6938. The mean annual EVI values range from 0.2645 to 0.5280. The mean annual LSWI values range from −0.0662 to 0.2251. All paddock types had low vegetation index values throughout most of 2007, with BN and BW performing slightly better than the rest.

Sensitivities Analysis of NDVI, EVI, and LSWI Indices
Among the sensitivity analyses, we first analyzed the linear relationships between the unique pairings of the three indices (e.g., NDVI vs. EVI, EVI vs. LSWI, and NDVI vs. LSWI) across climate conditions (dry, wet, and normal; Figure 4). NDVI and EVI had strong linear relationships in all climate conditions (normal, wet, or dry). LSWI had linear relationships with NDVI and EVI in normal and wet years. However, the linear relationship between LSWI and NDVI/EVI was not strong in dry years. There are also many dispersed points in the LSWI versus NDVI/EVI relationships. We first analyzed sensitivities between NDVI/EVI because they are highly correlated, and then compared the more sensitive option of the two to LSWI. It is important to note that we scaled LSWI to positive values by adding the greatest negative value to all LSWI values. This was done because the sensitivity calculations do not accept negative values. According to the sensitivity analysis ( Figure 5), EVI was more sensitive than NDVI in wet years, but LSWI was more sensitive across a wider range of values. There was little deviation between each of the indices in normal years, but LSWI showed sensitivity across a wider range than both NDVI and EVI. NDVI showed much greater sensitivity than EVI in dry years, but LSWI showed sensitivity across a value range nearly triple that of NDVI. In agriculture, we are generally more interested in the extreme values (i.e., to identify healthy and non-healthy portions of the paddock) than central values. NDVI showed higher peak sensitivity values overall, but its sensitivity carried across a shorter range in nearly every scenario. We first analyzed sensitivities between NDVI/EVI because they are highly correlated, and then compared the more sensitive option of the two to LSWI. It is important to note that we scaled LSWI to positive values by adding the greatest negative value to all LSWI values. This was done because the sensitivity calculations do not accept negative values. According to the sensitivity analysis ( Figure 5), EVI was more sensitive than NDVI in wet years, but LSWI was more sensitive across a wider range of values. There was little deviation between each of the indices in normal years, but LSWI showed sensitivity across a wider range than both NDVI and EVI. NDVI showed much greater sensitivity than EVI in dry years, but LSWI showed sensitivity across a value range nearly triple that of NDVI. In agriculture, we are generally more interested in the extreme values (i.e., to identify healthy and non-healthy portions of the paddock) than central values. NDVI showed higher peak sensitivity values overall, but its sensitivity carried across a shorter range in nearly every scenario.

Time Series of LSWI Values Across Paddocks
The values of the LSWI index consistently followed the pattern of the precipitation data, with minor variations appearing due to missing dates of imagery. Dry months were followed by lower index values, while wet months, during every season except winter, were followed by higher index values. The time series charts ( Figure 6; see Supplemental  Table S1 for in-depth summary statistics) show discrepancies between paddocks of the same type. There are sharp declines showing in different months within paddock types in the 2006 charts. Reflecting the number of images shown previously in Table 2

Time Series of LSWI Values Across Paddocks
The values of the LSWI index consistently followed the pattern of the precipitation data, with minor variations appearing due to missing dates of imagery. Dry months were followed by lower index values, while wet months, during every season except winter, were followed by higher index values. The time series charts ( Figure 6; see Supplemental  Table S1 for in-depth summary statistics) show discrepancies between paddocks of the same type. There are sharp declines showing in different months within paddock types in the 2006 charts. Reflecting the number of images shown previously in Table 2

Boxplot Comparison of Fertilizer Type with Endophyte Association
Separated by paddock type, the boxplots (Figure 7) show which paddock types had higher mean LSWI values during the dry, normal, and wet years. During years with normal precipitation, the inorganic fertilizer combined with the wild-type endophyte had a slightly higher median LSWI value than the others, but the maximum values for all of the broiler litter paddock types were higher than all of the inorganic paddock types. During the dry year, the endophyte-free type fertilized with broiler litter performed the worst. However, all paddock types performed nearly the same during the wet year. Remote Sens. 2021, 13, x FOR PEER REVIEW 13 of 20

Boxplot Comparison of Fertilizer Type with Endophyte Association
Separated by paddock type, the boxplots ( Figure 7) show which paddock types had higher mean LSWI values during the dry, normal, and wet years. During years with normal precipitation, the inorganic fertilizer combined with the wild-type endophyte had a slightly higher median LSWI value than the others, but the maximum values for all of the broiler litter paddock types were higher than all of the inorganic paddock types. During the dry year, the endophyte-free type fertilized with broiler litter performed the worst. However, all paddock types performed nearly the same during the wet year.   Table S1. Each paddock was named based on if broiler litter (B) or inorganic (I) fertilizer was applied and if the tall fescue type was free of endophytes (F), novel-tall fescue MaxQ type with novel endophytes (N), or wild type endophytes (W), with two paddocks of each combination type within the study area (i.e., BF, IF, BN, IN, etc.).

Boxplot Comparison of Fertilizer Type with Endophyte Association
Separated by paddock type, the boxplots (Figure 7) show which paddock types had higher mean LSWI values during the dry, normal, and wet years. During years with normal precipitation, the inorganic fertilizer combined with the wild-type endophyte had a slightly higher median LSWI value than the others, but the maximum values for all of the broiler litter paddock types were higher than all of the inorganic paddock types. During the dry year, the endophyte-free type fertilized with broiler litter performed the worst. However, all paddock types performed nearly the same during the wet year. (1 year) years for each endophyte association (EA; free, novel, and wild) and fertilizer type (broiler litter or inorganic). The large circles represent the mean value for the respective paddocks. It is worth noting that inorganic wild and broiler litter free were significantly different (p < 0.05) from other pairs of endophyte and fertilizer types.

Statistical Significance of Fertilizer Type and Endophyte Association on LSWI
The interaction ANOVA test showed statistical significance (p = 0.05) by fertilizer type within the grazing period during only the dry year. Every other test resulted in relatively high p-values (Tables 4 and 5). Table 4. Two-way ANOVA results for fertilizer (broiler litter or inorganic) and endophyte associations (EA; free, novel, or wild). The median LSWI values within the grazing period during the normal years were similar across all fertilizer/endophyte association types, with broiler litter/free resulting in the lowest and broiler litter/wild resulting in the highest. The wet year resulted in nearly identical results. However, the dry year showed broiler litter/free performing the worst and inorganic/wild performing slightly above all other types.

Satellite Options
Landsat 7 ETM+ imagery was used for this study because the timeframe did not allow for the use of satellites with higher temporal and spatial resolutions, which presented some potential issues with data availability. While Landsat 7 ETM+ has a temporal resolution of 16 days, the SLC failure and cloud cover removal caused many of the images to have missing data within some of the paddocks. This resulted in an uneven number of images between years and paddocks within the same scan. Particularly, some paddocks had only one or two images during the grazing period of one of the years (2009). For more recent years, we now have the option of Sentinel-2 Multispectral Instrument (MSI) imagery with spatial resolutions of up to 10 m and a temporal resolution of five days. Without the SLC failure issue with Landsat 7 ETM+, Sentinel-2 MSI has potential to still provide a better, more complete dataset to work with.

Vegetation Index Sensitivities
While we found that LSWI had the greatest sensitivity, there are some overlapping results that show NDVI or EVI would be better for particular vegetation index value ranges. For years with normal precipitation, NDVI was more sensitive than LSWI between the value range of~0.4 and~0.8. Similarly, EVI was more sensitive than LSWI between the value range of~0.37 and~0.6 during wet years. Nevertheless, LSWI was more sensitive across a greater range of values, especially those of the extreme values that farm managers are generally more interested in identifying within their paddocks.

Tall Fescue Grassland Performance
Tall fescue supports many cattle operations in Georgia and it is important to understand what combination of endophyte association paired with fertilizer type performs best across varying climatic conditions. A recent study concerned with environmental conservation (e.g., hydrology) found that tall fescue fertilized with inorganic fertilizer has 30% greater runoff than with poultry litter [24]. From the same study, Franzluebbers et al. [23] found results suggesting that broiler litter was as effective as inorganic fertilizer for pasture growth and cattle production. They also found that cattle body weight gain from novel endophyte tall fescue paddocks were as good as those from endophyte free tall fescue paddocks.
In this study, the mean annual LSWI values showed a clear decline with precipitation, but the two years following the dry year did not follow the same pattern ( Figure 8). The year directly following the dry year showed the highest mean annual LSWI values in all types except for the combination of broiler litter and endophyte free (BF) tall fescue paddocks. However, the BF paddocks were the only paddocks that showed an increase during the following year with nearly twice as much precipitation. It is also important to note that broiler litter and novel endophyte (BN) paddocks showed the highest LSWI values during the dry year. The combination of these two paddocks' performance, combined with the lesser runoff of broiler litter when compared to inorganic fertilizers (Endale et al. 2013), suggests that the use of broiler litter on tall fescue paddocks could benefit both farm managers and the environment.

Manuresheds
There are plentiful poultry-based manureshed sources in Georgia, which can make broiler litter more readily available at comparable or cheaper prices than inorganic fertilizer (see [23]) in the same area ( Figure 9). The results here show that there is a significant difference in LSWI values among paddock types in certain climatic conditions (especially dry years). The results of this research combined with the other studies mentioned here support increased use of LSWI as a proxy for tall fescue plant health, as well as consideration for increased adoption of broiler litter for fertilization purposes around the Georgia Piedmont.

Manuresheds
There are plentiful poultry-based manureshed sources in Georgia, which can make broiler litter more readily available at comparable or cheaper prices than inorganic fertilizer (see [23]) in the same area ( Figure 9). The results here show that there is a significant difference in LSWI values among paddock types in certain climatic conditions (especially dry years). The results of this research combined with the other studies mentioned here support increased use of LSWI as a proxy for tall fescue plant health, as well as consideration for increased adoption of broiler litter for fertilization purposes around the Georgia Piedmont.

Climate Trends
Climate data show an increase in the variability of dry and wet periods worldwide [32]. Following climate prediction models, farm managers should choose the fertilizer type and endophyte association that performs best in their predicted climate conditions [23]. The findings here provide a remote sensing based method incorporating the use of LSWI to determine fescue plant health. Moreover, employing these methods suggests that the plant health of those fields using broiler litter versus those that incorporated inorganic fertilizers did not have significant differences. Thus, looking forward, grazingland research/managers can continue to explore the use of the more sustainable process of using manure-based nutrients.

Climate Trends
Climate data show an increase in the variability of dry and wet periods worldwide [32]. Following climate prediction models, farm managers should choose the fertilizer type and endophyte association that performs best in their predicted climate conditions [23]. The findings here provide a remote sensing based method incorporating the use of LSWI to determine fescue plant health. Moreover, employing these methods suggests that the plant health of those fields using broiler litter versus those that incorporated inorganic fertilizers did not have significant differences. Thus, looking forward, grazingland research/managers can continue to explore the use of the more sustainable process of using manure-based nutrients.

Conclusions
With all available options for vegetation indices to choose from, it can be difficult to decide which to use for a given application. Each research question has its own data needs and availability. A sensitivity analysis can identify the best option when comparing two indices, which in this case was LSWI. However, it is worth mentioning that for years with normal precipitation, NDVI was more sensitive than LSWI between the value range of ~0.4 and ~0.8. Similarly, EVI was more sensitive than LSWI between the value range of ~0.37 and ~0.6 during wet years. Nevertheless, LSWI was more sensitive across a greater range of values, especially for those of the extreme values that farm managers are generally more interested in identifying within their paddocks; thus, LSWI was utilized in the remainder of the research. The only significant difference in LSWI values between fertilizer type and endophyte combinations occurred during the year with low precipitation. The uneven temporal distribution of imagery available for the study area could have impacted some comparisons across years, so it may be worthwhile to repeat this study and focus on recent years with better available satellite data. Regardless of this limitation, the findings of this study support previous recommendations for using broiler litter to fertilize tall fescue in locations like the Georgia Piedmont. The data and findings also point to

Conclusions
With all available options for vegetation indices to choose from, it can be difficult to decide which to use for a given application. Each research question has its own data needs and availability. A sensitivity analysis can identify the best option when comparing two indices, which in this case was LSWI. However, it is worth mentioning that for years with normal precipitation, NDVI was more sensitive than LSWI between the value range of~0.4 and~0.8. Similarly, EVI was more sensitive than LSWI between the value range of~0.37 and~0.6 during wet years. Nevertheless, LSWI was more sensitive across a greater range of values, especially for those of the extreme values that farm managers are generally more interested in identifying within their paddocks; thus, LSWI was utilized in the remainder of the research. The only significant difference in LSWI values between fertilizer type and endophyte combinations occurred during the year with low precipitation. The uneven temporal distribution of imagery available for the study area could have impacted some comparisons across years, so it may be worthwhile to repeat this study and focus on recent years with better available satellite data. Regardless of this limitation, the findings of this study support previous recommendations for using broiler litter to fertilize tall fescue in locations like the Georgia Piedmont. The data and findings also point to the potential use of satellite imagery in assessing grazingland tall fescue health and advancing the concept of poultry manureshed in the region or elsewhere where poultry manure production is extensive.