Extended-Range Forecast of Regional Persistent Extreme Cold Events Based on Deep Learning

: Regional persistent extreme cold events are meteorological disasters that cause serious harm to people’s lives and production; however, they are very difﬁcult to predict. Low-temperature weather systems and their effects have a signiﬁcant low-frequency oscillation period (10–20 d and 30–60 d). This paper uses deep learning to analyze the extended-range time scale and predict regional persistent extreme cold events. The dominant low-frequency oscillation components of cold events are obtained via wavelet transform and Butterworth ﬁltering. The low-frequency oscillation component is decomposed via empirical orthogonal function decomposition to extract the main spatial mode and time coefﬁcient. A convolutional neural network is used to establish the correlation between large-scale circulations and the time coefﬁcient of the low-frequency oscillation component of the lowest temperature. The proposed deep learning model exhibits good prediction accuracy for regional persistent extreme cold events with low-frequency oscillations.


Introduction
With the intensification of global warming, persistent extreme cold events are occurring more frequently [1][2][3].The persistent cold rain and snow disasters in southern China in 2008 [4], persistent cold events in northern China in 2009, and continuous extreme cold weather in Eurasia in 2012 have caused substantial harm to transportation and people's production and life.
Cold weather events can be categorized as cold air, cold wave, and extreme cold events.The National Climate Center of the China Meteorological Administration (CMA) uses the following method to define extreme cold events.The lowest temperature of a specific day and the temperatures 5 days before and after that day are related.The temperatures are sorted in ascending order, and the 10th percentile is considered the extreme threshold of the station.An extreme cold event is characterized by a temperature lower than the extreme threshold in the same period of the year within the normal climate period.A cold wave event is defined as a 24 h drop in the daily minimum temperature greater than or equal to 8 • C, a 48 h drop in temperature greater than or equal to 10 • C, or a 72 h drop in temperature greater than or equal to 12 • C and a daily minimum temperature drop below 4 • C. Regional persistent extreme cold events (RPECEs) are primarily cold wave events.According to the statistics of the National Climate Center, most cold wave events (approximately 89%) are not RPECEs because the minimum temperature of a cold wave event may be greater than the extreme threshold in the same period of the year; thus, a cold wave does not constitute a persistent extreme cold event.
Cold waves occur more frequently in winter.It is generally believed that they are caused by the collapse of the high-pressure area in the Ural region or the re-establishment of the East Asian trough, causing cold air to invade northern China from the southward low-pressure trough [5].The key circulation system that causes different RPECEs compared to cold waves is the large oblique ridge over the Eurasian continent [6] that extends from the Caspian Sea to the Bering Strait.A cold wave occurs when the high-latitude two-wave structure is transformed into a three-wave structure, and the wave train of the Okhotsk Sea declines; the duration is usually 2-3 d.However, the oblique ridge does not form three waves but slowly changes to a larger scale; therefore, its duration is longer [7].
On the other hand, the slow change in large-scale oblique ridges with longer time scales corresponding to RPECEs is the reason for the predictability of the extended range of these events.Therefore, the changes in the circulation of the lower tropospheric atmosphere can be analyzed dynamically to predict RPECEs.Although the cold wave formation can be theoretically predicted similarly, the predictable upper limit is often only 7-8 d, much shorter than the extended-range forecast.Therefore, many cold waves are difficult to predict for longer periods of time [8].
Existing studies on the extended period prediction of temperature and precipitation primarily used statistical or machine learning methods.Yang [9,10] used a principal component-complex autoregressive model and an extended complex autoregressive model to establish a relationship between V850 hPa, T850 hPa, and the forecast quantity for the daily low-frequency component of precipitation and average temperature in the lower reaches of the Yangtze River.Hsu et al. [11,12] established a space-time projection method (STPM) model, based on which Zhu et al. [13] used augmented singular value decomposition (SVD) to predict the extended period of winter surface air temperature and extreme cold weather.They used seven predictors (T850, temperature advection, sea level pressure, 500 hPa vertical motion, H200, SH700, and the time coefficient of the surface air temperature space mode) and calculated their weighted average to obtain one predictor.Lee et al. [14] used historical data to estimate the probability of extremely heavy precipitation for different phases in different months of the Boreal summer intraseasonal oscillation (BSISO) index in the Asian summer monsoon region.Lei et al. [15] utilized a convolutional neural network (CNN) to obtain the future hourly time coefficients of the previous forecast factors (H200, H500, and outward long-wave radiation (OLR)).They used the first K spatial modes of the daily maximum temperature data and synthesized the spatial modes and time coefficient.Wu et al. [16] developed a dynamical-statistical forecasting model (DPSM) to forecast precipitation in an extended period.Augmented SVD was used to establish an equation describing the relationship between the SVD time coefficient and the predicted measurements for the predictors of different historical pentads (OLR, U850, U200, H850, H500, H200, SH700, and SLP).The dynamic and SVD modes were coupled, and the projected time coefficient was used in the linear equation to obtain the predictions.Zhang et al. [17] used SVD to obtain a coupled spatial mode, the time coefficient of the precipitation field, and the predictors of tropical OLR and mid-high latitude (H500).They replaced the time coefficient of the precipitation field with that of the predictor.The subseasonal to seasonal S2S model forecast field of the predictor was projected onto the coupling mode of the actual forecast factor.The result was substituted into the trained model to obtain the predicted precipitation field.In addition, many scholars have proposed methods using dynamical statistics to predict meteorological elements for an extended range [18][19][20].
It should be pointed out here that predicting meteorological factors, such as precipitation and temperature (maximum temperature and minimum temperature), requires statistical models for averaging or fusion.Commonly used statistical methods include linear regression, probability models, or neural networks.However, relatively few studies performed extended-range forecasts for RPECEs.The objectives of this study are to determine the probability of RPECEs using deep learning methods, develop a method for analyzing them, and use deep learning to create extended-range forecasts for RPECEs.

Materials
In this paper, winter refers to October to March; for example, the winter in 1995/1996 is denoted as 1995/1996.The JRA55 reanalysis data of the Japan Meteorological Agency were used.These data are global daily average grid data with a horizontal resolution of 1.25 • × 1.25 • and a temporal range of 1960/1961-2019/2020.The sea-level pressure field, 850 hPa zonal wind field, and meridional wind field were selected.Other meteorological datasets included NOAA's OLR data with a range of 1979/1980-2019/2020 and the daily dataset of meteorological elements of national ground weather stations from 1960/1961 to 2019/2020 in southern China (latitude below 36 • N and longitude below 108 • E) (Figure 1).Following Chen et al. [21], the research area was divided into an area north of the Yangtze River (blue dots in Figure 1) and another south of the Yangtze River (pink dots in Figure 1), with a dividing line at 30 Atmosphere 2023, 14, x FOR PEER REVIEW 3 of 11

Materials
In this paper, winter refers to October to March; for example, the winter in 1995/1996 is denoted as 1995/1996.The JRA55 reanalysis data of the Japan Meteorological Agency were used.These data are global daily average grid data with a horizontal resolution of 1.25° × 1.25° and a temporal range of 1960/1961-2019/2020.The sea-level pressure field, 850 hPa zonal wind field, and meridional wind field were selected.Other meteorological datasets included NOAA's OLR data with a range of 1979/1980-2019/2020 and the daily dataset of meteorological elements of national ground weather stations from 1960/1961 to 2019/2020 in southern China (latitude below 36° N and longitude below 108° E) (Figure 1).Following Chen et al. [21], the research area was divided into an area north of the Yangtze River (blue dots in Figure 1) and another south of the Yangtze River (pink dots in

Methods
The RPECEs were extracted as follows.First, the extreme threshold was determined.We selected the 10th percentile of the lowest temperatures for each day within a 10 d window (5 before and after) from 1960/1961 to 2019/2020, resulting in 682 d (62 years × 11 d).The criteria for stations with persistent low-temperature events were as follows: the lowest temperature remained below the extreme threshold continuously for 5 d, and the temperature drop during this period exceeded 5 °C [22].The following conditions had to be met for RPECEs: (1) The 10 neighboring stations of 1 station recorded persistent low-temperature events for at least 5 d.(2) The number of stations meeting condition 1 exceeded 20 on at least one day during this period [21].
EOF maximizes the concentration of most spatiotemporal variables based on the principal components of a few independent variables.It extracts the main spatial modes and their time coefficients.In this study, the EOF decomposition of the low-frequency oscillation component of the 2 m minimum temperature in the training set was conducted from 1979 to 2020 to obtain the first spatial mode, which was used as the main feature to characterize the minimum temperature oscillation in the winter.A projection formula was used to project the test set onto the first spatial mode to obtain the time coefficient.The reason was to prevent using information from the test set in the training set.

Methods
The RPECEs were extracted as follows.First, the extreme threshold was determined.We selected the 10th percentile of the lowest temperatures for each day within a 10 d window (5 before and after) from 1960/1961 to 2019/2020, resulting in 682 d (62 years × 11 d).The criteria for stations with persistent low-temperature events were as follows: the lowest temperature remained below the extreme threshold continuously for 5 d, and the temperature drop during this period exceeded 5 • C [22].The following conditions had to be met for RPECEs: (1) The 10 neighboring stations of 1 station recorded persistent low-temperature events for at least 5 d.(2) The number of stations meeting condition 1 exceeded 20 on at least one day during this period [21].
EOF maximizes the concentration of most spatiotemporal variables based on the principal components of a few independent variables.It extracts the main spatial modes and their time coefficients.In this study, the EOF decomposition of the low-frequency oscillation component of the 2 m minimum temperature in the training set was conducted from 1979 to 2020 to obtain the first spatial mode, which was used as the main feature to characterize the minimum temperature oscillation in the winter.A projection formula was used to project the test set onto the first spatial mode to obtain the time coefficient.The reason was to prevent using information from the test set in the training set.
We employed two signal processing methods.The first one was the Butterworth bandpass filter, which separates the climate cycle components from the observation sequences of stations or grid points.The upper and lower cut-off angle frequencies of the bandpass were used to separate the climate signal in a specific frequency interval from the original signal.The second method was wavelet analysis, which utilizes a convolution kernel to transform the time domain signal of the time series x(t) into the frequency domain ω to decompose the x corresponding to t at different times, obtaining the signal to be analyzed.Morlet wavelet analysis was used to transform the minimum temperature time series data to obtain the wavelet power spectrum corresponding to the period and time in order to determine the statistically significant low-frequency change period.
An FCNN uses a fully connected layer to integrate the inputs from different dimensions; however, it has shortcomings.First, when the data have high dimensionality, the computational efficiency of the FCNN is relatively low.Second, the direct processing of two-dimensional data by the FCNN results in information loss.The CNN was first proposed in 1980 [24].Due to its structural characteristics, such as local area connection, weight sharing, and downsampling, it is highly suitable for signal processing and extraction from 2D/3D data; and for grid data analysis, such as reanalysis data and modeling.In this study, a CNN was used to extract large-scale signal information (42.5 • -150 • E, 0 • -67.5 • N) from the 850 hPa, 500 hPa geopotential height field, and OLR field.The extracted features were used to train the FCNN to predict low-frequency oscillation components.

Signal Processing of Intraseasonal Oscillation and RPECEs
According to the definition of RPECEs in this paper, we sorted the data by the start and end dates and obtained the duration (days) and peak dates of RPECEs in Table 1.We analyzed the 850 hPa temperature field, 500 hPa geopotential height field, 850 hPa wind field, and sea-level pressure field.As shown in Figure 2, the low-level warm anomaly in Siberia corresponds to the positive geopotential anomaly in the middle and upper atmosphere, and the low-level cold anomaly in southern China corresponds to the negative geopotential anomaly in the middle and upper atmosphere.The result shows that we can use the geopotential height at 500 hPa in the middle layer and the geopotential height at 200 hPa in the upper layer to indicate changes in the atmospheric temperature in the lower layer.Figure 3 shows that the extremely low temperatures in southern China were caused by the significant northerly wind component extending southward to southern China via an anticyclone.After removing the annual cycle and annual average from the daily minimum temperature series from 1960 to 2020, Morlet wavelet analysis was performed.Figures 4b and 5b show the wavelet power spectra.The spectral values passing the 0.05 significance level are surrounded by thick solid lines, indicating that the period is significant.Figures 4c and 5c show the global wavelet power spectrum.The dotted line represents the critical value of the 0.05 significance level in the cycle.The spectral value on the right side of the dotted line indicates a 0.05 significance level.Figure 4 shows a low-frequency change period of about 20 d in the low-temperature events in the northern Yangtze River region, and Figure 5      After removing the annual cycle and annual average from the daily minimum perature series from 1995 to 2020, the Butterworth filtering of the low-frequency o tion series was conducted after station-by-station filtering for 10-20 d and 30-60 d decomposition was performed on the low-frequency oscillation series at multiple st As shown in Figure 6, the explained variances of the first mode of the minimum te ature on the two time scales (10-20 d and 30-60 d) for the northern Yangtze River were 78.9% and 80%, respectively.The physical meaning of the first mode is that the ern region experienced uniform cold conditions.The time coefficient represents the d of coldness, and the higher the positive phase, the greater the degree of coldness region.Similarly, Figure 7 shows the explained variance for the southern Yangtze region (79%).The time coefficient indicates uniform coldness.After removing the annual cycle and annual average from the daily minimum temperature series from 1995 to 2020, the Butterworth filtering of the low-frequency oscillation series was conducted after station-by-station filtering for 10-20 d and 30-60 d.EOF decomposition was performed on the low-frequency oscillation series at multiple stations.As shown in Figure 6, the explained variances of the first mode of the minimum temperature on the two time scales (10-20 d and 30-60 d) for the northern Yangtze River region were 78.9% and 80%, respectively.The physical meaning of the first mode is that the northern region experienced uniform cold conditions.The time coefficient represents the degree of coldness, and the higher the positive phase, the greater the degree of coldness in this region.Similarly, Figure 7 shows the explained variance for the southern Yangtze River region (79%).The time coefficient indicates uniform coldness.After removing the annual cycle and annual average from the daily minimum temperature series from 1995 to 2020, the Butterworth filtering of the low-frequency oscillation series was conducted after station-by-station filtering for 10-20 d and 30-60 d.EOF decomposition was performed on the low-frequency oscillation series at multiple stations.As shown in Figure 6, the explained variances of the first mode of the minimum temperature on the two time scales (10-20 d and 30-60 d) for the northern Yangtze River region were 78.9% and 80%, respectively.The physical meaning of the first mode is that the northern region experienced uniform cold conditions.The time coefficient represents the degree of coldness, and the higher the positive phase, the greater the degree of coldness in this region.Similarly, Figure 7 shows the explained variance for the southern Yangtze River region (79%).The time coefficient indicates uniform coldness.We calculated the proportion of the 46 RPECEs occurring since 1979 in the phases of the low-frequency components, and their time coefficient amplitude was higher than one standard deviation.The proportion was 90% in the northern Yangtze River region and 80% in the southern Yangtze River region, indicating that most RPECEs were strongly correlated with low-frequency oscillations.Thus, we can predict the occurrence of RPECE events by predicting the time coefficient of low-frequency oscillations.The standard deviation of the positive phase of the time coefficient of the low-frequency oscillation was calculated.If the time coefficient of the low-frequency oscillation was higher than the positive phase of one standard deviation in a certain period, it was highly likely that it was a lowfrequency oscillation of RPECEs.Our goal was to predict the time coefficient of low-frequency oscillations.However, it is worth noting that not every high positive phase indicates a large-scale low-temperature event.

Extended-Range Forecast
The CNN deep learning method commonly used in climate prediction [26,27] was used to obtain an extended-range forecast of RPECEs.Specifically, we predicted the time coefficient of low-frequency oscillations.

Deep Learning Model for Extended-Range Forecast
The CNN was used to determine the spatial-temporal characteristics of the low-frequency field and the OLR field in the current time (t) and the previous time (t-25 d, t-20 d, t-15 d, t-10 d, and t-5 d).The correlation between the feature set of large-scale circulation factors and the time coefficient of the low-frequency oscillation component of the minimum air temperature was calculated.The structure of the CNN model is shown in Figure 8.It has three convolutional layers, two pooling layers, and one fully connected layer.The convolution operation extracts features from the two-dimensional longitude and latitude data, and the pooling operation reduces the feature dimensions of adjacent data at certain locations.The fully connected layer improves the nonlinear fitting performance of the CNN by stacking multiple neurons.
The first layer (18) represents 18 channels, i.e., 6 times H200, H500, and OLR, and the two-dimensional vector of each channel is 28×44, because the spatial range of each element is (0°-67.5°N, 42.5°-150° E), and the spatial resolution is 2.5° × 2.5°.The output is 1 × 1, which is the time coefficient of the low-frequency oscillation component of the minimum temperature.It is predicted for the following 15 d.
The convolution kernel size during network training was 3 × 3, and the size of the pooling kernel was 2 × 2. Maximum pooling was adopted; the learning rate was 0.005, the sample batch size was 32, and the loss function was  -loss.The ReLU activation function and Kaiming initialization were used.A dropout layer was inserted before the fully connected layer to prevent overfitting.The features extracted by the CNN were discarded or kept based on their probability.We calculated the proportion of the 46 RPECEs occurring since 1979 in the phases of the low-frequency components, and their time coefficient amplitude was higher than one standard deviation.The proportion was 90% in the northern Yangtze River region and 80% in the southern Yangtze River region, indicating that most RPECEs were strongly correlated with low-frequency oscillations.Thus, we can predict the occurrence of RPECE events by predicting the time coefficient of low-frequency oscillations.The standard deviation of the positive phase of the time coefficient of the low-frequency oscillation was calculated.If the time coefficient of the low-frequency oscillation was higher than the positive phase of one standard deviation in a certain period, it was highly likely that it was a low-frequency oscillation of RPECEs.Our goal was to predict the time coefficient of low-frequency oscillations.However, it is worth noting that not every high positive phase indicates a large-scale low-temperature event.

Extended-Range Forecast
The CNN deep learning method commonly used in climate prediction [26,27] was used to obtain an extended-range forecast of RPECEs.Specifically, we predicted the time coefficient of low-frequency oscillations.

Model Validation
We As shown in Figure 9, three of the four events were correctly predicted.Similar trends The first layer (18) represents 18 channels, i.e., 6 times H200, H500, and OLR, and the two-dimensional vector of each channel is 28×44, because the spatial range of each element is (0 • -67.5 • N, 42.5 • -150 • E), and the spatial resolution is 2.5 • × 2.5 • .The output is 1 × 1, which is the time coefficient of the low-frequency oscillation component of the minimum temperature.It is predicted for the following 15 d.
The convolution kernel size during network training was 3 × 3, and the size of the pooling kernel was 2 × 2. Maximum pooling was adopted; the learning rate was 0.005, the sample batch size was 32, and the loss function was l 1 -loss.The ReLU activation function and Kaiming initialization were used.A dropout layer was inserted before the fully connected layer to prevent overfitting.The features extracted by the CNN were discarded or kept based on their probability.As shown in Figure 9, three of the four events were correctly predicted.Similar trends were observed for the predicted and observed values.The prediction for the period in Figure 9a exhibits a phase deviation.Two criteria were used to assess the prediction accuracy: where NA represents hits, NB represents false alarms, and NC represents misses.POD denotes the probability of detection, and FAR is the false alarm rate.POD is 0.75, and FAR is 0.76.A positive phase higher than one standard deviation corresponds to low-frequency oscillation.The average correlation coefficient between the observed and predicted lowfrequency oscillation component for the test set is 0.41 for the deep learning method and 0.1 for the linear method.

Model Validation
We used the 2010/2011-2019/2020 dataset for model validation.Four RPECEs occurred during the test period: 28 October to 2 November in 2010, 10 February to 15 February in 2014, 23 January to 27 January in 2016, and 25 January to 30 January in 2018.
As shown in Figure 9, three of the four events were correctly predicted.Similar trends were observed for the predicted and observed values.The prediction for the period in Figure 9a exhibits a phase deviation.Two criteria were used to assess the prediction accuracy: where NA represents hits, NB represents false alarms, and NC represents misses.POD denotes the probability of detection, and FAR is the false alarm rate.POD is 0.75, and FAR is 0.76.A positive phase higher than one standard deviation corresponds to low-frequency oscillation.The average correlation coefficient between the observed and predicted lowfrequency oscillation component for the test set is 0.41 for the deep learning method and 0.1 for the linear method.(1979/1980-1988/1989) had more extreme cold events.We obtained predictions for 15 RPECEs.POD is 0.6, and FAR is 0.55.POD and FAR are lower than those in the previous experiment.The results indicate that the model has a reasonable predictive ability for the low-frequency oscillation time coefficient in this period and can describe the trend of the low-frequency oscillation curve.

Conclusions
The innovation of this article lies in the following points: • N. Since the OLR data started in 1979, our prediction test was conducted for the data within the range of 1979/1980-2019/2020.First, 10 years out of 31 years (1979/1980-2009/2010) were randomly selected as the validation set, the remaining 21 years represented the training set, and 2010/2011-2019/2020 was the test set (1670 d).
Figure 1), with a dividing line at 30° N. Since the OLR data started in 1979, our prediction test was conducted for the data within the range of 1979/1980-2019/2020.First, 10 years out of 31 years (1979/1980-2009/2010) were randomly selected as the validation set, the remaining 21 years represented the training set, and 2010/2011-2019/2020 was the test set (1670 d).

Figure 1 .
Figure 1.Map of weather stations in southern China.

Figure 1 .
Figure 1.Map of weather stations in southern China.

Figure 2 .
Figure 2. Composite of temperature anomalies in 850 hPa (colors, unit: °C) and geopotential height anomalies in the 500 hPa (contours, unit: gpm) of RPECEs (the dots indicate that the relationship does not exceed the 0.05 significance level).

Figure 2 .
Figure 2. Composite of temperature anomalies in 850 hPa (colors, unit: • C) and geopotential height anomalies in the 500 hPa (contours, unit: gpm) of RPECEs (the dots indicate that the relationship does not exceed the 0.05 significance level).

Figure 3 .
Figure 3. Composite of wind anomalies in 850 hPa (vectors, unit: m/s) and sea level p alies (contours, unit: hPa) of RPECEs (the dots indicate that the relationship does not significance level).

Figure 3 .
Figure 3. Composite of wind anomalies in 850 hPa (vectors, unit: m/s) and sea level pressure anomalies (contours, unit: hPa) of RPECEs (the dots indicate that the relationship does not exceed the 0.05 significance level).

Figure 3 .
Figure 3. Composite of wind anomalies in 850 hPa (vectors, unit: m/s) and sea level p alies (contours, unit: hPa) of RPECEs (the dots indicate that the relationship does not e significance level).

Figure 4 .
Figure 4. Minimum temperature index (a), Morlet wavelet power spectrum (b), global trum (c), and wavelet coefficient (d) from 1960 to 2020 in the northern Yangtze River r

Figure 6 .
Figure 6.The explained variance of the first mode of the minimum temperature in the no

Figure 5 .
Figure 5. Minimum temperature index (a), Morlet wavelet power spectrum (b), global wavelet spectrum (c), and wavelet coefficient (d) from 1960 to 2020 in the southern Yangtze River region.

Figure 6 .
Figure 6.The explained variance of the first mode of the minimum temperature in the northern Yangtze River region for (a) 10-20 days and (b) 30-60 days.

Figure 6 .
Figure 6.The explained variance of the first mode of the minimum temperature in the northern Yangtze River region for (a) 10-20 days and (b) 30-60 days.

Figure 7 .
Figure 7.The same as Figure 6 but for the southern Yangtze River region.

Figure 7 .
Figure 7.The same as Figure 6 but for the southern Yangtze River region.
used the 2010/2011-2019/2020 dataset for model validation.Four RPECEs occurred during the test period: 28 October to 2 November in 2010, 10 February to 15 February in 2014, 23 January to 27 January in 2016, and 25 January to 30 January in 2018.

4. 2 .
Model Validation We used the 2010/2011-2019/2020 dataset for model validation.Four RPECEs occurred during the test period: 28 October to 2 November in 2010, 10 February to 15 February in 2014, 23 January to 27 January in 2016, and 25 January to 30 January in 2018.

Figure 9 .
Figure 9. Predicted and observed RPECEs: (a) 10 February 2014-15 February 2014 (30-60 d), (b) 23 January 2016-27 January 2016 (30-60 d), and (c) 25 January 2018-30 January 2018 (30-60 d).We conducted another numerical experiment using test data in the range of 1979/1980-1988/1989; 10 years out of 31 years (1989/1990-2009/2010) were randomly selected as the validation set, and the remaining 21 years were used as the training set.Previous periods (1979/1980-1988/1989) had more extreme cold events.We obtained predictions for 15 RPECEs.POD is 0.6, and FAR is 0.55.POD and FAR are lower than those in the previous experiment.The results indicate that the model has a reasonable predictive ability for the low-frequency oscillation time coefficient in this period and can describe the trend of the low-frequency oscillation curve.

Figure 9 .
Figure 9. Predicted and observed RPECEs: (a) 10 February 2014-15 February 2014 (30-60 d), (b) 23 January 2016-27 January 2016 (30-60 d), and (c) 25 January 2018-30 January 2018 (30-60 d).We conducted another numerical experiment using test data in the range of 1979/1980-1988/1989; 10 years out of 31 years (1989/1990-2009/2010) were randomly selected as the validation set, and the remaining 21 years were used as the training set.Previous periods (1979/1980-1988/1989) had more extreme cold events.We obtained predictions for 15 RPECEs.POD is 0.6, and FAR is 0.55.POD and FAR are lower than those in the previous experiment.The results indicate that the model has a reasonable predictive ability for the low-frequency oscillation time coefficient in this period and can describe the trend of the low-frequency oscillation curve.

Table 1 .
The start and end dates, durations, and peak dates of the RPECEs.