The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, whic...The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, which is a satisfactory approximation to the observed troposphere and lower stratosphere standard atmosphere, is introduced into the equations of the model to permit a more accurate calculation of the vertical transport terms, especially near the tropopause; (2) The vertical levels of the model are carefully selected to guarantee a smooth dependence of layer thickness upon pressure in order to reduce the truncation error involved in the unequal interval vertical finite-differencing. For testing the model, two kinds of linear baroclinic Rossby-Haurwitz waves, one of which has a dynamically stable vertical structure and the other has a relatively unstable one, are constructed to provide initial conditions for numerical experiments. The two waves have been integrated for more than 300 days and 100 days respectively by using the model and both of them are propagating westward with almost identical phase-speed during the time period of the integrations. No obvious change of the wave patterns is found at the levels in the model's troposphere. The amplitudes of both two waves at the uppermost level, however, exhibit rather significant oscillation with time, of which the periods are exactly 20 days and 25 days espectively.The explanation of this interesting phenomena is still under investigation.展开更多
Features of structure and propagation of the 30 to SO day atmospheric oscillations are investigated using the ECMWF analysis of 1980-1983. Evidence is provided to confirm the characteristics of the oscillation in the ...Features of structure and propagation of the 30 to SO day atmospheric oscillations are investigated using the ECMWF analysis of 1980-1983. Evidence is provided to confirm the characteristics of the oscillation in the equatorial region. Those in the mid-high latitudes, however, are revealed to be very different from the tropics and pose a strong barotropic structure. Horizontal coherence shows teleconnection patterns which can be identified as EAP and PNA. The wind field of the specified time scale of the oscillation appears as long-lived vortices and vortex pairs. Mid-latitude perturbations propagate clearly westwards, especially during the winter season. In the high latitudes, they propagate westwards in the winter but eastwards in the summer. Meridional propagations are rather different from region to region.展开更多
Problems of instability of rotating atmospheric motions are investigated by using nonlinear governing equations and the variational principle. The method suggested in this paper is universal for obtaining criteria of ...Problems of instability of rotating atmospheric motions are investigated by using nonlinear governing equations and the variational principle. The method suggested in this paper is universal for obtaining criteria of instability in all models with all possible basic flows. For example, the model can be barotropic or baroclinic, layer or continuous, quasi-geostrophic or primitive equations; the basic flow can be zonal or nonzonal, steady or unsteady.Although the basic flows possess a great deal of variety, they all are the stationary points in the functional space determined by an appropriate invariant functional. The basic flow is an unsteady one if the conservation of angular momentum is included in the associated functional.The second variation, linear or nonlinear, gives the criteria of instability. Especially, the general criteria of instability for unsteady basic flow, orographically disturbed flow as well as nongeostrophic flow are first obtained by the method described in this paper.It is also shown that the difference between the criteria of instability obtained by the linear theory and our variational principle clearly indicates the importance of using nonlinear governing equations.In the appendix the theory is extended to cases such as in a β-plane where the fluid does not possess finite total energy, hence the variational principle can not be directly applied. However, a generalized Liapbunoff norm can still be obtained on the basis of variational consideration.展开更多
In this part, the temporal evolution and interaction across the equator of 30-50 day oscillation in the atmosphere are investigated further. The annual variation of 30-50 day oscillation is quite obvious in the mid-hi...In this part, the temporal evolution and interaction across the equator of 30-50 day oscillation in the atmosphere are investigated further. The annual variation of 30-50 day oscillation is quite obvious in the mid-high latitudes. In the tropical atmosphere, the obvious interannual variation is an important property for temporal evolution of 30-50 day oscillation. The low-frequency wavetrain across the equator over the central Pacific and central Atlantic area, the movement of the long-lived low-frequency system across the equator and the meridional wind component across the equator will obviously show the interaction of 30-50 day oscillation in the atmosphere across the equator.展开更多
The CO_2-seawater system and the method for calculating the partial pressure of CO2 (pCO2) in seawater are studied. The buffer capability of the ocean to increasing atmospheric CO2 is expressed in terms of the differe...The CO_2-seawater system and the method for calculating the partial pressure of CO2 (pCO2) in seawater are studied. The buffer capability of the ocean to increasing atmospheric CO2 is expressed in terms of the differential buffer factor and buffer index. Dissolutions of aragonite and calcite have a significant inffluence on the differential buffer factor. The trend of change in the buffer factor is obtained by a box model.展开更多
In this paper, the principle to determine the atmospheric columnar Mie optical depth from downward total solar radiative flux is theoretically studied, and the effect on Mie optical depth solution of the errors in sur...In this paper, the principle to determine the atmospheric columnar Mie optical depth from downward total solar radiative flux is theoretically studied, and the effect on Mie optical depth solution of the errors in surface albedo, sin-gle scattering albedo, asymmetrical factor of scattering phase function, instrumental constant and the approximate expression of diffusion flux is analy/ed, and then a method for determining surface albedo in shorter wavelength range is presented.展开更多
A new remote sensing method is described to determine the vertical distribution and total content of atmospheric ozone. The method combines surface infrared, satellite infrared and ultraviolet channels. The width of t...A new remote sensing method is described to determine the vertical distribution and total content of atmospheric ozone. The method combines surface infrared, satellite infrared and ultraviolet channels. The width of the infrared channels is 0.01 cm-1, less than Lorentz half-width at the earth's surface, rather than the present width, because these channels can obtain information about variations in the ozone profile below the profile main-peak. The numerical experiments show that the method has a satisfactory precision in determining total ozone content, just about I percent error, and vertical distribution from the earth to 65 km space. In addition, some semi-analysis functions lor calculating backscattered ultraviolet and a relaxation equation are described in this paper.展开更多
The seasonal cycle of the climate of 9000 years before present was simulated with the IAP two-level atmospheric general circulation model.The incoming solar radiation was specified from the orbital parameters for 9000...The seasonal cycle of the climate of 9000 years before present was simulated with the IAP two-level atmospheric general circulation model.The incoming solar radiation was specified from the orbital parameters for 9000 years ago.The boundary conditions of that time were prescribed to the present value because of the small differences between the two.The change in radiation makes temperature to be higher in summer and lower in winter over large areas of the land;and the increased temperature contrast between the land and the ocean strengthens the summer monsoon circulation and increases the precipitation over there.The asymmetry of temperature change between the Northern Hemisphere and the Southern Hemisphere and between summer and winter still exists,which agrees with that get from the previous perpetual experiments.展开更多
In this paper, the winter atmospheric circulation, the convection along the equator and their variations of 1982 and 1983 are investigated. It is suggested that there was a well organized three dimensional structure o...In this paper, the winter atmospheric circulation, the convection along the equator and their variations of 1982 and 1983 are investigated. It is suggested that there was a well organized three dimensional structure of anomalies of the atmospheric circulations during 1982 winter which may be related to the variations of the convection in the equatorial region.展开更多
After my paper (Zeng, 1986b) was published and another (Zeng, 1989) was submitted to the journal, I found two papers written by Arnold (1966) and McIntyre et al. (1987) and received some reprints of Ripa’s papers (19...After my paper (Zeng, 1986b) was published and another (Zeng, 1989) was submitted to the journal, I found two papers written by Arnold (1966) and McIntyre et al. (1987) and received some reprints of Ripa’s papers (1983; 1984; 1987; 1988) in the same field. I thank Drs. Mu Mu and Pedro Ripa very much for showing and sending me these interesting papers.展开更多
The influence of ENSO on the summer climate change in China and its mechanism from the observed data is discussed. It is discovered that in the developing stage of ENSO, the SST in the western tropical Pacific is cold...The influence of ENSO on the summer climate change in China and its mechanism from the observed data is discussed. It is discovered that in the developing stage of ENSO, the SST in the western tropical Pacific is colder in summer, the convective activities may be weak around the South China Sea and the Philippines. As a consequence, the subtropical high shifted southward. Therefore, a drought may be caused in the Indo-China peninsula and in the South China. Moreover, in midsummer the subtropical high is weak over the Yangtze River valley and Huaihe River valley, and the flood may be caused in the area from the Yangtze River valley to Huaihe River valley. On the contrary, in the decaying stage of ENSO. the convective activities may be strong around the Philippines, and the subtropical high shifted northward, a drought may be caused in the Yangtze River valley and Huaihe River valley.展开更多
In this paper, a close relationship between the intraseasonal variation of subtropical high over East Asia and the convective activities around the South China Sea and the Philippines is analysed from OLR data.This re...In this paper, a close relationship between the intraseasonal variation of subtropical high over East Asia and the convective activities around the South China Sea and the Philippines is analysed from OLR data.This relationship is studied by using the theory of wave propagating in a slowly varying medium and by using a quasi-geoslrophic, linear, spherical model and the IAP-GCM, respectively. The results show that when the SST is warming around the western tropical Pacific or the Philippines, the convective activities are intensified around the Philippines. As a consequence, the subtropical high will be intensified over East Asia. The computed results also show that when the anomaly of convective activities are caused around the Philippines, a teleconnection pattern of circulation anomalies will be caused over South Asia, East Asia and North America.展开更多
Non-acceleration theorem in a primitive equation system is developed to investigate the influences of waves on the mean flow variation against external forcing. Numerical results show that mechanical forcing overwhelm...Non-acceleration theorem in a primitive equation system is developed to investigate the influences of waves on the mean flow variation against external forcing. Numerical results show that mechanical forcing overwhelms thermal forcing in maintaining the mean flow in which the internal mechanical forcing associated with horizontal eddy flux of momentum plays the most important roles. Both internal forcing and external forcing are shown to be active and at the first place for the mean flow variations, whereas the forcing-induced mean meridional circulation is passive and secondary. It is also shown that the effects on mean flow of external mechanical forcing are concentrated in the lower troposphere, whereas those due to wave-mean flow interaction are more important in the upper troposphere. These act together and result in the vertically easterly shear in low latitudes and westerly shear in mid-latitudes. This vertical shear of mean flow is to some extent weakened by thermal forcing.展开更多
A brief description of the dynamic framework of the IAP 2-L AGCM (referred as 'IAP DF') is presented in contrast with the corresponding 'OSU DF' and preliminary comparison and examination of the IAP DF...A brief description of the dynamic framework of the IAP 2-L AGCM (referred as 'IAP DF') is presented in contrast with the corresponding 'OSU DF' and preliminary comparison and examination of the IAP DF and the OSU DF with the aid of numerical experiments using a set of baroclinic Rossby-Haurwitz wave initial conditions is performed.The results of the numerical experiments show that both the IAP DF arid the OSU DF are long-term computational stable and are able to simulate the fundamental behavior of the Rossby-Haurwitz wave under the limitation of 7.5 minutes time interval. In the respect of the prediction of the large-scale wave, the IAP DF is better than the OSU DF, however, the OSU DF gives better predictions of zonal mean status than the IAP DF does. Another important contrast is that the level of small-scale noise of the IAP DF is higher than the OSU DF, on the contrary, the damping effect of the OSU DF on large-scale vortex field is stronger than the IAP DF.展开更多
A dynamic clustering method based on multispectral satellite imagery to identify the different features is described. The channel combinations selected are for the different purposes in classification. Several cases a...A dynamic clustering method based on multispectral satellite imagery to identify the different features is described. The channel combinations selected are for the different purposes in classification. Several cases are presented using the polar-orbiting satellite imageries.展开更多
This paper introduces some advanced subjects on lidar remote sensing of the atmosphere, emphasizing recent studies and developments in lidar application for measuring ozone, cloud, aerosol, atmospheric temperature, mo...This paper introduces some advanced subjects on lidar remote sensing of the atmosphere, emphasizing recent studies and developments in lidar application for measuring ozone, cloud, aerosol, atmospheric temperature, moisture, pressure and wind.展开更多
Some more proper criteria for the nonlinear stability of three-dimensional quasi-geostrophic motions are given by combining variational principle with a prior estimates method. The criteria are suitable for perturbati...Some more proper criteria for the nonlinear stability of three-dimensional quasi-geostrophic motions are given by combining variational principle with a prior estimates method. The criteria are suitable for perturbations of initial condition as well as parameters in the model. The basic flow can be steady or unsteady. Particularly the difficulty due to the nonlinear boundary condition is completely overcome by the use of our method.展开更多
The gravity wave breaking is crucial to the large-scale circulation of middle atmosphere. In this paper, we follow Lindzen (1981) to draw out the parameterization of two-dimensional gravity wave breaking including ine...The gravity wave breaking is crucial to the large-scale circulation of middle atmosphere. In this paper, we follow Lindzen (1981) to draw out the parameterization of two-dimensional gravity wave breaking including inertial effect. Also we present some properties of critical levels and inertial critical levels.展开更多
The introduction of 'hydrostatic extraction' scheme, or 'standard stratification approximation', into spectral model gained some advantages compared with commonly used schemes. However, computational i...The introduction of 'hydrostatic extraction' scheme, or 'standard stratification approximation', into spectral model gained some advantages compared with commonly used schemes. However, computational instability may occur for high vertical resolution versions if the stratification parameter C0 taken as a constant. In this paper, the possible cause leading to the instability is discussed and an improved scheme presented where C0 is generalized to be a function of both height and latitudes. Hence the reference atmosphere gets closer to the real atmosphere and the temperature deviation field to be expanded becomes smoother everywhere. Test by real case forecasts shows good computational stability of the new scheme and better prediction performance than-usual schemes of spectral model.展开更多
The Bayes Decision method (BD) has been implemented on GMS Stretched VISSR (S.VISSR) images. Several objects are identified. These are earth's surface, stratiform, convective and cirrus clouds. Overlaid 6-hour sur...The Bayes Decision method (BD) has been implemented on GMS Stretched VISSR (S.VISSR) images. Several objects are identified. These are earth's surface, stratiform, convective and cirrus clouds. Overlaid 6-hour surface rainfall observations show that the rain area on the ground matches with precipitable cloud area identified by BD method. Finally the data of classified image are compressed to about one-eighth amount by use of Line-Run-Coding technique. It is helpful for data archive and remote transmission.展开更多
文摘The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, which is a satisfactory approximation to the observed troposphere and lower stratosphere standard atmosphere, is introduced into the equations of the model to permit a more accurate calculation of the vertical transport terms, especially near the tropopause; (2) The vertical levels of the model are carefully selected to guarantee a smooth dependence of layer thickness upon pressure in order to reduce the truncation error involved in the unequal interval vertical finite-differencing. For testing the model, two kinds of linear baroclinic Rossby-Haurwitz waves, one of which has a dynamically stable vertical structure and the other has a relatively unstable one, are constructed to provide initial conditions for numerical experiments. The two waves have been integrated for more than 300 days and 100 days respectively by using the model and both of them are propagating westward with almost identical phase-speed during the time period of the integrations. No obvious change of the wave patterns is found at the levels in the model's troposphere. The amplitudes of both two waves at the uppermost level, however, exhibit rather significant oscillation with time, of which the periods are exactly 20 days and 25 days espectively.The explanation of this interesting phenomena is still under investigation.
文摘Features of structure and propagation of the 30 to SO day atmospheric oscillations are investigated using the ECMWF analysis of 1980-1983. Evidence is provided to confirm the characteristics of the oscillation in the equatorial region. Those in the mid-high latitudes, however, are revealed to be very different from the tropics and pose a strong barotropic structure. Horizontal coherence shows teleconnection patterns which can be identified as EAP and PNA. The wind field of the specified time scale of the oscillation appears as long-lived vortices and vortex pairs. Mid-latitude perturbations propagate clearly westwards, especially during the winter season. In the high latitudes, they propagate westwards in the winter but eastwards in the summer. Meridional propagations are rather different from region to region.
文摘Problems of instability of rotating atmospheric motions are investigated by using nonlinear governing equations and the variational principle. The method suggested in this paper is universal for obtaining criteria of instability in all models with all possible basic flows. For example, the model can be barotropic or baroclinic, layer or continuous, quasi-geostrophic or primitive equations; the basic flow can be zonal or nonzonal, steady or unsteady.Although the basic flows possess a great deal of variety, they all are the stationary points in the functional space determined by an appropriate invariant functional. The basic flow is an unsteady one if the conservation of angular momentum is included in the associated functional.The second variation, linear or nonlinear, gives the criteria of instability. Especially, the general criteria of instability for unsteady basic flow, orographically disturbed flow as well as nongeostrophic flow are first obtained by the method described in this paper.It is also shown that the difference between the criteria of instability obtained by the linear theory and our variational principle clearly indicates the importance of using nonlinear governing equations.In the appendix the theory is extended to cases such as in a β-plane where the fluid does not possess finite total energy, hence the variational principle can not be directly applied. However, a generalized Liapbunoff norm can still be obtained on the basis of variational consideration.
基金This study was supported in part by National Natural Science Foundation of China
文摘In this part, the temporal evolution and interaction across the equator of 30-50 day oscillation in the atmosphere are investigated further. The annual variation of 30-50 day oscillation is quite obvious in the mid-high latitudes. In the tropical atmosphere, the obvious interannual variation is an important property for temporal evolution of 30-50 day oscillation. The low-frequency wavetrain across the equator over the central Pacific and central Atlantic area, the movement of the long-lived low-frequency system across the equator and the meridional wind component across the equator will obviously show the interaction of 30-50 day oscillation in the atmosphere across the equator.
文摘The CO_2-seawater system and the method for calculating the partial pressure of CO2 (pCO2) in seawater are studied. The buffer capability of the ocean to increasing atmospheric CO2 is expressed in terms of the differential buffer factor and buffer index. Dissolutions of aragonite and calcite have a significant inffluence on the differential buffer factor. The trend of change in the buffer factor is obtained by a box model.
文摘In this paper, the principle to determine the atmospheric columnar Mie optical depth from downward total solar radiative flux is theoretically studied, and the effect on Mie optical depth solution of the errors in surface albedo, sin-gle scattering albedo, asymmetrical factor of scattering phase function, instrumental constant and the approximate expression of diffusion flux is analy/ed, and then a method for determining surface albedo in shorter wavelength range is presented.
文摘A new remote sensing method is described to determine the vertical distribution and total content of atmospheric ozone. The method combines surface infrared, satellite infrared and ultraviolet channels. The width of the infrared channels is 0.01 cm-1, less than Lorentz half-width at the earth's surface, rather than the present width, because these channels can obtain information about variations in the ozone profile below the profile main-peak. The numerical experiments show that the method has a satisfactory precision in determining total ozone content, just about I percent error, and vertical distribution from the earth to 65 km space. In addition, some semi-analysis functions lor calculating backscattered ultraviolet and a relaxation equation are described in this paper.
文摘The seasonal cycle of the climate of 9000 years before present was simulated with the IAP two-level atmospheric general circulation model.The incoming solar radiation was specified from the orbital parameters for 9000 years ago.The boundary conditions of that time were prescribed to the present value because of the small differences between the two.The change in radiation makes temperature to be higher in summer and lower in winter over large areas of the land;and the increased temperature contrast between the land and the ocean strengthens the summer monsoon circulation and increases the precipitation over there.The asymmetry of temperature change between the Northern Hemisphere and the Southern Hemisphere and between summer and winter still exists,which agrees with that get from the previous perpetual experiments.
文摘In this paper, the winter atmospheric circulation, the convection along the equator and their variations of 1982 and 1983 are investigated. It is suggested that there was a well organized three dimensional structure of anomalies of the atmospheric circulations during 1982 winter which may be related to the variations of the convection in the equatorial region.
文摘After my paper (Zeng, 1986b) was published and another (Zeng, 1989) was submitted to the journal, I found two papers written by Arnold (1966) and McIntyre et al. (1987) and received some reprints of Ripa’s papers (1983; 1984; 1987; 1988) in the same field. I thank Drs. Mu Mu and Pedro Ripa very much for showing and sending me these interesting papers.
文摘The influence of ENSO on the summer climate change in China and its mechanism from the observed data is discussed. It is discovered that in the developing stage of ENSO, the SST in the western tropical Pacific is colder in summer, the convective activities may be weak around the South China Sea and the Philippines. As a consequence, the subtropical high shifted southward. Therefore, a drought may be caused in the Indo-China peninsula and in the South China. Moreover, in midsummer the subtropical high is weak over the Yangtze River valley and Huaihe River valley, and the flood may be caused in the area from the Yangtze River valley to Huaihe River valley. On the contrary, in the decaying stage of ENSO. the convective activities may be strong around the Philippines, and the subtropical high shifted northward, a drought may be caused in the Yangtze River valley and Huaihe River valley.
文摘In this paper, a close relationship between the intraseasonal variation of subtropical high over East Asia and the convective activities around the South China Sea and the Philippines is analysed from OLR data.This relationship is studied by using the theory of wave propagating in a slowly varying medium and by using a quasi-geoslrophic, linear, spherical model and the IAP-GCM, respectively. The results show that when the SST is warming around the western tropical Pacific or the Philippines, the convective activities are intensified around the Philippines. As a consequence, the subtropical high will be intensified over East Asia. The computed results also show that when the anomaly of convective activities are caused around the Philippines, a teleconnection pattern of circulation anomalies will be caused over South Asia, East Asia and North America.
基金Research Project No.[75-09-01] on medium-range numerical weather forecasts.
文摘Non-acceleration theorem in a primitive equation system is developed to investigate the influences of waves on the mean flow variation against external forcing. Numerical results show that mechanical forcing overwhelms thermal forcing in maintaining the mean flow in which the internal mechanical forcing associated with horizontal eddy flux of momentum plays the most important roles. Both internal forcing and external forcing are shown to be active and at the first place for the mean flow variations, whereas the forcing-induced mean meridional circulation is passive and secondary. It is also shown that the effects on mean flow of external mechanical forcing are concentrated in the lower troposphere, whereas those due to wave-mean flow interaction are more important in the upper troposphere. These act together and result in the vertically easterly shear in low latitudes and westerly shear in mid-latitudes. This vertical shear of mean flow is to some extent weakened by thermal forcing.
文摘A brief description of the dynamic framework of the IAP 2-L AGCM (referred as 'IAP DF') is presented in contrast with the corresponding 'OSU DF' and preliminary comparison and examination of the IAP DF and the OSU DF with the aid of numerical experiments using a set of baroclinic Rossby-Haurwitz wave initial conditions is performed.The results of the numerical experiments show that both the IAP DF arid the OSU DF are long-term computational stable and are able to simulate the fundamental behavior of the Rossby-Haurwitz wave under the limitation of 7.5 minutes time interval. In the respect of the prediction of the large-scale wave, the IAP DF is better than the OSU DF, however, the OSU DF gives better predictions of zonal mean status than the IAP DF does. Another important contrast is that the level of small-scale noise of the IAP DF is higher than the OSU DF, on the contrary, the damping effect of the OSU DF on large-scale vortex field is stronger than the IAP DF.
文摘A dynamic clustering method based on multispectral satellite imagery to identify the different features is described. The channel combinations selected are for the different purposes in classification. Several cases are presented using the polar-orbiting satellite imageries.
文摘This paper introduces some advanced subjects on lidar remote sensing of the atmosphere, emphasizing recent studies and developments in lidar application for measuring ozone, cloud, aerosol, atmospheric temperature, moisture, pressure and wind.
文摘Some more proper criteria for the nonlinear stability of three-dimensional quasi-geostrophic motions are given by combining variational principle with a prior estimates method. The criteria are suitable for perturbations of initial condition as well as parameters in the model. The basic flow can be steady or unsteady. Particularly the difficulty due to the nonlinear boundary condition is completely overcome by the use of our method.
文摘The gravity wave breaking is crucial to the large-scale circulation of middle atmosphere. In this paper, we follow Lindzen (1981) to draw out the parameterization of two-dimensional gravity wave breaking including inertial effect. Also we present some properties of critical levels and inertial critical levels.
基金This work has been carried out under the support of the Medium-range Numerical Weather Forecast research project
文摘The introduction of 'hydrostatic extraction' scheme, or 'standard stratification approximation', into spectral model gained some advantages compared with commonly used schemes. However, computational instability may occur for high vertical resolution versions if the stratification parameter C0 taken as a constant. In this paper, the possible cause leading to the instability is discussed and an improved scheme presented where C0 is generalized to be a function of both height and latitudes. Hence the reference atmosphere gets closer to the real atmosphere and the temperature deviation field to be expanded becomes smoother everywhere. Test by real case forecasts shows good computational stability of the new scheme and better prediction performance than-usual schemes of spectral model.
文摘The Bayes Decision method (BD) has been implemented on GMS Stretched VISSR (S.VISSR) images. Several objects are identified. These are earth's surface, stratiform, convective and cirrus clouds. Overlaid 6-hour surface rainfall observations show that the rain area on the ground matches with precipitable cloud area identified by BD method. Finally the data of classified image are compressed to about one-eighth amount by use of Line-Run-Coding technique. It is helpful for data archive and remote transmission.