Detecting Recent Crop Phenology Dynamics in Corn and Soybean Cropping Systems of Kentucky

Accurate phenological information is essential for monitoring crop development, predicting crop yield, and enhancing resilience to cope with climate change. This study employed a curve-change-based dynamic threshold approach on NDVI (Normalized Differential Vegetation Index) time series to detect the planting and harvesting dates for corn and soybean in Kentucky, a typical climatic transition zone, from 2000 to 2018. We compared satellite-based estimates with ground observations and performed trend analyses of crop phenological stages over the study period to analyze their relationships with climate change and crop yields. Our results showed that corn and soybean planting dates were delayed by 0.01 and 0.07 days/year, respectively. Corn harvesting dates were also delayed at a rate of 0.67 days/year, while advanced soybean harvesting occurred at a rate of 0.05 days/year. The growing season length has increased considerably at a rate of 0.66 days/year for corn and was shortened by 0.12 days/year for soybean. Sensitivity analysis showed that planting dates were more sensitive to the early season temperature, while harvesting dates were significantly correlated with temperature over the entire growing season. In terms of the changing climatic factors, only the increased summer precipitation was statistically related to the delayed corn harvesting dates in Kentucky. Further analysis showed that the increased corn yield was significantly correlated with the delayed harvesting dates (1.37 Bu/acre per day) and extended growing season length (1.67 Bu/acre per day). Our results suggested that seasonal climate change (e.g., summer precipitation) was the main factor influencing crop phenological trends, particularly corn harvesting in Kentucky over the study period. We also highlighted the critical role of changing crop phenology in constraining crop production, which needs further efforts for optimizing crop management practices.


Introduction
Vegetation phenology is defined as the development, differentiation, and initiation of plant organs [1]. Accurate retrieval of crop phenology information is a prerequisite for evaluating crop adaptation to climate change, modeling agricultural ecosystem carbon exchange, and predicting future agricultural production [2][3][4][5]. The Intergovernmental Panel on Climate Change has reported a change in global mean temperature of 1.5 • C above pre-industrial levels, along with changes in precipitation and an increased frequency of extreme climate events (IPCC, 2018). This shift in climate may bring varying degrees of impacts on agricultural ecosystems at different temporal and spatial scales. Crop phenology is closely related to climate change and is a critical indicator of optimum yield [6,7]. Therefore, it is essential to consider changes in crop phenology when assessing climate impacts on agricultural productivity, carbon cycling, and land-atmosphere feedbacks [8,9].
Many studies have shown that the climate impacts on agricultural ecosystems are reflected by variations in crop phenology, such as the advanced or delayed planting and harvesting dates [10][11][12]. For example, He et al. [13] reported that soybean planting dates were delayed by an average of 1.78 days/decade, and the growing season length was shortened by an average of 1.16 days/decade during 1981-2010 across the major soybeanproducing areas in China. Climate warming is a primary factor that drives phenological shifts [14], with temperature responses varying with crop types, locations, and study periods [15,16]. Many studies have investigated the responses of crop phenology to historical climate warming at regional to global scales. For example, Estrella et al. [17] reported that corn and oats sowing dates in Germany advanced in response to increases in March-May temperature at a rate of 0.60 days/ • C and 4.15 days/ • C, respectively. Based on corn phenology observations collected from agro-meteorological stations in China, Tao et al. [18] reported that the growing season lengthened during 1981-2009 due to combined effects of warming temperature, changing field practices, and shifting varieties. Model simulation results from Tubiello et al. [19] have shown that warmer temperatures accelerated plant phenology and further shortened the crop growing period, which resulted in crop yield reduction and potential food insecurity. Nevertheless, other climatic factors such as precipitation could determine the planting date more directly than the temperature in some regions [20,21]; however, few studies have explored the crop phenological changes and their relations with precipitation.
Remote sensing imagery can be considered an essential tool that complements fieldbased data collection approaches [22]. Numerous studies have reported the use of satellitebased Normalized Differential Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) for detecting crop phenology [23][24][25]. Some studies have shown good performance in identifying phenological stages of specific crop types using pre-defined VI thresholds [26]. For example, Sakamoto et al. [27] used a two-step filtering approach to detect the phenological stages of corn and soybean and achieved high accuracies at the site and regional levels. Huang et al. [28] applied dynamic thresholds of VI time series to detect the start and end of the season of different crop types and obtained higher accuracy than the results from fixed thresholds.
Kentucky is a traditional agricultural state, with corn and soybean being leading field crops. As a typical climatic transition zone, agriculture in Kentucky faces mixed climates that blend northern and southern weather patterns. Over the past 100 years, this region has not seen significant seasonal changes in temperature, especially during the crop growing season in this region [29]. Although crop phenological changes such as earlier planting dates have been widely reported under a warming climate [17,18,30], the associated spatial patterns are highly varied [31,32]. Uncertainties remain regarding how crop phenology has changed over areas like Kentucky, where temperature trends were generally flat over the past decades.
This study adopted a curve-change-based dynamic threshold approach along with MODIS NDVI time series and ground observations to detect the planting and harvesting dates for corn and soybean in Kentucky from 2000 to 2018. We also quantified the temporal trends of crop phenology and its responses to climatic factors (i.e., temperature and precipitation) and the correlations with crop yields. The objectives of the study are (1) to identify phenological dates of corn and soybean using MODIS NDVI time series in Kentucky from 2000 to 2018; (2) to evaluate the accuracy of estimated crop phenological stages using ground data at the state and county levels; (3) to characterize the spatialtemporal trends of crop phenological stages for corn and soybean in Kentucky during the study period; (4) to examine the correlations between crop planting/harvesting dates and temperature/precipitation variations; and (5) to analyze the effects of crop phenological change on crop yields.  (Figure 1). In general, Kentucky has a humid subtropical climate characterized by hot summers and cold to mild winters, with an oceanic climate found in the highlands of the southeast. The mean annual temperatures in Kentucky range from 11.67 • C in the northeast to 15 • C in the southwest. The state-wide annual precipitation is 1143 mm, with 965.2 mm for the northern part and 1270 mm for the southern part). Crops in Kentucky are predominantly corn and soybean, which account for more than 90% of total cropland in the state.
Remote Sens. 2021, 13, x FOR PEER REVIEW 3 period; 4) to examine the correlations between crop planting/harvesting dates and temperature/precipitation variations; and 5) to analyze the effects of crop phenological change on crop yields.

Study Area
In this study, we focused on the Commonwealth of Kentucky (36° 30′ N to 39° 9′ N and 81° 58′ W to 89° 34′ W) ( Figure 1). In general, Kentucky has a humid subtropical climate characterized by hot summers and cold to mild winters, with an oceanic climate found in the highlands of the southeast. The mean annual temperatures in Kentucky range from 11.67 °C in the northeast to 15 °C in the southwest. The state-wide annual precipitation is 1143 mm, with 965.2 mm for the northern part and 1270 mm for the southern part). Crops in Kentucky are predominantly corn and soybean, which account for more than 90% of total cropland in the state.

Ground Data
We acquired crop planting and harvesting dates of corn and soybean in Kentucky at both the state and county levels. Crop reports released by USDA National Agricultural Statistics Service (NASS) provided the state-level progress of crop phenology information of Kentucky from 2004 to 2018. We extracted the dates of 80% progress of planting and harvesting stages of corn and soybean from the crop reports. We also obtained 5-year averaged crop planting and harvesting dates from the same data source. The state-level crop yields were obtained from the USDA survey. The crop phenology datasets were derived from the Kentucky Hybrid Corn Performance Tests and Kentucky Soybean Variety Performance Tests. These tests offered annual planting and harvesting dates of corn and soybean from 2000 to 2018. The details of the ground data are shown in Table 1.

Ground Data
We acquired crop planting and harvesting dates of corn and soybean in Kentucky at both the state and county levels. Crop reports released by USDA National Agricultural Statistics Service (NASS) provided the state-level progress of crop phenology information of Kentucky from 2004 to 2018. We extracted the dates of 80% progress of planting and harvesting stages of corn and soybean from the crop reports. We also obtained 5-year averaged crop planting and harvesting dates from the same data source. The state-level crop yields were obtained from the USDA survey. The crop phenology datasets were derived from the Kentucky Hybrid Corn Performance Tests and Kentucky Soybean Variety Performance Tests. These tests offered annual planting and harvesting dates of corn and soybean from 2000 to 2018. The details of the ground data are shown in Table 1.

Satellite and Ancillary Data
In this study, the MODIS NDVI time series calculated from the MCD43A4 product (version 6, https://lpdaac.usgs.gov/products/mcd43a4v006/; accessed on 20 February 2021) was used to detect the planting and harvesting dates of corn and soybean in Kentucky from 2000 to 2018 [33]. MCD43A4 provides 500-m and daily surface reflectance of seven bands in a Sinusoidal projection system, available from February 2000 to the present.
The crop maps from NASS Cropland Data Layers (NASS-CDL) (https://nassgeodata. gmu.edu/CropScape/; accessed on 20 February 2021) were used to identify specific locations of corn and soybean fields. The NASS-CDL classifies specific crop types and provides multi-year crop classification maps at 30 m resolution for the conterminous United States. This classification map is available from 2008 to 2018 for Kentucky.
We used gridded monthly air temperature and precipitation from Daymet to examine the relationships between climate change and crop phenological development [34]. These climate datasets include minimum/maximum temperature and precipitation at a 1km spatial resolution (https://daymet.ornl.gov/; accessed on 20 February 2021). We calculated the monthly average air temperature based on the maximum and minimum temperatures.

Time Series Data Processing
We processed MODIS daily reflectance data on the Google Earth Engine Platform. The NDVI was calculated from the reflectances of the RED and NIR bands as follows [35]: where ρ RED and ρ N IR are band 1 (0.620-0.670 µm) and band 2 (0.841-0.876 µm) reflectances from the MODIS product, respectively. It was necessary to smooth the time series data using smoothing functions before extracting phenological dates. The smoothing methods consider the noise bias caused by snow or clouds and can handle missing data. Here, we used the Harmonic analysis method to smooth the NDVI time series. This Harmonic algorithm can smooth and reconstruct remotely sensed VI time series while reducing the influence of clouds at the pixel level [36].

Detection of Crop Planting Dates, Harvesting Dates, and Crop Growth Period
In this study, the definitions of crop phenological stages were from USDA NASS ( https://www.nass.usda.gov/Publications/National_Crop_Progress/terms_definitions; accessed on 20 February 2021). We considered the silking stage of corn and the blooming stage of soybean as heading dates, respectively. We used a curve-change-based dynamic threshold approach on NDVI time series to identify crop planting and harvesting dates for corn and soybean in Kentucky from 2000 to 2018.
The corn and soybean areas were extracted using the NASS-CDL maps from 2008 to 2018. The original 30 m CDL maps were aggregated into 500 m maps with the percentages of corn or soybean areas being calculated in each 500 m pixel, respectively, to match the size of the MODIS pixel. Pixels with corn or soybean > 50% were retained for crop phenology detection. Previous studies have shown that the NDVI increases with leaf green-up during the spring season and decreases with leaf senescence in the fall [37,38]. As VI values in croplands generally exceed 0.4 at peak growth [39], spurious peaks were discarded if the corresponding NDVI values were less than 0.35. We then set a threshold of 0.35 to limit the cropland, i.e., the pixels with the maximum NDVI values less than 0.35 were excluded as non-cropland cover types [40].
For each crop pixel in a given year, the first and the second derivatives of the NDVI curve were defined by the following equations: Remote Sens. 2021, 13, 1615 where f ' and f " are the first-and second-order derivatives of the smoothed NDVI time series (f), i is the time sequence number of values in the smoothed NDVI time series (2, 3 . . . 365), 1 is the time step of NDVI time series, and f is the smoothed NDVI time series. We then identified crop phenological dates based on the characteristics of the derivatives: Heading dates: Previous studies have shown that the maximum NDVI occurs around the heading dates [41]. We, therefore, used the point at the NDVI peak to capture crop heading dates and constrained the valid range according to the five-year averaged planting dates from the crop reports dataset ( Table 2).
where f ' is the first-order derivative of the NDVI curve; f is smoothed NDVI curve; i means the ith of NDVI/NDVI' values in the time series (1, 2, 3 . . . 365), a and b are the upper and lower boundaries of the valid time range for NDVI peak, respectively. Planting dates: The NDVI curve shows lower values before crop planting when agricultural lands are plowed or cultivated ( Figure 2). After the crop planting, photosynthetic activity starts with plant expanded leaves, and thereby, the NDVI curve begins to increase. It is reasonable to expect the NDVI value of the planting date is located at the low point at the early stage of the NDVI curve. We, therefore, applied the peak of the second-order derivative of the NDVI curve (before the heading date) to detect the crop planting date. The crop planting dates were constrained within 40-120 days before the heading dates based on the Corn and Soybean Production Calendar in Kentucky (https://simpson.ca.uky.edu/files/corn_ and_soybean_production_calendar.pdf; accessed on 20 February 2021). Additionally, we also used more accurate ranges to filter out all possible outlier estimates according to the 5-year averaged phenology derived from the crop reports dataset ( Table 2). heading date) of the second-order derivative of the NDVI curve can catch the lowest value of NDVI at the last period of the NDVI curve ( Figure 2). Here we used this transition point to detect the crop harvesting date. Similarly, the harvesting dates were constrained to occur within the time range of 30-110 days after the heading date, according to the crop calendar in Kentucky. Similarly, we retained estimates that fall into the valid time range as determined by the 5-year averaged harvesting dates (Table 2).
Subsequently, crop growing season length was calculated for each pixel using the time difference between planting and harvesting dates.

Harvesting dates:
Plant leaves continue to wither and die during the harvesting season. Crop canopy can be harvested in this stage. Correspondingly, the NDVI value decreases to the lowest point when the crop is harvested from fields. After the heading date, the peak (after the heading date) of the second-order derivative of the NDVI curve can catch the lowest value of NDVI at the last period of the NDVI curve ( Figure 2). Here we used this transition point to detect the crop harvesting date. Similarly, the harvesting dates were constrained to occur within the time range of 30-110 days after the heading date, according to the crop calendar in Kentucky. Similarly, we retained estimates that fall into the valid time range as determined by the 5-year averaged harvesting dates (Table 2).
Subsequently, crop growing season length was calculated for each pixel using the time difference between planting and harvesting dates.

Evaluation and Trend Analysis
At the state level, we calculated the dates when the areas of estimated phenological dates occupied 80% of the total planting areas across the whole state for corn and soybean. For county-level evaluation, the mean values of the estimations were calculated for corn (68 counties) and soybean (74 counties) in top producer counties. The coefficient of determination (R 2 ) and root mean square error (RMSE) were used to evaluate the estimated crop phenology against the ground data at both the state and county levels.
where n represents the number of samples. y i and x i are the ground data and remote sensing estimates, respectively. Linear regression analysis was applied for generating the changing trends of the phenological estimations at the state level over the study period. We also used the Mann-Kendall test [42,43] and the Sen's slope estimator [44] to analyze the temporal trends of phenological stages at the pixel scale. During the process, pixels with more than 12 years being identified as an individual crop (corn or soybean) were included in the Mann-Kendall test. The analytical method was implemented using the R computing environment [45].
We used linear regression analysis to examine temporal patterns of climatic factors and crop yields and their relationships with crop phenology. The Pearson correlation coefficient was adopted to describe the sensitivity of crop phenology to climate change. Climatic factors include minimum, maximum, average temperatures, and accumulated precipitation during three seasons (spring: March-May, summer: June-August, fall: September-November) and the whole crop growing period.

State-Level Evaluation
The state-level evaluation results showed that crop phenology estimated by remote sensing was at a high level of agreement with the crop reports from the survey data ( Figure 3). The estimated harvesting dates closely matched those from the crop reports, with R 2 of 0.92 and 0.90 for corn and soybean, respectively ( Figure 3b). The R 2 of the estimated planting dates of corn and soybean against survey data was 0.87 and 0.79, respectively. The accuracy of the estimated harvesting dates of soybean was the highest, with an RMSE of 3.34 days. The RMSE value of corn harvesting dates was 3.82 days. The accuracies of the estimated planting dates of corn and soybean were 6.05 and 3.70 days, respectively ( Figure 3a).

County-Level Evaluation
The county-level assessment appeared to show lower accuracies compared to the statelevel assessment ( Figure 4). Generally, the estimated crop phenological dates were later than those observed from field tests. Overestimations were larger in estimated planting dates than harvesting dates for both corn and soybean. The RMSE values of corn planting and harvesting dates were 10.84 and 10.93 days, respectively. For soybean, the RMSE of harvesting dates was 9.17 days, and the RMSE value of planting dates was 12.26 days.

County-Level Evaluation
The county-level assessment appeared to show lower accuracies compared to the state-level assessment ( Figure 4). Generally, the estimated crop phenological dates were later than those observed from field tests. Overestimations were larger in estimated planting dates than harvesting dates for both corn and soybean. The RMSE values of corn planting and harvesting dates were 10.84 and 10.93 days, respectively. For soybean, the RMSE of harvesting dates was 9.17 days, and the RMSE value of planting dates was 12.26 days.

Spatial Distribution of Crop Phenology
The corn planting and harvesting dates were mapped every six years from 2000 to 2018 in Kentucky ( Figure 5). The corn cultivation areas were mainly distributed in western Kentucky. The corn planting dates showed obvious spatial differences in 2000 and 2006,

Spatial Distribution of Crop Phenology
The corn planting and harvesting dates were mapped every six years from 2000 to 2018 in Kentucky ( Figure 5). The corn cultivation areas were mainly distributed in western Kentucky. The corn planting dates showed obvious spatial differences in

Changing Trends of Crop Phenology
Phenological trends were analyzed for corn and soybean at the state level over the study period (Figure 7). The crop planting dates were slightly delayed by 0.01 days/year for corn and 0.07 days/year for soybean. Corn harvesting dates were delayed by an average rate of 0.67 days/year, while a slightly advanced pattern (0.05 days/year) in the soybean harvesting dates was detected. The inter-annual variation in the crop growing season length was related to the changing planting and harvesting dates. For soybean, a slightly shortening trend was found at a rate of 0.12 days/year, i.e., 2.28 days over the entire study period. However, the corn growing season experienced an increasing tendency by an average rate of 0.66 days/year, i.e., 12.54 days over the study period.

Changing Trends of Crop Phenology
Phenological trends were analyzed for corn and soybean at the state level over the study period (Figure 7). The crop planting dates were slightly delayed by 0.01 days/year for corn and 0.07 days/year for soybean. Corn harvesting dates were delayed by an average rate of 0.67 days/year, while a slightly advanced pattern (0.05 days/year) in the soybean harvesting dates was detected. The inter-annual variation in the crop growing season length was related to the changing planting and harvesting dates. For soybean, a slightly shortening trend was found at a rate of 0.12 days/year, i.e., 2.28 days over the entire study period. However, the corn growing season experienced an increasing tendency by an average rate of 0.66 days/year, i.e., 12.54 days over the study period. In addition, widespread negative tendencies were detected for the phenological estimations of corn and soybean from pixel to pixel in Kentucky from 2000 to 2018 (Figures 8  and 9). For corn, the slope values in Figure 8a,d,g showed that more than 40% of corn production areas for planting and harvesting dates, about a third of the area for growing season length, experienced phenological changes. Pixels with unchanged slopes (slope = 0) accounted for more than half of the total pixels for corn planting dates, harvesting dates, as well as growing season length. The pixels with significant changing trends were scattered across the corn production areas (Figure 8b,e,h). The statistical histograms in Figures  8 and 9 displayed all changing slopes that were significant or not in crop planting dates, harvesting dates, and growing season length. According to the statistics (Figure 8c,f,i), the proportions of negative trends (slope < 0) were much larger than those of positive trends (slope > 0), indicating the advanced corn planting dates, harvesting dates, and shortened growing season length over the study area. Figure 8c showed an evident advanced trend over the significant values for corn planting dates. However, for corn harvesting dates (Figure 8f) and the growing season length (Figure 8i), proportions of significant trends with negative and positive slope values were roughly equivalent. In addition, widespread negative tendencies were detected for the phenological estimations of corn and soybean from pixel to pixel in Kentucky from 2000 to 2018 (Figures 8 and 9). For corn, the slope values in Figure 8a,d,g showed that more than 40% of corn production areas for planting and harvesting dates, about a third of the area for growing season length, experienced phenological changes. Pixels with unchanged slopes (slope = 0) accounted for more than half of the total pixels for corn planting dates, harvesting dates, as well as growing season length. The pixels with significant changing trends were scattered across the corn production areas (Figure 8b,e,h). The statistical histograms in Figures 8 and 9 displayed all changing slopes that were significant or not in crop planting dates, harvesting dates, and growing season length. According to the statistics (Figure 8c,f,i), the proportions of negative trends (slope < 0) were much larger than those of positive trends (slope > 0), indicating the advanced corn planting dates, harvesting dates, and shortened growing season length over the study area. Figure 8c showed an evident advanced trend over the significant values for corn planting dates. However, for corn harvesting dates (Figure 8f) and the growing season length (Figure 8i), proportions of significant trends with negative and positive slope values were roughly equivalent.   In Figure 9a,d,g, more than a third of pixels with positive or negative trends were observed in soybean planting dates, harvesting dates, and growing season length. Similarly, pixels with unchanged trends (slope = 0) made up to more than half of the total values for three soybean phenological variables. Larger proportions with significant trends occurred in soybean planting and harvesting dates compared with those of corn (Figure 9b,e). Similar significant proportions were found in the growing season length of soybean and corn (Figures 8h and 9h). Advanced in soybean planting and harvesting dates were detected with high proportions of negative values (Figure 9c,f). However, Figure 9i showed that pixels with extended soybean growing season length (slope > 0) accounted for more than 60%, and comparable percentages of positive and negative slopes showed significant trends.

Trends of Climatic Factors and Its Correlation with Crop Phenology
Maximum temperatures decreased in three seasons and ranged from −0.001 to −0.01 • C/year in Kentucky from 2000 to 2018 (Table 3). Warming trends in minimum and average temperatures were observed, ranging from 0.03 to 0.05 • C/year and from 0.01 to 0.03 • C/year, respectively. Specifically, the minimum temperature during the growing season showed a significant increasing trend with a rate of 0.05 • C/year. Accumulated precipitation increased over time in all seasons in Kentucky. Historical climate records showed that summers from 2014 to 2018 are among the ten wettest summers over the last 30 years in Kentucky (http://kyclimate.org/climtrends.html; accessed on 20 February 2021). Thus, over the years tested, the summer climate trended wetter in Kentucky. The crop planting/harvesting dates were negatively correlated with three temperature variables but positively correlated with the precipitation for both crops (Tables 4 and 5). Corn planting dates showed significant correlations with the accumulated precipitation in spring (r = 0.56 for corn, r = 0.49 for soybean). Compared with soybean, corn planting dates were more sensitive to spring temperature, i.e., significant responses were −3.95 days/ • C and −3.70 days/ • C to Tmax and Tavg, respectively. For harvesting dates, higher correlation coefficients with temperature and precipitation were observed for corn and soybean. In particular, the correlations between crop harvesting dates and Tmax in summer/April-October were higher than 0.7 for corn and 0.65 for soybean. Significant relationships were also found between harvesting and the accumulated precipitation in summer/April-October for corn, and in fall/April-October for soybean, respectively. Corn growing season length exhibited negative sensitivities to temperature variables. Apart from a negative correlation in spring, positive relationships were detected between corn growing season length and the accumulated precipitation. Soybean growing season length was negatively correlated with all climatic factors except with the accumulated precipitation in fall and April-October. Significant correlations between growing season length and precipitation were mainly concentrated in summer/April-October for corn and summer/fall for soybean, respectively (Tables 4 and 5).

Trends of Crop Yield and Its Correlation with Crop Phenology
Crop yields showed significant increases in corn (2.19 Bu/acres per year, p < 0.05) and soybean (0.75 Bu/acres per year, p < 0.05), respectively, in Kentucky over the study period. A more noticeable increment was found in corn yield. However, we observed that corn yield consistently increased over time except for the sharp decrease in 2012 (68 Bu/acre), dramatically lower than the average corn yield (143 Bu/acre) of the study period. The reduced crop production was relevant to extreme heatwaves and drought during the summer [46].
We further investigated the relationships between the crop phenological dates and crop yields of corn and soybean using the linear regression analysis (Table 6). Over the 2000-2018 period, a significant positive correlation was found between corn growing season length and corn yield, suggesting that a one-day extension of the growing period increased 1.67 Bu/acres (p < 0.01) in corn yield. Furthermore, significant responses of harvesting dates to crop yields were detected for corn (trend = 1.37 Bu/acre per day, p < 0.01) and soybean (trend = 0.39 Bu/acre per day, p < 0.05), respectively.

Comparisons of Remote Sensing-Based Crop Phenology with Other Studies
Various vegetation phenology detecting methods have been developed in previous studies. However, mapping crop phenology is still challenging because the land surface vegetation dynamics or remote sensing phenology is different from crop physiological growth stages [47]. In this study, we detected the crop planting dates, harvesting dates, and growing season length for corn and soybean using an NDVI curve-change-based dynamic threshold approach. This approach linked characteristics of remote sensing vegetation index to crop physiological growth stages. NDVI time series at a daily time step allows for a high degree of coupling between remote sensing and crop growth stages. The threshold setting for each crop phenological stage based on survey data can further improve the reliability of the approach.
Our evaluation results suggested that our estimated crop phenological stages were favorably comparable with the results in previous studies. For example, using a remote sensing approach, Sakamoto et al. [27] reported the RMSEs of estimated phenological dates ranged from 0.7 to 8.6 days for corn and 1.9 to 14.5 days for soybean. In our study, RMSEs of estimated crop phenological dates were between 3.34 and 6.05 days at the state level and between 9.17 and 12.26 days at the county level. The lower accuracy in the county-level estimates was mainly due to limited available site-level field observations for evaluation. However, the state-level evaluation was based on the USDA crop report; its favorable performance illustrates the potential of the NDVI curve-change-based dynamic threshold approach in realistically estimating phenology for corn and soybean.

Spatial-Temporal Trends of Crop Phenology
We analyzed the state-level linear trends of three estimated crop phenological variables and built their spatial-temporal patterns by the Mann-Kendall test. Many studies have reported that earlier crop planting and extended growing seasons occurred during recent decades, but the changing trends vary depending on the study period [48,49].
Menzel et al. [50] showed that phenological trends were weaker for the most recent 30-year period (1989-2018) compared to the 1976-2005 period for both agricultural and wild plants. Kucharik [30] found that the planting date in approximately 75% across the 12 Corn Belt states was advanced by 0.37 days/year from 1979 to 2005. Notably, Kentucky was among the states with significant changes, with advanced corn planting dates at a rate of 0.8 days/year [30]. Sacks and Kucharik [12] reported that soybean planting dates advanced by 0.49 days/year averaged across the U.S. from 1981 to 2005. However, our study showed a slight delay in crop planting for both corn (0.01 days/year) and soybean (0.07 days/year) at the state level in Kentucky from 2000 to 2018. The changing patterns over different study periods implied that the earlier trend of crop planting season slowed down during the last two decades over the study area. In the Midwest U.S., soybean is usually planted after completing corn planting. Therefore, delayed corn planting dates might cause delayed soybean planting as well [51]. Sacks and Kucharik [12] also showed that the growing season length of corn and soybean was significantly extended by 0.67 and 0.30 days/year, respectively, in the U.S. during 1981-2005. Their findings of the prolonged corn growing season were similar to the results in our study (0.66 days/year).
Meanwhile, we found that soybean experienced a shorter growing season (0.12 days/year) in Kentucky during 2000-2018. According to Sacks and Kucharik [12], both corn (1 day/year) and soybean (0.83 days/year) experienced a trend of earlier harvesting in Kentucky over 1981-2015. However, our study showed largely postponed harvesting for corn (0.67 days/year) and slightly advanced harvesting for soybean (0.05 days/year) during 2000-2018. We also found that the longer corn growing season length could significantly benefit corn yield. The shortened soybean growth period may have undesired consequences for yield but allow more intercropping or earlier sowing of winter cereals [50]. The advanced harvesting dates and shortened growing season of soybean were probably related to the increasing double cropping system in Kentucky [52].

Effects of Climate Change and Other Factors on Crop Phenology
Temperature is often considered the most critical factor that influences crop phenological change. Many studies suggested advanced trends of crop phenology across the northern hemisphere due to the rising temperatures [53,54]. A long-term study showed that nearly all earlier planting events occurred in warmer years, and more than 80% of them were related to seasonal spring and summer temperatures [17]. In this study, only two climatic variables showed significant changing trends (minimum temperature in April-October and precipitation from June to August) in Kentucky from 2000 to 2018. Sensitivity analysis showed that crop phenology responded negatively to temperature and positively to precipitation, but no significant response was found with the growing season temperature. The overall analysis revealed that the phenology shifts (crop planting and harvesting) were not related to the increasing temperature during April-October.
Climate change raises the concern about how field management could be optimized to adapt to the changes in crop phenological development. The trend analysis showed that temperature did not have distinct warming trends in Kentucky over the study period. However, the crop phenology was observed with significant changes for both corn and soybean. As we discussed, the sensitivity analysis found that only the summer precipitation was significantly related to the delayed corn harvesting dates. The weak linkage between crop phenology and climatic variables implied that changing phenology is relatively dominated by human factors (e.g., management practices and variety improvement). This finding agrees with previous studies showing non-climatic factors (e.g., crop varieties, farmer decisions, cropping systems, and agronomic practices) may lead to changes in crop phenology. For example, Kucharik [30] suggested no strong evidence to support that the warming temperature was the most important factor driving the corn planting trends across the majority of the Corn Belt from 1979 to 2005.
Previous studies have also suggested the use of more advanced equipment and improved field management may be major contributors to planting trends in the spring season [55]. Lithourgidis et al. [56] suggested that farmer decisions were becoming much more efficient at field operations, especially adopting management practices that can facilitate earlier planting. Kucharik [30] presented that a trend toward performing tillage immediately after fall harvests might be one of the most significant changes in agronomic practices. Additionally, crop insurance policy restricts the final planting dates for different crop types in different regions, resulting in shifts in crop phenology [57].

Effects of Crop Phenological Shift on Crop Yield
Agricultural crop production is closely related to crop phenological change. Previous research presented that the optimum range of crop phenological stages can lead to high crop production [58]. Some studies suggested that warming climate advanced phenological phases and, consequently, shortened crop growth duration, potentially reducing crop yield [2,19]. However, this study found no significant effects of earlier planting on crop yield in Kentucky. Our result is consistent with Sacks and Kucharik [12], which verified that earlier planting did not show significant effects on crop yields across the U.S. Corn Belt. In addition, we found that planting dates did not show significant correlations with crop yield for both corn and soybean in Kentucky. In contrast, corn growing season length and harvesting dates contributed to the increased yield during the last two decades. This result is in agreement with Wu et al. [7], who suggested that a longer growth duration might increase agricultural production. These findings can serve as a benchmark by farmers to access crop phenology and its associated impacts on crop yield in Kentucky.

Uncertainty and Expectations
This study detected crop phenological stages using the remote sensing-based approach. The state-and county-level evaluations against ground-based datasets illustrated the robustness of this approach. However, this study still involves some uncertainties. First, agriculture in Kentucky is mainly concentrated in the north and west regions, with highly fragmented cropland areas in the rest of the state. In this case, MODIS products at a 500 m spatial resolution may not accurately capture the crop phenological stages in fragmented areas due to the effects of mixed and perimeter pixels. In the next step, data fusion of high-resolution satellite imagery is needed for reducing uncertainties in areas where mixed cropland pixels are dominant. Second, our county-level evaluation was based on very limited field observations; further effort is needed to collect more crop phenology data for robust validation and evaluation. Third, crop phenology is affected by many non-climatic factors. For example, the growing degree days for different crops are species-level characteristics. They may vary among varieties and are highly independent of the circumstances. At present, we cannot identify the impacts of crop variety on phenological shifts due to unavailable associated data not available. This knowledge gap might be filled as high spatial-temporal-spectrum resolution remote sensing images and site-specific variety information are available. Understanding the microclimate impacts on crop phenology will also require more available observations, such as data recorded from flux towers. Fourth, other crop phenological stages, such as the flowering and grain filling stages, also play critical roles in affecting crop development. Future work should endeavor to cover more phenological stages to analyze the effects of climate change and improve the capability of crop yield prediction. Additionally, future studies should explore the application of machine learning techniques on crop phenology detection, which might be a solution to improve the estimation accuracy [2,3].

Conclusions
In this study, using MODIS NDVI time series and ground datasets, we detected the planting dates, harvesting dates, and growing season length of corn and soybean in Kentucky from 2000 to 2018. We also investigated their temporal patterns and correlations with climate change and yields. Trend analysis showed that corn experienced delayed planting/harvesting dates and extended growing season length over the study period.
However, soybean was found to have delayed planting dates, an advanced harvesting season, and a shortened growing season length. Sensitivity analysis showed that increased seasonal climate temperature could significantly advance the planting and harvesting dates for both corn and soybean. Combining the climate variables and crop phenological patterns revealed that increasing accumulated precipitation in summer was substantially related to the delayed harvesting dates of corn in Kentucky over the study period. This study also suggested that the increasing corn yield had a significant correlation with the delayed harvesting dates and prolonged growing season. No significant correlation was found between climate change and soybean changing phenology. Moreover, changing phenological stages did not contribute to soybean yield. Our findings highlight the future needs to explore the impacts of non-climate-related factors on soybean phenology. The quantitative crop phenology responding to climate change and crop yields may provide farmers and local policy-makers guidelines for optimizing the field operations.