ABSTRACT

Cloud detection and classification form a basis in weather analysis. Split window algorithm (SWA) is one of the simple and matured algorithms used to detect and classify water and ice clouds in the atmosphere using satellite data. The recent availability of Himawari-8 data has considerably strengthened the possibility of better cloud classification owing to its enhanced multi-band configuration as well as high temporal resolution. In SWA, cloud classification is attained by considering the spatial distributions of the brightness temperature (BT) and brightness temperature difference (BTD) of thermal infrared bands. In this study, we compare unsupervised classification results of SWA using the band pair of band 13 and 15 (SWA13-15, 10 and 12 µm bands), versus that of band 15 and 16 (SWA15-16, 12 and 13 µm bands) over the Japan area. Different threshold values of BT and BTD are chosen in winter and summer seasons to categorize cloud regions into nine different types. The accuracy of classification is verified by using the cloud-top height information derived from the data of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). For this purpose, six different paths of the space-borne lidar are selected in both summer and winter seasons, on the condition that the time span of overpass falls within the time ranges between 01:00 and 05:00 UTC, which corresponds to the local time around noon. The result of verification indicates that the classification based on SWA13-15 can detect more cloud types as compared with that based on SWA15-16 in both summer and winter seasons, though the latter combination is useful for delineating cumulonimbus underneath dense cirrus regions.

Keywords:

Cloud Type Detection, Himawari-8, Split Window Algorithm, Brightness Temperature

1. Introduction

Precise classification of cloud types in satellite images is indispensable for the better understanding of cloud formation and development processes. While low Earth orbit satellites provide cloud images only once or twice a day, much more frequent images are available from meteorological satellites in the geostationary orbit. Himawari-8 satellite, launched on October 7, 2014 and operated by the Japan Meteorological Agency (JMA), started the official data dissemination from July 7, 2015. One of the advantages of Himawari-8 satellite data is the very high frequency of data acquisition―10 min for the full disk and 2.5 min for the area around Japan. Himawari-8 is equipped with a sensor called Advanced Himawari Imager (AHI), which has a total of 16 bands in the visible (VIS, 3 bands), near-infrared (NIR, 3 bands) and infrared (IR, 10 bands) portion of the electromagnetic spectrum [1]. As compared with the former meteorological satellite of JMA (MTSAT, Himawari-7) that had a single, panchromatic channel in the visible part, a true-color image can be formed by assigning blue, green, and red to band 1, 2 and 3 of AHI, respectively.

Detection of clouds and estimation of cloud physical parameters are important topics in the data analysis of meteorological satellites. Previous works have revealed, for instance, that thick cloud can be detected using 0.64 µm channel reflectance (R0.64) using a certain threshold [2], though this technique is not valid over the bright surface [3]. Cloud parameters, on the other hand, can be estimated by means of the information provided from both NIR and IR bands. For example, optical depth and ice particle size of cirrus clouds were retrieved using 1.38 and 1.88 µm channels [4], which correspond to band 5 (centered at 1.61 µm) of Himawari-8. During the daytime, the reflectance of visible and NIR band are more informative than during the night time. According to Bessho et al. (2016) [1], band 7, a mid-IR band centered at 3.88 µm, has the capability of detecting low-level cloud and fog. During the daytime, low-level water clouds are effectively detected by using the brightness temperature difference (BTD) between the 11 and 3.7 µm bands (BTD11-3.7) [2]. Bessho et al. (2016) [1]indicated that water vapor bands located at 6.2, 6.9 and 7.3 µm are useful for cloud phase detection. In relation to the BTD and cloud properties, the use of TIR bands were proposed for other meteorological satellites, such as the cloud top height estimation using MODIS [5]and Meteosat [6], cirrus and water vapor detection using Meteosat [7], cloud phase and cloud optical thickness using Meteosat [8], and rainfall detection [9]. Also, the movement of water vapor detected by means of 6.2 and 7.3 µm bands (BTD 6.2-7.3) was used to study wind movement and convection [6]. Strabala et al. (1994) [10]reported cloud type detection using BTD values from three spectral bands at 8, 11, and 12 µm. In their work, cloud phase was determined by comparing BTD8-11 and BTD11-12, because of the observation that the scattering of water and ice clouds is different near 12 µm. Also, seven types of clouds were discriminated on the basis of the comparison between BT11.2 and BTplev11.2 of Himawari-8 IR band [11]: here “plev” denotes a constant pressure value between surface and tropopause. Shang et al. (2017) [12]modelled a formula using BT11.2, BTD11.2-3.9 and R0.64 to delineate cold cloud over inland water, liquid/ice cloud, and cloud over land, respectively.

Split window algorithm (SWA) has widely been used for the analysis and classification of satellite imagery. Recent researches used SWA for the estimation of land surface temperature [13][14][15][16][17]as well as the estimation of cirrus cloud top height [18]. In the region of cloud classification, SWA was first introduced by Inoue (1987) [19]by utilizing two-dimensional threshold values of BT and BTD. Subsequently, Inoue [19]used the 11 and 12 µm channels of NOAA 7 to detect the cumulonimbus (Cb), cumulus (Cu), dense cirrus (Dc) and cirrus (Ci) clouds over the tropical ocean. It was found that distinguishing Ci cloud from the land area was difficult, especially in the night time, though the threshold value was not explicitly shown. Hamada et al. (2004) [20]tried to make cloud classification using GMS-5 VISSR by modifying Inoue’s approach, and identified cloud types by using ground-based observation. The threshold values of BTD11-12 and the BT11 were determined to be 1.8 and 250 K, respectively. Inoue continued his research in cloud classification using SWA over the eastern sub-tropical Pacific area [21]. They used 11 and 12 µm bands of NOAA-9 data in combination with coincident and collocated Earth Radiation Budget Experiment (ERBE) S-8 data. As a result, generally three cloud types, namely, Cb, Cu, and Ci, were classified using BT and BTD thresholds. Three types of Ci were found as inferred from the thickness and cloud-top temperature. Subsequently, Lutz et al. (2003) [5]reported cloud type classification by applying the SWA using 11 and 12 µm bands in the Meteosat Second Generation (MSG) method. Although high-Cb, middle-Cb, low-Cb, Dc, thick cirrus, Ci, and thin cirrus were detected separately, their approach could not distinguish between ice and water clouds. A brief overview of the previous SWA researches is shown in Table 1.

In order to validate the classification results based on Himawari-8 imagery, we utilize the observation data of a space-borne lidar, namely, the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) [22]. CALIPSO was launched in 2006 by NASA and CNES to the low-Earth orbit with the altitude of approximately 690 km above surface. A Mie-scattering lidar provides the vertical profiles of aerosols and clouds at the two wavelengths generated from a Nd:YAG laser (532 and 1064 nm). Also, the information on water/ice cloud

Note: Cb = Cumulonimbus; Cu = Cumulus; DCi = Dense Cirrus; Ci = Cirrus; Hi = High; Mid = Middle; Lo = Low.

phase is available using the depolarization ratio of the lidar backscatter signal [22].

The remaining part of this paper is organized as follows. In Section 2, the methodology of present work will be described, with the pertinent information on data processing and validation. In Section 3, the results are shown for particular cases chosen with the criteria of concurrent observation with the Himawari-8 satellite and CALIPSO, and comparison is made between the two pairs of IR bands for implementing SWA. Finally, conclusions will be given in Section 4.

2. Method

2.1. Band Selection of Himawari-8

Himawari-8 data are obtained from the archive of the Center for Environmental Remote Sensing (CEReS), Chiba University (http://www.cr.chiba-u.jp/japanese/database.html), which provides the original and gridded data to the science community under the contract with JMA. In terms of energy transfer process between the surface and atmosphere, the atmospheric window wavelength is ranged between 8 - 12 µm, which is largely or partially transparent to most of atmospheric gases [23]. As seen from Table 1, most of the previous researches employed the combination of 10 and 11 µm bands for implementing SWA.

Figure 1 shows an example of BT histograms of a Himawari-8 image (observed on April 1, 2016 around Japan area) for ten IR bands covering the wavelength range of 3.9 - 13.3 µm. Similar comparison of histograms has also been made for images captured in winter and summer seasons. Based on the shapes of

the histograms of atmospheric window bands (band 13, 14, 15 and 16) seen in Figure 1, here we select the band pairs of 13 - 15 and 15 - 16 for implementing cloud classification. Among the four bands, band 13 exhibits the narrowest histogram. The wavelength difference of the combination of band 13 - 14 or band 14 - 15 is too small, which may result in the loss of accuracy in the calculation of BTD. In the present work, we compare the classification accuracy for the two cases of SWA13-15 (inside the atmospheric window range) and SWA15-16 (slightly outside the atmospheric window range). The area of interest is the Japan area, at which the Himawari-8 imagery is available every 2.5 min. The Himawari-8 bands utilized in this research are band 1 (centered at 0.47 µm), 4 (0.86 µm) 13 (10.4 µm), 15 (12.4 µm) and 16 (13.3 µm).

2.2. CALIPSO Data

CALIPSO data are obtained from the Atmospheric Science Data Center of NASA Langley Research Center. The BT information from satellite observation is pertinent to the temperature of cloud top, and CALIPSO lidar data are useful for validating the cloud-top height. Concurrent data sets of Himawari-8 and CALIPSO have been selected during winter and summer seasons in 2016. Among the six cases summarized in Table 2, January 1, January 6, and March 19 represent the winter season, while June 1, July 2, August 29 the summer season. These dates were selected considering the availability of the data of CALIPSO that passes over the Chiba area in the time range of 01:00 - 05:00 UTC or 10:00 - 14:00 Japan Standard Time (JST). Attenuated backscattering, depolarization ratio, and ice/water phase parameter of CALIPSO data are used in this research. The Himawari-8 data that coincide with the overpass time are downloaded in the geo-corrected format cropped in the Japan area (22.02˚N - 47.74˚N and 120.11˚E - 156.99˚E). In addition to the high frequency observation, ten minutes interval cloud image movie (hereafter image movie) is available by using band 4

for the time range of around 00:00 - 06:00 UTC (09:00 - 15:00 JST). This tool helps to facilitate the visual interpretation process during the cloud interpretation and threshold adjustment.

2.3. Albedo Calculation

Band 1 data of Himawari-8 centered at 0.47 µm are used for the simple detection of cloud areas characterized with high reflectance (albedo) values. In this band, most of bright objects are considered to be cloud, though snow cover in land areas exhibits similarly high reflectance especially in the middle latitude zone. The calculation of albedo, A, is carried out by using the following formula:

$I=\text{Gain}\ast \text{Count}+\text{Constant},A=c^{\prime }I.$ (1)

Here I is the radiance, which is calculated for each pixel from the digital count of Himawari-8 band 1. Then, the albedo value, A, is calculated by multiplying the radiance with the transformation coefficient, c’ [24]. In Equation (1), the values of Gain, Constant and c' are derived from the header information file of Himawari-8. The resulting value of A is between 0 and 1, and the present examination of a number of pixels has led to a threshold value of A = 0.2, i.e., if A > 0.2, it is considered as a cloud pixel for daytime images.

2.4. BT and BTD Calculation

The BT value is derived for band 13 (hereafter BT13), band 15 (BT15) and 16 (BT16) of Himawari-8 using the principle of Planck’s radiation law. The BT values of these three bands are employed also for calculating BTD13-15 and BTD15-16. For an IR band with a center wavelength of λ, the value of BT can be calculated as [24]:

$BT=c_0+c_1{T}_{\text{e}}+c_2{T}_{\text{e}}^2$ (2)

and

$T_{\text{e}}=\frac{hc}{k\lambda }{\left[\ln \left(\frac{2h{c}^2}{{\lambda }^{5}I}+1\right)\right]}^{-1}.$ (3)

Here, c₀, c₁, and c₂ are constants that can be obtained from the header information file of Himawari-8 for each band. The radiance, I, is calculated for each pixel from the digital count of each band using Equation (1); h is the Planck constant, k is the Boltzmann constant, and c is the speed of light. The BTD13-15 is calculated as the difference between BT13 and BT15 while BTD15-16 between BT15 and BT16.

2.5. SWA Implementation

Here we examine two schemes of SWA, namely SWA13-15 that is based on TB13 (x axis) and BTD13-15 (y axis) and SWA15-16 based on TB15 and BTD15-16. Figure 2 shows the nine types of cloud classes resulting from the present unsupervised classification. In this matrix, lower values of BT correspond to higher clouds, and smaller values of BTD correspond to thicker clouds. Cirrus (Type #8) and thin cirrus (type #9) are exceptional in that they exhibit higher values of BT, which generally suggest middle or low-level clouds. This is due to the influence of the surface below the cirrus or thin cirrus clouds that have semi-transparent features. The threshold values (BT-1, BT-2, BTD-1, and BTD-2) are tabulated in Table 3 for the two SWA schemes for both winter and summer cases: see Section 3 below for the determination procedure of these threshold values.

2.6. Verification of Cloud Phase and Altitude

The depolarization ratio (ρ) provided in the CALIPSO data is employed here for

checking the cloud phase. The water and ice cloud typically have ρ less than 0.1 and more than 0.3, respectively. The value of ρ between 0.1 and 0.3 indicates mixed-phase cloud [25]. Generally high-level clouds tend to exhibit large ρ values, and low-level cloud small ρ values.

Clouds can be roughly classified into three levels on the basis of their altitudes, namely, low-level (<3 km), middle-level (3 - 6 km), and high-level cloud (>6 km). Because of the average lapse rate of troposphere (~6.5 K/km), the difference in altitude leads to the difference in BT. Thus, the cloud top height, h, can be estimated from the value of BT observed for the cloud pixel as [26]

$h=\frac{1}{\delta }\left(T_{\text{s}}-T_{\text{CT}}\right).$ (4)

Here δ = 6.5 K/km is the standard value of the vertical lapse rate, T_S is the surface temperature, and T_CT is the cloud top temperature. In the present work, we apply Equation (4) to check the cloud top height over the Chiba area by assuming that T_CT is equal to BT13 with T_S value observed at Chiba University.

3. Results and Discussion

Six datasets of Himawari-8 imagery listed in Table 2 are processed and verified during the test. Figure 3(a) shows the albedo distribution generated from the band 1 image observed at 03:40 UTC (12:40 JST) on July 2, 2017, while Figure 3(b) the result of cloud extraction using the threshold value of A > 0.2. As compared with other visible bands, band 1 centered at 0.47 µm gives the best discrimination of cloud and non-cloud (surface) areas. Figures 3(c)-(h) show albedo images for each of the dataset listed in Table 2. In Figure 3(c), streak patterns due to the formation of cumulus clouds are seen over the Japan Sea area. These clouds often develop into cumulonimbus due to the orographic effect, resulting in heavy snowfall in the northern part of central Japan. Both in panels (d) and (e), frontal clouds are seen in the southern part of the Honshu Island. In panels (f) and (g), clouds due to the Baiu front are seen in addition to the cloud system related to a low pressure in the northern part. A developed cloud system in association with a typhoon is seen in panel (h). In addition, it is noted that in winter season such as in panel (c), some land areas in northern Japan are covered with snow, which tend to be misclassified as cloud areas. The examination of image movie, as mentioned above, is useful for discriminating clouds from snow-covered surface areas.

The threshold values of BT and BTD listed in Table 3 have been determined through the examination of a number of images (not limited to the six datasets in Table 2) with the help of the image movie. The visual observation of the movie greatly facilitated the assignment of cloud types based on their movement and appearance. In the following, some explanations are given on the recognition of cloud characteristics.

Cumulonimbus clouds (Cb, type #1) always exhibit high brightness in the albedo image, indicating the very thick nature as a result of vertical development.

Sometimes they are overlapped with a large area of dense and thick cirrus clouds (Dc, type #4) that move together [27], as actually seen in the case of Figure 3(e). Also, they are detected around a typhoon, as seen in Figure 3(h). Sometimes it appears as middle-cloud without anvil as middle cumulonimbus (Mi-Cb, type #2) and ice cloud (type #5). In the pre-development stage of cumulonimbus, cumulus clouds (type #3) are detected as low-level clouds.

Low-level cloud is mostly dominated by water cloud (type #6) with a very large spatial coverage and non-uniform appearance. Cumulus clouds (Cu, type #3) form as small clusters. Since they are low-level clouds, their movement is often affected by the existence of high mountains. Stratocumulus (Sc) and nimbostratus (Ns) exhibit similar appearance as cumulus, but with larger spatial coverage. Nimbostratus (Ns) is rain cloud with substantial thickness. In some cases, we found difficulty in delineating low-level clouds only from their appearance in image movie because of their similarity and lack of spatial resolution.

It is found that the movement of cirrus cloud (Ci, type #8) is mostly on top of other cloud types. Cirrus clouds are characterized with semi-transparent and smooth texture. Thin cirrus (type #9) and cirrus are high clouds with cloud top height of more than 6 km.

Type #4 is dense cirrus (Dc): they appear as high clouds, sometimes forming an anvil during the development process of cumulonimbus (type #1). The appearance is always very dense with wide spatial coverage. Sometimes cirrus clouds exhibit non-uniform shapes and located in the surrounding of the dense cirrus as ice clouds (type #5).

The determination of threshold values in BT and BTD is carried out by constructing scatter diagrams. Figure 4 shows examples of scatter diagram for dataset #1 (Figure 4(a) and Figure 4(b)) and dataset #5 (Figure 4(c) and Figure 4(d)) analyzed by means of SWA13-15 (Figure 4(a) and Figure 4(c)) and SWA15-16 (Figure 4(b) and Figure 4(d)). In the case of dataset #5, the mean value of BTD13-15 is 4.32 K, which is substantially smaller than that of BTD15-16 (11.68 K). The maximum values, on the other hand, are similar, 27.86 K for BTD13-15 and 27.91 K for BTD15-16.

Wider spread is seen for the scatter plot of SWA13-15 than that for SWA15-16. From the scatter plot of SWA15-16 (Figure 4(d)), it is understood that cumulus clouds with high BT and low BTD cannot be detected: this is in clear contrast with SWA13-15, where cumulus is successfully detected. The detection of thin cirrus (with high BT and BTD), on the other hand, can be facilitated using the scheme of SWA15-16 (Figure 4(d)). The similar trend of scatter plot diagram is found in other datasets, e.g. dataset #1. The shape of the BTD13-15 tends to spread more widely than BTD15-16. When compared with the winter diagram (Figure 4(a) and Figure 4(b)), summer scatter plots (Figure 4(c) and Figure 4(d)) tends to spread more widely on both BT (x axis) and BTD (y axis).

The results of final cloud classification, made by combining the results of cloud area extraction (Figure 3(b)) and split window implementation, are shown in Figure 5 for the two cases of (A) SWA13-15 and (B) SWA15-16. In Figure 5(A) based on SWA13-15, the most noticeable cloud type is dense cirrus

(type #4 in red color), followed by water cloud that includes stratocumulus and stratus (type #6, light blue). Also, cumulus clouds (type #3, dark green) are found in most of the panels, especially in the northern part of panel (c). Ice clouds (type #5, pink) are found in the northern part of panel (b), in addition to the rim of dense cirrus. In panel (d), a region of thick cirrus (type #7, light green) appears over that of thin cirrus clouds (type #9, grey).

In Figure 5(B) based on SWA15-16, on the other hand, the distribution of dense cirrus (type #4, red) is more or less similar to that in Figure 5(A), though central portions are sometimes replaced with cumulonimbus (type #1, blue), especially in images in summer season, i.e., panels (d)-(f). Since these cases represent rainfall events in relation to the frontal system ((d) and (e)) or a typhoon (f), the result in Figure 5(B) is more reasonable than that in Figure 5(A). The spread of water cloud (type #6, light blue) is also similar to the case in Figure 5(A), but most of cumulus clouds (type #3, green) that have been seen in Figure 5(A) are not found in Figure 5(B). Similarly, thick cirrus regions (type #7, light green) are not seen in the case of Figure 5(B).

All the results of cloud type detection are summarized in Table 4. From this table, it is apparent that generally SWA13-15 exhibits more sensitivity for different cloud types. SWA13-15 can detect by 50%, 43%, 29%, 29%, 38%, and 17% more cloud types than SWA15-16 for dataset #1 through #6, respectively. However, SWA15-16 gives more appropriate discrimination between dense cirrus and cumulonimbus. Although SWA15-16 exhibits rather limited capability in discriminating different cloud types, the results are not contradictory to the cloud types inferred from SWA13-15. For example, cumulus clouds (type #3, dark green) that have been detected on the northern part of panels (a), (b) and (c) using SWA13-15 are detected as water clouds (type #6, light blue) using SWA15-16. This is not an incompatible situation, since the cumulus type is part of water cloud. Similarly, thick cirrus clouds (type #7, light green) detected on SWA13-15 in panels (d) and (e) are detected as dense cirrus (type #4, red color) using SWA15-16, but both types belong to the cirrus cloud family.

In both Figure 5(A) and Figure 5(B), snow-covered land areas have been

Note: D = Detected; ND = Not detected; SWA = Split Window Algorithm.

detected as water cloud because of the similarity in albedo values. Usually such situations occur during winter season in the northern area of Japan. Practically speaking, however, the snow areas can easily be delineated using the image movie of Himawari-8. The other difficulty encountered in the present SWA is the precise detection of very thin cirrus clouds, which are mostly detected as clear area.

The verification of cloud classification results is implemented using the CALIPSO data. Figure 6(a) shows an example, where the CALIPSO flight track is superposed on the map of cloud classification based on SWA13-15 for the case of 03:40 UTC on July 2, 2016. Two ranges (A and B) are defined along the flight track to examine the temporal change of parameters, namely, the depolarization ratio (ρ), ice/water phase, and the cloud type. The time-height indicator (THI) representation of the depolarization ratio is shown in Figure 6(b), while that of the ice/water phase in Figure 6(c), both of which have been downloaded from

the CALIPSO website. The level of depolarization ratio is indicated in the track using different colors showing no value (yellow), ρ < 0.3 (brown), and ρ > 0.3 (dark brown). This last category stands for ice cloud, which corresponds to dense cirrus (type #4, red) or cumulonimbus (type #1, blue). Thus, the present result of cloud classification around region A (dense cirrus) can be confirmed using the large ρ value from the concurrent CALIPSO data. In region B, on the other hand, the ρ value is larger than 0.3, but the classification result suggests clear area. More close examination of the image movie, however, indicates that the actually, the region is covered with very thin cirrus. Thus, this discrepancy is ascribable to the difficulty in detecting very thin cirrus clouds by using the SWA based on BT and BTD. Figure 6(c) shows the THI representation of the cloud phases. The dominance of ice particles (white color) is consistent with the presence of high-level clouds such as dense cirrus or cumulonimbus. Similar validation using the concurrent CALIPSO data has led to the conclusion that the classification based on both SWA is acceptable.

The altitude verification is performed at Chiba, indicated with a blue circle in Figure 5(a). The value of surface temperature, T_S, is derived from the ground observation at Chiba University, while that of cloud top, T_CT, from the BT of the corresponding pixel of Himawari-8. The results of altitude verification are summarized at the bottom of Table 4. For dataset #1, the estimated cloud top height, h, is 0 km, which agrees with the classification result of clear (no cloud) area. For dataset #2, the resulting value of h = 7.1 km is reasonable for the cloud type of dense cirrus. Similarly, good agreements have been found for other datasets except dataset #4. The resulting value of h = 2.8 km is inconsistent with the classification result thin cirrus, a type of high clouds. Thin cirrus clouds exhibit the value of T_CT that is close to the surface temperature, T_S, because of their semi-transparent features.

The comparation of SWA13-15 and SWA15-16 has indicated that generally the first combination can exhibit more sensitivity for different cloud types than the latter, as manifested in the wider spread in the scatter plot diagrams (Figure 4). However the latter combination of SWA tends to discriminate rain cloud regions among cumulonimbus and dense cirrus clouds. Nevertheless, the results from both SWAs are not contradictory to each other, as verificated by using the CALIPSO depolarisation ratio data.

4. Conclusion

This paper has described the effectiveness of SWA for detecting the cloud types using the IR bands of Himawari-8 satellite. Verification has been made with the space-borne lidar (CALIPSO) data and lapse rate consideration. Through the implementation of SWA, the following nine types of clouds have been successfully distinguished around the Japan area: high cumulonimbus, middle cumulonimbus, cumulus, dense cirrus, ice cloud, water cloud, thick cirrus, cirrus, and thin cirrus as summarized in the form of a SWA matrix. The classification results have been compared for the two cases of SWA, namely, SWA13-15 and SWA15-16. Scatter plot diagrams between BT and BTD are examined for both SWA, showing that generally more cloud types are detected with SWA13-15 (10.4 vs. 12.4 µm) as compared with SWA15-16 (12.4 vs. 13.3 µm). This is ascribed to the fact that the SWA13-15 diagram tends to spread more widely than SWA15-16 diagram in both winter and summer seasons. The latter SWA, however, is found to be more effective in detecting the region of cumulonimbus (i.e. likely rainfall region) amidst that of dense cirrus. In accordance with the variation in climate conditions, different values of BT and BTD thresholds have been determined from the close examination of cloud types. It is difficult to detect thin cirrus and very thin cirrus clouds using BT and BTD because of their semi-transparency. The coupled information on BT and surface temperature, however, is useful for solving this problem, leading to better discrimination between thin cirrus and clear regions. Combining the ancillary data from either space-borne or ground-based lidar will also be useful in this context. Finally, it is emphasized that the use of high frequency data of Himawari-8 is available every 10 min at full disk and 2.5 min around Japan, as an image movie makes it possible to detect detailed movement in cloud systems. Snow-covered areas, which tend to be classified as water cloud or cirrus cloud, can easily be delineated in such a movie type examination.

Acknowledgements

One of the authors (B.P.) expresses his gratitude to the financial supports, Ministry of Research, Technology and Higher Education (RISTEKDIKTI), Republic of Indonesia through the Scholarship Program for Research and Innovation in Science and Technology (RISET-Pro).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Purbantoro, B., Aminuddin, J., Manago, N., Toyoshima, K., Lagrosas, N., Sumantyo, J.T.S. and Kuze, H. (2018) Comparison of Cloud Type Classification with Split Window Algorithm Based on Different Infrared Band Combinations of Himawari-8 Satellite. Advances in Remote Sensing, 7, 218-234. https://doi.org/10.4236/ars.2018.73015

References