A Review of the Detection Methods for Climate Regime Shifts

An abrupt climate changemeans that the climate system shifts from a steady state to another steady state. Study on the phenomenon and theory of the abrupt climate change is a new research field ofmodern climatology, and it is of great significance for the prediction of future climate change. The climate regime shift is one of the most common forms of abrupt climate change, which mainly refers to the statistical significant changes on the variable of climate system at one time scale. These detection methods can be roughly divided into five categories based on different types of abrupt changes, namely, abrupt mean value change, abrupt variance change, abrupt frequency change, abrupt probability density change, and the multivariable analysis. The main research progress of abrupt climate change detection methods is reviewed. What is more, some actual applications of those methods in observational data are provided. With the development of nonlinear science, many new methods have been presented for detecting an abrupt dynamic change in recent years, which is useful supplement for the abrupt change detection methods.


Introduction
The research on the phenomenon of abrupt climate change and the related theory is an important field in modern climatology, which is one of the core issues concerning the global climate change.It is also important for the prediction of climate change [1][2][3][4][5][6][7][8][9][10].Since Lorenz [11] and Charney and Devore [12] revealed abrupt climate change on a theoretical basis, the studies on abrupt climate change have been broadly conducted [13][14][15][16][17][18][19][20][21][22].Although the studies on abrupt climate change are relatively late in China, many Chinese climatologists have noticed this phenomenon [23][24][25][26].For example, Chinese scholars first noted the phenomenon of seasonal abrupt change in atmospheric circulation and climate [27][28][29][30][31][32][33][34].Abrupt climate change is a ubiquitous phenomenon in nature.Regarding its physical meaning, it is a transition of the climate system into a different mode on a time scale that is faster than the responsible forcing.Regarding its influence, the occurrence of abrupt change is so rapid and unexpected that human or natural system have difficulty adapting to it [35,36].These two definitions are complementary: the former provides what abrupt climate change is, and the latter explains why there are so many studies on abrupt climate change.
Because of the suddenness for the occurrence of abrupt climate change, it has a potentially significant impact on the sustainable development of the ecological environment, the social economy, or even the extinction of species.For example, historically, the drought in the region of Mesopotamia led to the decline of ancient civilizations, such as the ancient Greek civilization, the Egyptian civilization in the Nile Valley, and the ancient Indian civilization in the Indus Valley [37]; the Permian-Triassic dinosaur extinction event was also related to abrupt climate change [38].In 2013, a report by the United States National Research Council noted that when the climate system reaches a certain critical threshold, even a gradual and stable change in the climate system would affect the human infrastructure and the ecosystem.Therefore, the accurate detection and attribution of historical abrupt climate change events have significant practical meaning and potential socioeconomic value for monitoring and forecasting the transition or abrupt change that the future climate system may possibly face [39].
At present, there are mainly three aspects of studies on abrupt climate change: early-warning signals for abrupt climate change, detection methods for abrupt climate change, and attribution studies on abrupt climate change events.Scheffer et al. [40] proposed that many complex systems had a critical threshold, that is, a critical point, and that the system would change from one state to another state at this point.This phenomenon is related to the "critical slowing down" phenomenon in the theory of dynamic systems.During the process in which the system approaches the critical point, the critical slowing down phenomenon causes three possible early-warning signals in the dynamical system: a slowing recovery from disturbance, an amplification of autocorrelation, and an increase of variance.There are also numerous studies on the possible reasons that have led to all types of abrupt climate change throughout history.Alley et al. [14] noted that the influence of human activity might increase the probability of the occurrence of abrupt climate change, and once abrupt climate change occurred, it was likely to have a serious potential impact on the economy and the ecosystem.The change in Atlantic thermohaline circulation is likely one reason that caused abrupt climate change [41].That is, under certain conditions, because of a rapid thawing of ice and melting of snow, more precipitation, and larger discharges from rivers and lakes, much more freshwater is injected into the North Atlantic, therefore causing the salinity and density decrease in the surface water in the North Atlantic.This phenomenon blocks or cuts off the Atlantic thermohaline circulation, causing the temperature to decline rapidly in the surrounding areas of the North Atlantic, leading to abrupt climate change.
The accurate identification and detection of abrupt climate change is the prerequisite and basis for studies on earlywarning signals and attribution studies of abrupt climate change.Many scholars have conducted a lot of relevant studies, such as detection with respect to transition in the climate system state.Rodionov [42] divided these methods into four groups, provided a brief description of each method, as well as their strong and weak points, and pointed out that shifts in the mean are the most common type.Mantua [43] provided a brief review and assessments of the relative strengths and weaknesses of five analytical methods: principal components analysis, vector autoregressive process modeling, autoregressive moving average and intervention analysis modelling, compositing average standard deviates, and Fisher information.Ramanayake and Gupta [44] proposed three statistics, explored their distribution theories, and applied these tests to air traffic data and Stanford Heart Transplant Data.Transition in the climate state, which mainly refers to the notion that there is a statistically significant change in the statistical features of the system state variable over a certain time scale, is one of the most common types of abrupt change phenomenon.This type of abrupt climate change can be roughly divided into five categories: abrupt mean value change, abrupt variance change, abrupt frequency change, abrupt probability density change, and multivariate analysis.
In this paper, we mainly review the progress on the theory and methodological studies in China and other countries regarding the detection of regime shifts in the climate system.The discussions are focused on describing the climate regime shifts detection methods, and some applications on actual observation data.This paper is organised as follows.In Section 2, we briefly introduce the abrupt change detection methods in mean value.Section 3 shows the detection methods for the abrupt probability density change, and the other abrupt change detection methods are presented in Section 4. Discussions and conclusions follow in Section 5.

Abrupt Mean Value Change Detection Methods
The detection methods for the abrupt mean value changes can be divided into six categories, which mainly include parameter detection method, non-parameter detection method, cumulative sum method, Bayesian analysis method, sequential method, and detection methods based on regression analysis.The moving -test (MTT) [45,46] is a classical parameter detection method, which is a widely applied method in meteorology.It is often used to determine whether abrupt change will occur by testing the difference between the mean values of two sequences.The advantages of MTT include simple procedures, relatively high operational efficiency, and clear physical meaning, that is, strong theoretical basis.However, the assumption of normality needs to be satisfied in using MTT, which makes its scope of application limited.Li et al. [45] presented a successive MTT method and used the MTT method to conduct a comprehensive detection on global abrupt climate change.They found that there was the phenomenon of abrupt climate change in the intensity of atmospheric circulation, South Asian summer monsoon, and global temperatures in the last century.Xiao et al. [16,18] also employed the successive MTT method to detect the decadal abrupt change in the time series of several atmosphere-ocean variables and found that the warming of the northern midlatitude SST might be responsible for the northern atmospheric decadal abrupt changes and caused the climate change events in the 1980s.The nonparametric detection methods mainly include the Mann-Kendall test (M-K test) [47,48], the Mann-Whitney  test [49], the Lepage test [50], and wavelet analysis [51].The M-K test method presumes that the sequence of elements to be analysed is a stationary time series, which obeys normal distribution, and individual samples are relatively independent.M-K method is mainly used to test the trend of abrupt change in a time series.For a time series of {  ,  = 1, 2, 3, . . ., }, the cumulative number for the records which are less than   in subsequences {  ,  = 1, 2, 3, . . ., } is represented by   , and the total cumulative number is denoted as follows: The mean and variance of the newly established sequence {  ,  = 1, 2, 3, . . ., } can be, respectively, defined as The standardization of the newly established sequence {  } can be obtained: where  (  ) obeys the standard normal distribution.If the probability density distribution is for a given significance level  0 , the original hypothesis is valid for  1 >  0 .The hypothesis is null for  1 <  0 , which indicates that there is a significant trend change in this series.
The standardized new sequence  (  ) consists of a curve of temporal variation.If the curve falls in the confidence interval, there is no significant variation trend in the original sequence; otherwise, the variation trend is significant.The M-K method also has wide application in detecting the abrupt trend change in a time series.For example, Ma and Fu [52] used the M-K test to analyse the time-space structure of drywet change in the northern China over the past 54 years and found that there was only one dry-wet conversion on the decadal scale in the eastern area of northwest China and northern China at the end of 1970s.He et al. [53] found that the M-K results were completely different for an artificial time series when the M-K method was applied to analyse different numbers of samples, which indicated that the M-K results heavily depend on sample sizes.There were also failures which occurred in the application of MK method in detecting the changes points, such as no change point in an obvious step time series [54].The Mann-Whitney  test, also called the rank-sum test, was proposed jointly by Mann and Whitney in 1947.This method assumes two samples that derive from two different populations for which all of the statistical characteristics are completely identical except for the mean value.On this basis, it tests whether there is a significant difference in the mean values of these two populations, that is, whether there is an occurrence of abrupt change.The detailed procedures of this algorithm are as follows: Original hypothesis  0 : There is no significant difference between the two samples.
1 : There is a significant difference between the two samples.
Step 1. Combine two sets of data, and assign numeric ranks to all the records, beginning with 1 for the smallest value (where there are groups of tied values, assign a rank equal to the midpoint of unadjusted rankings [e.g., the ranks of (3,5,5,9) are (1, 2.5, 2.5, 4)]).
Step 2. Calculate the sum of the ranks for the two samples, and  1 and  2 are the sums of the ranks in the first sample and another one, respectively.
Step 3. Calculate the statistical quantity  of the Mann-Whitney  test: where  1 is the sample size for the first sample, and  2 is the sample size for the second sample.The smaller value of  1 and  2 is the one used when consulting significance tables and comparing with the given critical value   .If the smaller value of  1 and  2 is smaller than   , the hypothesis  1 is accepted; that is, there is the occurrence of an abrupt change with a significance level   .
If the original hypothesis is true, then the mean value and variance of random variable  are, respectively, ( If  1 and  2 are both not smaller than 10, the random variable  approximately meets the normal distribution.
The Mann-Whitney  test is more widely applicable than independent samples Student's -test, and the question is, which one of them should be preferred.
The Lepage test is a nonparametric statistical method for testing whether there is a significant difference between the mean values of two independent samples, even if the distributions of parent populations are unknown.If there is a significant difference between the two subsequences at a significant level, it is concluded that there is an abrupt change.
Suppose that the number of samples for the subsequence before the reference point is  1 , the number of samples for the subsequence thereafter is  2 , and  12 is the sum of  1 and  2 .And then, we can calculate the rank sequence {  ,  = 1, 2, 3, . . .,  12 }.If the minimum appears before the reference point,   = 1; otherwise,   = 0. We then construct a quantity of rank statistics as : is the cumulative number of the two subsequences.Its mean value and variance are, respectively, Additionally, we then construct another quantity of rank statistics as : is the sum of the cumulative number of the two subsequences.The first half is the cumulative number of the subsequence before the reference point, and the second half is the cumulative number of the subsequence after the reference point.If  12 is even, the mean and variance of  are as follows: Otherwise, A combine of the squares of the standardized Wilcoxon and Ansariy-Braley's statistics is the Lepage test: The Lepage test has more powerful functions than the moving -test and the rank-sum test.Not only can it be used to test abrupt changes but it can also be applied to study the linear trend and periodic variation of climate change.Yonetani [55,56] provided a detailed description of the procedures for using the Lepage method in climate change and applied this method to test abrupt climate change in Japan since 1900.The wavelet analysis is also another abrupt change detection method based on the non-parametric regression function.Antoniadis and Gijbels [57] applied the wavelet analysis method to detect abrupt change and found that this method could well detect an abrupt change in the discharge of the Nile River that occurred in 1898.
The cumulative sum method is the third type of detection method for the abrupt mean value change, and it mainly includes the cumulative sum test and the cumulative frequency test.The cumulative sum method was proposed by Page [58] and then it was subsequently continuously improved [59][60][61].The main idea of this method is to obtain the evolutionary trend of a time series by calculating the cumulative sum of the differences between the values and the average in the time series, and then we can identify whether there is a significant trend change in the time series.If the values are all above average during a period of time, the amounts added to the cumulative sum will be positive and the sum will steadily increase, and vice versa.The algorithm of this method is relatively simple, and the cumulative anomaly,   (), for a time series can be obtained as follows: where   is the sample anomaly.When the absolute value of   () reaches the maximum, the corresponding time is likely the time for the occurrence of the abrupt change.Janssen [62] applied the model to test the cumulative sum method and used this method to test the change in GPS data.Elsner et al. [63] used the cumulative frequency method to test the occurrence frequency of hurricanes in two areas where hurricane often happens and found that there was a significant change in the frequency of hurricanes.
The result was consistent with climatic warming in North Africa.
The Bayesian analysis is the fourth detection method type for the abrupt mean value change, which is developed based on Bayesian theory.Firstly, we can integrate the a priori information and sample information of unknown parameters and then derive the a posteriori information according to Bayesian theory.The unknown parameters can be inferred according to the derived posteriori information.The Markov Chain Monte Carlo (MCMC) approach [64] is one of Bayesian analysis.MCMC was developed in the early 1950s.Within the framework of Bayesian theory, the Markov process was introduced into the Monte Carlo simulation, which compensated for the defect in traditional Monte Carlo integration; that is, it only can be used to conduct a static simulation.The procedures of the MCMC algorithm are shown as follows.
Step 1. Construct the Markov chain and make it converge with the stationary distribution.
Step 2. Draw the sample from a target space, apply the chain constructed in Step 1 to sample and simulate, and then generate the sequence.
Step 3. Conduct the Monte Carlo integration, and the expected estimate of any function is as follows: Chu and Zhao [65] used the Bayesian analysis to test the variation in the annual number of tropical cyclones in the central North Pacific during 1966-2002 and found that it was extremely possible that there was an abrupt mean change in the number of tropical cyclones in 1982 (Figure 1).Elsner et al. [66] applied MCMC to test the variation in the annual number of hurricanes in North Africa during the 20th century and disclosed that there were three abrupt climate changes in the annual number of hurricanes in the middle of the 1940s, the 1960s, and 1995, respectively.
The sequential method can find the maximum of the mean value from a number of normal distributions with different mean values and then can be used to conveniently and rapidly detect an abrupt mean change.Rodionov [67] noted that the sequential method was better than MTT.In using the sequential method, it is not necessary for making an a priori assumption for a time series and can be used to detect the change point in advance.The sequential -test analysis method (STARS) has several advantages, including analysing the real-time data, timely monitoring of the signal of the change point, and the evolutionary trend of the signal strength.STARS does not need to make an a priori assumption on the change point in a time series and also can be used to deal with the time series with multiple change points.Rodionov and Overland [68] used STARS to test the marine ecosystem index of the Bering Strait and found that the climate transition in 1977 was the strongest in all sea-air index records related to the Bering Strait since the beginning of the 20th century.Reid et al. [69] had employed the Rodionov method to perform a comprehensive detecting on the 1980s regime shift.The shift in temperature that occurred at this time is hypothesized as the main forcing for a cascade of abrupt environmental changes.
Bernaola-Galvan et al. [70] introduced a heuristic segmentation algorithm to probe the temporal organization of heterogeneities in human heartbeat interval time series.
They found that the lengths of segments with different local mean heart rates followed a power-law distribution.Feng et al. [71] verified the performances of the segmentation algorithm in dealing with nonlinear time series by using a number of artificial time series and found that the method could be effectively used to detect an abrupt mean change.Gong et al. [72] used the heuristic segmentation algorithm to analyse northern hemisphere tree rings and Beijing stalagmite based on high and low frequency series and distinguished abrupt changes in different scales and disclosed its physical mechanism.
The two-phase regression technique (TPR) is a typical method for detecting abrupt change based on regression analysis.Hinkley [73,74] proposed the TPR model and presumed that the regression response is continuous at the point of abrupt change.This assumption is unrealistic, and the influence of the actual environment and instruments may result in the sequence discontinuous.Lund and Reeves [75] found this problem and made an improvement (i.e., the LR method), but the improved method can cause several false points of abrupt change.Subsequently, Wang [76] also found that although the LR method could make the presumed climatic conditions unrealistic because abrupt change might not affect the mean value and trend simultaneously.In order to deal with this problem, Wang proposed a trend model, and the advantage of this model is that it presumes the presence of a trend in a time series and allows the trend to change at the point of abrupt change.
The composite analysis [77,78] is a detection method for a regime shift, and a composite analysis of climatic records is likely to find them in Gaussian, red noise with stationary statistics.Rudnick and Davis [77] examined the composite analysis using random independent time series generated to have the same frequency content as the Pacific Decadal Oscillation (PDO) and found that it was likely to find ostensibly significant steps even if the input time series were red noise.Rodionov [78] proposed a prewhitening procedure which was designed to remove the red noise component from time series prior to applying a regime shift detection method.The application of this method to analyze PDO shows that PDO appears to be more than just a manifestation of a red noise process.The size of subsamples should be as small as possible to minimize the effect of regime shifts on the parameter in the first order autoregressive model.On the other hand, it should be as large as possible to minimize the sampling variability of the parameter.These opposite requirements to the size of subsamples make it difficult to apply this procedure to short time series with decadal or shorter regimes.

Detection Methods for the Abrupt Probability Density Change
The methods that detect the change in the probability density distribution of elements include the abrupt density change, skewness and kurtosis coefficient, and sliding skewness index method based on the Box-Cox transformation.For a stable dynamic system, the probability density of the system variables has a relatively stable distribution pattern, although, when the dynamic structure of the system changes, it may cause different degrees of change in the distribution patterns of the system variables.From the perspective of identifying small changes in the probability density distribution of the system variables, Cheng et al. [79] applied the coefficient of skewness and kurtosis, which described the characteristics of the probability distribution for a time series, to detect an abrupt change in time series.The numerical test indicates that the coefficients of skewness and kurtosis have a very good capability to identify an abrupt probability density change in artificial signals.The detection results have a relatively small dependence on the sample sizes.
Skewness is a measure of the asymmetry of the probability distribution of a real-valued random variable about its mean.Coefficient of skewness can be positive or negative, or even undefined.For a time series {  ,  = 1, 2, 3, . . ., }, coefficient of skewness  is defined as The parameter  is the sample standard deviation.The coefficient of skewness  = 0 when the probability density curve is normal distribution.For  < 0, the distribution has a negative deviation and the left tail is longer than the right one, and vice versa.The skewness is not strictly connected with the relationship between the mean and median: a distribution with negative skew can have the mean greater than or less than the median, and likewise for positive skew.
Kurtosis is a descriptor of the shape of a probability distribution and, just as for skewness, there are different ways of quantifying it for a theoretical distribution and corresponding ways of estimating it from a sample from a population.Coefficient of kurtosis is mainly used to describe the steep degree of the probability density distribution curve.The standard measure of kurtosis is based on a scaled version of the fourth moment of the data or population, and the coefficient of kurtosis  is defined as The coefficient of kurtosis of a normal distribution is  = 3.
For  > 3, the probability distribution is steep compared to that for normal distribution, and the distribution is flat compared to normal distribution for  < 3.
It should be noted that when the difference in the steepness or skewness of the probability density curve before and after an abrupt change is relatively small, it might need more samples to identify the abrupt change to ensure the statistical reliability of the detection results.Consequently, under certain conditions, the coefficients of skewness and kurtosis most likely have a relatively high requirement regarding the sample size and are likely not applicable to those time series with relatively a small sample size.
He et al. [80] constructed 1000 groups of daily precipitation records with two change points which generated by a weather generator.In these artificial precipitation data, the probability density distributions are obviously differences before and after the change points.The testing results showed that the coefficients of skewness and kurtosis both cannot detect the abrupt change in most of the 1000 groups of artificial data.The coefficient of skewness can only detect two groups, and the coefficient of kurtosis cannot detect these preset change points.The main causes for this were that the daily rainfall data had very strong skewness and the coefficient of skewness and kurtosis might not be sensitive enough to the changes of the probability density distribution of daily precipitation data.In view of this, He et al. [80] put forward the sliding transformation parameter (STP) based on the Box-Cox transformation.A large amount of numerical experiments showed that STP could identify the abrupt probability density change in the artificial series, and it performed much better than the kurtosis coefficient and the skewness coefficient.
Box-Cox transformation [81] is a mathematical transformation tool which can transform a time series of skewed distribution into a new time series with approximative normal distribution.For a time series {  ,  = 1, 2, 3, . . ., }, Box-Cox transform is defined as where   is the response variable by Box-Cox transformation, and  is the transformation parameter.This formula only can be used when observational data   > 0. Thus, Box and Cox [82] proposed a modified version which could be used for all kinds of data, and the revised Box-Cox transform is defined as where  is a positive number and   + > 0. The parameter  can be estimated by the maximum likelihood method.
In order to test the performance of the STP method, the 1000 artificial daily precipitation series have been analyzed, and the results indicate that STP can almost recognize all of the first change point and the accuracy for the second change point is 99.8%.It shows the ability of STP to identify abrupt probability density changes in daily precipitation series.To examine the performances of STP in the observed daily precipitation records, He et al. randomly select six observational stations in China and find that all the STP results for six daily precipitation records display almost identical evolutionary trend.The results indicate that the transformation parameters undergo a rapid shift from one stable state to an alternative stable state with obviously different means between 1979 and 1980 (Figure 2), which can be detected by MTT with a significance level  = 0.01 with the sample size of 10 years.
Moreover, Ning and Gupta [83] studied the performance of generalized lambda distributions (GLD) for detecting multiple change points of parameters.Compared to the other distributions, the advantage of GLD is that the GLD family is broad and flexible.The method was applied on fibroblast cancer cell line data, and the change points are successfully located.

Other Types of Abrupt Detection Methods
The abrupt variance change can be determined by the change in the second-order statistics of elements.There are few methods to test this type of abrupt change.The Downton-Katz test [84] does not need assumptions for the distribution frequency, but it needs a reference time series without abrupt change.The Nikiforov method [85] can be used to detect an abrupt change in the frequency structure of elements.In using the Nikiforov method, the time series before and after an abrupt change are first introduced into a differential autoregression and moving average model, which can be used to approximately describe the original time series.And then it can detect an abrupt frequency change in the time series by combining with a likelihood ratio test.The Nikiforov test has a strong theoretical basis, but it is only suitable to detect single abrupt change point.Multivariate analytical methods, such as principal component analysis [43], average standard deviation, and vector autoregression are often used to identify a shift of the spatial pattern of a system.Principal component analysis is a statistical procedure that uses an orthogonal transformation to convert a set of observations of possibly correlated variables into a set of values of linearly uncorrelated variables called principal components.The number of principal components is less than or equal to the number of original variables.This transformation is defined in such a way that the first principal component has the largest possible variance.Principal component analysis is the simplest of the true eigenvector-based multivariate analyses and is widely used to identify the coherent patterns of spatial variability.Hare and Mantua [86] applied principal component analysis to study the marine ecosystem in the North Pacific and Bering Strait and found that the abrupt climate change in 1976-1977 had a certain influence on the ecosystem in the North Pacific.

Discussion and Conclusions
In this paper, we briefly review some advances in the abrupt climate change detection methods in recent years.According to the types of abrupt change, the methods for detecting a state transition in the climate system can be divided into five categories: abrupt mean value change, abrupt variance change, abrupt frequency change, abrupt probability density change, and multivariate analysis.There are relatively numerous research contents on abrupt mean value change and abrupt probability density change.Based on this, we provide a considerably detailed summary.The detection methods for abrupt mean value change mainly have six subcategories, which include the parametric method, the non-parametric method, the cumulative sum method, the Bayesian analysis method, the sequential method, and test algorithms based on regression analysis.The methods for detecting abrupt change in probability density distribution of elements include the coefficient of kurtosis and skewness coefficient and the STP method based on the Box-Cox transformation.There are few testing methods based on an abrupt change in variance, frequency, and multivariate analysis, and the details on the algorithmic methods are not described in the present paper.
In this paper, we mainly describe the methods for climate regime shift detection.These methods have one common shortcoming; that is, the detection results heavily depend on the used sample size.In other words, the detected abrupt change point possesses the characteristics of multiple time scales.The reason for this is that an abrupt change event is often determined according to a statistical significant change of a state variable.So, these detection methods for regime shift are not applicable to the abrupt change in the climate dynamic structure.With the development of nonlinear science, many new methods and new techniques have been presented in recent years, especially for the detection methods for abrupt dynamical change.These methods and techniques are mainly based on the long-range correlation in climate system, the complexity property of time series, and the reconstruction theory of phase space.Different from the detection methods for regime shift, the detection results of these new methods have less dependence on sample size, which are applicable to identify whether there is an abrupt dynamic change in the climate system from single time series.

Figure 1 :
Figure 1: The time series of tropical cyclones number in the central of North Pacific during 1966-2002 and its MCMC detection results [65].

Figure 2 :
Figure 2: The STP results for daily precipitation records (from 1961 to 2010) in six meteorological observational stations in China, and the sample size of sliding window is 1 year: (a) Tonghe in Heilongjiang Province; (b) Hailun in Heilongjiang Province; (c) Qiqihaer in Heilongjiang Province; (d) Zhalantun in Inner Mongolia; (e) Sunwu in Heilongjiang Province; (f) Nenjiang in Heilongjiang Province [80].