1. Introduction
Considerable progress has been made in numerical weather prediction (NWP) during the past decades due partly to a steady accumulation of scientific knowledge and partly to technological advances in utilizing a variety of observations and gaining computing power (Bauer et al. 2015). Despite the steady progress, significant forecast errors in today’s NWP models are still present, especially in processing atmospheric statistical properties that are not directly available from the NWP outputs, such as tropical cyclone (TC) intensity, and other mesoscale and smaller-scale features. Thus, it is highly desirable to develop alternative approaches to reducing NWP errors. To our knowledge, Koo (1958a) was among the first to make such an attempt by treating weather forecasts as both an “initial-value” and “evolution” problem. He hypothesized that the three-dimensional flow structures of a large-scale baroclinic system at a given moment could be reconstructed from the evolution of flow contours and temperature fields before and after. In a subsequent study, Koo (1958b) stressed the importance of using historical data in NWP models to improve weather forecasts. Later, Charney et al. (1969) pointed out that historical data, even incomplete sometimes, are useful for inferring the present atmospheric state in NWP models. By transforming a solution of differential equations problem into an extreme-value of functional-variational problem, Chou (1974) expanded the notion of the solution of differential equations by introducing the “generalized solution” concept through incorporating past observations into the temporal evolution of meteorological fields in NWP. Then, Chou (1986) explained why a dynamical model and statistical methods can be combined and how the combined dynamical-statistical model could improve NWP. Briefly, such a dynamical-statistical model can retain the physical significance of the atmospheric processes described by the primitive equations, while it adds the statistical properties of certain quantities that cannot be accurately determined dynamically.
In general, NWP errors come from the following two types of sources: one from the model’s initial conditions and the other from the model’s dynamical and physical deficiencies. The former can be significantly reduced through the application of sophisticated data assimilation and ensemble forecasts, while the latter still remains a challenging research subject (Chou 1974; Lorenc 1986; Kalnay 2002). To reduce prediction errors from model deficiencies, an important alternative is to develop empirical strategies, in addition to increasing model resolution and improving physical parameterizations.
Empirical strategies to reduce forecast errors may be classified into (i) tendency error corrections, which can be made online by nudging the model tendencies toward observations during model’s integrations, and (ii) direct forecast-error corrections, which can be made offline after obtaining model forecasts (Danforth and Kalnay 2008a). The first empirical strategy is NWP model dependent and has been widely used with various data assimilation techniques in both the operational and research communities (e.g., Kalnay 2002). Leith (1978) developed a statistical method in which a state-dependent empirical correction is applied to the dynamical model through minimizing the model tendency errors at every time step. Subsequently, numerous studies have shown improvements of NWP with different models by following Leith’s methodology (e.g., DelSole and Hou 1999; Danforth et al. 2007; DelSole et al. 2008; Danforth and Kalnay 2008a,b; Yang et al. 2008).
The second strategy involves the diagnostic or prognostic use of analogs with the assumption that two similar flow patterns will likely remain similar for some time to follow (Van Den Dool 1994). For the former, Roebber and Bosart (1998) found a strong dependence of observed rainfall structures on mesoscale details after exploring their sensitivities to small variations in the synoptic circulation. The latter has been applied to predicting short-term climate fluctuations (Barnett and Preisendorfer 1978; Livezey et al. 1994), and estimating atmospheric predictability (Lorenz 1969). Further, Chou (1979) suggested that if the forecast state could be regarded as a small perturbation superimposed on a historical analog field, an analog-statistical technique can be applied to the dynamical forecast field. By estimating the model’s tendency error with respect to the historical analog state, Qiu and Chou (1989) obtained a deviation equation in a quasigeostrophic model and then illustrated the improvements of model forecast performance. With the same philosophy, improvements have been reported in medium-range to seasonal forecasts by some studies (Huang et al. 1993; Bao et al. 2004; Yu et al. 2014). A similar methodology, applied to four-dimensional variational assimilation to estimate the initial tendency errors of analog samples, was proposed independently by D’Andrea and Vautard (2000). In addition, some other studies have developed various strategies to reduce model tendency errors (Johansson and Saha 1989; Kaas et al. 1999; Klinker and Sardeshmukh 1992; Saha 1992; Yang and Anderson 2000).
It is worthwhile to mention some recent studies on the hybrid analog-ensemble-based forecasts. Diomede et al. (2014) investigated the calibration of limited-area ensemble rainfall forecasts for driving discharge prediction, and results revealed that the analog-based method has the best performance among three different methods in capably correcting position errors and spreading deficiencies. Some other studies (e.g., Delle Monache et al. 2013; Nagarajan et al. 2015; Junk et al. 2015; Eckel and Delle Monache 2016) focused on developing analog-based ensemble schemes for ensemble model output statistics (EMOS), and their results suggested that the analog-based ensemble forecasts are better than traditional EMOS or MOS forecasts. Zhou and Zhai (2016) developed an analog prediction system for persistent extreme precipitation (PEP), and showed that it is capable of identifying approaching PEP events earlier and yielding a more accurate forecast for the location and intensity of PEP than direct model outputs at 3-day and longer lead times in the Yangtze–Huai River valley. Frediani et al. (2017) used an object-analog technique with 15 ensemble members to forecast 89 storm events in the northeastern United States during the period of 2004–13, and found that the object-analog forecasts are competitively skillful compared to simpler analog techniques. Elsberry and Tsai (2014) developed a situation-dependent intensity prediction (SDIP) technique for western North Pacific TCs, based on the mean intensity changes from the 10 best historical track analogs. This SDIP technique was then successfully applied by Tsai and Elsberry (2015, 2016, 2017) to the other TC cases.
In the present study, we will present a new analog-ensemble-based approach with more attention given to the full use of historical observational data. Thus, the purposes of this study are to (i) improve the ACE approach by replacing historical model forecasts with historical observations; (ii) develop a dynamical-statistical-analog ensemble forecast (DSAEF) model by introducing a concept of the perfect model; and (iii) demonstrate the forecast credibility of the DSAEF model through its application to the heavy rainfall forecasts associated with 21 landfalling TCs.
The next section provides a review of the theory and the current status of the ACE approach. Section 3 shows the forecast improvements of the ACE approach by taking an arithmetic mean of analogs from observations, and how it is used as a backbone for the development of the DSAEF model including the functionals of finding the maximum or minimum values of ensemble members from observations. Some fundamental differences between the DSAEF model and the previous analog-ensemble-based studies will be discussed. Section 4 presents an example of applying the DSAEF model to heavy rainfall forecasts associated with 21 landfalling TCs. A summary and conclusions are given in the final section.
2. Current status of the analog correction of errors
As presented by Ren and Chou (2006, 2007) and Chou and Ren (2006), the ACE approach consists of two steps: an inverse problem of estimating model errors, and the ACE approach to correcting the errors. They are described below separately.
a. An inverse problem of estimating model errors
Since available historical data are essentially the special solutions to (3) or (4), methodologically, the initial-value problem set of PDEs can be transformed into an inverse problem that solves the PDEs by rewriting the formulation of the forecast problem. In doing so, the PDE solutions are known (i.e., from the historical data), whereas the unknowns in (4) can be determined by estimating the model error operator E(ψ). This can be achieved by solving the inverse problem through the ACE as described in the next subsection.
b. The analog correction of errors (ACE)
A simple approach can be used to estimate E(ψ) in a numerical model (Ren and Chou 2006). That is, given an initial value ψ0, the predictand P(ψ0) and “(unknown) observation” P0(ψ0) can be obtained from the model (2) and perfect model (3), respectively. Then, the forecast error is simply E′(ψ0) = P0(ψ0) − P(ψ0). Turning this equation around gives P0(ψ0) = P(ψ0) + E′(ψ0), implying that the model prediction problem boils down to finding E′(ψ0).
Mathematically, the more observational data used for a linear model system, the better are the forecast results to be obtained, just like a pure statistical model. However, this is not the case for a nonlinear model system (such as the atmosphere and oceans) (Chou and Ren 2006). Instead, the use of some optimally selected data with certain relations (e.g., weather patterns or storm tracks) is better than just the use of more observational data. Thus, for a specific forecast, an analog forecast approach could be adopted, in which analogs of the initial value ψ0 from historical observations are first found and then used to make future forecasts (i.e., the ACE approach). Note that in this approach, NWP models play an important role in producing realistic model initial conditions, but are not needed to make forecasts because of adding E′(ψ0), valid at T = 0, to the forecast variables P(ψi). In this sense, the ACE approach should provide more realistic forecasts than those classical dynamical-statistical forecasts, in which all possible historical data are typically used.
However, in estimating the forecast error vector E′(ψ0) associated with the model initial conditions, it is often not possible to obtain a long history of model forecasts for the predictands corresponding to the top m analogs (e.g., due to the continued updates of various model physics and initialization procedures). This makes it difficult to implement the ACE approach in these cases, which represents one of the limitations of the ACE approach.
3. An improved analog prediction approach
a. Improving the ACE-based prediction by taking an ensemble mean
That is, when the analog forecast error is taken as an arithmetic mean solution [i.e.,
It follows from the above analysis that the improved ACE-based prediction is equivalent to the arithmetic mean of an ensemble of m observations (i.e., the exact solutions of a perfect model) that correspond closely to the top m analogs in the historical data. In the next subsection, we will show another form of the ensemble function that can be used to improve the ACE-based prediction with historical observational data by introducing a concept of the perfect model.
b. A dynamical-statistical-analog ensemble forecast (DSAEF) model
As mentioned before, the atmospheric motion fully satisfies the perfect model operator A(ψ) in the dynamical model (3) or (4). Then, a forecast problem is to determine how to obtain the predicted field vector PO(ψO) using the perfect model A from the initial value ψO.

Schematic diagram of the DSAEF model or perfect model, in which ψ is the model state vector, and ψ0, which is mainly composed of observations, is the generalized initial value. It is desirable to combine dynamical forecast products, which show acceptable forecast performance in some key physical variables (e.g., TC track as described in section 4), with initial condition analogs from historical observations. In addition, all the m best analogs to ψ0 and the corresponding m predictands are all from observations.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

Schematic diagram of the DSAEF model or perfect model, in which ψ is the model state vector, and ψ0, which is mainly composed of observations, is the generalized initial value. It is desirable to combine dynamical forecast products, which show acceptable forecast performance in some key physical variables (e.g., TC track as described in section 4), with initial condition analogs from historical observations. In addition, all the m best analogs to ψ0 and the corresponding m predictands are all from observations.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Schematic diagram of the DSAEF model or perfect model, in which ψ is the model state vector, and ψ0, which is mainly composed of observations, is the generalized initial value. It is desirable to combine dynamical forecast products, which show acceptable forecast performance in some key physical variables (e.g., TC track as described in section 4), with initial condition analogs from historical observations. In addition, all the m best analogs to ψ0 and the corresponding m predictands are all from observations.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
For some forecast problems, the ensemble function, F(X1, X2, …, Xm), may also take the functionals of finding the maximum or minimum values of the ensemble members. As an example, we will illustrate in section 4, through an application of the DSAEF model to rainfall forecasts associated with a landfalling TC over South China, that the performance of the ensemble function of “finding the maximum values” is better than that of taking an arithmetic mean in the improved ACE approach.
Figure 2 shows the flowchart of applying the DSAEF model, which consists of three steps: (i) construct the generalized initial value (GIV); (ii) identify analogs; and (iii) produce an ensemble of predictands. In the first step, the GIV includes not only initial conditions at T = 0 but also, if necessary, the subsequent temporal evolution of main variables affecting the predictands in the context of dynamical models, including the perfect model (i.e., observations) and typical NWP models, and statistical methods. Thus, the GIV in four dimensions (4D; i.e., three dimensions at the initial state plus the subsequent temporal evolution) may be composed of all possible variables associated with the predictands, which can be a prognostic meteorological variable such as pressure and horizontal winds, an index such as the East Asian monsoon intensity index or a weather phenomenon such as TC tracks. Both the second and third steps involve the application of statistical methods. Identifying analogs consist of three substeps: (i) Identify analogs for each variable, which can be done after applying a similarity index, following Ren et al. (2018). (ii) Identify analogs for the GIV. Obviously, when there are more than one variable involved in the GIV, a certain order of variables for identifying analogs needs to be considered. (iii) Determine the m best analogs. Generally, the order of variables considered in step (ii) may influence the final solution of the DSAEF. Thus, it is desirable to obtain the best possible order after comparing different ones, and then the m (m ≥ 1) best analogs can be obtained. The third step consists of two parts: identify the observations of predictands for the m top analogs; and select a functional form for an ensemble of the m observations of predictands. After completing the above three steps, a forecast using the DSAEF model can be obtained.

Flowchart of the dynamical-statistical-analog ensemble forecast (DSAEF) model.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

Flowchart of the dynamical-statistical-analog ensemble forecast (DSAEF) model.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Flowchart of the dynamical-statistical-analog ensemble forecast (DSAEF) model.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
One can see from the above description that in the DSAEF model, “dynamics” implies the theoretical basis of the perfect model, while “statistics” includes the generalized initial value, identifying analogs, and optimizing an ensemble of predictands.
c. Comparison with some existing analog-based statistical methods
Given many analog-ensemble-based statistical studies in the literature, as reviewed in section 1, it is of interest to discuss what unique features of the DSAEF model are. First, its theoretical framework originates from the GIV problem of a perfect model, which is the ultimate goal of NWP model development, whereas almost all the existing analog-ensemble approaches do not have their own theoretical frameworks. Moreover, the DSAEF model can be used to not only make analog ensemble forecasts but also show what physical implications of each variable involved are.
In addition, the DSAEF model differs in methodology from some existing analog-ensemble approaches such as the hybrid analog ensemble (AnEn) method (Delle Monache et al. 2013; Nagarajan et al. 2015; Junk et al. 2015; Eckel and Delle Monache 2016). That is, (i) unlike the AnEn method, the DSAEF model requires the construction of a suite of appropriate initial conditions; (ii) the AnEn method needs a long history of deterministic NWP model forecast data from a frozen modeling system, which may not be easy to obtain whereas the DSAEF model searches optimal analogs directly from historical observations (i.e., the so-called generalized initial value); (iii) in the AnEn method analogs are sought independently at each location and for each lead time, whereas in the DSAEF model analogs are sought solely depending on the generalized initial value and for all lead times; and (iv) although both the AnEn method and the DSAEF model obtain ensemble forecasts from historical observations (of predictands) with different schemes, detailed procedures to find those optimal analogs are quite different, as can be seen from an example shown in the next section.
4. An application to landfalling TC precipitation forecasts
a. How to apply the DSAEF model to landfalling TC precipitation forecasts?
In this section, we illustrate the applicability of the DSAEF model to predicting accumulated precipitation associated with landfalling TCs. Figure 3 shows a flowchart on how to apply the DSAEF model to the prediction of landfalling TC precipitation (LTP), hereafter referred to as DSAEF_LTP model, following the procedures listed in Fig. 2. The DSAEF_LTP model involves four steps: obtaining TC track forecast, constructing generalized initial value, identifying analogs, and finding ensemble LTP of the analogs.

Flowchart of the DSAEF model for predicting landfalling TC precipitation (LTP) (abbreviated as DSAEF_LTP model).
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

Flowchart of the DSAEF model for predicting landfalling TC precipitation (LTP) (abbreviated as DSAEF_LTP model).
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Flowchart of the DSAEF model for predicting landfalling TC precipitation (LTP) (abbreviated as DSAEF_LTP model).
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Obtaining accurate TC track forecasts is the first important step toward reasonable forecasts of LTP, since it is a zero-order variable to determine the LTP distribution. After reviewing the previous studies of heavy LTP and LTP forecasts, Ren and Xiang (2017) found that a TC track with respect to different geographical locations is also the key variable in determining the intensity of LTP. Given considerable progress made during the past decades in TC track forecasts (Rappaport et al. 2009; Langmack et al. 2012; Cangialosi and Franklin 2016), it is reliable to use the predicted TC tracks by today’s NWP models as the first step in the DSAEF_LTP model.
In the second step of constructing generalized initial value, it is highly desirable to include all the physical variables that may influence LTP, including both TC’s internal variables (e.g., intensity, size) and environmental variables (e.g., vertical wind shear, subtropical high, summer monsoon), many of which remain to be examined. Their GIV includes not only their initial state but also during the subsequent evolution period. For example, both the observed and predicted tracks for a target TC are treated as the generalized initial value. At this stage of the model development, only two physical variables, TC track and landfalling date, are incorporated into the GIV.
The third step involves identifying m top analogs for the GIV. It consists of three substeps: identifying analogs for each variable, identifying GIV’s analogs, and determining the m best analogs, which are similar to those in the second step of the DSAEF model (Fig. 2). At present, the GIV includes both TC tracks and landfalling dates. The former can be done by applying the TC track similarity area index (TSAI) (Ren et al. 2018), which is an objective index representing an enclosed area bounded by two TC tracks with two lines connecting their initiating and ending points, while the latter is performed by taking seasonal similarities of the TCs landfalling dates. The smaller the TSAI value is, the greater similarities are the two TC tracks, where a null value indicates the complete overlap of the two tracks. Based on the TSAI obtained within a region of concern that is referred to as the similarity region, it is straightforward to obtain the historical landfalling TCs in a similarity order to the target TC. Then, seasonal similarities can be divided into three landfalling types according to the TC landfalling date as detailed in Table 1.
Parameters and the total number of forecast schemes of the DSAEF_LTP model, and the best forecast scheme for experiment 1 and experiment 2, respectively.


In the fourth step, we identify accumulated rainfall amounts for each of the top m analogs, and then select a functional form for an ensemble of the m analogs. The former can be obtained by performing the objective synoptic analysis technique (OSAT) (Ren et al. 2006, 2007), whose function is to partition TC precipitation objectively. The latter can be obtained by taking the two functional forms: ensemble mean and maximum values. Inputting the above information into the DSAEF model would yield the targeted TC’s accumulated rainfall field.
b. An application to LTP forecasts
With the above qualitative description, we may apply the DSAEF_LTP model to the prediction of accumulated rainfall associated with TCs that produced more than 100 mm of daily precipitation at the minimum one rain gauge station during 2012–16 in South China. There are a total of 27 TCs, 17 of which occurring during 2012–14 are used as a training sample, and while the remaining 10 occurring during 2015–16 are used as an independent set of forecast samples for validation. In doing so, historical data (1960–2016) used herein include the best track data from the TC database (http://tcdata.typhoon.org.cn/en/zjljsjj_sm.html) of the China Meteorological Administration (CMA) and the daily rainfall data of 191 rain gauge stations in South China (Fig. 4). Despite the above long historical archival, real-time TC track forecast data of CMA/National Meteorological Center (NMC) for the 27 TCs are also used. In addition, the rainfall forecast data from the European Centre for Medium-Range Weather Forecasts (ECMWF) model, the Global Forecast System (GFS) of the National Centers for Environmental Prediction, and the global spectral model (T639) of CMA/NMC are used herein for the purpose of comparison with the predicted rainfall by the DSAEF_LTP model.

Horizontal distribution of 191 rain gauge stations over South China and tracks of the 21 targeted tropical cyclones adopted for the study. Note that all these TCs have produced the maximum daily rainfall of more than 100 mm over South China, although a few TC tracks do not go through the region.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

Horizontal distribution of 191 rain gauge stations over South China and tracks of the 21 targeted tropical cyclones adopted for the study. Note that all these TCs have produced the maximum daily rainfall of more than 100 mm over South China, although a few TC tracks do not go through the region.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Horizontal distribution of 191 rain gauge stations over South China and tracks of the 21 targeted tropical cyclones adopted for the study. Note that all these TCs have produced the maximum daily rainfall of more than 100 mm over South China, although a few TC tracks do not go through the region.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Figure 5 presents scatterplots of TSs for the training sample of the accumulated rainfall amounts of ≥100 and ≥250 mm from the 5184 forecast schemes associated with experiment 1. Obviously, there are two groups of TSs: one is the ensemble mean (gray dots) and the other is the maximum value ensemble (blue dots) according to the ensemble scheme selected at each station. It is evident from Fig. 5 that the maximum value ensemble is much better than the ensemble mean for intense precipitation forecasts. By comparison, the two highest TS values from the three NWP models, denoted by hollow circles (○), have the values of TS250= 0.0237 by the ECMWF model and TS100= 0.1312 by the GFS model. Clearly, a total of 1048 out of 5814 schemes are better than the NWP forecasts, as they appear in the first quadrant encompassed by the two inner dashed axes. Moreover, the red dot indicates the highest value of TS250 + TS100 that gives the best forecast scheme. See Table 1 for the values of the six parameters associated with the best scheme in this experiment.

Scatterplots of threat scores (denoted by TS250 and TS100, respectively, where the subscripts, “250” and “100”, refer to as the accumulated precipitation of ≥250 and ≥100 mm, respectively) for the 1048 forecast schemes associated with forecast experiment 1 for the training sample of years 2012–13. Gray and blue dots indicate the ensemble “mean” and “maximum,” respectively. The hollow circle (○) is used to denote the LTP forecast TS values by the three NWP models (i.e., ECMWF, GFS, and T639), and the intersect of two inner dashed axes indicates the highest values of TS250 = 0.0237 from the ECMWF model forecasts and TS100 = 0.1312 from the GFS model forecasts.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

Scatterplots of threat scores (denoted by TS250 and TS100, respectively, where the subscripts, “250” and “100”, refer to as the accumulated precipitation of ≥250 and ≥100 mm, respectively) for the 1048 forecast schemes associated with forecast experiment 1 for the training sample of years 2012–13. Gray and blue dots indicate the ensemble “mean” and “maximum,” respectively. The hollow circle (○) is used to denote the LTP forecast TS values by the three NWP models (i.e., ECMWF, GFS, and T639), and the intersect of two inner dashed axes indicates the highest values of TS250 = 0.0237 from the ECMWF model forecasts and TS100 = 0.1312 from the GFS model forecasts.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Scatterplots of threat scores (denoted by TS250 and TS100, respectively, where the subscripts, “250” and “100”, refer to as the accumulated precipitation of ≥250 and ≥100 mm, respectively) for the 1048 forecast schemes associated with forecast experiment 1 for the training sample of years 2012–13. Gray and blue dots indicate the ensemble “mean” and “maximum,” respectively. The hollow circle (○) is used to denote the LTP forecast TS values by the three NWP models (i.e., ECMWF, GFS, and T639), and the intersect of two inner dashed axes indicates the highest values of TS250 = 0.0237 from the ECMWF model forecasts and TS100 = 0.1312 from the GFS model forecasts.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Following the same procedures as described above, we can find the best forecast scheme for experiment 2. It is evident from Fig. 6 that a total of 3087 schemes out of 15 552 ones perform better than the NWP forecasts in predicting the accumulated precipitation of ≥100 and ≥250 mm. Similarly, the best scheme should have the highest TS value of TS250 + TS100 (i.e., as indicated by a red dot). See Table 1 for the values of the seven parameters associated with the best forecast scheme for experiment 2.

As in Fig. 5, but for the 15 552 forecast schemes associated with forecast experiment 2.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

As in Fig. 5, but for the 15 552 forecast schemes associated with forecast experiment 2.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
As in Fig. 5, but for the 15 552 forecast schemes associated with forecast experiment 2.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
A comparison of the parameter values used for the two best schemes between experiments 1 and 2 reveals that (i) five of the six parameters in experiment 1 (i.e., P1, P2, P3, P4, and P7) are the same as those in experiment 2; and (ii) P6, the number of analogs, is 6 and 7 in experiment 1 and 2 (see Table 1), respectively. Furthermore, a comparison of the TSs by the two best schemes of the DSAEF_LTP model between experiments 1 and 2 shows that adding P5 (i.e., seasonal similarity) tends to improve the overall performance of the best forecast scheme. That is, experiment 2 has the summed TS of TS250 +TS100 = 0.2941, where TS250 = 0.0833, and TS250 = 0.2108, while experiment 1 has the summed TS of TS250 +TS100 = 0.26, where TS250 = 0.0662, and TS100 = 0.1938.
Using the best scheme of experiment 2 to represent the DSAEF_LTP model, Fig. 7 compares its mean forecast skills to those of the three NWP models for the independent sample of years 2014–16. Results show that the summed TS (i.e., TS250 + TS100) of the DSAEF_LTP model is comparable to that of the NWP models, with the GFS being the best one (0.3094), followed by the DSAEF_LTP (0.2406), the ECMWF (0.2182), and the T639 (0.2065). In terms of the individual TS250 and TS100, the DSAEF_LTP model is ranked the second and third, respectively. Similar forecast skills are true for the best scheme of experiment 1 (not shown).

Comparison of the threat scores for predicting accumulated intense precipitation (≥250 and ≥100 mm) by the DSAEF_LTP model and the three NWP models (i.e., ECMWF, GFS, and T639) for the independent sample of years 2014–16.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

Comparison of the threat scores for predicting accumulated intense precipitation (≥250 and ≥100 mm) by the DSAEF_LTP model and the three NWP models (i.e., ECMWF, GFS, and T639) for the independent sample of years 2014–16.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Comparison of the threat scores for predicting accumulated intense precipitation (≥250 and ≥100 mm) by the DSAEF_LTP model and the three NWP models (i.e., ECMWF, GFS, and T639) for the independent sample of years 2014–16.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Figure 8 compares the individual TS250 associated with the 9 targeted TCs during 2014–16 from the four different models. Note that the three NWP models have TS250 ≥ 0 only for the most intense rain-producing TC(1410), indicating that they tend to underpredict intense rainfall due partly to the use of coarse grid resolution. Meanwhile, the DSAEF_LTP model has TS250 ≥ 0 for three cases: TC(1410), TC(1418), and TC(1522). In this case, why is the DSAEF_LTP model not ranked the first? A further analysis of the results reveals that our DSAEF_LTP model tends to overpredict the coverage of more intense precipitation (i.e., ≥250 mm), namely giving higher false alarms (i.e., 9) than those of the NWP models (i.e., 5), as shown by thin color bars on the horizontal axis in Fig. 8. As a result, the mean TS for the DSAEF_LTP is smaller than those of the GFS and T639 shown in Fig. 7. The false alarms with TS = 0, causing the lower ranking of the DSAEF_LTP model as calculated by (12), can be attributed partly to the ensemble prediction scheme P7 being “the maximum”.

Comparison of the threat scores of the accumulated precipitation of ≥250 mm associated with 9 targeted TCs during 2014–16 by the DSAEF_LTP model (red column) to the ECMWF (green column), GFS (purple column), and T639 (blue column). The horizontal axis stands for TS = 0, black solid line represents the observed maximum accumulated precipitation for individual TCs, and thin-dashed line denotes the threshold of 250 mm. Note the notations used for the chronological order of TCs: for example, 1410 denotes the 10th TC in 2014, similarly for the other 8 TCs. Note also that false alarms with TS = 0 are represented with thin color bars on the horizontal axis.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

Comparison of the threat scores of the accumulated precipitation of ≥250 mm associated with 9 targeted TCs during 2014–16 by the DSAEF_LTP model (red column) to the ECMWF (green column), GFS (purple column), and T639 (blue column). The horizontal axis stands for TS = 0, black solid line represents the observed maximum accumulated precipitation for individual TCs, and thin-dashed line denotes the threshold of 250 mm. Note the notations used for the chronological order of TCs: for example, 1410 denotes the 10th TC in 2014, similarly for the other 8 TCs. Note also that false alarms with TS = 0 are represented with thin color bars on the horizontal axis.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
Comparison of the threat scores of the accumulated precipitation of ≥250 mm associated with 9 targeted TCs during 2014–16 by the DSAEF_LTP model (red column) to the ECMWF (green column), GFS (purple column), and T639 (blue column). The horizontal axis stands for TS = 0, black solid line represents the observed maximum accumulated precipitation for individual TCs, and thin-dashed line denotes the threshold of 250 mm. Note the notations used for the chronological order of TCs: for example, 1410 denotes the 10th TC in 2014, similarly for the other 8 TCs. Note also that false alarms with TS = 0 are represented with thin color bars on the horizontal axis.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
The TS100 from the four models are compared in Fig. 9, revealing that all the four models have good average TS100 values (i.e., 0.23–0.40) for intense rain-producing TCs (i.e., with precipitation ≥250 mm), whereas much smaller values (i.e., 0.056–0.094) are obtained for TCs with precipitation between 100 and 250 mm.

As in Fig. 8, but for the accumulated precipitation of ≥100 mm.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1

As in Fig. 8, but for the accumulated precipitation of ≥100 mm.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
As in Fig. 8, but for the accumulated precipitation of ≥100 mm.
Citation: Monthly Weather Review 148, 4; 10.1175/MWR-D-19-0174.1
It is apparent from the above analyses that the DSAEF_LTP model can provide numerous possible schemes for a given set of parameters and available historical data, and that the best scheme could be obtained through a training sample. This differs from previous statistical or analog studies, in which only one forecast scheme could be obtained. In particular, in most climate prediction studies and the model output statistics (MOS) for objective weather forecasts, a training test needs to be conducted in order to determine a statistical relationship in the form of a regression equation between a predictand and certain variables. In this sense, the best scheme in our model is equivalent to “the regression equation.”
5. Summary and conclusions
In this study, we have (i) presented an innovative way to improve the prediction of certain meteorological variables with the ACE approach by taking an arithmetic mean of the observations corresponding to the m best historical analogs; (ii) developed the DSAEF model based on the concepts of a perfect model and its initial perturbations, after considering different ensemble functional forms, and (iii) demonstrated the applicability of the DSAEF model to LTP forecasts of 21 TCs occurring in South China.
Our work suggests that although a perfect forecast with today’s dynamical models is impossible to obtain, combining dynamical forecasts with initial condition analogs from historical observations to a maximum possible extent constitutes an innovative way to reduce forecast errors. The DSAEF model developed herein is such a promising approach along this direction, in which analogs from historical observations play an important role in determining forecast accuracy.
Despite the successful application of the DSAEF model, there are several issues on its general applications that are worth discussing. First, it is necessary for the initial perturbations to be as small as possible. This implies that having sufficient analogs from historical observations is critical, and that the more analogs, the better forecasts are. Second, appropriate construction of the GIVs is the key to forecast success with the DSAEF model. A useful strategy is to consider the physical implications of the generalized initial value associated with predictands and the easiness to find their corresponding analogs. Our experience suggests that (i) as many key physical factors associated with the predictands as possible be included as the main body of the GIV; (ii) the GIV be mainly composed of historical observations, and well-predicted key physical variables such as TC tracks; and (iii) the GIV be limited to a (small) key region rather than too large an area. In this regard, Van Den Dool (1994) also pointed out that finding an analog over a small area is much easier than that over a large one. Third, since the quality of identified analogs depend critically on the degree of analogs, a loose analog criterion for one variable may yield a few more analogs than those with a strict one. Thus, an effective strategy is to design different analog criteria for different variables. Fourth, it is always a good idea to test different functional forms for an ensemble of predictands and then choose the best one. Some common functional forms may include an ensemble mean, an extreme (maximum or minimum) value, or a certain percentile.
As demonstrated in section 4, a successful strategy for applying the DSAEF model to LTP forecasts is to include as many relevant variables as possible to the GIV in order to improve forecast skill. For example, the DSAEF_ LTP model with the two variables of “TC track” and “landfalling date” shows better performance than the model with only one variable of “TC track.” In addition, an important step to obtain a successful application of the DSAEF model is to find the best forecast scheme. Thus, it is necessary to figure out how many basic schemes should be and how they could be optimized to yield the best forecast scheme. Consider a case with k “similarity” parameters plus the other ones (e.g., the initial time and the number of analogs). If the DSAEF model has a total of K parameters (K > k) and if the number of values for the jth parameter pj(j = 1, 2, …, K) is nj, then the number of schemes is n1 × n2 × … × nK. To find the best scheme, a reasonable way is to use one of the most important indices to evaluate the forecast performance for each scheme, and then determine the best one based on the evaluations.
It should be mentioned that the present study is just a good start to the application of the DSAEF model to heavy LTP forecasts. There is considerable room for improvements with the DSAEF_LTP model. This could be achieved mainly by two ways: one is incorporating more variables associated with TC internal characteristics (e.g., intensity and size) and environmental conditions (e.g., vertical wind shear, relative humidity, and summer monsoon), another one is improving values of the parameters in the model such as similarity region (P2) and ensemble prediction scheme (P7). The DSAEF model could also be applied to the other meteorological problems (e.g., intense horizontal winds in TCs, regional weather conditions, some of which will be explored in our future studies). In particular, we wish to mention that with its encouraging performance shown herein and our other recent studies (i.e., Ren et al. 2018; Ding et al. 2019; Jia et al. 2020), the DSAEF_LTP model has now been put under operational heavy precipitation foresting tests for landfalling TCs at the CMA/NMC and the Hainan Meteorological Observatory of China. The related results will be reported in future studies.
Acknowledgments
We would like to express sincere thanks to Profs. Russ Elsberry, Jifan Chou, Jianping Huang, Weijing Li, and Guolin Feng for their helpful discussion and comments. This work was supported by the National Natural Science Foundation of China (Grant 41675042), the Hainan Provincial Key R&D Program of China (SQ2019KJHZ0028), the project “dynamical–statistical ensemble technique for predicting landfalling tropical cyclones precipitation,” the National Basic Research Program of China (973 Program: 2014CB441402), and the Jiangsu Collaborative Innovation Center for Climate Change. DLZ was also funded by the U.S. Office of Naval Research Grants N000141410143 and N000141712210.
REFERENCES
Arpe, K., and E. Klinker, 1986: Systematic errors of the ECMWF operational forecasting model in the mid-latitudes. Quart. J. Roy. Meteor. Soc., 112, 181–202, https://doi.org/10.1002/qj.49711247111.
Bao, M., Y. Ni, and J. Chou, 2004: The experiment of monthly mean circulation prediction using the analogue-dynamical model (in Chinese). Chin. Sci. Bull., 49, 1296–1300, https://doi.org/10.1360/03wd0497.
Barnett, T. P., and R. W. Preisendorfer, 1978: Multifield analog prediction of short-term climate fluctuations using a climate state vector. J. Atmos. Sci., 35, 1771–1787, https://doi.org/10.1175/1520-0469(1978)035<1771:MAPOST>2.0.CO;2.
Bauer, P., A. Thorpe, and G. Brunet, 2015: The quiet revolution of numerical weather prediction. Nature, 525, 47–55, https://doi.org/10.1038/nature14956.
Cangialosi, J., and J. Franklin, 2016: National Hurricane Center forecast verification report: 2015 hurricane season. NOAA/NWS/NCEP/National Hurricane Center, 69 pp., https://www.nhc.noaa.gov/verification/pdfs/Verification_2015.pdf.
Charney, J., M. Halem, and R. Jastrow, 1969: Use of incomplete historical data to infer the present state of the atmosphere. J. Atmos. Sci., 26, 1160–1163, https://doi.org/10.1175/1520-0469(1969)026<1160:UOIHDT>2.0.CO;2.
Chen, H., and Z. Lin, 2006: A correction method suitable for dynamical seasonal prediction. Adv. Atmos. Sci., 23, 425–430, https://doi.org/10.1007/s00376-006-0425-3.
Cheng, Y., H. Ren, and G. Tan, 2016: Empirical orthogonal function-analogue correction of extra-seasonal dynamical prediction of East-Asian Summer Monsoon (in Chinese). J. Appl. Meteor. Sci., 27, 285–292.
Chou, J., 1974: A problem of using past data in numerical weather forecasting. Sci. China Ser., 17A, 108–119.
Chou, J., 1979: Some Problems of Long-Term Numerical Weather Prediction. Water Resources and Electric Power Press, 277 pp.
Chou, J., 1986: Why do the dynamical models and statistical methods need to combine?—Also on how to combine (in Chinese). Plateau Meteor., 5, 77–82.
Chou, J., and H. Ren, 2006: Numerical weather prediction—Necessity and feasibility of an alternative methodology (in Chinese). J. Appl. Meteor., 27, 240–244.
Danforth, C. M., and E. Kalnay, 2008a: Impact of online empirical model correction on nonlinear error growth. Geophys. Res. Lett., 35, L24805, https://doi.org/10.1029/2008GL036239.
Danforth, C. M., and E. Kalnay, 2008b: Using singular value decomposition to parameterize state-dependent model errors. J. Atmos. Sci., 65, 1467–1478, https://doi.org/10.1175/2007JAS2419.1.
Danforth, C. M., E. Kalnay, and T. Miyoshi, 2007: Estimating and correcting global weather model error. Mon. Wea. Rev., 135, 281–299, https://doi.org/10.1175/MWR3289.1.
Delle Monache, L., F. A. Eckel, D. L. Rife, and B. Nagarajan, 2013: Probabilistic weather prediction with an analog ensemble. Mon. Wea. Rev., 141, 3498–3516, https://doi.org/10.1175/MWR-D-12-00281.1.
DelSole, T., and A. Y. Hou, 1999: Empirical correction of a dynamical model. Part I: Fundamental issues. Mon. Wea. Rev., 127, 2533–2545, https://doi.org/10.1175/1520-0493(1999)127<2533:ECOADM>2.0.CO;2.
DelSole, T., M. Zhao, P. A. Dirmeyer, and B. P. Kirtman, 2008: Empirical correction of a coupled land–atmosphere model. Mon. Wea. Rev., 136, 4063–4076, https://doi.org/10.1175/2008MWR2344.1.
Ding, C., and Coauthors, 2019: A dynamical statistical ensemble forecast model for landfalling tropical cyclone precipitation based on track similarity (in Chinese). Meteor. Mon., 45, 29–37.
Diomede, T., C. Marsigli, A. Montani, F. Nerozzi, and T. Paccagnella, 2014: Calibration of limited-area ensemble precipitation forecasts for hydrological predictions. Mon. Wea. Rev., 142, 2176–2197, https://doi.org/10.1175/MWR-D-13-00071.1.
Duan, H., G. Tan, and H.-L. Ren, 2011: Statistical correction for model prediction on winter circulation in the extra-tropical of northern hemisphere. Int. Conf. on Information Science and Technology, Nanjing, China, IEEE, 1367–1370.
D’Andrea, F., and R. Vautard, 2000: Reducing systematic errors by empirically correcting model errors. Tellus, 52A, 21–41, https://doi.org/10.1034/j.1600-0870.2000.520103.x.
Eckel, F. A., and L. Delle Monache, 2016: A hybrid NWP–analog ensemble. Mon. Wea. Rev., 144, 897–911, https://doi.org/10.1175/MWR-D-15-0096.1.
Elsberry, R. L., and H. C. Tsai, 2014: Situation-dependent intensity skill metric and intensity spread guidance for western North Pacific tropical cyclones. Asia-Pac. J. Atmos. Sci., 50, 297–306, https://doi.org/10.1007/s13143-014-0018-5.
Feddersen, H., A. Navarra, and M. N. Ward, 1999: Reduction of model systematic error by statistical correction for dynamical seasonal predictions. J. Climate, 12, 1974–1989, https://doi.org/10.1175/1520-0442(1999)012<1974:ROMSEB>2.0.CO;2.
Frediani, M. E. B., T. M. Hopson, J. P. Hacker, E. N. Anagnostou, L. Delle Monache, and F. Vandenberghe, 2017: Object-based analog forecasts for surface wind speed. Mon. Wea. Rev., 145, 5083–5102, https://doi.org/10.1175/MWR-D-17-0012.1.
Gao, L., H. Ren, J. Li, and J. Chou, 2006: Analogue correction method of errors and its application to numerical weather prediction. Chin. Phys., 15, 882–889, https://doi.org/10.1088/1009-1963/15/4/038.
Glahn, H. R., and D. A. Lowry, 1972: The use of model output statistics (MOS) in objective weather forecasting. J. Appl. Meteor., 11, 1203–1211, https://doi.org/10.1175/1520-0450(1972)011<1203:TUOMOS>2.0.CO;2.
Huang, J., Y. Yi, S. Wang, and J. Chou, 1993: An analogue-dynamical long-range numerical weather prediction system incorporating historical evolution. Quart. J. Roy. Meteor. Soc., 119, 547–565, https://doi.org/10.1002/qj.49711951111.
Jia, Z., F. Ren, D.-L. Zhang, C. Ding, M.-J. Yang, T. Feng, B. Chen, and H. Yang, 2020: An application of the LTP_DSEF model to heavy precipitation forecasts of landfalling tropical cyclones over China in 2018. Sci. China Earth Sci., 63, 27–36, https://doi.org/10.1007/s11430-019-9390-6.
Johansson, A., and S. Saha, 1989: Simulation of systematic error effects and their reduction in a simple model of the atmosphere. Mon. Wea. Rev., 117, 1658–1675, https://doi.org/10.1175/1520-0493(1989)117<1658:SOSEEA>2.0.CO;2.
Junk, C., L. Delle Monache, and S. Alessandrini, 2015: Analog-based ensemble model output statistics. Mon. Wea. Rev., 143, 2909–2917, https://doi.org/10.1175/MWR-D-15-0095.1.
Kaas, E., A. Guldberg, W. May, and M. Déqué, 1999: Using tendency errors to tune the parameterization of unresolved dynamical scale interactions in atmospheric general circulation models. Tellus, 51A, 612–629, https://doi.org/10.3402/tellusa.v51i5.14481.
Kalnay, E., 2002: Atmospheric Modeling, Data Assimilation and Predictability. Cambridge University Press, 364 pp.
Klein, W. H., and H. R. Glahn, 1974: Forecasting local weather by means of model output statistics. Bull. Amer. Meteor. Soc., 55, 1217–1227, https://doi.org/10.1175/1520-0477(1974)055<1217:FLWBMO>2.0.CO;2.
Klinker, E., and P. D. Sardeshmukh, 1992: The diagnosis of mechanical dissipation in the atmosphere from large-scale balance requirements. J. Atmos. Sci., 49, 608–627, https://doi.org/10.1175/1520-0469(1992)049<0608:TDOMDI>2.0.CO;2.
Koo, C.-C., 1958a: On the equivalency of formulations of weather forecasting as an initial value problem and as an “evolution” problem (in Chinese). Acta Meteor. Sin., 29, 93–98.
Koo, C.-C., 1958b: On the utilization of past data in numerical weather forecasting (in Chinese). Acta Meteor. Sin., 29, 176–184.
Langmack, H., K. Fraedrich, and F. Sielmann, 2012: Tropical cyclone track analog ensemble forecasting in the extended Australian basin: NWP combinations. Quart. J. Roy. Meteor. Soc., 138, 1828–1838, https://doi.org/10.1002/qj.1915.
Leith, C., 1978: Objective methods for weather prediction. Annu. Rev. Fluid. Mech., 10, 107–128, https://doi.org/10.1146/annurev.fl.10.010178.000543.
Liu, Y., and H.-L. Ren, 2017: Improving ENSO prediction in CFSv2 with an analogue-based correction method. Int. J. Climatol., 37, 5035–5046, https://doi.org/10.1002/joc.5142.
Livezey, R. E., A. G. Barnston, G. V. Gruza, and E. Y. Ran’kova, 1994: Comparative skill of two analog seasonal temperature prediction systems: Objective selection of predictors. J. Climate, 7, 608–615, https://doi.org/10.1175/1520-0442(1994)007<0608:CSOTAS>2.0.CO;2.
Lorenc, A., 1986: Analysis methods for numerical weather prediction. Quart. J. Roy. Meteor. Soc., 112, 1177–1194, https://doi.org/10.1002/qj.49711247414.
Lorenz, E. N., 1969: Atmospheric predictability as revealed by naturally occurring analogues. J. Atmos. Sci., 26, 636–646, https://doi.org/10.1175/1520-0469(1969)26<636:APARBN>2.0.CO;2.
Lorenz, E. N., 1977: An experiment in nonlinear statistical weather forecasting. Mon. Wea. Rev., 105, 590–602, https://doi.org/10.1175/1520-0493(1977)105<0590:AEINSW>2.0.CO;2.
Nagarajan, B., L. Delle Monache, J. P. Hacker, D. L. Rife, K. Searight, J. C. Knievel, and T. N. Nipen, 2015: An evaluation of analog-based postprocessing methods across several variables and forecast models. Wea. Forecasting, 30, 1623–1643, https://doi.org/10.1175/WAF-D-14-00081.1.
Qiu, C., and J. Chou, 1989: An analogue-dynamical method of weather forecasting (in Chinese). Chin. J. Atmos. Sci., 13, 22–28.
Rappaport, E. N., and Coauthors, 2009: Advances and challenges at the National Hurricane Center. Wea. Forecasting, 24, 395–419, https://doi.org/10.1175/2008WAF2222128.1.
Ren, F., and C. Xiang, 2017: Review and prospect of researches on the causes of intense precipitation and the prediction of precipitation associated with landfalling tropical cyclones (in Chinese). J. Mar. Meteor., 37, 8–18.
Ren, F., G. Wu, W. Dong, X. Wang, Y. Wang, W. Ai, and W. Li, 2006: Changes in tropical cyclone precipitation over China. Geophys. Res. Lett., 33, L20702, https://doi.org/10.1029/2006GL027951.
Ren, F., Y. Wang, X. Wang, and W. Li, 2007: Estimating tropical cyclone precipitation from station observations. Adv. Atmos. Sci., 24, 700–711, https://doi.org/10.1007/s00376-007-0700-y.
Ren, F., W. Qiu, C. Ding, X. Jiang, L. Wu, Y. Xu, and Y. Duan, 2018: An objective track similarity index and its preliminary application to predicting precipitation of landfalling tropical cyclones. Wea. Forecasting, 33, 1725–1742, https://doi.org/10.1175/WAF-D-18-0007.1.
Ren, H., and J. Chou, 2006: Analogue correction method of errors by combining statistical and dynamical methods. Acta Meteor., 20, 367–373.
Ren, H., and J. Chou, 2007: Strategy and methodology of dynamical analogue prediction. Sci. China Ser. D Earth Sci., 50, 1589–1599, https://doi.org/10.1007/s11430-007-0109-6.
Ren, H., P. Zhang, W. Li, and J. Chou, 2006: A new method of dynamical analogue prediction based on multi-reference-state updating and its application. Acta Phys. Sin., 55, 4388–4396.
Ren, H., J. Chou, J. Huang, and P. Zhang, 2009: Theoretical basis and application of analogue-dynamical model in the Lorenz system. Adv. Atmos. Sci., 26, 67–77, https://doi.org/10.1007/s00376-009-0067-3.
Ren, H., L. Ying, J. Fei-Fei, Y. Yu-Ping, and L. Xiang-Wen, 2014a: Application of the analogue-based correction of errors method in ENSO prediction. Atmos. Oceanic Sci. Lett., 7, 157–161, https://doi.org/10.3878/j.issn.1674-2834.13.0080.
Ren, H., P. Zhang, W. Li, and L. Chen, 2014b: The dynamical-analogue ensemble method for improving operational monthly forecasting. Acta Meteor. Sin., 72, 723–730.
Roebber, P. J., and L. F. Bosart, 1998: The sensitivity of precipitation to circulation details. Part I: An analysis of regional analogues. Mon. Wea. Rev., 126, 437–455, https://doi.org/10.1175/1520-0493(1998)126<0437:TSOPTC>2.0.CO;2.
Saha, S., 1992: Response of the NMC MRF model to systematic error correction within integration. Mon. Wea. Rev., 120, 345–360, https://doi.org/10.1175/1520-0493(1992)120<0345:ROTNMM>2.0.CO;2.
Sun, C., L. Wei-Jing, R. Hong-Li, Z. Pei-Qun, W. Dong-Yan, 2006: A dynamic analogue error correction model for ENSO prediction and its initial hindcast verification (in Chinese). Chinese J. Atmos. Sci., 30, 965–976.
Tan, G., H. Duan, and H. Ren, 2012: Statistical correction for dynamical prediction of 500 hPa height field in mid-high latitudes. J. Appl. Meteor. Sci., 23, 304–311.
Tsai, H. C., and R. L. Elsberry, 2015: Seven-day intensity and intensity spread predictions for western North Pacific tropical cyclones. Asia-Pac. J. Atmos. Sci., 51, 331–342, https://doi.org/10.1007/s13143-015-0082-5.
Tsai, H. C., and R. L. Elsberry, 2016: Skill of western North Pacific tropical cyclone intensity forecast guidance relative to weighted-analog technique. Asia-Pac. J. Atmos. Sci., 52, 281–290, https://doi.org/10.1007/s13143-016-0001-4.
Tsai, H. C., and R. L. Elsberry, 2017: Seven-day intensity and intensity spread predictions for Atlantic tropical cyclones. Wea. Forecasting, 32, 141–147, https://doi.org/10.1175/WAF-D-16-0165.1.
Van Den Dool, H. M., 1994: Searching for analogues, how long must we wait? Tellus, 46A, 314–324, https://doi.org/10.3402/tellusa.v46i3.15481.
Yang, X., T. DelSole, and H.-L. Pan, 2008: Empirical correction of the NCEP global forecast system. Mon. Wea. Rev., 136, 5224–5233, https://doi.org/10.1175/2008MWR2527.1.
Yang, X.-Q., and J. L. Anderson, 2000: Correction of systematic errors in coupled GCM forecasts. J. Climate, 13, 2072–2085, https://doi.org/10.1175/1520-0442(2000)013<2072:COSEIC>2.0.CO;2.
Yu, H., J. Huang, and J. Chou, 2014: Improvement of medium-range forecasts using the analog-dynamical mthod. Mon. Wea. Rev., 142, 1570–1587, https://doi.org/10.1175/MWR-D-13-00250.1.
Zheng, Z., H. Ren, and J. Huang, 2009: Analogue correction of errors based on seasonal climatic predictable components and numerical experiments (in Chinese). Acta Phys. Sin., 58, 7359–7367.
Zhou, B., and P. Zhai, 2016: A new forecast model based on the analog method for persistent extreme precipitation. Wea. Forecasting, 31, 1325–1341, https://doi.org/10.1175/WAF-D-15-0174.1.