J. Meteor. Res.  2017, Vol. 31 Issue (6): 1007-1017   PDF    
The Chinese Meteorological Society

Article Information

XU, Huanbin, and Jinfang YIN, 2017.
Key Issues in Developing Numerical Models for Artificial Weather Modification. 2017.
J. Meteor. Res., 31(6): 1007-1017

Article History

Received July 7, 2017
in final form September 4, 2017
Key Issues in Developing Numerical Models for Artificial Weather Modification
Huanbin XU, Jinfang YIN     
State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081
ABSTRACT: The scientific foundation of artificial weather modification is meso- and small-scale dynamics and cloud–precipitation microphysics. Artificial weather modification requires the realistic coupling of weather patterns, dynamical processes, and microphysical processes. Now that numerical models with weather dynamical characteristics have been widely applied to artificial weather modification, several key points that should not be neglected when developing numerical models for artificial weather modification are proposed in this paper, including the dynamical equations, model resolution, cloud–precipitation microphysical processes, numerical computation method, and initial and boundary conditions. Based on several examples, approaches are offered to deal with the problems that exist in these areas.
Key words: artificial weather modification     numerical model     dynamical processes     cloud–precipitation microphy-sical processes    
1 Introduction

Compared with other countries worldwide, China has carried out the largest amount of artificial weather modification (AWM) activities. Various approaches have been employed, and the practice of AWM has become a part of China’s operational weather service. The key factor for the rise in weather modification activities is societal demand. It is also a representation of the progress that has been made in AWM science and technology, as well as operational application. However, AWM is still unable to fully meet the urgent needs of society in terms of water resources, disaster (flood, drought, hail, etc.) mitigation, and environmental protection, because of the substantial limitation and lack of understanding of the basic scientific issues involved in AWM that remain at present. Therefore, it is necessary to thoroughly investigate the macro- and microscale processes of natural clouds, as well as adopt scientific approaches to modify the evolution of clouds. In fact, observational or theoretical studies can only provide clues to the structure and evolution of clouds. To reproduce these aspects in their entirety, a numerical model is required. This is because a numerical model with a full set of equations will enable the reproduction of various weather phenomena through the integration of observations into theoretical and experimental results. Accordingly, it is necessary to continuously develop the next generation of numerical models for AWM on the basis of the characteristics of AWM activities.

Zheng and Guo (2012) proposed that the development of detailed numerical models that can be used in AWM operational work is driven by developments in computing technology. They emphasized that numerical models can be used to conceptually design and study, as well as offer guidance for, AWM operations; plus, they can help in analyzing the results following AWM activities. Although considerable progress in AWM has been made by combining observations from aircraft, radar, satellites, and others with meso- and small-scale numerical model results, there are still many scientific issues and tech-nical problems to be addressed (Liu and Lou, 2002).

In fact, a great deal of effort has already been devoted to numerical model applications in AWM (e.g., Mao and Hu, 1993; Huang and Xu, 1994; Hong, 1998; Hong and Lei, 2012; Guo et al., 2013, 2015). The National Meteorological Center of China developed an AWM prediction system based on the GRAPES model in which a detailed cloud microphysical scheme and a seeding module were included (Chen et al., 2007; Ma et al., 2009). Based on the Beijing rapid update cycle of assimilation forecast system (BJ-RUC), a cold cloud seeding potential recognition platform of weather modification in Beijing (BJ-CCSPR platform) was designed by the Beijing Weather Modification Office (Ji et al., 2014). Moreover, the Weather Modification Center of the China Meteorologi-cal Administration developed two operational systems called MM5-CAMS (Chinese Academy of Meteorologi-cal Sciences) and GRAPES-CAMS, respectively. In weather modification work, using an AWM numerical model can obtain information on the timing, spatial range, amount of seed, and other factors related to AWM. Additionally, numerical modeling is also a useful approach for verifying the results of a weather modification operation.

Presently, numerical models used in weather modification are developed on the basis of small- and/or mesoscale numerical weather prediction (NWP) models, which include cloud microphysical schemes and cloud seeding modules, and these models show potential as tools in theoretical studies and operational applications in the field of AWM. However, it remains a key point as to whether these models and their component microphysi-cal schemes are suitable for AWM, and to what extent they can reproduce various types of weather phenomena. Indeed, thus far, certain problems have emerged from the application of such models and, in view of this, updates or improvements to the dynamical equations, numerical computation method, and cloud–precipitation microphy-sical processes, as well as the techniques for coupling the dynamical and physical processes, are required.

Much attention has been paid in China over the past two decades to developing AWM numerical models. For example, several three-dimensional (3D) complete elas-tic nonhydrostatic models have been built, e.g., by Kong et al. (1990), Wang and Hu (1990), Xu and Wang (1990), and Zhou et al. (1996). Besides, Guo et al. (2006) coupled a cloud microphysical scheme with discrete mass categories of hail/graupel into a hail storm model.

Alongside these developments, achievements have also been made in mesoscale weather dynamics and large-scale computational mathematics, and these successes have been employed in developing AWM numerical models. In doing so, however, some shortcomings have also been introduced, and these shortcomings in certain functions have appeared gradually. If these issues are not comprehensively investigated, they could potentially cause problems. Therefore, it is necessary for AWM numerical model developers to understand the characteristics of AWM, thus enabling them to design suitable models for this field of study and its associated operational applications. One specific point is that, when attempting to build a reliable AWM numerical model, it would be meaningful to include a seeding module and verify its results.

In this study, we investigate some of the basic issues in the development and application of AWM numerical models. Specific questions related to the building of the next generation of AWM numerical models are discussed, and the answers to which might provide useful guidance for designing improved AWM numerical models in the future.

2 Characteristics of AWM numerical models

The characteristics of an AWM numerical model are determined by the functions of the model. Users expect to gain the original state, evolution processes, and final state of cloud systems from AWM models. The AWM activities can then be applied based on these conditions, to achieve an evolutionary direction of the cloud system, to reduce the probability of natural disasters, and/or to gain some other public benefits. Cloud structure and its evolution are the result of the interaction between cloud dynamics and cloud–precipitation physics. Therefore, AWM models should have the ability to describe macro- and microphysical processes as closely as possible to capture how they occur in nature, as well as the feedbacks between them. It is important to emphasize that the cloud structure and evolution mentioned here are directly described in detail in the model, rather than the overall effects of the cloud–precipitation processes. To capture the overall effects of clouds in climate models, the most important aspect—in addition to their vertical distribution—is the effects of clouds on radiation and condensation phase-change related latent heat release; while the specific properties, such as liquid water content, cloud top height, and others, of clouds are not necessarily required. However, most weather models focus only on the precipitation amount; the detailed cloud attributes are not considered. In other words, the most commonly used weather models usually enable the capturing of rainfall intensity and spatial distribution, but pay little attention to cloud structure and evolution.

In summary, AWM numerical models have their own characteristics, and have significant differences to other (weather, climate, environment, single-column cloud–precipitation microphysical models) numerical models. Other models are usually only concerned with the cloud field, ground precipitation, cloud microphysics, latent heat profiles, or similar. However, AWM numerical models need to cover the structure of the cloud fields, cloud–precipitation particle properties, and the cloud structure and evolution of cloud particle fields. Therefore, AWM numerical models need to explicitly describe cloud processes. However, cloud structure and evolution are the result of the interaction between cloud macroscopic fields and cloud microphysical properties, which cannot be described by simple or fixed parameterization schemes. Therefore, more stringent functions are required for the dynamical–thermal framework design and the explicit description of cloud–precipitation microphysics in AWM numerical models.

3 General requirements of an AWM numerical model

A set of accurate and reliable dynamical equations is required to describe fine-scale cloud structures. Therefore, complete three-dimension (3D) dynamical equations, conservative numerical computation methods, physical schemes applicable on high temporal and special resolutions, and suitable initial and boundary conditions, are required. Moreover, cloud microphysical schemes that can determine the evolution of cloud processes at a bifurcation point are particularly necessary. These are the key points when designing an AWM numerical model.

Additionally, the phase of hydrometeors and their spatial distributions are also key parameters. Consequently, it is necessary for an AWM numerical model to include cloud microphysical schemes that enable the capture of main evolutionary properties of cloud microphysical processes and provide the best quantitative conditions for seeding and tracing of cloud evolution. Moreover, models should be able to provide a large number of diagnos-tic variables to provide guidance for seeding verification and assessment.

4 Specific aspects of the general requirements 4.1 Complete 3D dynamical equations

Given that atmospheric movement is 3D, a numerical model should contain a 3D description of such motion. Owing to limitations in computing capability, one-dimensional (z) and two-dimensional (x, z) models are often used in the early stages of numerical model development. However, the reduced dimensionality causes a series of distortions that make it difficult to correctly consider the wind field, and thus misrepresents the convergence of horizontal airflow. This is particularly true for one-dimensional models. The two-dimensional models can only describe horizontal movement in the positive and negative directions, and cannot describe the horizontal rotation of flow, thus limiting the energy exchange between the ambient wind and the convective circulation. Therefore, to reasonably describe various phenomena of the atmosphere, a numerical model should take into account the 3D structure of atmospheric movement. Since atmospheric phenomena are usually 3D, as mentioned above, it is useful to consider as much of the atmospheric compressibility as possible in cloud-scale convective events, and it is necessary to develop a 3D cloud-scale model for compressible atmosphere, in which the equations should be simplified as little as possible (or little transformation with assumptions, except for improving some functions), and should contain three-phase hydrometeors.

The atmospheric equations can be expressed by using Eulerian or Lagrangian coordinate. But which one performs better? The representation of natural clouds should involve a group of hydrometeor particles growing as they move under certain dynamic and thermodynamic conditions. The movement of these particles is determined by both the dynamic field and the particle movement characteristics such as terminal velocity, mass, shape, and surface roughness; and their growth is determined by water vapor, hydrometeor interactions, and ambient thermal conditions. Although particles occur in groups, the movement of each particle is controlled by the physical properties of the particle itself. The movement of a large particle is mainly affected by its own attributes and the environmental field, but small particles are also affected by the movement of large particles. However, the distance of interaction between large and small particles is difficult to obtain through computational physics and experiments. Generally, it is assumed that the movement of a single particle is not constrained by surrounding particles of the same type. Although the hydrometeor field is a discontinuous field, it is not a dynamic rupture field. It does not block the interactions of cloud macro fields with microphysical ones. To achieve numerical model results that are as close as possible to the physical picture of the real atmospheric movement, the following options are available: the macroscopic field of cloud can be described by the Eulerian method; the background field of cloud, which can be treated as a continuous field, can be described by the semi-Lagrangian method; and the precipitation field can be dealt with via the total Lagrangian method (Xu and Duan, 2001; Xu, 2012).

4.2 High spatial resolution

Since AWM models require a detailed understanding of the structure and evolution of cloud systems, one direct approach is to improve the temporal and spatial resolution of models. First, enabling clouds to be resolved is necessary, but to achieve this, what resolution is needed? Importantly, the resolution of radar observations should be matched, which is currently about 200 m. Therefore, model resolution should be comparable to this level.

There are several advantages to increasing model resolution, including: (1) a reduction in calculation errors; (2) good representation of small-scale information; (3) accommodating and describing the excitation of short waves; and (4) expanding model simulation capabilities for multi-wave (components of a spatial wave group at a given time) and varying-wave (component variation of a wave group at different times) moments. However, meso- and small-scale movements are usually 3D, non-hydrostatic, and nonlinear, which is difficult to control during model integration. Therefore, improving model resolution is not only a single technical strategy; it also involves several relevant scientific and technological issues.

4.3 Adaptive explicit cloud microphysical scheme

During AWM operations, it is necessary to obtain cloud structure and its evolution in detail from model outputs. Therefore, the use of an explicit cloud microphysics scheme is required to directly describe the cloud–precipitation microphysical processes at play. The study byZhang et al. (1988) showed that a cloud microphysics (explicit) scheme has advantages for stratiform cloud precipitation, while a convection parameterization (implicit) scheme is beneficial for convective precipitation. A direct (explicit) description of cloud–precipitation processes is suitable for stratiform cloud with a simple dynamic structure, while an indirect (implicit) parameterization scheme is beneficial for convective cloud with complex dynamic structures under coarse-resolution conditions. Recent studies have shown that the description of small-scale movement is poor at high resolutions (grid distance of less than 2 km), implying that some flaws may exist when explicitly describing subgrid-scale movement over a fine grid. This is contrary to the original intention of developing explicit schemes, meaning that it is important to resolve why explicit schemes have advantages for large-scale and uniform cloud systems but weaker ability than implicit schemes when describing medium- and small-scale and non-uniform cloud systems.

One possible and important reason is that explicit schemes are used to describe the clouds directly. Although it is possible to describe the details of cloud processes with an explicit scheme, it is the dynamical processes that support the occurrence, development, and decay of the cloud, and therefore constitute the key factor. When a dynamical framework is appropriate, the general structure of a cloud would become realistic. Only under these conditions does a detailed description make sense. For a dynamical framework for large-scale and uniform clouds, the average airflow in the grid area is close to that of the grid itself, and thus the cloud structure supported by the average airflow is similar to that of real cloud. Therefore, when the cloud structure is generally suitable, describing the fine-scale characterization with an explicit scheme works well. As far medium- and small-scale non-uniform clouds, considerably large differences exist between the average airflow and the real one. As a result, the average airflow over the grids is significantly different to the actual airflow. The difference might be between convective clouds and stratiform clouds, or even between cloud and cloudless conditions. If the model cannot capture the basic characteristics of no-cloud or cloud area, the detailed explicit scheme is meaningless.

Regarding the resolution of models, currently, the subgrid-scale processes are still fairly complex and important, such as the subgrid cloud and its interaction with the relevant scale, the subgrid gravity-wave drag, turbulent diffusion, and so on. Therefore, an explicit cloud microphysical scheme needs to be placed in coordination with subgrid physical processes to describe these physical processes reasonably.

Next, we use an example to illustrate this issue. As shown in Fig. 1a, at the beginning, the average ascending velocity in the grid is zero. This means that there is no dynamical process for cloud formation, and the condition should be no-cloud. However, the average being equal to zero does not mean that there is no ascending movement in the grid. In fact, it can be an average result over a distribution with multiple vertical movements. Since there is indeed ascending movement, there should be cloud formed (Fig. 1). In the same way, if the grid has some uniform ascending movement, stratiform cloud will form. Depending on the intensity of non-uniform ascending movement, mixed cloud or convective cloud should be generated (Fig. 2). Therefore, correctly describing the cloud field depends primarily on a suitably represented airflow field. That is, the simulated airflow field should has the same pattern as that of the actual airflow field.

Figure 1 Potential patterns of vertical motion and cloud distribution under the condition that the gridbox-averaged vertical velocity is 0 m s–1: (a) no cloud; (b) one cloud cell; (c) two cloud cells; and (d) three cloud cells.
Figure 2 Potential patterns of vertical motion and cloud types under the condition that the gridbox-averaged vertical velocity is greater than 0 m s–1: (a) stratiform cloud; (b) stratiform and cumulus mixed cloud; and (c) cumulus (convective) cloud.

Subgrid movement is contained in implicit schemes, which do not simply estimate the formation of cloud only by ascending movement. In fact, a variety of convective parameterization schemes have taken into account the effects of subgrid movement on clouds and precipitation, as well as the influence of subgrid vertical transfer on temperature and humidity structure. In view of this point, although the cumulus convective parameterization is described in an implicit, indirect, and approximate manner, its representation of physical processes is relatively complete. On the other hand, despite explicit cloud microphysical schemes describing cloud processes in detail, they still need to be based on a reasonable and adaptive dynamic framework. Otherwise, fake phenomena might appear, with fine-scale cloud microphysics described but important physical mechanisms missed. Therefore, mixed-description schemes have been introduced—that is, both explicit and implicit descriptions are applied to a scheme (Tiedtke, 1993; Tan et al., 2013). However, this approach should heed the possibility of “repetition.” That is, a cloud process might be considered in both the explicit and implicit schemes.

Experimental studies using hybrid schemes have also shown that the importance of implicit cumulus convection parameterization is significantly reduced when the horizontal resolution of a model is below 5 km (Zhang et al., 1988; Yu and Lee, 2010). Therefore, some researchers have suggested that cumulus convection parameterization schemes should be avoided as much as possible at horizontal resolutions below 5 km, and all cloud–precipitation microphysical processes should be described by the explicit scheme (Zhang, 1998).

From the above analysis, it can be concluded that when the difference between the average grid movement of a model and the actual movement is very small, the model provides a movement style that is close to that of the real atmosphere. This indicates that the dynamic framework of the model is correct and reliable. In such a case, a good explicit cloud–precipitation scheme may achieve decent results (both macro- and micro-fields are accurate). When the macroscopic dynamical field is inappropriate, even though the implicit convective parameterization schemes enable approximate representation of subgrid movement, they are still not capable to produce the specific subgrid movement. This requires increasing the model resolution as much as possible to modify the subgrid movement into a grid-distinguishable movement, and to reduce the effects of subgrid movement.

However, many studies have indicated that increasing the model resolution—especially up to around 1 km—does not definitively improve simulation results (Schwartz et al., 2014; Seiki et al., 2015). The reason for this is that there are problems in the existing dynamic framework design and numerical computation method, which cannot correctly describe the movement episodes with a spatial resolution of less than 2 km. This is possible because the equations are nonlinear, and accurate numerical solution is difficult to obtain. Therefore, a suitable numerical integration scheme, which has some conservation requirements such as mass conservation and energy conservation, is required. For weather and climate models, such a requirement may be satisfied with a resolution of around 2 km. For AWM models, however, there may be other requirements at a resolution of less than 2 km, such as whether the components of these conservative quantities are appropriate. Momentum conservation means that the momentum of all scale movements is constant. However, in the case when the momentum of large-scale movement is marginally large, and that of small-scale movement is marginally small, or vice versa, the total momentum might still be conserved. In such a case, the composition of the momentum component of the movement is different, or the characteristics of the movement or the dominant scale have been changed; that is, the nature of the movement has been changed, which might result in a distortion of movement. Consequently, “fidelity” is required in a numerical model.

In terms of model resolution, aside from considering the coordination between model dynamic processes and cloud microphysical processes, the effects of model resolution on cloud microphysical and parameterization schemes should also be considered. For instance, the atmosphere may be saturated with a relative humidity of 70% for a horizontal grid distance of 10 km, and the part greater than 70% can be condensed into clouds. However, saturation can only take place under relative humidity up to 100% in the case of the resolution being around 1 km. In addition, the effect of turbulence on cloud microphysical processes should also be considered at resolutions of around 1 km or higher.

4.4 Conservative numerical computational schemes

In the early stages of numerical model development, the first author of this paper has noted that choosing a differential scheme with short waves attenuated as small as possible, and paying attention to controlling the effects of the waves with wave length of doublet grid distance, play a slight positive role in reducing computational errors when the existing differential schemes are used. However, the idea of “fidelity” had not been proposed for numerical computation at that time.

In fact, a better simulation of wind field relies on a better dynamical framework of numerical models. Generally, it seems difficult to improve the dynamical framework, and only slight benefits can be obtained by improving the dynamical framework. This might be an excuse for shrinking back from difficulties given by the developers of large- and medium-scale models. As for the dynamic requirements of AWM models, regardless of the benefits being large or small in improving the dynamical framework, it is worthwhile improving the dynamical framework when facing such difficulties. At the very least, it can be seen that only a set of complete equations is not enough in terms of the dynamical framework, and a conservative numerical computational scheme for differential calculation of the equations is required too. Zhong (1992a) and Zhong et al. (2002) proposed a constructive principle for spatial discretization in which the idea and method for all-energy conservation are contained. This method deserves attention, and it might be worth trying to apply it to the improvement of the dynamical framework in AWM models.

A conservative numerical computation scheme for spatial discretization was employed in the fidelity model developed by Zhong (1992b). This new method, proposed in the 1990s but received scarce attention, is based on the characteristics of energy conversion in the actual atmosphere. In view of energy balance, the energy conversion and constraint mechanism among the internal energy, gravity energy, internal energy, and latent heating, as well as the energy distribution mechanism between these elements, are reconstructed, so as to maintain the energy and mass characteristics of the element. Compared to the existing dynamic framework of nonhydrostatic models, the newly proposed scheme can effectively suppress the “noise” due to false sources and sinks caused by the damage of energy balance after discretization. Although the “conservation” scheme keeps the shortwaves of atmosphere, the development of false noise may introduce computational problems leading to instability (Fig. 3). Figure 3 compares the analytic solution and the general numerical solution. It can be seen that there is considerable noise in the conventional numerical method. These false and small noises not only distort the physical picture, but also cause problems of calculation instability.

Figure 3 Comparison of (a) a conservative scheme with (b) a conventional scheme (Qing Zhong, personal communication).

It appears that local energy conversion conservation is required at the grid point (i, j, k), rather than only in the full computational domain (Fig. 4). This idea is quite reasonable. The local energy conversion and conservation (fidelity) in a grid means the absence of obvious distortion in each grid, and thus it is possible to retain conservation in all fields in the whole domain.

Figure 4 Flowchart of a conservative scheme for local energy transfer at the grid (i, j, k) (Qing Zhong, personal communication).

As for the fidelity model discretization design, it is proposed that the local energy conversion should be satisfied at the grid point (i, j, k) as follows:


where K is kinetic energy, Φ is the gravity potential, T is temperature, F is the external force, V is wind speed, Q is specific humidity, cυ is specific heat at constant volume, cp is specific heat at constant pressure, and ρ is air density.

Guaranteeing the stability of the calculation and maintaining the energy characteristics in a grid enable the accurate description of small-scale components and multi-scale structure. It also helps to avoid the weakness of component variability due to excessive numerical smoothing and additional dissipation caused by artificial treatments. The preliminary tests of Zhong (personal communication) show that the stability of discrete integration was improved because of the improvement in the local energy conversion and conservation in a grid, with an effective resolution of 3–4 grids. In this way, the new scheme makes the discrete integration stable through satisfying the necessary and sufficient condition. On the other hand, it is noted that the effective resolution of the WRF model is 5–7 grids, the MM5 model is 8–9 grids, and the ECMWF model is 10 grids. Therefore, this approach should be taken into account in further research and experiments by AWM numerical model developers.

4.5 Shortcomings from component distortion, excessive numerical smoothing, and additional dissipation

Atmospheric movement alternates between quiet and drastic change, and secondary small-scale movement is often accompanied by transitions between these states. However, small-scale movement is difficult to deal with in numerical models, and has been removed because it seems that small-scale movement makes no sense for weather or dynamical processes as viewed in the past. But is this really the case? In reality, it does not seem so. Next, we provide an example to explain this issue.

Figure 5 is an example of a closed streamline, with its possible distorted patterns. Figure 5a shows that when the air rises due to buoyancy, there is a “closed convection circulation,” which is a new component on the move-ment scale and is a small-scale movement. If the component is eliminated as noise in the calculation, the “embryonic convection circulation” initializing convections might be choked. Even though it is not eliminated, the movement of the closed convection might be replaced by open convective airflow or wavy patterns due to unconsciously handling of the small-scale movement or its distortion (Figs. 5b, c). From the structural characteristics of the flow patterns, it is seen that there are three emergent flow patterns that are essentially different in nature. For example, differences exist in their maintenance and development mechanism, as well as advantageous conditions. Besides, the range for unstable energy collection, the evolution of the final state, and the accompanying weather phenomena, are significantly different. If this were true in natural clouds, it would be a highly meaningful dynamics topic for numerical model development.

Figure 5 Patterns of (a) closed convective airflow, (b) open convective airflow, and (c) wavy upstream airflow.
4.6 Explicit microphysical scheme verification and development of intelligent schemes for determining the bifurcation of weather system evolution

In addition to improving the dynamic framework, existing explicit microphysical schemes also need to be verified and improved. As far as high-resolution model results are concerned, there are three possible reasons for simulation errors. First, the model dynamical processes fail to match the cloud microphysical processes. Second, the parameter values in a cloud microphysical scheme mismatch the model resolution. Third, even though the coupled cloud microphysical scheme is complete, there are several steps to link dynamical processes with cloud microphysical processes. Therefore, are the rules employed in a model the same as those in natural clouds? How much of an influence does inconsistency in the model have? It is necessary to develop natural cloud–precipitation processes by combining field observations, radar measurements, satellite retrievals, and other data, and then compare these with numerical simulations for investigating the differences between simulations and observations and propose new methods to match dynami-cal processes with cloud microphysical processes. Additionally, model results and associated diagnostic information may be used to predict the evolution of cloud microphysical processes.

One important question is how to evaluate existing cloud microphysical schemes. Comparing the modeled precipitation with ground observations of precipitation is the most simple and common method. However, this method only understands the terminal state of the cloud microphysical processes, and is unable to investigate the evolution of cloud and precipitation processes. The radar and satellite observations make up the defect in terms of temporal and spatial resolution of ground observations of precipitation to some extent. Therefore, AWM numeri-cal models require more in-depth explorations and applications of the observations from radar, satellite, conventional ground instruments, and aircraft data, than those of weather and climate models.

Medium- and small-scale weather systems are the product of macro- and micro-field interactions in the atmosphere. Consequently, there are feedbacks between dynamical processes and cloud microphysical processes and, in fact, the feedback from cloud microphysical processes to dynamical fields is stronger. In addition, such a feedback plays a major role in the direction of evolution of medium- and small-scale weather systems, and the feedback pathways vary in various ways. As for the description of cloud–precipitation processes, explicit schemes offer a variety of coherent pathways between the various links. It varies as to which link plays the major role and which downstream link it evolves into. Natural evolution in clouds is one of several possible paths, and there are significant differences in the final state among the various evolutionary pathways. Therefore, it is difficult to use a single-cloud microphysical scheme to describe all types of cloud–precipitation microphysical processes. This might be the main reason for the model uncertainties caused by cloud microphysical processes and their feedbacks.

More recently, Yin et al. (2015) pointed out that it is necessary to reveal the evolution of the dominant cloud microphysical process forming severe weather systems to improve numerical simulations of heavy rainfall occurrence and development. The development of “intelligent” cloud microphysical schemes that make transformation of cloud microphysical processes from “hard rules” to ones close to “natural transformation” is required so as to make the transformation of cloud microphysical processes follow their own natural tendency in numerical models. It is well-known that the evolution of cloud–precipitation microphysical processes is variable, with many possible pathways. However, defining the evolution of cloud–precipitation microphysical processes in numerical models is based on a few cases of observations, theoretical research, or experimental results, which greatly limits the development of cloud–precipitation in numerical models. Therefore, it is necessary to combine multi-source observations for various weather systems to analyze the evolution of cloud–precipitation processes under different macro- and microphysical conditions. Based on these results, a variety of options could then be designed in a cloud microphysical scheme. By doing so, cloud microphysical processes in a numerical model can be selected properly among the various approaches according to the cloud–precipitation microphysical properties and macroscopic conditions provided by the numerical model.

4.7 Classifying hydrometeor particles

Hydrometeor particles comprise cloud droplets, raindrops, ice crystals, snow, graupel, and hail. It is difficult to consider the effects of particle phase, shape, and density by classifying them only according to their mass or scale, if not sorting the hydrometeor particles into seve-ral types. To distinguish the phase, is it only necessary to divide the particles into two types (liquid and ice particles)? The shape of liquid particles can be assumed to be quasi-spherical, and the density is equivalent to water. The shape of ice particles is too diverse, making it difficult to scrutinize. However, the conditions for dry or wet growth can be used to determine whether ice particles are dry or wet (water), and thus to diagnose the density. Can solid particles be divided into ice crystals, snow, graupel or hail according to scale, density, and other characteristics? In fact, it is not necessary to divide hydrometeor particles into ice crystals, frozen drops, snow, graupel, and hail by the artificial hard rules given by the scheme designer, since such a classification increases the tasks that involve equations describing the interactions among the ice particles in a cloud microphysical scheme, and thus hinders the natural conversion between the ice particles.

4.8 Suitable initial and lateral boundary conditions

AWM numerical models require a better scheme for initial and lateral boundary conditions than those of weather and climate models. For the initial conditions, objective analysis should pay attention to the structural characteristics of medium- and small-scale systems, and data at short distances should be used to fill the grid field as far as possible to retain medium- and small-scale information. Figure 6 shows the fitting curves using data at long distances and short distances. It is indicated that the fitting with data at short distances highlights the characteristics of medium- and small-scale systems.

Figure 6 Schematic diagram of the natural distribution of data (solid line), and the fitting curves using data at long distances (dashed line) and at short distances (dot-dashed line), respectively.

Regional models involve issues related to boundary conditions. In particular, the requirements of a wide movement scale spectrum and high fidelity are needed for weather models, and thus false wave information should be prevented. Therefore, thoughtful consideration is required for the processing of boundary conditions. Our previous tests indicate that several commonly used boundary treatment methods are far from perfect, and thus the suspension boundary and the hierarchical boundary are designed for this purpose. Large-scale movement from boundary inflow may be blocked at the boundary. Therefore, the suspended boundary method was utilized (Xu and Wang, 1988). At the outflow boundary, parasi-tic short waves, wave reflection, or standing waves can often be found. Consequently, the hierarchical boundary scheme was used to “escort” false waves (Xu, 2014).

5 Conclusions and discussion

There are certain difficulties in building a numerical model with comprehensive functions. Special functions are required for AWM numerical models, so it is hard to employ the previous achievements of weather and climate models. Consequently, scientists involved in AWM should start some work for improving AWM numerical models at first. Based on different types of issues in AWM numerical models, special AWM numerical models should be designed and developed first—for example, building models separately for theoretical or mechanism studies, and detailed simulation (reproduction, integrated observations, conceptual model design, etc.) for real cases, as well as results verification of AWM. The following characteristics are associated with an AWM numerical model.

(1) The 3D and compressible atmospheric equations are required for AWM numerical models, as well as equations with three-phase hydrometeors included. The Eularian method describes the cloud field; the background field of cloud water particles, which is treated as a continuous field, is described by the semi-Lagrangian method; and the precipitation particles are described by the full Lagrangian method. Boundary inflow is a large-scale movement, which may be blocked at the boundary, and thus the suspended boundary can be used to solve the problem. As for the outflow boundary, at which parasitic short waves, wave reflection, or standing waves often occur, the hierarchical boundary scheme can be used to “escort” false waves.

(2) The horizontal resolution of AWM numerical models should be comparable to the resolution of radar observations (around 200 m), which enables the description of the inhomogeneity of the atmospheric movement. Under the high-resolution condition, developing a conservative differential method ensures the conservation of momentum, and thus generates a reasonable evolution for each component of the momentum.

(3) Explicit cloud microphysical schemes are necessary to directly describe the cloud–precipitation microphysical processes in AWM numerical models, and several types of hydrometeors can be dealt with by the bin method to meet certain needs. Dynamic processes and model resolution should be taken into account when an explicit cloud microphysics scheme is developed.

If an AWM numerical model is built with the above-mentioned functions, it would be an effective instrument for simulating or predicting the structure and evolution of natural cloud. The evolution of real cases can be reproduced by the newly built model and this helps in making decisions for what should be done with respect to AWM to achieve an evolution of the clouds toward our expectation. With the newly built model, it is possible to do some work for AWM with special purposes “according to the actual circumstance” and the “structure” of a real case from observation. Then, the strategies for AWM are made based on the actual “structure” of a weather phenomenon. To obtain reliable results, AWM should be launched “according to the actual circumstance” of clouds, rather than “according to rough plans.”

In fact, the functions of data assimilation and cloud analysis are also necessary in an AWM numerical model. Besides, complete modules for seeding processes and technologies are required. Indeed, improving the model is also helpful for data assimilation and cloud analysis. Based on those functions of a model, it is possible to make available macro- and micro-field cloud products, as well as modules with key parameters for cloud microphysics diagnosis and verification.

Acknowledgments. The authors thank the three anonymous reviewers for providing constructive comments and suggestions, which have greatly improved the quality of this paper.

Chen, X. M., Q. J. Liu, and J. C. Zhang, 2007: A numerical simulation study on microphysical structure and cloud seeding in cloud system of Qilian Mountain region. Meteor. Mon., 33, 33–43.
Guo, X. L., G. G. Zheng, and D. Z. Jin, 2006: A numerical comparison study of cloud seeding by silver iodide and liquid carbon dioxide. Atmos. Res., 79, 183–226. DOI:10.1016/j.atmosres.2005.04.005
Guo, X. L., D. H. Fu, and Z. X. Hu, 2013: Progress in cloud physics, precipitation, and weather modification during 2008–2012. Chinese J. Atmos. Sci., 37, 351–363. DOI:10.3878/j.issn.1006-9895.2012.12321
Guo, X. L., D. H. Fu, X. Y. Li, et al., 2015: Advances in cloud physics and weather modification in China. Adv. Atmos. Sci., 32, 230–249. DOI:10.1007/s00376-014-0006-9
Hong, Y. C., 1998: A 3D hail cloud numerical seeding model. Acta Meteor. Sinica, 56, 641–653. DOI:10.11676/qxxb1998.060
Hong, Y. C., and H. C. Lei, 2012: Research advance and thinking of the cloud precipitation physics and weather modification. Climatic Environ. Res., 17, 951–967. DOI:10.3878/j.issn.1006-9585.2012.06.32
Huang, Y., and H. Y. Xu, 1994: Numerical experiments on hail suppression by AgI seeding. Chinese J. Atmos. Sci., 18, 612–622.
Ji, L., J. Zhou, H. Y. Li, et al., 2014: Cold cloud seeding potentials recognition platform of weather modification based on the BJ-RUC system. Meteor. Mon., 40, 981–991.
Kong, F. Y., M. Y. Huang, and H. Y. Xu, 1990: Three-dimen-sional numerical simulation of ice phase microphysics in cumulus clouds. Part I: Model establishment and ice phase parameterization. Chinese J. Atmos. Sci., 14, 441–453.
Liu, Q. J., and X. F. Lou, 2002: Demands, status and scientific and technical problems on rain enhancement. Sci. Technol. Rev., 53–55, 64.
Ma, Z. S., Q. J. Liu, Y. Y. Qin, et al., 2009: Verification of forecasting efficiency to cloud microphysical characters of mesoscale numerical model for artificial rainfall enhancement by using TRMM satellite data. Acta Meteor. Sinica, 67, 260–271. DOI:10.11676/qxxb2009.026
Mao, Y. H., and Z. J. Hu, 1993: The 2-D numerical study of rain-enhancement and hail-suppression principles on convective clouds. Acta Meteor. Sinica, 51, 184–194. DOI:10.11676/qxxb1993.023
Schwartz, C. S., G. S. Romine, K. R. Smith, et al., 2014: Characterizing and optimizing precipitation forecasts from a convection-permitting ensemble initialized by a mesoscale ensemble Kalman filter. Wea. Forecasting, 29, 1295–1318. DOI:10.1175/WAF-D-13-00145.1
Seiki, T., C. Kodama, M. Satoh, et al., 2015: Vertical grid spacing necessary for simulating tropical cirrus clouds with a high-resolution atmospheric general circulation model. Geophys. Res. Lett., 42, 4150–4157. DOI:10.1002/2015GL064282
Tan, C., Q. J. Liu, and Z. S. Ma, 2013: Influences of sub-grid convective processes on cloud forecast in the GRAPES global model. Acta Meteor. Sinica, 71, 867–878. DOI:10.11676/qxxb2013.067
Tiedtke, M., 1993: Representation of clouds in large-scale models. Mon. Wea. Rev., 121, 3040–3061. DOI:10.1175/1520-0493(1993)121<3040:ROCILS>2.0.CO;2
Wang, Q., and Z. J. Hu, 1990: Three-dimensional elastic atmospheric numerical model and the simulations of a severe storm case. Acta Meteor. Sinica, 48, 91–101. DOI:10.11676/qxxb1990.011
Xu, H. B., 2012: The Physics of Severe Convective Storms and Associated Application. China Meteorological Press, Beijing, 340 pp. (in Chinese)
Xu, H. B., 2014: The Studies of Dynamics in Weather Modification. China Meteorological Press, Beijing, 145 pp. (in Chinese)
Xu, H. B., and S. W. Wang, 1988: Two-dimensional hailcloud modeling. Acta Meteor. Sinica, 46, 227–236. DOI:10.11676/qxxb1988.028
Xu, H. B., and S. W. Wang, 1990: A three-dimensional cloud-scale model suitable for compressible atmosphere. Acta Meteor. Sinica, 48, 80–90. DOI:10.11676/qxxb1990.010
Xu, H. B., and Y. Duan, 2001: The mechanism of hailstone formation and the hail-suppression hypothesis: " Beneficial Competition”. Chinese J. Atmos. Sci., 25, 277–288.
Yin, J. F., D. H. Wang, and G. Q. Zhai, 2015: An attempt to improve Kessler-type parameterization of warm cloud microphysical conversion processes using Cloudsat observations. J. Meteor. Res., 29, 82–92. DOI:10.1007/s13351-015-4091-1
Yu, X., and T.-Y. Lee, 2010: Role of convective parameterization in simulations of a convection band at grey-zone resolutions. Tellus A, 62, 617–632. DOI:10.1111/j.1600-0870.2010.00470.x
Zhang, 1998: Roles of various diabatic physical processes in mesoscale models. Chinese J. Atmos. Sci., 22, 548–561.
Zhang, D.-L., E.-Y. Hsie, and M. W. Moncrieff, 1988: A comparison of explicit and implicit predictions of convective and stratiform precipitating weather systems with a meso-β-scale numerical model . Quart. J. Roy. Meteor. Soc., 114, 31–60. DOI:10.1002/qj.49711447903
Zheng, G. G., and X. L. Guo, 2012: Status and development of sciences and technology for weather modification. Eng. Sci., 14, 20–27.
Zhong, Q., 1992a: A general inverse formulation principle of perfect conservative scheme and its applications. Chinese J. Comput. Phys., 9, 758–764.
Zhong, Q., 1992b: The formulation of a perfect square conservative semi-implicit time difference scheme and its preliminary check. Chinese Sci. Bull., 37, 497–502.
Zhong, Q., J. T. Chen, and Z. L. Sun, 2002: Elimination of computational systematic errors and improvements of weather and climate system models in relation to baroclinic primitive equations. Adv. Atmos. Sci., 19, 1103–1112. DOI:10.1007/s00376-002-0068-y
Zhou, X. P., and D. H. Wang, 1996: Study on the very short-range storm numerical forecasting model. Part I: Theoretical frame. Chinese J. Atmos. Sci., 20, 1–11.