Research article 13 Sep 2018
Research article | 13 Sep 2018
In this paper, a novel approach is proposed for solving conditional nonlinear optimal perturbations (CNOPs), called the adaptive cooperative coevolution of parallel particle swarm optimization (PSO) and the Wolf Search algorithm (WSA) based on principal component analysis (ACPW). Taking Fitow (2013) and Matmo (2014), two tropical cyclone (TC) cases, CNOPs solved by the ACPW algorithm are used to investigate the sensitive regions identified by TC adaptive observations with the fifth-generation Mesoscale Model (MM5). Meanwhile, the 60 and 120 km resolutions are adopted. The adjoint-based method (short for the ADJ method) is also applied to solve CNOPs, and the result is used as a benchmark. To evaluate the advantages of the ACPW algorithm, we run the PSO, WSA and ACPW programs 10 times and then compare the maximum, minimum and mean objective values as well as the RMSEs. The analysis results prove that the hybrid strategy and cooperative coevolution are useful and effective. To validate the ACPW algorithm, the CNOPs obtained from the different methods are compared in terms of the patterns, energies, similarities and simulated TC tracks with perturbations. The results of our study may be summarized as follows:
The ACPW algorithm can capture similar CNOP patterns as the ADJ method, and the patterns of TC Fitow are more similar than TC Matmo.
At the 120 km resolution, similarities between the CNOPs of the ADJ method and the ACPW algorithm are more than those at the 60 km resolution.
Compared to the ADJ method, although the CNOPs of the ACPW method produce lower energies, they can have improved benefits gained from the reduction of the CNOPs not only across the entire domain but also in the identified sensitive regions.
The sensitive regions identified by the CNOPs from the ACPW algorithm have the same influence on the improvements of the skill of TC-track forecasting as those identified by the CNOPs from the ADJ method.
The ACPW method is more efficient than the ADJ method. All conclusions prove that the ACPW algorithm is a meaningful and effective method for solving CNOPs and can be used to identify sensitive regions of TC adaptive observations.
Tropical cyclones (TCs) are one of the most frequent and influential natural hazards in the world. An accurate forecast of TCs is conducive to the response of the government and people. Thus, it is essential to improve the skill of TC forecasting. One effective way is to identify the sensitive regions of TC adaptive observations (TCAOs) (Franklin and Demaria, 1992; Bergot, 1999; Aberson, 2003). Once observations in sensitive regions are identified and added to reduce initial errors, better forecasts will be expected (Bender et al., 1993; Zhu and Thorpe, 2006; Froude et al., 2007). Conditional nonlinear optimal perturbations (CNOPs) proposed by Mu and Duan (2003) are a nonlinear extension of the linear singular vector (SV) method and have been applied to study the successful identification of sensitive regions by TCAOs (Mu et al., 2009; Qin, 2010; Zhou and Mu, 2011, 2012a, b; Zhou and Zhang, 2014; Qin and Mu, 2012; Qin et al., 2013; Qin and Mu, 2014; Wang et al., 2010, 2013).
Comparing between the sensitive regions identified from CNOPs and those identified through SVs, Qin (2010) concludes that the former is more appropriate for TCAOs. Zhou and Mu (2011) use the CNOP method to investigate different verification areas and how to affect the identification of sensitive regions. They also studied the influence of different horizontal resolutions (2012a). Moreover, different times and regime dependency were also researched (2012b). These research results directed further research. Zhou and Zhang (2014) propose three schemes for identifying sensitive regions based on the CNOP method and recommend the vertically integrated energy scheme. Moreover, some researchers analyze the sensitivity of dropwindsonde observations on TC predictions, which can be used in the CNOP method, and conclude that the sensitive regions identified by CNOPs have a positive impact on TC-track predictions (Qin and Mu, 2012; Qin et al., 2013). In studies of improving the sensitivity of CNOPs in TC intensity forecasts, Qin and Mu (2014) suggest that the use of an ocean-coupled model needs to be considered as well as the better initialization of the TC vortex. Wang et al. (2013) use the CNOP method to study the mutual effects of binary typhoons. Previous studies have shown that the CNOP method is a useful and meaningful method for studying the aforementioned phenomenon (Zhou et al., 2013; Mu and Zhou, 2015).
There are generally two types of methods for solving CNOPs, one based on adjoint models (ADJ method) and one without adjoint models. As useful and effective methods for solving CNOPs without adjoint models, some modified intelligent algorithms (IAs) based on dimension reduction have been successfully proposed and applied to solve CNOPs in the Zebiak–Cane (ZC) (Zebiak and Cane, 1987) model, such as SAEP (simulated annealing based ensemble projecting method) (Wen et al., 2014), PPSO (principal component analysis-based particle swarm optimization – Mu et al., 2015a; principal component analysis, PCA – Jolliffe, 1986), PCGD (principal component-based great deluge) (Wen et al., 2015a), RGA (robust PCA-based genetic algorithm) (Wen et al., 2015b), CTS-SS (continuous Tabu search algorithm with sine maps and staged strategy) (Yuan et al., 2015) and PCAGA (principal component analysis-based genetic algorithm) (Mu et al., 2015b). Compared to the ADJ method, these methods all obtain CNOPs with similar spatial patterns and acceptable objective function values. Several of them have been parallelized with the message passing interface (MPI), reducing the computation time. In TC adaptive observations, such adjoint-free methods are also required because the lack of adjoint models and solution spaces with too many dimensions have become obstacles for solving CNOPs; this is a focal point of this study.
We have adopted the PCAGA method to solve CNOPs for the sensitive regions identified by TCAOs with the fifth-generation Mesoscale Model (MM5) and obtained meaningful results (Zhang et al., 2017). However, we used a resolution of 120 km, which is the lowest in such research. When using a higher resolution, information on a smaller scale can be predicted and more accurate sensitive regions can be expected. It is necessary to use a higher resolution. Moreover, although the PCAGA method achieves meaningful results, its performance is not sufficient because it is based on a genetic algorithm, which has a good global searching ability but a slow convergence rate. In addition, the PCAGA method was not parallelized in the previous study.
Therefore, in this paper, we propose a novel approach, the adaptive cooperative coevolution of parallel particle swarm optimization (PSO) and Wolf Search algorithm (WSA) (ACPW) based on the PCA to solve CNOPs for the sensitive regions identified by TCAOs. We take two tropical cyclones as study cases, Fitow (2013) and Matmo (2014), and simulate them with the MM5 model using two different resolutions, 60 and 120 km. According to the study of Zhou and Zhang (2014), we adopt the total dry energy as the objective function. The CNOPs from the ADJ method are referred to as a benchmark. Specific details of the ADJ method can be found in Zhou (2009). To validate the ACPW method, the CNOPs from the ACPW method are compared with the benchmark in terms of the patterns, energies, similarities and benefits from the CNOPs reduced in the entire domain and in sensitive regions. Further, the CNOPs with different resolutions are also compared in terms of these aspects. To evaluate the sensitive regions located by the ACPW algorithm, we simulate TC tracks with the initial states perturbed by the amended CNOPs in the location of the sensitive regions from the ACPW algorithm and ADJ method. Moreover, we design two schemes to amend the CNOPs using the same points and the equivalent proportional points. In addition, we evaluate the efficiency of the ACPW algorithm. All experimental results show that the ACPW method is a meaningful and effective method to solve CNOPs for selecting the sensitive regions of TCAOs.
The organization of the paper is as follows. Section 2 describes the formalized definition of CNOPs and the ACPW method. In Sect. 3, we give the design of the experiments in this study. Section 4 presents the experimental analysis and results. Summaries and conclusions are provided in Sect. 5.
The mathematical formalism of CNOPs is described in Eq. (1). Under the constraint condition ${\u2225{u}_{\mathrm{0}}\u2225}^{\mathrm{2}}\le \mathit{\delta}$, an initial perturbation $\mathit{\delta}{u}_{\mathrm{0}}^{\ast}$ of vector U_{0} (initial basic state) is called a CNOP if and only if
where
and P represents a local projection operator and the value within the verification region is 1 and 0 elsewhere. In addition,
where M expresses a nonlinear propagation operator, and U_{t} is the development of U_{0} at time t.
In this paper, we propose the ACPW method to solve CNOPs for identifying sensitive regions of TCAOs. The core of this approach is the cooperative coevolution of two intelligent algorithms, the PSO and WSA, and the adaptive number of two sub-swarms. PSO is a classic population-based stochastic optimization technique developed by Kennedy and Eberhart (1995) and inspired by the social behaviors of bird flocking or fish schooling. The technique has been successfully and effectively applied to solve CNOPs in the ZC model for studying El Niño–Southern Oscillation (ENSO) predictions (Mu et al., 2015a). The WSA is a new bio-inspired heuristic optimization algorithm based on wolf preying behaviors, which was proposed by Tang et al. (2012) and has been applied to studying the traveling salesman problem with test functions. Their experiments showed that the WSA is an effective global optimizing algorithm but requires long computation times.
We have adopted the PSO and WSA methods to solve CNOPs in the MM5 model, although the results exhibit slow convergence or premature convergence. Hence, we combine the advantages of these two algorithms. We use the WSA to explore the global space due to its independence and use PSO to examine the local space and ensure the convergence of the ACPW algorithm. Moreover, we design the adaptive sub-swarms of the PSO and WSA for cooperative coevolution. The ACPW framework is shown in Fig. 1.
In Fig. 1, the most important part of the ACPW algorithm is inside the dotted box. We divide the entire initial swarm into two sub-swarms with the same number of individuals; one updates the individuals with the PSO's rule and the other with the WSA's rule. Then, the two sub-swarms are adaptively varied along with the convergence state of the ACPW algorithm. When the change in the objective function adaptive value is less than a threshold value, the number of individuals in the sub-swarm belonging to the WSA is increased and the other sub-swarm belonging to PSO is decreased by an equal number of individuals to keep the same number for the entire swarm. A more specific analysis of the ACPW algorithm is discussed in Sect. 4.
The process of solving CNOPs with the ACPW algorithm is described as follows:
Randomly generate an initial swarm with N individuals. An individual u_{i} needs to satisfy the boundary constraint in the terms of Eq. (4). Once u_{i} goes beyond the boundary, it must be pulled back, i.e.,
Divide the entire initial swarm into two sub-swarms with an adaptive coefficient α. One sub-swarm updates individuals with the PSO's rule and the other with the WSA's rule.
Calculate the adaptive value of the objective function in parallel, i.e.,= J(u_{i}) in Eq. (1).
Update individuals by the PSO (Eq. 5) or the WSA (Eq. 6). When
the superscript k or k+1 is the iterative step, ${v}_{i}^{k+\mathrm{1}}$ is the velocity of the individual ${u}_{i}^{k}$ calculated by the first sub-formula, ω is the inertia coefficient, c_{1} and c_{2} are the learning factors, α and β are the random numbers uniformly distributed on the interval from 0 to 1, ${o}_{i}^{k}$ is the local optimum, ${o}_{\mathrm{g}}^{k}$ is the global optimum in the kth iteration, γ is the restraint factor to control the speed and ${u}_{i}^{k+\mathrm{1}}$ is the updated individual based on PSO.
There are two ways for updating individuals in the WSA, prey and escape, which represent the functions of searching in a local region and escaping from a local optimum. These are represented as
where the superscript k or k+1 is also the iterative step, θ is the velocity, r is the local optimizing radius that is smaller than the global constraint radius δ, rand( ) is the random function whose mean value is distributed in $[-\mathrm{1},\mathrm{1}]$, escape( ) is the function for calculating a random position that is 3 times larger than r and s is the step size of the updating individual.
As described in Eq. (6), the wolf has two behaviors, i.e., prey and escape. The prey behavior uses the first sub-formula, and the second one is for the escape function that happens in every iteration when the condition p>p_{a} is satisfied, where p is a random number in [0,1] and p_{a} is the probability of an individual escaping from the current position.
Judge whether the change in the adaptive value of the objective function is smaller than ε. If so, set a new value for the adaptive sub-swarm coefficient α. If not, continue running the process. The detailed updating procedure for α is described as
In this paper, before we update the individuals, α is calculated and we divide the entire initial swarm into two sub-swarms according to the α value i.e., the number of individuals depending on the PSO's rule is α×N and the other number is $(\mathrm{1}-\mathit{\alpha})\times N$. We set the initial value of ε and α to 0.1 and 0.5, respectively.
Judge whether the termination condition is satisfied. If so, terminate the iteration. Otherwise, go to step 2.
All of the above processes are based on the dimension reduction within the PCA, a procedure that has been described in the study of Mu et al. (2015a). After many experiments, the parameters of the ACPW algorithm can be set, as shown in Table 1.
Although there are more parameters than demanded for each single algorithm, most retain the empirical value of each algorithm and do not require adjustments. The reason for using a different number of individuals is that the internal storage memory was not sufficient when using more than 200 individuals, resulting in the premature termination of the ACPW algorithm.
All the experiments are run on a Lenove Thinkserver RD430 with two Intel Xeon E5-2450 2.10 GHz CPUs, 32 logical cores and 132G RAM. The operating system is CentOS 6.5. All the codes are written in the FORTRAN language and compiled by the PGI Compiler 10.2.
In this paper, we adopt the MM5 model to study the sensitive region identification of TCAOs and the corresponding adjoint system of the MM5 model (Zou et al., 1997) is used to obtain the benchmark. The ERA interim daily analysis data (1^{∘} × 1^{∘}) (Dee et al., 2011) from the European Centre for Medium range Weather Forecasts (ECMWF) are used to generate the initial and boundary conditions. The physical parameterization schemes are defined as dry convective adjustment, the high-resolution planetary boundary-layer scheme, grid-resolved large-scale precipitation, and the Kuo cumulus parameterization scheme.
We also utilize the best TC-track data (Ying et al., 2014) from the China Meteorological Administration–Shanghai Typhoon Institute (CMA–SHTI) as TC tracks observed for evaluating the simulation TC tracks of the MM5 model.
TCs Fitow (2013) and Matom (2014) are taken as the study cases and introduced below. Fitow was the 23rd TC in 2013 and developed to the east of the Philippines on 29 September, striking China at Fuding in Fujian Province on 6 October. Matom was the 10th named typhoon in 2014. It formed on 17 July and reached land in Taiwan on 22 July. In these two cases, 24 h control forecasts are set as background fields based on integration from 00:00 UTC 5 October 2013 to 00:00 UTC 6 October 2013 (TC Fitow) and from 18:00 UTC 21 July 2014 to 18:00 UTC 22 July 2014 (TC Matom). After the 24 h period, TC Fitow had a maximum sustained wind of 162 km h^{−1} whereas TC Matmo had a maximum wind speed of 151.2 km h^{−1}. In addition, the forecasts were executed at the 60 km and 120 km resolutions with 11 vertical levels, and the model domain covered 55 × 55 and 21 × 26 grids, respectively.
The simulated TC tracks from the MM5 model for these two cases are acceptable, as has been shown in our previous study (Zhang et al., 2017). The following analysis is based on those simulations.
Because slight changes in the verification area never hurts the results (Zhou and Mu, 2011), we design the verification areas as rectangles covering the potential typhoon tracks at the forecast time.
The initial perturbation sample δu_{0} is composed of the perturbed zonal wind u_{0}^{′}, meridional wind v_{0}^{′}, temperature ${{T}^{\prime}}_{\mathrm{0}}$ and surface pressure ${{p}^{\prime}}_{{\mathrm{s}}_{\mathrm{0}}}$. Each component can be represented as a matrix $a\phantom{\rule{0.125em}{0ex}}m\times \phantom{\rule{0.125em}{0ex}}n\times \phantom{\rule{0.125em}{0ex}}l$, where m×n is the distribution of the horizontal grid, and l denotes the number of vertical levels. To extract features for reducing the dimensions and solving CNOPs, the $m\times n\times l$ matrix is reshaped to a k×1 vector, where $k=m\times n\times l\times S$ (S is the number of the components). Assuming we have R vectors to represent the features of the solution space, we recombine the R vectors to a k×R matrix and use the PCA to capture the feature space with lower dimensions. Then, the CNOP is solved in the space of the feature until we obtain the global CNOP, which will be projected to the original solution space. When using the ACPW algorithm to solve CNOPs, its initial inputs are produced randomly in the feature space, and the CNOP has the largest nonlinear evolution at the prediction time, i.e., the largest adaptive value of the objective function in Eq. (9). The objective function is measured by the total dry energy (Zhou and Zhang, 2014) since it has been proven that the sensitive regions gained by the dry energy are more beneficial than those obtained from the moist energy (Zhou, 2009).
The following is defined as
where ${E}_{\mathrm{T}}\left(i,j,\mathit{\sigma}\right)$ denotes the total dry energy of the CNOP at the MM5 grid point $\left(i,j,\mathit{\sigma}\right)$.
Corresponding to Formulas (1) and (2), we have
where ${{u}^{\prime}}_{\mathrm{t}}$, ${v}_{\mathrm{t}{}^{\prime}}$, ${{T}^{\prime}}_{\mathrm{t}}$ and $p{{}^{\prime}}_{\mathrm{st}}$ are the components of u_{NT}, which is the nonlinear development of the perturbed U_{0} (i.e., U_{0}+δu_{0}) from the initial time t_{0} to the prediction time t, and σ is the vertical coordinate. Table 2 illustrates the other reference parameters.
For the convenience of optimization, solving CNOPs can be transformed into a minimized problem as
To facilitate understanding, all symbols are listed in Table 2, and their meanings are explained.
To evaluate the advantages of the ACPW algorithm, we run the PSO, WSA and ACPW programs 10 times and then compare the maximum, minimum and mean objective values as well as the RMSE. We also exhibit the objective value scope after the first iteration to analyze the effect of initial objective values on the different algorithms. Meanwhile, to illustrate the performance of the algorithms, we compare the degree of change of the objective function value for the three algorithms.
Because the statistical analysis results are similar for the two TCs with two resolutions, we only describe the analysis of Fitow at a resolution of 60 km. Table 3 presents the maximum objective value, the minimum objective value, the mean objective value and the RMSE of the 10 results.
In Table 3, the maximum objective value is gained from the ACPW algorithm, and its mean value is also more than the other two algorithms. However, the RMSE of the PSO is the smallest, which shows the most stability.
For additional analysis, we draw a box plot of the 10 results for the PSO, WSA and ACPW algorithms, as shown in Fig. 2. PSO has the narrowest range of values, although the objective values are smaller than the other two algorithms. The WSA has the widest range of values, although the objective values are also smaller than the ACPW algorithm. The ACPW algorithm has the second best stability, although it has the best objective values. The experiments display the stability of the PSO and the exploitation of the WSA. We combine the advantages of the PSO and WSA methods and use them to develop the ACPW algorithm to solve CNOPs. The analysis results demonstrate that the hybrid strategy and cooperative coevolution is both useful and effective.
Since these three algorithms are all heuristic algorithms generated randomly and the initial inputs are also generated by random way, the initial objective value is different for every run. To analyze the effect of initial objective values on the different algorithms, we exhibit the objective value scope of the PSO, WSA and ACPW algorithms after the first iteration in Fig. 3.
In Fig. 3, for convenience, only the integer is indicated in the coordinate system. In 10 experiments, the PSO has the narrowest scope, from 467.1719 to 781.6482. The WSA and ACPW algorithms have similar value spans that are wider than the PSO, but the objective values of the ACPW are higher. And the value scope is reasonable according to the characteristics of these three algorithms. The WSA is the most random, the PSO is the most stable and the ACPW combines the advantages of the two. From the results, we cannot find the direct relationship between the initial objective value and the final results, but a better first objective value is beneficial in finding the optimal value.
To illustrate the improved performance of the ACPW algorithm, we calculate the average objective value of every step in 10 program results and obtain the change degree between the two iterations. We draw them in Fig. 4. If the objective value is continuously changing, then the algorithm has better global searching ability. Otherwise, the algorithm tends to experience a drop in local optimization.
In Fig. 4, the degree of change is calculated from the subtraction of two objective values. For example, the objective value of the second iteration minus the first objective value is the first degree of change m has better performance than the PSO and WSA, because we combine their strengths using hybrid strategy and cooperative coevolution.
To validate the ACPW algorithm for solving CNOPs and to identify the sensitive regions, we compare the ADJ method and the ACPW algorithm results in terms of the CNOP patterns, energies, similarities, benefits from reduction of the CNOPs and simulated TC tracks with perturbations.
In this subsection, we compare the CNOPs obtained from the ADJ method and the ACPW algorithm in terms of the patterns of temperature and wind. Experimental results show that TC Fitow has more similar CNOP patterns than TC Matmo. The CNOP patterns are described in Fig. 5.
At the 120 km resolution for TC Fitow (Fig. 5a, b), the two methods have nearly the same major warm locations and similar cold regions, while the wind vectors have opposite directions. The ADJ method captures the CNOP with two major locations. The red (warm) location is distributed to the west of the initial cyclone (IC), while the green (cold) location is distributed to the north of the IC. The ACPW algorithm also captures the CNOP with two main locations. The warm one is distributed to the west and the cold one is located to the northwest of the IC. In this subsection, the spatial orientation is relative to the position of the IC. Therefore, in the following discussion, we explain the spatial orientation in the figures without repeating the IC.
For the TC Fitow analysis with a 60 km resolution (Fig. 5c, d), the CNOP spatial distribution based on the ACPW algorithm is very similar to the ADJ method's results. In the northwest of the verification area, the two CNOPs have two similar major parts, a warm area and a cold area. The difference between these two patterns is that the ADJ method has another major warm area located in the northwest, while the ACPW method produces another major warm area in the east. The distribution of the secondary parts exhibits only a slight difference.
For the same method with different resolutions (Fig. 5a, c and b, d), the CNOP patterns have similar major distributions in the northwest, although these occur within a different region. The reason is that when using a higher resolution, more small-scale phenomena can be resolved (Zhou and Mu, 2012a).
For the analysis of TC Matmo with a 120 km resolution (Fig. 6a, b), the ADJ method and the ACPW algorithm obtain CNOPs with different spatial patterns in terms of temperature and wind. The ADJ method has two major parts, with the warm part located in the west and the cold one in the east. The ACPW algorithm results in two main parts distributed in the northeast, with a warm area near the IC and a cold one far from the IC. For the analysis of TC Matmo with a 60 km resolution (Fig. 6c, d), in the verification area, the two CNOP patterns have similar spatial distributions, with two warm areas located at nearly the same positions. However, the parts outside the verification area are distributed in different locations. Moreover, the CNOP of the ADJ method has more regular distributions than the ACPW's distributions. For the same method with a different resolution (Fig. 6a, c and b, d), the CNOP patterns cover similar areas but with different ranges and details.
Based on the above analysis regarding the patterns of temperature and wind, we can conclude that when using a resolution of 60 km, the CNOPs predicted by the ADJ method and the ACPW algorithm have more similar major patterns than those predicted at a resolution of 120 km. In addition, the ACPW algorithm can obtain CNOPs with more similar patterns in TC Fitow than in TC Matmo.
The vertically integrated energies of the CNOPs for TC Fitow are displayed in Fig. 7. Compared to the ADJ method, at the 120 km resolution, the CNOPs of the ACPW method have much lower energy and differing positions. However, when using a resolution of 60 km, similar energies and positions are obtained. Moreover, the energy of the CNOPs obtained from the ACPW algorithm has a larger range in the center.
Vertically integrated energies of the CNOPs for TC Matmo are displayed in Fig. 8. Compared with the ADJ method, at the 120 km resolution, the CNOPs of the ACPW algorithm have a lower energy and cover larger areas. However, when using a resolution of 60 km, although the energy is still lower, the positions are more similar.
When we evaluate the CNOPs, in addition to the characteristics and distributions of the CNOP patterns, consideration should also be given to the numerical similarities and the benefits of the CNOPs. Therefore, we calculate the similarity between the CNOPs determined from the ADJ method and the ACPW algorithm and use X and Y to represent them in the following formula, in which
The results are shown in Table 4. The similarity values can reflect the similarities among the CNOP patterns (Figs. 5 and 6).
In Table 4, for TC Fitow, the similarity at 120 km is −0.83, whereas the similarity with a resolution of 60 km is 0.43. For the analysis of TC Matmo, the similarity at 120 km is 0.42, whereas that with a resolution of 60 km is 0.37. The negative sign indicates that portions of the CNOPs from these two methods have opposite wind–vector directions, which is shown in Fig. 5. We also find that when using a higher resolution, the similarity is lower. The reason for this finding is that although the major patterns of the CNOPs are similar, the secondary parts differ and they cover larger areas. When using a higher resolution, we can achieve information on a smaller scale, and the identification of sensitive regions becomes more accurate. Regarding the analysis of the CNOP patterns, we obtain more similar major patterns when for a resolution of 60 km. However, compared with the different parts, the similar parts are very small. The similarities decreased do not affect the identification of the sensitive regions because the adaptive observations only focus on the points with larger influences, which will be demonstrated Sect. 4.4.
We also compare the energy for 24 h of nonlinear development under the initial states perturbed by different CNOPs, i.e., $J\left(M\right({\mathit{U}}_{\mathrm{0}}+\mathit{\delta}{u}_{\mathrm{0}}^{\ast}\left)\right)$. The results are shown in Table 5. All CNOPs obtained using the ACPW produce lower energies than the those of the ADJ method. However, when reducing the CNOPs to W×CNOPs in the entire domain and reducing the CNOPs by a factor of 0.5 in the sensitive regions, the ACPW algorithm has better results, which will be discussed in following subsection.
In this subsection, we design two groups of idealized experiments to investigate the validity of the sensitive regions identified using CNOPs based on two assumptions.
First, when adding adaptive observations in sensitive regions, the surrounding environment is idealized, and the improvements from adding observations reduce the original errors by a factor of 0.5.
Second, the obtained CNOPs can be seen as the optimal initial perturbations. Once we reduce them in the sensitive regions, the benefits are the highest.
Under these assumptions, by reducing the CNOPs to W×CNOPs and inserting them into the initial states, we can investigate how the reductions in the CNOPs influence the skill of TC forecasting. Moreover, reducing the CNOPs by a factor of 0.5 in the identified sensitive regions by vertically integrating the energies can be used investigate how the addition of adaptive observations in the sensitive regions can impact the skill of TC forecasting.
First, because CNOPs can be seen as the optimal initial perturbations in the TCAOs, we reduce the CNOPs to W×CNOPs, where W is a coefficient in (0, 1), insert the reduced CNOPs into the initial state and allow for 24 h of evolution of the MM5 model. Then, we calculate the forecast error using Formula (14) to determine the benefits of the reductions. Second, we determine the sensitive regions via vertically integrated energies using two schemes, namely the same points in the different resolutions and the equivalent percentage of points from the different grids. Then, we reduce the CNOPs by a factor of 0.5 in only the sensitive regions and insert the amended CNOPs into the initial states. The model is run for 24 h. The experimental results are described below.
We explore the forecast improvements induced by reducing the CNOPs to W × CNOPs for the entire domain. The approach requires using the reduced CNOPs in their initial state for a 24 h simulation of the MM5 model. The prediction error is computed by Formula (12), where
and the definitions of u_{NT}, P, M and U_{0} are the same as in Eqs. (1), (2), and (3).
The prediction error after reducing the CNOPs for the entire domain is computed by Formula (13), where
and W is the weighting coefficient, which is set to 0.25, 0.5 or 0.75 for decreasing error. The benefit from such reductions is calculated by Formula (14), represented as
The prediction benefit increases for decreasing W. Figures 9 and 10 also show that the ACPW algorithm can obtain CNOPs with better benefits from reducing the CNOPs to W×CNOPs for the entire domain than the ADJ method except for when W is 0.25 for TC Fitow at a resolution of 60 km. This is because the ACPW algorithm optimizes a low-dimensional feature space due to the PCA and focuses on more effective points in the entire domain, which has positive effects on improving the forecast.
We explore the forecast improvement caused by reducing the CNOPs by a factor of 0.5 in the sensitive regions. We determine the sensitive regions based on vertically integrated energies using two schemes, the 20 points with the highest energy at the different resolutions and 1∕100 points of the different grids, which is 30 points at the 60 km resolution (55 × 55) and 6 points at the 120 km resolution (21 × 26). The sensitive regions with the 20 points having the highest energy are denoted in Figs. 11 and 12.
In Figs. 11 and 12, when the equivalent points approach is adopted, a larger scope is covered with the 120 km resolution than with the 60 km resolution. When using the 20 points from the ADJ method and the ACPW algorithm and reducing the CNOPs by a factor of 0.5, the benefits are displayed in Table 6.
In Table 6, for TC Fitow, compared to the ADJ method, i.e., 5.93 % at the 120 km resolution and 3 % at the 60 km resolution, the ACPW algorithm obtains a higher benefit (8.05 %) for a resolution of 120 km and a lower benefit (−0.84) for a resolution of 60 km. Here, −0.84 % means that a reduction in the CNOPs results in no benefit and narrows the quality of the initial state. For the analysis of TC Matmo, the ACPW algorithm achieves a much higher benefit (20.48 %) than the ADJ method (6.12 %) at the 60 km resolution and a lower benefit (16.26 %) than the ADJ method (20.90 %) at the 120 km resolution. In addition, when using the same number of energy points, the benefits from using the 120 km resolution are nearly as high as those for the 60 km resolution except for the ACPW algorithm at 60 km resolution for TC Matmo.
The sensitive regions with 1∕100 points from the different grids are denoted in Figs. 13 and 14.
Figures 13 and 14 show that when using different resolutions, the sensitive regions identified by the same method are different. The sensitive regions identified by the ACPW algorithm are more dispersive than those identified by the ADJ method, which is attributed to the randomness of the intelligent algorithms. Table 7 shows the benefits gained from reducing the CNOPs by a factor of 0.5 in the sensitive regions identified by the ADJ method and the ACPW algorithm with different points in the different resolutions.
According to Table 7, for TC Fitow, the ACPW algorithm achieves a 4.23 % benefit, which is higher than the ADJ method (3.9 %) at the 60 km resolution and a lower benefit 0.01 % than the ADJ method (1.72 %) at the 120 km resolution. For the analysis of TC Matmo, the ACPW algorithm also has a higher benefit (9.75 %) and a lower benefit (6.86 %) than the ADJ method (1.21 % and 13.24 %, respectively).
Combined with Tables 6 and 7, we can conclude that the sensitive regions cover a larger scope and higher benefits are obtained. When using the same proportion of grids with the different resolutions, the sensitive regions under higher resolution achieve higher benefits. These results also demonstrate that the CNOPs obtained from the ACPW algorithm can identify sensitive regions with higher benefits at the 60 km resolution.
We further investigate the validity of the sensitive regions identified by the CNOPs through using a comparison of simulated TC tracks predicted by the MM5 model for each case by inserting the CNOPs or W × CNOPs into the initial states. We also simulate the TC tracks through the insertion of the amended CNOPs in the different sensitive regions (20 or 30 points). Because 120 km is the lowest resolution in this research and the tracks cannot be drawn under this resolution in our study, we only analyze the simulated TC tracks at the 60 km resolution. We draw two tracks in a sub-figure, which are represented by the observed TC tracks from the CMA–SHTI and the simulated TC track from the MM5 model, and the different perturbations are overlayed onto the same initial states. According to the experimental results, when overlaying the CNOPs or amended CNOPs onto the same initial states, although the CNOPs are obtained from different methods, the simulated tracks are the same. Therefore, we only discuss one group of figures for each case. The results are presented in Figs. 15 and 16.
Figure 15 demonstrates the simulated TC tracks of the MM5 by inserting the CNOPs or W × CNOP into the initial state for TC Fitow; the four sub-figures are the same. The reason is that the deviations of the simulated TC track and the observed TC track are very small. Therefore, it is not easy to make improvements. Hence, when inserting different CNOPs into identical initial states to simulate TC tracks, a change is not evident. Moreover, the resolution we used was 60 km, which is not high enough to show more details about changing tracks.
Figure 16 demonstrates the simulated TC tracks from the MM5 model by inserting the CNOPs or W × CNOP into the initial state for TC Matmo. Figure 16a and b are the same, and from Fig. 16b to d, the simulated positions after 24 h become closer to the observed positions. These results illustrate that when the CNOPs obtained by the ACPW algorithm and ADJ method are used as the optimal initial perturbations, reducing the CNOPs has a positive effect on the skill of the forecasting of the simulated tracks. Moreover, the ACPW algorithm is a meaningful and effective method for solving the approximate CNOPs of the ADJ method.
We also simulate TC tracks by inserting the amended CNOPs, which are reduced by a factor of 0.5 in only the sensitive regions. We use 20 and 30 points as the sensitive regions to study how the number of points affects the skill of forecasting. The results are shown in Figs. 17 and 18.
In Figs. 17 and 18, the simulated TC tracks are the same not only for different methods but also for different sensitive regions. We can conclude that the ACPW algorithm, an adjoint-free method, is a meaningful and effective method for solving the approximate CNOPs of the ADJ method. According to these results, we can also conclude that using 20 or 30 points as the sensitive regions results in the same improvement in the TC tracks in terms of forecasting. Thus, fewer points can be used in real adaptive observations to reduce costs.
To promote the efficiency of the ACPW algorithm, we parallelize it with MPI technology. The time consumption of each case is nearly the same. Hence, we can use a group of experimental results to elucidate the efficiency of the ACPW algorithm. Because the ADJ method cannot be parallelized because each input depends on the output of the previous step, its time consumption is not changed. Moreover, because this method generally uses 4 ∼ 8 initial guess fields to obtain the optimal value, we use 1 and 4 initial guess fields to determine the CNOPs. The time consumption of the ADJ method and the ACPW algorithm are shown in Table 8.
^{*} ADJ method (1) means using 1 initial guess field and ADJ method (4) means using 4 initial guess fields.
At the 120 km resolution, the time consumption of the ADJ method using 1 and 4 initial guess fields is 12.4 and 49.7 min, respectively. At the 60 km resolution, the time consumption is 79.9 and 321.1 min, respectively. Unlike the ADJ method, the ACPW algorithm can be parallelized. When using 22 cores, the ACPW method requires much less time, i.e., 2.74 min at the 120 km resolution and 20.8 min at the 60 km resolution. Obviously, the ACPW is more efficient. Compared to the ADJ method (1), the speedup reaches 4.53 and 3.84 for the different resolutions. Compared to the ADJ method (4), the speedup reaches 18.14 and 15.44. Although the different initial guess fields are calculated in parallel, the time consumption must be higher than that of the ADJ method (1), since the ACPW algorithm is also faster than the ADJ method.
In this study, we present a novel approach, the adaptive cooperative coevolution of the parallelized PSO and WSA (ACPW), to solve CNOPs. The CNOPs based on the ACPW algorithm are applied to study the identification of sensitive regions by TCAOs in the MM5 model without using an adjoint model. We study two TC cases, Fitow (2013) and Matmo (2014), with 60 and 120 km resolutions. The objective function is set as the total dry energy based on 24 h simulations starting with initial perturbations at the prediction time within the verification area. We also calculate CNOPs with the ADJ method, using results as a benchmark. To validate the ACPW algorithm, the CNOPs obtained from the different methods are compared in terms of the patterns, energies, similarities, benefits in reducing the CNOPs and simulated TC tracks with perturbations. To evaluate the advantages of the ACPW algorithm, we run the PSO, WSA and ACPW programs 10 times and compare the maximum, minimum and mean objective values as well as the RMSE. We also exhibit the objective value scope after the first iteration to analyze the effect of initial objective values on the different algorithms. To illustrate the performance of the algorithms, we compare the degree of change of the objective function value for the three algorithms. The analysis results demonstrate that the hybrid strategy and cooperative coevolution are useful and effective.
According to all of the experiments, the following five conclusions are obtained:
Compared with the ADJ method, the ACPW algorithm can obtain CNOPs with more similar patterns of temperature and wind for TC Fitow than those for TC Matmo.
At the 120 km resolution, the similarities in the CNOPs achieved by the ADJ method and the ACPW algorithm are higher than those at the 60 km. The reason is that although the major patterns of the CNOPs are similar, the other parts differ and cover larger areas. At a higher resolution, we can find information on a smaller scale. Moreover, sensitive region identification becomes more accurate. Regarding the CNOP patterns, more similar major patterns are obtained at the 60 km resolution, although the similar parts are very small compared with the other differing parts. However, the decreased similarities do not affect identifying sensitive regions because the adaptive observations only focus on the points with a larger influence.
When adding adaptive observations in the sensitive regions for a surrounding environment that is idealized, the original errors are reduced by a factor of 0.5. Thus, the CNOPs can be seen as the optimal initial perturbations. Once they are reduced in the sensitive regions, the benefits are highest. We design two groups of idealized experiments to investigate the validity of the sensitive regions identified by the CNOPs for the skill of TC-track forecasting. This involves reducing CNOPs to W × CNOPs and reducing the CNOPs by a factor of 0.5 in the sensitive regions identified using vertically integrated energies. The experimental results show that the CNOPs of the ACPW algorithm produce lower energies than the ADJ method but can obtain better benefits when reducing the CNOPs.
The ACPW algorithm can be effective for identifying the sensitive regions, which have the same influence on the forecast improvements of the simulated TC tracks as the ADJ method. We compare the different forecast improvements of the TC tracks with the different reduced perturbations, including reducing the CNOPs to W × CNOPs for the entire domain and reducing the CNOPs by a factor of 0.5 in the sensitive regions. The experimental results all support our conclusions.
The ACPW algorithm is more efficient than the ADJ method. Compared to the ADJ method using 1 initial guess field, the speedup reaches 4.53 at the 120 km resolution and 3.84 at the 60 km resolution. Compared to the ADJ method using 4 initial guess fields, the speedup reaches 18.14 and 15.44 for the 120 km and 60 km resolutions, respectively.
All of the conclusions demonstrate that the ACPW algorithm is a meaningful and effective method for solving approximate CNOPs and identifying the sensitive regions of TCAOs. In addition, as we reduce the dimensions within the PCA, the CNOPs lose some energy. Compared to the CNOPs form the ADJ method, the CNOPs from the ACPW algorithm are all local CNOPs. However, for the ACPW algorithm, they are global CNOPs. Because the PCA makes our optimization focus on more effective points with higher energies, the ACPW algorithm can achieve the CNOPs with better benefits and the same improvements in the skill of TC-track forecasting.
We are restricted to computation sources for the time being. We are also limited by the parallelization of the ACPW algorithm. We will improve the conditions of computation and use the parallel ACPW algorithm to solve CNOPs in the Weather Research and Forecasting (WRF) model with a finer grid and higher resolution. In addition, we will apply this type of method to solve CNOPs in the Community Earth System Model (CESM) model, which does not have an adjoint model.
The data in our paper are all obtained by ourselves and all data uploaded can be accessed.
The supplement related to this article is available online at: https://doi.org/10.5194/npg-25-693-2018-supplement.
LZ proposed the ACPW algorithm and designed the experimental scheme. BM and SY helped apply the work to the studying of the sensitive areas' identification of tropical cyclone adaptive observations. FZ provided the theoretical guidance about the atmosphere science. All four authors contributed to the writing of the paper.
The authors declare that they have no conflict of interest.
In this paper, research was sponsored by the Foundation of National
Natural Science Fund of China (No. 41405097). Edited
by: Christian Franzke
Reviewed by: two anonymous referees
Aberson, S. D.: Targeted Observations to Improve Operational Tropical Cyclone Track Forecast Guidance, Mon. Weather Rev., 131, 1631–1628, 2003.
Bender, M. A., Ross, R. J., and Tuleya, R. E., and Kurihara, Y.: Improvements in tropical cyclone track and intensity forecasts using the GFDL initialization system, Mon. Weather Rev., 121, 2046–2061, 1993.
Bergot T.: Adaptive observations during FASTEX: A systematic survey of upstream flights, Q. J. Roy. Meteorol. Soc., 125, 3271–3298, 1999.
Dee, D. P., Uppala, S. M., Simmons, A. J., and 33 co-authors: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597, 2011.
Franklin, J. L. and Demaria, M.: The Impact of Omega Dropwindsonde Observations on Barotropic Hurricane Track Forecasts, Mon. Weather Rev., 120, 381–391, 1992.
Froude, L. S. R., Bengtsson, L., and Hodges, K. I.: The Predictability of Extratropical Storm Tracks and the Sensitivity of Their Prediction to the Observing System, Mon. Weather Rev., 135, 315–333, 2007.
Jolliffe, I. T.: Principal Component Analysis, Springer Berlin, 87, 41–64, 1986.
Kennedy, J. and Eberhart, R.: Particle swarm optimization, in: Proc. of IEEE Int. Conf. Neural Networks, 1942–1948, 1995.
Mu, B., Wen, S., Yuan, S., and Li, H.: PPSO: PCA based particle swarm optimization for solving conditional nonlinear optimal perturbation, Comput. Geosci., 83, 65–71, 2015a.
Mu, B., Zhang, L., Yuan, S., and Li, H.: PCAGA: principal component analysis based genetic algorithm for solving conditional nonlinear optimal perturbation, in: 2015 International Joint Conference on Neural Networks (IJCNN), IEEE, Ireland, 12–17 July 2015, 1–8, 2015b.
Mu, M. and Duan, W. S.: A new approach to studying ENSO predictability: Conditional nonlinear optimal perturbation, Chinese Sci. Bull., 48, 1045–1047, 2003.
Mu, M. and Zhou, F. F.: The Research Progress of the Typhoon Targeted Observations Based on CNOP Method, Adv. Meteor. Sci. Technol., 3, 6–17, 2015.
Mu, M., Zhou, F., and Wang, H.: A Method for Identifying the Sensitive Areas in Targeted Observations for Tropical Cyclone Prediction: Conditional Nonlinear Optimal Perturbation, Mon. Weather Rev., 137, 1623–1639, 2009.
Qin, X. and Mu, M.: Can Adaptive Observations Improve Tropical Cyclone Intensity Forecasts?, Adv. Atmos. Sci., 31, 252–262, 2014.
Qin, X., Duan, W., and Mu, M.: Conditions under which CNOP sensitivity is valid for tropical cyclone adaptive observations, Q. J. Roy. Meteorol. Soc., 139, 1544–1554, 2013.
Qin, X. H.: A Comparison Study of the Contributions of Additional Observations in the Sensitive regionss Identified by CNOP and FSV to Reducing Forecast Error Variance for the Typhoon Morakot, Atmos. Ocean. Sci. Lett., 3, 258–262, 2010.
Qin, X. H. and Mu, M.: Influence of conditional nonlinear optimal perturbations sensitivity on typhoon track forecasts, Q. J. Roy. Meteorol. Soc., 662, 185–97, 2012.
Tang, R., Fong, S., Yang, X. S., and Deb, S.: Wolf search algorithm with ephemeral memory, in: Seventh International Conference on Digital Information Management, University of Macau, Macau, June 2012, 165–172, 2012.
Wang, X. L., Zhou, F. F., and Zhu, K. Y.: The application of conditional nonlinear optimal perturbation to the interaction between two binary typhoons fengshen and fung-wong, J. Trop. Meteorol., 29, 235–244, 2013.
Wang, X., Zhum, K., and Zhou, F.: The Study of Conditional Nonlinear Optimal Perturbation's Application in Typhoon Over the South China Sea, Tournal of Chengdu University of Information Technology, 6, 640–646, 2010.
Wen, S., Yuan, S., Mu, B., Li, H., and Chen, L.: SAEP: Simulated Annealing Based Ensemble Projecting Method for Solving Conditional Nonlinear Optimal Perturbation, in: Algorithms and Architectures for Parallel Processing, 14th international conference, ICA3PP 2014, Dalian, China, 24–27 August 2014, 655–668, 2014.
Wen, S., Yuan, S., Mu, B., Li, H., and Ren, J.: PCGD: Principal components-based great deluge method for solving CNOP, in: Evolutionary Computation (CEC), 2015 IEEE Congress on. IEEE, 25–28 May 2015, 1513–1520, 2015a.
Wen, S., Yuan, S., Mu, B., and Li, H.: Robust PCA-Based Genetic Algorithm for Solving CNOP, in: Intelligent Computing Theories and Methodologies, 11th International Conference, ICIC 2015, Fuzhou, China, 20–23 August 2015, 597-606, 2015b.
Ying, M., Zhang, W., Yu, H., Lu, X., Feng, J., Fan, Y., Zhu, Y., and Chen, D.: An overview of the China Meteorological Administration tropical cyclone database, J. Atmos. Ocean. Technol., 31, 287–301, 2014.
Yuan, S., Qian, Y., and Mu, B.: Paralleled Continuous Tabu Search Algorithm with Sine Maps and Staged Strategy for Solving CNOP, in: Algorithms and Architectures for Parallel Processing, 15th International Conference, ICA3PP 2015, Zhangjiajie, China, 18–20 November 2015, 281–294, 2015.
Zebiak, S. E. and Cane, M. A.: A Model El Niño-Southern Oscillation, Mon. Weather Rev., 10, 2262–2278, 1987.
Zhang, L. L., Yuan, S. J., Mu, B., and Zhou, F. F.: CNOP-based sensitive areas identification for tropical cyclone adaptive observations with PCAGA method, Asia-Pac. J. Atmos. Sci., 53, 63–73, 2017.
Zhou, F.: Application of conditional nonlinear optimal perturbation method to typhoon target observation, Institute of Atmospheric Physics, Chinese Academy of Sciences, 43–55, 2009.
Zhou, F. and Mu, M.: The impact of verification area design on tropical cyclone targeted observations based on the CNOP method, Adv. Atmos. Sci., 28, 997–1010, 2011.
Zhou, F. and Mu, M.: The Impact of Horizontal Resolution on the CNOP and on Its Identified Sensitive Areas for Tropical Cyclone Predictability, Adv. Atmos. Sci., 29, 36–46, 2012a.
Zhou, F. and Mu, M.: The Time and Regime Dependencies of Sensitive Areas for Tropical Cyclone Prediction Using the CNOP Method, Adv. Atmos. Sci., 29, 705–716, 2012b.
Zhou, F. and Zhang, H.: Study of the Schemes Based on CNOP Method to Identify Sensitive Areas for Typhoon Targeted Observations, Chinese J. Atmos. Sci., 2, 261–272, 2014.
Zhou, F., Qin, X., Chen, B., and Mu, M.: The Advances in Targeted Observations for Tropical Cyclone Prediction Based on Conditional Nonlinear Optimal Perturbation (CNOP) Method, in: Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications (Vol. II), Springer Berlin Heidelberg, 577–607, 2013.
Zhu, H. and Thorpe, A.: Predictability of Extratropical Cyclones: The Influence of Initial Condition and Model Uncertainties, J. Atmos. Sci., 63, 1483–1497, 2006.
Zou, X., Vandenberghe, F., Pondeca, M., and Kuo, Y.-H.: Introduction to adjoint techniques and theMM5 adjoint modeling system, NCAR Tech. Note, NCAR/TN-435-STR, 107, 49–66, 1997.
An interactive open-access journal of the European Geosciences Union