INTRODUCTION
Climate change has been a major interest to various fields of science over the past two decades, including agriculture [1,2]. Particular interest needs to be focused on changes in the frequency and intensity of abnormal climate events that can greatly affect crop production. In Korea, elevating winter-spring temperatures and summer heavy rainfall are considered as significant weather events caused by climate change [3-5]. The quantity of rainfall in the Korean monsoon season accounts for 40%–60% of total average annual rainfall [6], with a high likelihood that the prevalence of torrential rainfall increases with the increasing occurrence of heavy rainfall events. Furthermore, during the monsoon season, the event of typhoons also needs to be looked at because various summer crops have been seriously damaged by typhoon ahead of harvests [7]. In contrast to the relatively short-term rainfall-related weather events, Korean spring droughts can be fatal to forage crops due to their long-term impact. Fortunately, the sorghum-sudangrass hybrid (Sorghum bicolor L.) has high adaptability against drought (rainfall range of 200–700 mm) due to its good heat tolerance [8]. Therefore, this study considered, heavy rainfall and typhoon occurrence as the weather events in the monsoon season.
Forage corn (Zea mays L.) and sorghum-sudangrass hybrid are typical summer annual forage crops in Korea. In particular, the sorghum-sudangrass hybrid can be harvested multiple times a year. In Korea, the headless type sorghum-sudangrass hybrid is preferred over heading type sorghum-sudangrass hybrid because it is popularly used as whole crops due to the slow hardening of the stem [9]. Thus, to achieve a stable yield, the sorghum-sudangrass hybrid has been widely cultivated as a silage crop where forage crop is hard to produce. Furthermore, due to the reliable regrowth, sorghum-sudangrass hybrid can produce more forage production through additional harvest [10]. For sorghum-sudangrass hybrid that grows at a warmer temperature than silage maize, livestock farmers prefer headless type to heading type in Korea [11]. Cultivation of sorghum-sudangrass hybrid is advantageous due to flexibility of management options. sorghum-sudangrass hybrid can be a useful forage crop due to high forage production, tolerance to environmental stresses, and ratoon capability [12]. Furthermore, sorghum-sudangrass hybrid is tolerant to drought condition. The crop accumulates great amount of lignocellulosic biomass with fermentable sugars and starch [13]. Although forage quality was negatively related to cellulosic biomass production, some headless sorghum-sudangrass hybrids remains longer in a vegetative growth stage and maintain its quality [11]. Furthermore, in terms of production economics, cost per ton of sorghum forage production was the lowest comparing with switchgrass (Panicum virgatum L.), reed canarygrass (Phalaris arundinacea L.), big bluestem (Andropogon gerardii Vitman var. gerardii) and alfalfa (Medicago sativa L.) [14]. Therefore, sorghum-sudangrass hybrid has been a valuable forage crop.
Typhoon accommodating heavy rainfall damage Korean agriculture during summer, in particular, property damage due to heavy rainfall accounted for about 70% of all agricultural meteorological disasters [15]. In spite of gradual disaster mitigation and prevention measures, meteorological disasters occur every year under various conditions and intensities in the case of heavy rainfall and typhoons [16]. Forage production has been also damaged severely by lodging [17]. According to the international convention [18], rainfall defined by the amount per hour such as drizzle (slow-falling water droplets between 0.2 and 0.5 mm in diameter), light rain (accumulations of 2.5 mm per hour), moderate rain (accumulations of 2.5 to 7.6 mm per hour), and heavy rain (accumulations of over 7.6 mm per hour). Meanwhile, the heavy rainfall also can be defined as the rainfall amount over a given time. According to Lee et al. [19], the daily mean intensity of heavy rainfall was approximately 150–200 mm day−1 for 10 continuous-day in Korea. Furthermore, the total heavy rainfall in July was more than 1,000 mm month−1. In particular, it was found that heavy rainfall of 1,000 mm during monsoon season reduces corn production [20]. Hence, it was worth considering both the frequency and the total amount of heavy rainfall during monsoon season in this study. Meanwhile, typhoon is a tropical cyclone having a wind speed near 17 m/s and 200–600 km in diameter [21]. Since the regional impact of typhoons is changeable depending on their path and radius [22], we were interested in local typhoon warnings in this study. The effects of climatic, soil and cultivar on forage production of sorghum-sudangrass hybrid have been assessed by field trials in Korea [23]. In general, the first harvest of sorghum-sudangrass hybrid in Korea is closely related to the monsoon season. Since the estimated forage production decreased after 5 days from the expected optimal harvesting time [24], the length of monsoon season determined by the frequency of heavy rainfall is also an important factor in forage production. Thus, it was recommended to seed and harvest within 5 days to secure a stable yield. Long-term sorghum-sudangrass hybrid yield trends affected by temperature, precipitation, and sunshine duration, the relationship was not clearly evaluated [25]. Hence, research is necessary to understand in the significant weather event impact on forage crop production during summer growing season. The effect of weather is meaningful not only in the form of events, but also in the independent variable itself. Dixon et al. [26], the mean temperature, accumulated solar radiation, mean soil moisture and accumulated precipitation were considered to predict corn yield. Furthermore, growing degree days (GDD), accumulated rainfall and irradiation were considered for modeling yield trend of summer forage crops [25,27]. Therefore, the prediction of future sorghum-sudangrass production was developed based on climate related sorghum-sudangrass production data appearing in the previous trials.
A time series analysis was applied to the data to estimate yearly climate impacts which is one of the most effective statistical methods [28]. The prediction models developed for vetch (Vicia sativa L.), silage maize, sainfoin (Onobrychis viciifolia Scop.) and alfalfa in Turkey [29] forecasted forage amount until 2025. However, there was no consideration of climate and weather. For the forage production trend of whole crop maize in Korea [11], heavy rainfall during the monsoon season lead to substantial production reduction. Furthermore, due to the increasing frequency and intensity of heavy rain, forecasted corn production can be reduced by up to 70% of the past 30-years’ average. Meanwhile, review on sorghum-sudangrass hybrid without considering the effect of climatic variables confirmed the trend is closely related to the past two years and had no significant relation to the previous three years’ production [27]. Thus, to develop more robust prediction, it will be necessary to consider forage production trend of sorghum-sudangrass hybrid along with weather event.
Therefore, this study analyzed the impact of heavy rainfall and typhoon events on the forage sorghum production during monsoon season in the central inland regions of Korea.
MATERIALS AND METHODS
The raw data (n = 1,076), including the forage production (kg/ha), number of harvests, year, address, seeding and harvesting dates, etc., were collected from the reports on the variety trials conducted by the Rural Development Administration (RDA) of Korea. Usually, the damage caused by heavy rainfall and typhoon in the monsoon season varies considerably depending on local and topographical conditions; therefore, it was necessary to select a single location for homogeneity of data. Based on the regional classification chart of the Republic of Korea [30], sorghum-sudangrass hybrid is predominantly cultivated in central inland and southern inland regions. Furthermore, it was important role to record long-term data to allow for trends to be detected. Considering these two conditions, Suwon (latitude: 37° 15′ N, longitudinal: 127° 04′ E, n = 388) was selected as a representative central inland region of Korea (1988–2013). Maximum likelihood algorithm was applied to adjust missing year data in the modeling of forage production [31].
The weather data, including temperature, rainfall and sunshine duration, were collected from the Korean Meteorological Administration [KMA] via open-application programming interface (API). For inclusion of the weather event during the monsoon season, rainfall per hour was considered to differentiate heavy rainfall from ordinary rainfalls [11]. Where, the frequency and total amount were considered based on the hourly heavy rainfall and accumulated rainfall in July, respectively. GDD (°C) with the base temperature of 10°C, accumulated rainfall (mm), sum of sunshine duration (hr), heavy rainfall 1 (frequency of over 7.6 mm/hr rainfall amount), heavy rainfall 2 (dummy variable, 1: over 1,000 mm/July, 0: else), and the frequency of typhoon were generated from seeding to harvesting date as a variable in the study.
The forage production and weather variables were compiled by the cultivation year to make the data in a time series format. However, the sorting based on field type (upland field or paddy field) and cultivars was impossible because some annual values were undetectable. Since time series analysis can decompose data into seasonality and circularity [32], an autoregressive integrated moving average with exogenous inputs (ARIMAX) model was used in this study to detect yearly production trends of sorghum-sudangrass hybrid according to the following equation [27]:
where, v is intercept; Yt is forage production at year t as the dependent variable; For the independent variables, Yt−p is forage production at year lag p, where the lag is the difference in time. For example, if p is 1, Yt−1 indicates forage production at time t−1, which is a year before time t. Xt are weather variables (GDD, accumulated rainfall and sunshine duration) and weather event (heavy rainfall 1, heavy rainfall 2 and typhoon); ϕ and θ are the coefficient of time lags in the autoregressive and moving average parts, respectively; ɛq is the white noise at year lag q, which indicates a sequence of random numbers as an unpredictable error. If the white noise is small, other model is no longer needed; βr is the coefficient of X at year lag r. Firstly, autoregressive of lag p, AR (p) is described by the association between forage production over the past p-year. Secondly, the integration of lag d, as I (d), serves to confirm the stationarity that enable a stable estimation of the forage production trend. Thirdly, the moving average of lag q, as MA (q), is formed from the linear combination of forage production’s white noise. It indicates a simple and common type of smoothing used in forage production forecasting. Finally, X is used to explain the impacts of weather variables and/or events on current forage production, such as GDD, accumulated rainfall, sunshine duration, heavy rainfall 1, heavy rainfall 2 and typhoon. Sorghum-sudangrass forage production trend was described by past year, lag time, and moving average of lag q of ARIMA (p, d, q), and the impacts of weather variables and/or events were estimated using independent variables. Outlier data were eliminated using box-plot for the application of consistent forage production trends.
For optimum ARIMAX model selection, root mean square error (RMSE), R2, maximum absolute percentage error (MAPE) and a Ljung-Box Q test were considered. Higher R2 value and lower error values of RMSE and MAPE were referred for a good fit [28]. Furthermore, the p-value greater than 0.05 of the Ljung-Box Q test was applied to confirm insignificance of residuals from the ARIMA model [33]. The proc ARIMA of SAS 9.4 (SAS Institute, Cary, NC) was used for forage production modeling.
RESULTS
To consider the impact of weather variables and weather events on sorghum-sudangrass hybrid forage production via ARIMAX model, it was necessary to estimate the forage production trend and weather. The GDD, accumulated rainfall and sunshine duration on forage production were significant as independent variables in the regression analysis (Table 1); however, the variance inflation factor of sunshine duration was greater than 4, meaning that this independent variable’s contribution to sorghum-sudangrass forage production was distorted probably due to multicollinearity. Therefore, sunshine duration should be ignored in the model.
Weather variables | Estimate | p-value | Variance inflation factor |
---|---|---|---|
Constant | 2,200.65 | 0.28 | |
Growing degree days | 13.72 | < 0.01 | 3.62 |
Rainfall amount | −7.49 | < 0.01 | 1.88 |
Sunshine duration | −16.82 | < 0.01 | 4.53 |
A yearly trend for the GDD and accumulated rainfall can be seen in Fig. 1A. The trends in forage production and GDD looked similar over the whole time period (1989–2013), especially between 1992–2007 and 2010–2013. Similar trends between forage production and accumulated rainfall were observed in some cases, but conflicting trends were suspected in most. For example, the accumulated rainfall trend in 2002–2009 showing rainfall of below 1,000 mm was similar to the forage production trend. In contrast, when accumulated rainfall was over 1,000 mm, forage production was found to decrease, with the years of 1990, 1995, 1998, 2010, 2011 and 2012 standing out. Therefore, it was hypothesized that the GDD’s trend was associated to the overall forage production trend and the accumulated rainfall trend was also correlated to the forage production trend depending on the specific range of rainfall.
The forage production trend was evaluated by ARIMA to determine the impact of weather events (Table 2). Unlike the GDD, rainfall amount, and sunshine duration, the three weather events did not contribute to the variance of the forage production (p > 0.05). This insignificance of the weather event probably reasoned from insufficient data of the extreme weather events. Similar to the heavy rainfall event, typhoon frequency did not contribute to the modeling of forage production (p = 0.56). The yearly trends of heavy rainfall and typhoon can be seen in Fig. 1B. heavy rainfall 2 was not considered because it was dummy scale. No meaningful association could be found between the trend in forage production influenced by heavy rainfall 1 and typhoon. The long-term field trial data of forage sorghum-sudangrass hybrid demonstrated some production peaks followed a similar pattern with GDD from 1992 to 2007, but demonstrated a reversed pattern with the accumulated rainfall (Fig. 1A).
From each trend in Fig. 1, it was hypothesized that the weather variables and events are related to the forage production trends. Considering the impacts of weather variables and extreme weather events on sorghum-sudangrass forage production, five models were developed (Table 3). The model 1 was developed based on time without considering climate variables and events. The forage production trends considering the impact of GDD or both GDD and accumulated rainfall presented in models 2 and 3, respectively. Inclusions of weather variables did not improve the sorghum-sudangrass forage production models. The estimates of the independent variables were significant only in models 2 and 3. To assess the potential sorghum-sudangrass forage production influenced by weather events, the impacts of heavy rainfall 2 (1: over 1,000 mm during monsoon season, 0: else) and typhoon (1: event occurrence, 0: else) were estimated in models 4 and 5, respectively. In models 2 and 3, estimates of GDD were positive while those of rainfall were negative, which support the yearly production trends appearing in Fig. 1A. Additionally, based on model fitness, other tested models adopting weather variables were not improved from model 1. The weather events during the monsoon season (heavy rainfall 2 in model 4 and typhoon in model 5) had no significant evidence implying potential impact on the current sorghum-sudangrass forage production (p > 0.05). Based on the Ljung-Box Q test, the model fitness of #2, 4 and 5 were inacceptable (p < 0.05). Although all the modeling attempts including GDD, accumulated rainfall, and extreme weather events did not substantially improve the model, a forage production model of sorghum-sudangrass hybrid may consider GDD and accumulated rainfall as presented in model 3 as following:
where Yt−i, GDDt−i and RAt−i are equivalent to forage production, GDD and RA (accumulated rainfall) at t−i time, respectively. i is time order; zero means the current point in time. The forage production trend was a first-order in the AR model with one order of non-seasonal differencing. Even though it looked like AR (2) because of the second-order term (Yt−2). Therefore, GDD (1) and RA (1) had the same order of AR in forage production.
With considering GDD and accumulated rainfall, the normality test demonstrated symmetrical data distribution probability density function of the kernel (Fig. 2A). To check the prediction performance based on the model 3, the expected forage production was plotted against the observed forage production (Fig. 2B). Compared to a diagonal line with a 1:1 incline, the majority of the data points were located within the 95% confidence range. Furthermore, in the residual diagnosis based on the autocorrelation function, the residual was close to the white noise. Therefore, the model 3 was judged as being suitable for forecasting future forage production trends of sorghum-sudangrass hybrid with considering weather variables. Based on the model 3, the forage production was forecasted from 2014 to 2045 (Fig. 3). The forecasted forage production (kg/ha) values were 14,908.06 ± 5,058.72, 14,904.43 ± 5,420.33, 14,909.66 ± 5,754.48, 14,915.13 ± 6,070.15, 14,920.60 ± 6,370.19 and 14,926.09 ± 6,656.72 for the years 2020, 2025, 2030, 2035, 2040 and 2045, respectively. The future trend of forage production was close to the flat shaped.
DISCUSSION
The impacts of weather events (heavy rainfall and typhoon), GDD and accumulated rainfall on sorghum-sudangrass hybrid forage production demonstrated potential in the development of yield prediction model. The extreme weather events tested in this study did not appear as a major sorghum-sudangrass hybrid forage production determination; however, this does not mean the irrelevancy of those weather factors with the forage production. Since the reduction of forage production caused by the weather factors are so erroneous, the accumulation of organized field data is required for future advanced model development. Perhaps, it is necessary to organize forage production damage caused by the type, which can be caused by heavy rainfall and strong winds. The correlation analysis indicated a weak relationship between heavy rainfall 1 and typhoon (Spearman’s r = 0.41*), which has been observed through heavy rainfall associated with typhoon. As previously reviewed, the effect of heavy rainfall and typhoon could not be confirmed through the long-term field trial data because of the insufficient sorghum-sudangrass hybrid cultivation data in the year when extreme summer weather conditions occurred. Extreme weather lasts only a period of time, which may impact crop growth temporarily. The Korean monsoon season lasts 20 to 30 days in general. Forage sorghum-sudangrass hybrid is harvested multiple times, therefore, the lodged sorghum-sudangrass hybrid can be harvested, and later re-growing forage recovers the whole season production [10,34]. Thus, even if the damage is caused by heavy rainfall accompanied by typhoon, it will be difficult to separate the damage from the impact on total forage production. Sorghum is known for strong drought tolerance than corn due to its deeper spreading root systems, and also tolerant to over 12.7 mm/hr rainfall [35]. According to Venuto and Kindiger [12], a single cut was more productive than the two-cut system for the photoperiod- and/or thermosensitive genotype sorghum-sudangrass hybrid. As reviewed through some forage sorghum-sudangrass field tests, single harvest can be efficient than multiple harvests for greater biomass accumulation. However, distribution of production risk during the monsoon season may justify multiple harvests throughout the sorghum-sudangrass growing season.
The effect of temperature was negative in a study that confirmed the forage production trend of sorghum-sudangrass hybrid in Korea [25]. They did not distinguish regions with different climates. It implies that impact of weather variables could be confounded due to poor data homogeneity. The optimum mean temperatures in the growing period for the highest yield are in the range of 23.89°C–26.67°C [36]. The range of mean temperature during growing period was in Suwon (37° 15′ N, 127° 04′ E) was 20.90°C–24.60°C, it indicates that it was at the lower limit of the optimum range. Sorghum-sudangrass hybrid is known to grow well from southern to central region in Korea [37]; however, we confirmed that regional classification is necessary to accurately estimate the impact of weather on forage production without the confounded effect. The mean accumulated rainfall in Suwon (37° 15′ N, 127° 04′ E) between 1988 and 2013 was 998.81 ± 264.90 mm, which was greater than the required rainfall amount (368.3–653.0 mm) for the growth and development of sorghum-sudangrass hybrid [35]. According to Chemere et al. [25], the annual precipitation in Korea was 956.9 mm, and it also found to be excessive for sorghum-sudangrass hybrid. The impact of accumulated rainfall on the forage production was estimated to be negative in both regression and time series modeling, indicating that some negative impact of rainfall on the sorghum-sudangrass hybrid forage production. Therefore, based on significant GDD and accumulated rainfall on forage production, but insignificant impact of the extreme weather events, the rainfall over longer period may be a greater impact factor than short term extreme weather events.
Finally, at any point of the forage production, only the past two years’ forage production was effective, and the impacts of forage production older than three years were not detected. Furthermore, the prediction of forage production with the relationship with GDD and accumulated rainfall, the variation range was small, however, there was a strong correlation between GDD and increasing forage production. In this study, the projected forage production in 2045 considers climate variability rather than weather events. Hence, this implies a stable forage production of sorghum-sudangrass hybrid in the future because of insignificant contribution of extreme weather events on sorghum-sudangrass hybrid, reduction of forage production risk through multi-harvesting is more promising in Korea. Based on the forage production trend considering the trends of GDD and accumulated rainfall (model 3), the forage production of sorghum-sudangrass hybrid in future was forecasted to be unchanged due to the simple pattern (single up-down). Furthermore, the positive impact of GDD and the negative impact of accumulated rainfall were offset in the forage production forecast. In future, it is necessary to interest in the impact of GDD in the long run than the short-term impact of accumulated rainfall. While the annual forage production was increasing in the historical data, the forecasted forage production was close to flat shaped trend. The future trend is determined by the trends of GDD and accumulated rainfall, and it implies that the negative impact of accumulated rainfall in particular was greater than expected. Meanwhile, if temperature of the climate change scenarios such as representative concentration pathways (RCP) scenarios were considered to GDD of model 3 in forecast, the forecasted forage production will be higher due to increasing temperature. However, in the RCP 8.5, the mean temperature is expected to increase by 4.8°C by the end of the 21st century, which is over to the upper limit of appropriate temperature for the growth and development of sorghum-sudangrass hybrid. Therefore, further study on yield forecasting and simulation using various RCP scenarios is currently underway.
In conclusion, the present study supported to hypothesis that the heavy rainfall and typhoons in Korean monsoon seasons was not enough to damage in forage production of sorghum-sudangrass hybrid. However, if the difference in forage production between frequencies exceeds 3,000 kg/ha, the damage caused by heavy rainfall and typhoons will be significant. No clear evidence of damage on annual forage production of sorghum-sudangrass hybrid due to extreme heavy rainfall and typhoons events in short term was confirmed in this study due to insufficient data about multiple harvesting management. Under the imminent climate change, development of harvest management and fertilization strategy for sorghum-sudangrass hybrid is necessary to secure consistent forage production.