RESEARCH ARTICLE

# Impact of abnormal climate events on the production of Italian ryegrass as a season in Korea

Moonju Kim1, Kyungil Sung2,*
1Institute of Animal Resources, Kangwon National University, Chuncheon 24341, Korea
2Department of Animal Industry Convergence, Kangwon National University, Chuncheon 24341, Korea
*Corresponding author: Kyungil Sung, Department of Animal Industry Convergence, Kangwon National University, Chuncheon 24341, Korea. Tel: +82-33-250-8635, E-mail: kisung@kangwon.ac.kr

© Copyright 2021 Korean Society of Animal Science and Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Aug 24, 2020; Revised: Oct 26, 2020; Accepted: Nov 02, 2020

Published Online: Jan 31, 2021

## Abstract

This study aimed to assess the impact of abnormal climate events on the production of Italian ryegrass (IRG), such as autumn low-temperature, severe winter cold and spring droughts in the central inland, southern inland and southern coastal regions. Seasonal climatic variables, including temperature, precipitation, wind speed, relative humidity, and sunshine duration, were used to set the abnormal climate events using principal component analysis, and the abnormal climate events were distinguished from normal using Euclidean-distance cluster analysis. Furthermore, to estimate the impact caused by abnormal climate events, the dry matter yield (DMY) of IRG between abnormal and normal climate events was compared using a t-test with 5% significance level. As a result, the impact to the DMY of IRG by abnormal climate events in the central inland of Korea was significantly large in order of severe winter cold, spring drought, and autumn low-temperature. In the southern inland regions, severe winter cold was also the most serious abnormal event. These results indicate that the severe cold is critical to IRG in inland regions. Meanwhile, in the southern coastal regions, where severe cold weather is rare, the spring drought was the most serious abnormal climate event. In particular, since 2005, the frequency of spring droughts has tended to increase. In consideration of the trend and frequency of spring drought events, it is likely that drought becomes a NEW NORMAL during spring in Korea. This study was carried out to assess the impact of seasonal abnormal climate events on the DMY of IRG, and it can be helpful to make a guideline for its vulnerability.

Keywords: Italian ryegrass; Autumn low-temperature; Severe winter cold; Spring drought; Dry matter yield

## INTRODUCTION

Some abnormal/extreme climate events are becoming more extreme in East Asia, depending on the season, due to global warming and rapid climate change [1]. It has been reported that abnormal climate events causes great damage to agricultural industry and crops [2,3]. In particular, climate events are important for forage crops that can be provided with high quality feed for livestock. Winter forage crops are vulnerable to abnormal climate events during winter, which is significant for survival. Therefore, on the Korean Peninsula, which belongs to the temperate climate zone and has four distinct seasons, the representative abnormal climate events which are fatal and prone to happen to winter crops are as follows: First, low-temperature events, which can occur after seeding in autumn, impede the growth and development of winter crops, making it difficult for them to obtain the energy needed for wintering [4]. According to the report on the climate of the Korean Peninsula [5], the damage caused by low-temperatures in autumn is expected to decrease gradually as global warming increases summer temperatures and extends the summer period. However, low-temperature event can still be fatal due to the increase in frequency and intensity of late monsoon rain (Kaul-Changma) and typhoon. Second, considering the characteristics of Italian ryegrass (Lolium multiflorum Lam., IRG), which has weak cold resistance, there is no doubt that the severe cold in winter is the most critical factor despite its good growth and development in autumn [6]. According to Min et al. [1], greenhouse warming has been assessed to lead to an increase in warm extremes and a decrease in cold extremes, but recent warming associated with Arctic sea ice-melting has tended to bring a cold snap to East Asia in winter. Therefore, the severe cold from warming is likely to become more deadly to IRG. Finally, if the spring drought persists for a long period, it causes a shortage of water resources following the winter dry season that damages the agricultural industry [7]. According to Kim et al. [8], there was a bimodal pattern of frequency between the short and long term for spring drought, which indicates the polarization of drought variation. Hence, as the dry season increases, spring drought will gradually emerge as an important abnormal climate event in the future. For abnormal/extreme climate events, the 90th percentile of temperature [9] and standard deviation [10,11] was mainly used as a standard. Although defining abnormal climate/weather by only considering temperature and precipitation was simple and efficient, the combined effects of several climate variables tend to be ignored. Furthermore, the selection by standard deviation with different multiples, such as one-fold, 1.5-fold, and two-fold, makes it difficult to handle consistent abnormal climate conditions. In this study, not only were temperature and precipitation considered to quantify the abnormal climate events, but also wind speed, humidity, and sunshine duration.

IRG has been generally cultivated in the southern inland and coastal regions of the Korean Peninsula, but this has gradually expanded from the southern to the central regions, and it is now cultivated in most regions except for the northern and central mountainous regions [12]. In general, IRG is vulnerable to the dry and cold conditions of the Korean winter due to lower cold and dry tolerance compared to other winter crops [13]. IRG is cultivated popularly in rotation system with rice, a major food crop, which is helpful to income of farmer in Korea. Furthermore, IRG is important winter forage crop in terms of feed value, such as feed nutritive and preference, on the southern and central regions of Korea [14]. Unfortunately, studies related to the effects of climate change and climate events on IRG have been seldom carried out in Korea, compared to the improvement of varieties and the development of cultivation techniques. Therefore, in this study, the impact of abnormal climate events on the IRG production was considered.

For studies on abnormal/extreme climate events on crop production, it is necessary to accumulate and analyze spatiotemporal climate information recorded in several places over a long period. Historical climate data has been used to assess climatic and soil factors on forage crops with various spatiotemporal characteristics. Climate information collected from the Korea Meteorological Administration (KMA) has been used for various purposes in the Korean agricultural field [15,16]. As part of a complex natural eco-system, the information on climate and soil factors was used to generate the structure of causality contains their direct and/or indirect pathways for forage crops, such as IRG [17,18], forage barley [19] and forage rye [20]. Climate data were also used to identify the trends in yields and climatic factors as a form of climate change research in Korea [2123]. The drawback of performing this study using historical climate information was that it could not reflect abnormal climate conditions because only normal climates conditions because only normal climate conditions were taken into account.

Therefore, this study aimed to assess the impact of seasonal abnormal climate events on IRG production considering various climate variables. Besides this, production trends by year were checked in consideration of significant abnormal climate events in different regions.

## MATERIALS AND METHODS

Data resources, variables and climate scenarios

The IRG metadata (n = 1,107) including location, year, cultivar, seeding and harvesting dates, height and dry matter yield (DMY, kg/ha) were collected from reports of the National Agricultural Cooperative Federation and the Rural Development Administration of Korea from 1988 to 2013 (26 years). In general, abnormal/extreme climates in some regions may be normal in other regions; thus, it is necessary to classify the regions. According to Ko et al. [24], the southern Korean Peninsula was divided into five regions by temperature distribution: northern inland, central inland, southern inland, coastal, and mountainous regions. In this study, the central inland, southern inland and coastal regions were selected because the northern inland and mountainous regions were not suitable for IRG cultivation [18]. Suwon (latitude: 37° 25′ N, longitude: 126° 98′ E, n = 321), Jeonju (latitude: 35° 84′ N, longitude: 127° 12′ E, n = 107) and Jeju (latitude: 33° 51′ N, longitude: 126° 53′ E, n = 135) were selected as representative locations for the central inland, southern inland and southern coastal regions, respectively. These locations are similar to the region classification of weather forecasts based on administrative district (KMA). The cultivar mainly consisted of ‘Hwasan 101’ and ‘Kowinearly’, bred for improved winter survival cold tolerance in the central inland regions, while the southern inland and coastal regions consisted of experiments on various cultivars (Greenfarm and 41 other cultivars).

Based on the location, the raw meteorological data contains daily mean temperature, wind speed, relative humidity, precipitation amount, snowfall amount and sunshine duration were collected by open-application programming interface from the weather information system of KMA (1988–2013). To quantify the climate events for IRG, the season was divided into autumn, winter and next spring based on the earliest seeding date (03/September) and the latest harvesting date (12/June). In Korea, IRG is generally sown between late September and early October and harvested in late May [25]. Thus, growing periods were divided into autumn (September to November), winter (December to February), and next spring (March to May), respectively. To quantify abnormal climate events, the variables were temperature (°C), relative humidity (RH, %), precipitation amount (PA, mm), snowfall amount (SA, cm), sunshine duration (SD, hr) and wind speed (WS, m/s). Where, the temperature was: accumulated temperature (AT, °C) over 5°C in autumn and next spring, and mean temperature (MT, °C) in winter. Where, the lower limitation of AT, 5°C (41°F) was the base temperature of growing degree days for IRG [26].

In general, abnormal climate can be defined by various criteria, including the type, intensity and frequency of climate variables. In this study, abnormal and normal climate scenarios for autumn, winter and next spring were categorized by seasonal climate variables using statistical analysis. Here, a category that satisfies climatic conditions favorable for survival, growth and development for the high IRG yield was defined as a normal climate scenario. Meanwhile, climatic characteristics unfavorable to the survival, growth and development of IRG at a specific season were reflected to the abnormal climate scenarios.

Statistical analysis

In this study, various climate variables were used to generate the abnormal climate events as a season. In general, temperature and precipitation are critical for most abnormal/extreme climate events [9], and wind, humidity and sunshine are also important factors for crops. Unfortunately, since all these variables are highly related to each other, setting events independently with a single variable might be effective, but could not accurately reflect climate conditions. Therefore, in order to generate the seasonal climate events (low temperature in autumn, severe cold in winter and drought in next spring) for aim to compare with normal condition, principal component (PC) analysis was carried out with climate variables as follows [27]:

$\begin{array}{c}{Z}_{1}={\lambda }_{11}{X}_{1}+{\lambda }_{12}{X}_{2}+\cdots +{\lambda }_{1p}{X}_{p}={\Lambda }_{1}^{T}X,\\ ⋮\\ {\stackrel{\to }{Z}}_{q}={\lambda }_{{q}_{1}}\stackrel{\to }{{X}_{1}}+{\lambda }_{{q}_{2}}\stackrel{\to }{{X}_{2}}+\cdots +{\lambda }_{qp}\stackrel{\to }{{X}_{p}}={\Lambda }_{q}^{T}X,\end{array}$

where Zq are climatic PC scores, Xp are climate variables, Λ is parameter of climate variable. p and q are the number of variable and climatic PC, respectively. These PC scores were used to characterize each climate scenarios in this study.

To consider various climate variables, distinct climatic characteristics were extracted via principal PC analysis with an eigenvalue over 1, and Euclidean-distance cluster analysis was used to classify abnormal and normal climate conditions by applying consistent criteria [27]. The year was classified by climatic PC scores using cluster analysis with Ward’s method; then, both were placed simultaneously to generate the climate scenario in scatter plots. Based on the year associated with the classified climate events, the DMY of IRG was compared between the events by t-test with a 5% significant level.

The PC analysis, cluster analysis, and t-test were performed by SPSS 24.0 (IBM, Armonk, NY, USA).

## RESULTS

Characteristics of climate variables by month at regions

For monthly climate variables, the temperature was remarkable between regions (Table 1). In January, the coldest month in Korea, MT only exceeded 5°C in the southern coastal region, while in the central inland and south inland regions, it recorded −2.30 and −0.20, respectively. This indicates that there may be no difficulty with overwintering in the south coastal region despite the low cold tolerance of IRG [28]. Hence, it was hypothesized that severe winter cold was fatal to the DMY of IRG in inland regions rather than coastal regions. After the Korean monsoon season, the PA trend rapidly decreased until winter, regardless of region. In the central inland and southern inland regions, the March PA (mm) was lowest during spring, with 44.20 and 55.60 respectively. On the other hand, the March PA in southern coastal regions tended to be high, at 91.70 mm. This implies that spring drought would be serious in inland regions. WS was generally high in the southern coastal regions, and its fluctuations were somewhat stable except during spring, regardless of region. SA was recorded in the inland regions until March, but in the southern coastal regions, it was concentrated in winter. These regional and seasonal climate variables were used to quantify climate events such as autumn low-temperature, severe winter cold and spring drought.

Table 1. Characteristics of climatic variables by season in central inland, southern inland, and southern coastal regions (1988–2013)
Regions Variables Autumn Winter Spring
9 10 11 12 1 2 3 4 5
Central inland Temperature1) (°C) 630.00 443.30 207.90 0.00 −2.30 0.10 164.30 351.70 536.30
RH (%) 74.00 71.00 68.00 66.00 64.00 63.00 63.00 62.00 66.00
PA (mm) 151.60 55.50 48.90 22.80 21.20 26.60 44.20 69.20 90.70
WS (m/s) 1.60 1.40 1.50 1.50 1.50 1.70 2.00 2.10 1.90
SD (h) 180.50 204.50 159.20 160.00 166.80 173.10 200.50 213.50 224.50
SA (cm) 6.70 16.60 16.00 6.20 5.00 0.40
Southern inland Temperature (°C) 648.00 471.20 255.00 2.10 −0.20 1.80 201.50 384.00 567.30
RH (%) 74.00 70.00 68.00 69.00 68.00 65.00 62.00 59.00 64.00
PA (mm) 133.20 52.30 52.60 33.00 32.30 38.80 55.60 74.10 85.40
WS (m/s) 1.40 1.40 1.50 1.50 1.50 1.60 1.90 2.00 1.80
SD (h) 167.50 197.10 155.20 140.30 149.60 158.20 188.80 211.90 220.80
SA (cm) 2.90 8.60 12.30 6.40 4.50
Southern coast Temperature (°C) 693.00 567.30 390.00 8.10 6.00 6.70 294.50 420.00 558.00
RH (%) 74.00 67.00 65.00 65.00 65.00 65.00 64.00 66.00 70.00
PA (mm) 209.70 78.60 68.60 51.90 64.60 62.90 91.70 88.00 103.10
WS (m/s) 3.00 3.20 3.60 4.30 4.30 4.00 3.80 3.40 3.00
SD (h) 164.30 178.40 125.20 82.00 70.80 104.10 160.60 196.00 212.10
SA (cm) 4.60 5.90 3.60 0.90

1) Temperature: accumulated temperature over 5°C in autumn and spring, mean temperature in winter.

RH, relative humidity; PA, precipitation amount; WS, wind speed; SD, sunshine duration; SA, snowfall amount.

Classification of abnormal climate events by seasons

For the regional and seasonal climate variables, two significant PCs were generated (p < 0.05), except for the southern coastal region in autumn (Table 2). The climatic characteristics of PCs were similar in both inland regions, regardless of the season: In autumn, PC 1 was associated with positive RH, positive PA and negative SD, respectively, and PC 2 with positive AT. In winter, PC 1 and PC 2 were strongly related to RH and WS, positive MT and negative SA, respectively. For spring, PC 1 and PC 2 were related to positive RH and PA, positive WS and negative SD, respectively; however, AT, PA, and SD did not noticeably contribute to only one PC. In the southern coastal regions, only PC 1 was strongly related to negative MT and positive WS in the winter. Finally, PC 1 and PC 2 were associated with positive RH, positive PA and negative SD, negative AT and positive WS in the spring, respectively.

Table 2. Characteristics of principal component score for climate variables to explain seasonal climate scenarios in central inland, southern inland, and southern coastal regions
Seasons Variables Central inland Southern inland Southern coast
PC 1 PC 2 Test PC 1 PC 2 Test PC 1 PC 2 Test
Autumn AT (°C) −0.06 0.91 χ2= 13.19* −0.36 0.83 χ2= 27.58* −0.38 −0.68 χ2= 6.01
RH (%) 0.86 0.22 0.88 0.29 −0.48 0.71
PA (mm) 0.68 −0.42 0.76 −0.35 0.68 0.28
SD (h) −0.67 −0.22 −0.71 −0.44 0.81 −0.14
Winter MT (°C) 0.22 0.76 χ2= 15.51* 0.19 0.79 χ2= 27.67* −0.70 0.55 χ2= 15.41*
RH (%) 0.89 0.08 −0.91 0.09 0.61 0.58
SA (cm) 0.28 −0.69 0.10 −0.78 0.52 0.51
WS m/s) −0.89 0.06 0.94 0.01 0.80 −0.29
Spring AT (°C) −0.36 0.38 χ2= 22.65* −0.28 0.54 χ2= 34.73* −0.29 −0.79 χ2= 23.06*
RH (%) 0.86 −0.27 0.93 −0.13 0.70 0.07
WS (m/s) −0.47 0.70 −0.53 0.71 0.33 0.80
PA (mm) 0.61 0.58 0.66 0.47 0.80 −0.32
SD (h) −0.54 −0.64 −0.51 −0.67 −0.70 0.41

* p < 0.05.

PC, principal component ; AT, accumulated temperature; RH, relative humidity; PA, precipitation amoun; SD, sunshine duration; MT, mean temperature; SA, snowfall amountt; WS, wind speed.

Based on the climatic PCs, the year and climate variables were plotted together to classify abnormal climate events from normal (Figs. 1, 2, and 3). In the central inland regions (Fig. 1), the autumn normal and low-temperature events could be divided into the second quadrant close to SD of PC 1 and AT of PC 2, and a fourth quadrant close to RH and PA of PC 1 and far from AT of PC 2, respectively. For winter, normal and severe cold events were placed in the first quadrant, where RH of PC 1 and MT of PC 2 were located, and the third quadrant, where WS of PC 1 and SA of PC 2 were located, respectively. For spring, normal and drought events were distributed on a right quadrant where RH and PA of PC 1 were mainly arranged, and a left quadrant where WS, SD, and AT of PC 2 were arranged, respectively. The characteristics of climatic PCs of the southern inland region were similar to those of the central inland region (Fig. 2). However, in winter, normal and severe cold events were located in the second and fourth quadrants, respectively; because the sign of climate variables in PC 1 was opposite between central and southern inland regions. In the southern coastal regions (Fig. 3), autumnal climate events were not classified as normal and low-temperature events. In the case of winter, it was difficult to distinguish severe cold from normal because of the ambiguous position of climate variables, although PCs were statistically distinct. The spring droughts were prominent in the third quadrant related to high AT, high SD and low PA.

Fig. 1. Scatter plots of year and climate variables to classify climate events in the central inland region. (A) autumn, (B) winter, (C) spring.
Fig. 2. Scatter plots of year and climate variables to classify climate events in the southern inland region. (A) autumn, (B) winter, (C) spring.
Fig. 3. Scatter plots of year and climate variables to classify climate events in the southern coastal region. (A) winter, (B) spring.

Years were also classified, based on the position of climate events. Years where the classification was not clear for each position were divided using individual cluster analysis based on the major climate variable related to PC 1 and PC 2. For example of the central inland region, 1988, 2002 and 1995, located on the third quadrant in autumn, were classified as low-temperature events using cluster analysis taking into consideration the PA of PC 1 and AT of PC 2. After classifying the year into each climate event, the impact of the climate event on IRG production was assessed.

Comparison of climate variables and dry matter yield between climate events

To assess the impact of abnormal climate events on IRG production, DMY was compared by t-test (Table 3). In the central inland, the difference in DMY was greater in order of severe winter cold (2,620 kg/ha), spring drought (2,048 kg/ha) and autumn low-temperature (1,452 kg/ha). In the southern inland region, the only severe winter cold reduced DMY by 4,127 kg/ha (p < 0.05). It was not expected that the damage was greater in the southern inland region than in the central inland region, where the MT was lower. In the southern coastal region, the low-temperature event was not outside of the normal range, and the severe winter cold event was not effective (p = 0.61). The difference in DMY was 2,061 kg/ha between normal and spring drought (p < 0.05).

Table 3. Comparison of dry matter yield between climate events for Italian ryegrass based on regions and seasons
Regions Seasons Events N Dry matter yield (kg/ha) t-test
Central inland Autumn Normal 162 10,698.40 ± 310.14 3.93*
Low-temperature 159 9,246.88 ± 201.20
Winter Normal 298 10,153.57 ± 195.39 3.66*
Severe cold 23 7,533.83 ± 478.68
Spring Normal 194 10,776.06 ± 253.86 5.79*
Drought 127 8,728.23 ± 237.99
Southern inland Autumn Normal 85 9,170.87 ± 361.30 −1.65 (p = 0.10)
Low-temperature 22 8,334.09 ± 355.37
Winter Normal 93 9,538.85 ± 291.65 −5.23*
Severe cold 14 5,411.50 ± 609.22
Spring Normal 70 9,140.99 ± 422.63 −0.78 (p = 0.44)
Drought 37 8,729.86 ± 319.28
Southern coast Autumn Normal Not classified
Low-temperature
Winter Normal 127 14,403.50 ± 1,769.99 0.53 (p = 0.61)
Severe cold 8 13,424.76 ± 483.68
Spring Normal 52 14,749.92 ± 515.27 5.89*
Drought 83 12.688.87 ± 673.40

* p<0.05.

The impacts of significant abnormal climate events on DMY of IRG were described by identifying differences in climate variables between normal and abnormal climate events as shown in Tables 4, 5, and 6. In the central inland regions (Table 4), the PA of the autumn low-temperature event was twice as high as the PA of the normal event (p < 0.05). In winter, the differences in RH, WS, and MT were clear between severe cold and normal events (p < 0.05). Finally, between normal and spring drought, the differences in RH and PA were noticeable (p < 0.05). In the southern inland regions (Table 5), only the severe winter cold event was effective, and the differences in MT, RH and SA were important between winter severe cold and normal events (p < 0.05). In the southern coastal regions (Table 6), for the significant spring drought event, the effective climate variables were AT, PA, and SD (p < 0.05). It was concluded that the difference in these climate variables led to the difference in the DMY of IRG between normal and abnormal climate events.

Table 4. Comparison of climate variables between climate events in central inland region
Seasons Variables Climate events t-test
Normal Abnormal1)
Autumn AT (°C) 1,327.68 ± 20.55 1,256.46 ± 26.53 2.15*
RH (%) 70.29 ± 1.01 71.56 ± 1.23 −0.80 (p = 0.43)
PA (mm) 56.85 ± 6.07 110.22 ± 12.56 −3.66*
SD (h) 193.16 ± 2.99 171.13 ± 3.49 4.79*
Winter MT (°C) −0.55 ± 0.21 °C1.60 ± 0.56 2.17*
RH (%) 66.00 ± 0.77 59.29 ± 1.90 3.28*
SA (cm) 8.40 ± 1.06 11.84 ± 2.27 −1.37 (p = 0.20)
WS (m/s) 1.56 ± 0.03 1.77 ± 0.09 −2.30*
Spring AT (°C) 1,074.10 ± 16.53 1,026.46 ± 22.96 1.68 (p = 0.10)
RH (%) 67.14 ± 0.92 61.13 ± 0.94 4.55*
WS (m/s) 1.91 ± 0.03 2.04 ± 0.07 −1.85 (p = 0.08)
PA (mm) 82.84 ± 4.48 55.10 ± 4.76 4.25*
SD (h) 219.77 ± 4.75 204.86 ± 4.88 2.19*

1) Abnormal: low-temperature in autumn, severe cold in winter, drought in spring.

* p < 0.05.

AT, accumulated temperature; RH, relative humidity; PA, precipitation amoun; SD, sunshine duration; MT, mean temperature; SA, snowfall amountt; WS, wind speed.

Table 5. Comparison of climate variables between climate events in southern inland region
Seasons Variables Climate events t-test
Normal Abnormal1)
Autumn Non-significant difference in DMY
Winter MT (°C) 1.32 ± 0.18 −0.65 ± 0.05 2.77*
RH (%) 67.28 ± 0.58 62.00 ± 1.00 2.33*
SA (cm) 12.26 ± 0.93 22.25 ± 3.45 −2.73*
WS (m/s) 1.49 ± 0.05 1.85 ± 0.15 −1.89 (p = 0.07)
Spring Non-significant difference in DMY

1) Abnormal: low-temperature in autumn, severe cold in winter, drought in spring.

* p < 0.05.

DMY, dry matter yield; MT, mean temperature; RH, relative humidity; SA, snowfall amountt; WS, wind speed.

Table 6. Comparison of climate variables between climate events in southern coastal region
Seasons Variables Climate events t-test
Normal Abnormal1)
Autumn Not classified
Winter Non-significant difference in DMY
Spring AT (°C) 1,255.56 ± 12.53 1,295.53 ± 18.77 −1.84*
RH (%) 68.09 ± 1.47 66.00 ± 1.36 1.04 (p = 0.31)
WS (m/s) 3.27 ± 0.10 3.41 ± 0.08 1.13 (p = 0.27)
PA (mm) 128.87 ± 10.24 74.25 ± 5.01 5.38*
SD (h) 175.40 ± 4.67 197.74 ± 2.71 −4.14*

1) Abnormal: low-temperature in autumn, severe cold in winter, drought in spring,

* p < 0.05

DMY, dry matter yield; AT, accumulated temperature; RH, relative humidity; WS, wind speed; PA, precipitation amoun; SD, sunshine duration.

The trend of abnormal climate events with climatic characteristics

The DMY trends of IRG by year with the happening of significant abnormal climate events were shown in Fig. 4, where the climatic characteristics were shown with heat-maps of PC scores to explain the intensity of abnormal climate events. First, in the case of the autumn low-temperature event in the central inland region (Fig. 4A), since 2000, the colors of PC 1 and 2 had tended to change to dark (low) and light (high), respectively. This indicates that the tendency of the low-temperature event was gradually weakening. The severe winter cold events in 1988, 1993, 1996, 1997, 2000, and 2005 were noticeable based on PC 1 related to WS, but PC 2 related to MT has been higher since 2000. Finally, for the spring drought, there were no specific patterns in PC 1, but the color of PC 2 was lighter since 2000. In the southern inland region (Fig. 4B), the high PC 1 in 2005 and low PC 2 in 2005, 2011, and 2012 were clear, indicating that the severe winter cold event based on WS and SA was serious. For the spring drought in the southern coastal region (Fig. 4C), the high PC 2 with low PC 1 stood out since 2005, which indicates low PA and strong WS. In these heat maps, the DMY fluctuation tended to change rapidly with the greater gap of color by year.

Fig. 4. Trend of dry matter yield of Italian ryegrass by year with significant abnormal climate event with heat-map of principal component scores. (A) central inland region, (B) southern inland region, (C) southern coastal region. PC, principal component.

## DISCUSSION

In this study, on abnormal climate events affecting the DMY of IRG, the abnormal climate events of severe winter cold and spring drought were considered in central inland, southern inland and southern coastal regions of Korea. All abnormal events were critical in the central inland region, while severe winter cold and spring drought events were fatal in the southern inland and coastal regions, respectively. The impact of severe winter cold event was greater in the southern inland region than in the central inland region, despite of lower winter MT in the central inland region. It concluded that despite frequent severe cold events in winter, yield damage was lower due to the cultivation of relatively strong cold-resistant varieties in the central inland region. In general, it is clear that wintering is difficult in the central inland region than in the southern inland region due to high temperature and humidity. Therefore, for the low yields in the southern inland region, it was concluded that other factors, such as soil physical properties and cultivation management for wintering, also contributed to the survival of IRG. According to Kim and Sung [18], IRG yield in paddy fields was lower than in upland fields due to the short growing period by the rice-IRG rotation system. Therefore, climate is the most important factor, but other factors in the severe cold below zero also affected the yield difference. Furthermore, the development and selection of new varieties have been intensively performed to expand to the central and northern regions of Korea [29]. Between normal and severe cold events in winter, the difference in MT was 1.05°C in the central inland region and 1.97°C in the southern inland region. This implies that not only frequency and intensity but the climatic fluctuation was important between normal and abnormal events. In both inland regions, increasing PC 2 related to winter MT was clear since 2000, but it was less important than PC 1 related to WS and RH. It is certain that the winters are getting warmer, but there is no scientific evidence that these changes lead to specific patterns for a severe cold event. Rather, it is likely that WS and RH were more involved in the severe winter cold event. Park et al. [30] reported a similar result that, in Korea, cold days tended to decrease with year, but there was no significant trend in cold surges.

In this study, the spring drought was fatal in central inland and southern coastal regions. In particular, the frequency of spring drought events in southern coastal regions was greater than normal; thus, the spring drought is abnormal in terms of intensity and production impact, but in terms of frequency, it seems to be a NEW NORMAL in Korea since 2000. There is no doubt that RH and PA were such important climate variables for the spring drought event that they were extracted on PC 1 regardless of region. Physically, the cooling of the air has the liquidizing properties of water vapor, so in inland regions where it is easily dried compared to coastal regions where the specific heat is low and sufficient moisture is maintained, the low RH in winter causes a severe cold [31]. The PA of spring drought in the southern coastal region was half as normal. For the spring drought in inland regions of the Korean Peninsula, spring PA in 2000 and 2001 was less than 50% of the average year, whereas spring PA was unusually high in 2003 [7,32]. The spring PA in the southern coastal regions is only 16 percent of the annual PA [33], and the Korean rainfall concentrated in the summer monsoon season will likely become even more severe, increasing the chances of a long dry season in spring [23]. As to the severity of the spring drought, Lee et al. [34] reported that if there was no rain for more than two weeks, concerns about agricultural drought were particularly high. However, it also reported that the spring drought response had improved and the damage had gradually decreased as the agricultural water facilities improved. Furthermore, the southern region was more susceptible to drought than the central region of the Korean Peninsula. Therefore, storage-oriented water management will be very important in the face of abnormal drought, such as regional torrential rainfall, increase in drought variability and incidence, and prolonged drought due to climate change. Finally, yield damage was the least for the low-temperature event in autumn, and this had significant impact only in the central inland region. However, given the noticeable trend of autumn low-temperature events in the central inland region since 2000, it was determined that there is a high risk of such events becoming prominent as a serious abnormal climate event in the future.

The information on climate events generated through this study will help to predict the production of winter crops under the climate event happening, and it can be expected to help to provide a guideline for cultivating IRG in abnormal climate events and to assess its vulnerability. Unfortunately, because of the small sample size, it is not possible to separate the impacts of abnormal climate events from two or more combinations of climate events by considering each seasonal climate event. Furthermore, the climate events were investigated based on monthly records without considering the growing stage of IRG; thus, it is difficult to reflect the duration and occurrence timing of climate events for the growth stages of winter crops. In the future, climate variables for the seasonal duration of abnormal climate events using daily records will be investigated to reflect the growth stage of IRG.

## CONCLUSION

The study was carried out to assess the impact of abnormal climate events such as autumn low-temperature, severe winter cold and spring drought on the production of Italian ryegrass in the central inland, southern inland, and southern coastal regions of Korea. The yield damage caused by severe winter cold event was serious in both inland regions; in particular, the damage was more fatal in the southern inland region, even though the frequency of events was lower. In consideration of the trend and frequency of spring drought events, it is likely that drought becomes a NEW NORMAL during spring in Korea. Finally, the damage caused by low-temperature events in autumn was lower than other abnormal events, although its tendency had changed since 2000. The greater the difference in climate characteristics by year, the more the production fluctuations tended to change rapidly. Therefore, for the abnormal climate events related to winter forage crops in Korea, a strategic approach focusing on prediction and quick judgment is required, since we cannot be routinely prepared for the characteristics of rare events.

## Competing interests

No potential conflict of interest relevant to this article was reported.

Not applicable.

## Acknowledgements

This study was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Science and ICT (NRF-2020R1C1C1004618).

## Availability of data and material

Upon reasonable request, the datasets of this study can be available from the corresponding author.

## Authors’ contributions

Conceptualization: Kim M.

Data curation: Kim M, Sung K.

Formal analysis: Kim M.

Methodology: Kim M.

Software: Kim M.

Validation: Sung K.

Investigation: Kim M.

Writing - original draft: Kim M.

Writing - review & editing: Sung K.

## Ethics approval and consent to participate

This article does not require IRB/IACUC approval because there are no human and animal participants.

## REFERENCES

1.

Min SK, Son SW, Seo KH, Kug JS, An SI, Choi YS, et al. Changes in weather and climate extremes over Korea and possible causes: a review. Asia Pac J Atmos Sci. 2015; 51:103-21

2.

Pachauri RK, Reisinger A. Climate change 2007: synthesis report. Contribution of working groups I, II and III to the fourth assessment report of the Intergovernmental Panel on Climate Change. Geneva, Switzerland: The Intergovernmental Panel on Climate Change [IPCC]. 2007; p p. 104

3.

Rosenzweig C, Iglesius A, Yang XB, Epstein PR, Chivian E. Climate change and extreme weather events-Implications for food production, plant diseases, and pests. Glob Chang Hum Health. 2:90-104

4.

Cooper JP. Short-day and low-temperature induction in Lolium. Ann Bot. 1960; 24:232-46

5.

Pachauri RK, Allen MR, Barros VR, Broome J, Cramer W, Christ R, et al. Climate change 2014: synthesis report. Contribution of working groups I, II and III to the fifth assessment report of the Intergovernmental Panel on Climate Change. Geneva, Switzerland: The Intergovernmental Panel on Climate Change [IPCC]. 2014; p p. 151

6.

Kobayashi H, Takahashi Y, Matsumoto T, Nishiguchi Y. Changes in nutritive value of italian ryegrass (Lolium multiflorum Lam.) during overwintering period. Plant prod sci. 2008; 11:228-31

7.

Kim S, Park CK, Kim MK. The regime shift of the northern hemispheric circulation responsible for the spring drought in Korea. Asia Pac J Atmos Sci. 2005; 41:571-85

8.

Kim DW, Byun HR, Choi KS, Oh SB. A spatiotemporal analysis of historical droughts in Korea. J Appl Meteorol Climatol. 2011; 50:1895-912

9.

Choi Y. Trends on temperature and precipitation extreme events in Korea. J Korean Geogr Soc. 2004; 39:711-21

10.

Sim GM, Kim YS, Jung MP, Kim JW, Park MS, Hong SH, et al. Recent changes in the frequency of occurrence of extreme weather events in South Korea. J Clim Change Res. 2018; 9:461-70

11.

World Meteorological Organization [WMO]. WMO guidelines on the calculation of climate normals. Geneva, Switzerland: WMO. 2017

12.

Choi GJ, Ji HC, Kim KY, Park HS, Seo S, Lee KW, et al. Growth characteristics and production of cold-tolerant “Kowinearly” Italian ryegrass in the northern part of South Korea. Afric J Biotech. 2011; 10:2676-82

13.

United States Department of Agriculture Natural Resources Conservation Service [USDANRCS]. Plant fact sheet: Italian ryegrass. Washington DC: USDANRCS. 2002

14.

Ko KH, Kim JK, Kim JD. Evaluation of agronomic characteristics and forage production of domestic and foreign italian ryegrass cultivar in Korea. J Korean Soc Grassl Forage Sci. 2015; 35:297-302

15.

Peng J, Kim M, Kim Y, Jo M, Kim B, Sung K, et al. Constructing Italian ryegrass yield prediction model based on climatic data by locations in South Korea. Grassl sci. 2017; 63:184-95

16.

Kim MJ, Peng JL, Kim BW, Kim YJ, Sung KI. Research on causal relationship between forage yield and various factors (climate and soil) as a way of forage yield prediction. In In: Proceeding of 6th Korea-Japan-China Grassland Conference; 2016 Aug 17-20, Jeju, Korea

17.

Kim M, Sung KI, Kim YJ. Analysis of climate effects on Italian ryegrass yield via structural equation model. Korean J Appl Stat. 2014; 27:1187-96

18.

Kim M, Sung K. Comparison of causality of temperature and precipitation on Italian ryegrass (Lolium multiflorum Lam.) yield between cultivation fields via multi-group structural equation model analysis in the Republic of Korea. Agriculture. 2019; 9:254

19.

Kim M, Jeon M, Sung KI, Kim YJ. Bayesian structural equation modeling for analysis of climate effect on whole crop barley yield. Korean J Appl Stat. 2016; 29:331-44

20.

Kim M, Peng J, Sung K. Causality of climate and soil factors affecting whole crop rye (Secale cereale L.) yield as part of natural ecosystem structure via longitudinal structural equation model in the Republic of Korea. Grassl Sci. 2020; 66:110-5

21.

Chemere B, Kim J, Lee B, Kim M, Kim B, Sung K. Detecting long-term dry matter yield trend of sorghum-sudangrass hybrid and climatic factors using time series analysis in the Republic of Korea. Agriculture. 2018; 8:197

22.

Chemere B, Kim M, Peng J, Kim B, Sung K. Detecting dry matter yield trend of whole crop maize considering the climatic factors in the Republic of Korea. Grassl Sci. 2019; 65:116-24

23.

Kim M, Chemere B, Sung K. Effect of heavy rainfall events on the dry matter yield trend of whole crop maize (Zea mays L.). Agriculture. 2019; 9:75

24.

Ko JW, Baek HJ, Kwon WT, Park JY. The characteristics of spatial distribution of temperature and regionalization in Korea. J clim res. 2006; 1:3-14

25.

Jang WG. Manual of new varieties and cultivation utilization for Italian ryegrass. Rural Development Administration [RDA]. Cheonan, Korea: Sangloksa. 2011

26.

Strand E. Base temperature levels in heat sum calculations. Acta Agric Scand. 1987; 37:279-86

27.

Johnson RA, Wichern DW. Applied multivariate statistical analysis. Upper Saddle River, NJ: Prentice hall. 2002

28.

Bélanger G, Rochette P, Castonguay Y, Bootsma A, Mongrain D, Ryan DAJ. Climate change and winter survival of perennial forage crops in eastern Canada. Agron J. 2002; 94:1120-30

29.

Kim KY, Lee SH, Ji HC, Hwang TY, Lee KW, Park SM. Growth characteristics and yield comparisons of domestic and foreign Italian ryegrass (Lolium multiflorum Lam.) varieties in Yecheon area of Korea. J Korean Soc Grassl Forage Sci. 2015; 35:43-9

30.

Park TW, Ho CH, Jeong SJ, Choi YS, Park SK, Song CK. Different characteristics of cold day and cold surge frequency over East Asia in a global warming situation. J Geophys Res Atmos. 2011; :116

31.

Elovitz KM. Understanding what humidity does and why. ASHRAE J. 1999; 41:84

32.

Jang MW, Yoo SH, Choi JY. Analysis of spring drought using NOAA/AVHRR NDVI for North Korea. J Korean Soc Agric Eng. 2007; 49:21-33

33.

Hong KO, Suh MS, Rha DK. Temporal and spatial variations of precipitation in South Korea for recent 30 years (1976–2005) and geographic environments. J Korean earth sci soc. 2006; 27:433-49

34.

Lee JH, Seo JW, Kim CJ. Analysis on trends, periodicities and frequencies of Korean drought using drought indices. J Korean Water Res Assoc. 2012; 45:75-89

##### 2019 JCR Impact Factor: 1.685

The 2019 Journal Citation Reports (JCR) was announced, and the impact factor of JAST was determined to be 1.685.

We would like to ask for your continued interest and support in our journal.

Thank you.

JAST Editorial Office

I don't want to open this window for a day.