Next Article in Journal
Deriving Water Quality Parameters Using Sentinel-2 Imagery: A Case Study in the Sado Estuary, Portugal
Previous Article in Journal
The Horizontal Distribution of Branch Biomass in European Beech: A Model Based on Measurements and TLS Based Proxies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced Dual-Satellite Method for Detection of Low Stratus and Fog near Japan at Dawn from FY-4A and Himawari-8

1
Department of Atmospheric Science and Engineering, Ewha Womans University, Seoul 120-750, Korea
2
Department of Science Education, Ewha Womans University, Seoul 120-750, Korea
*
Author to whom correspondence should be addressed.
Submission received: 10 February 2021 / Revised: 5 March 2021 / Accepted: 6 March 2021 / Published: 9 March 2021
(This article belongs to the Section Atmospheric Remote Sensing)

Abstract

:
The detection of low stratus and fog (LSF) at dawn remains limited because of their optical features and weak solar radiation. LSF could be better identified by simultaneous observations of two geostationary satellites from different viewing angles. The present study developed an advanced dual-satellite method (DSM) using FY-4A and Himawari-8 for LSF detection at dawn in terms of probability indices. Optimal thresholds for identifying the LSF from the spectral tests in DSM were determined by the comparison with ground observations of fog and clear sky in/around Japan between April to November of 2018. Then the validation of these thresholds was carried out for the same months of 2019. The DSM essentially used two traditional single-satellite tests for daytime such as the 0.65-μm reflectance (R0.65), and the brightness temperature difference between 3.7 μm and 11 μm (BTD3.7-11); in addition to four more tests such as Himawari-8 R0.65 and BTD13.5-8.5, the dual-satellite stereoscopic difference in BTD3.7-11 (ΔBTD3.7-11), and that in the Normalized Difference Snow Index (ΔNDSI). The four were found to show very high skill scores (POD: 0.82 ± 0.04; FAR, 0.10 ± 0.04). The radiative transfer simulation supported optical characteristics of LSF in observations. The LSF probability indices (average POD: 0.83, FAR: 0.10) were constructed by a statistical combination of the four to derive the five-class probability values of LSF occurrence in a grid. The indices provided more details and useful results in LSF spatial distribution, compared to the single satellite observations (i.e., R0.65 and/or BTD3.7-11) of either LSF or no LSF. The present DSM could apply for remote sensing of environmental phenomena if the stereoscopic viewing angle between two satellites is appropriate.

Graphical Abstract

1. Introduction

Reliable monitoring of low stratus and fog (LSF) is of vital importance because restricted visibility can cause hazards in transportation and navigation including LSF-related accidents (e.g., [1,2]). Thick LSF at dawn is a major hazard during the morning rush hours. The persistence and frequency of LSF on land can be associated with various factors: air pollution and quality [3,4,5] and climate variabilities such as Arctic Oscillation [6] and El Niño/Southern Oscillation [7,8].
According to Yoo et al. [9], 71–76% of fog phenomena in surface synoptic observations (SYNOP) tend to accompany the low stratus over the Korean Peninsula at summertime dawn. This indicates that the better the LSF detection, the higher the possibility for fog detection. Low stratus and fog are often decoupled from the upper part of fog at night. Thus, low stratus and fog are examined in terms of LSF without individual distinctions in this study. LSF is also an appropriate term in satellite remote sensing since the distinction between low stratus and fog is fundamentally difficult without additional ground information such as surface elevation and cloud base height, etc. [2,3,10,11,12]. The only difference between the two weather phenomena is whether their bases are in contact with the ground [13].
Since LSF has various scales in space and time [6,14], its detection from both ground observations and numerical forecast models has intrinsic weaknesses. The observations cannot reflect spatial distributions of LSFs because of a limited number of stations [15]. The models cannot simulate LSFs as fast as needed due to their long spin-up time [16], and with inaccuracy due to coarse spatiotemporal resolutions (e.g., [2,17]) and limited parameterizations of microphysical and turbulent mixing processes (e.g., [1,18]). The most important measure for continuous monitoring of LSF over a wide area is geostationary satellites (GEOs). However, it is also regarded that the sensing accuracy in LSF at dawn is in need of improvement.
Traditional detection methods for LSF from a single satellite utilize thresholds that come from either various satellite channels or their spectral differences [19]. During the nighttime, Eyre et al. [20] developed a fog detection method by the difference in brightness temperature at 3.7 μm and 11 μm (BTD3.7-11) from the Advanced Very High-Resolution Radiometer (AVHRR), using the principle that 3.7 μm has a lower emissivity for water droplets than 11 μm. The BTD3.7-11 threshold test for LSF from GEOs has been continuously used [21,22,23,24,25,26,27].
During the daytime, BTD3.7-11 is less useful because the brightness temperature at ~3.7 μm (BT3.7) includes both the solar and earth radiation. Thus, the visible 0.68 μm reflectance (R0.68) test has been attempted instead (e.g., [28]), based on that a water droplet has a higher reflectivity than the clear sky. Kim et al. [29] have carried out the daytime marine fog detection based on the decision tree technique from the Geostationary Ocean Color Imager (GOCI) and Himawari-8 Advanced Himawari Imager (AHI); although the hit rate (HR) is 0.66–1.00, there is room for improvement due to the high false alarm ratio (FAR) and the limited amount of data at three Korean islands (FAR = 0.31–0.33). The visible technique has obvious drawbacks for distinguishing LSF from mid/high-level clouds. Thereby, researchers have tried to overcome this hindrance using both visible and infrared thresholds [14,28,30].
Despite many studies on LSF remote sensing [31,32], the detection accuracy remains low, particularly at dawn [5,33]. The BTD3.7-11 value at dawn is nearly zero due to rapidly changing optical features with different signs (i.e., negative to positive) [34]. Like snow, the R0.67 value at dawn also increases in inaccuracy depending on SZA or surface conditions [33,35]. To address this issue, Yoo et al. [9] and Yang et al. [36] proposed the dual-satellite method (DSM) that utilizes the observed stereoscopic differences (Δ) in R0.67 and BTD3.7-11, between two GEOs (i.e., COMS and FY-2D), defining ΔR0.67 and ΔBTD3.7-11 in this study. Both studies applied DSM to LSF detection at dawn near the Korean Peninsula. While the ΔR0.67 test for LSF detection was useful in the summertime [9], the advantage of the ΔBTD3.7-11 test for the detection was emphasized in the springtime [36].
An example of the basic principle of DSM is that the detection of an object is better by two eyes (or satellites) than a single eye [37]. Since DSM for LSF detection is based on the bidirectional reflectance function (BRDF) of the two GEOs, it is most useful in summertime dawn when the stereoscopic differences are maximized under the Viewing Zenith Angle difference (ΔVZA) of ~46.5° across the Korean Peninsula [9]. However, the accuracy for LSF detection highly depended on ΔVZA, season, and satellite-observed channels. Therefore, new tests are required for the LSF detection improvement, utilizing the additional channels and high spatiotemporal resolution of the recently advanced imagers of GEOs. The purpose of this study was to improve the detection accuracy and spatial information of LSF at dawn near Japan in terms of probability index (PI) by devising additional tests from the advanced satellites, Feng-Yun-4A (FY-4A; [38,39]) and Himawari-8 [40] to the previous tests [9,36]. Compared to previous GEOs, these two GEOs have many channels in high spatiotemporal resolution, and they can provide additional independent information. The two GEOs have ΔVZA of ~40.4° in the region of interest near the nadir of Himawari-8 (Figure 1). The selection of the region is reasonable to enhance the signal-to-noise ratio in DSM. Fog in Japan occurs more frequently during the warm season than the cold season, with high regional variability [41].
The Radiative Transfer Model (RTM) simulations of this study which were arranged and stored in the Look-Up Table (LUT) in advance, were comprehensively used for an approximate diagnostic check of the LSF optical properties. In Section 2, the data and method were described. Here spectral properties of newly developed tests were explained. The results were shown in terms of the LSF and clear-sky classification using optimal thresholds with radiative transfer simulations in Section 3. The case study for validation was shown with the PI formulation. Discussion and conclusions were presented in Section 4 and Section 5, respectively.

2. Materials and Methods

To derive and validate the threshold values for LSF detection at dawn near Japan, the following three datasets were used: geostationary satellite observations from Himawari-8 and FY-4A, ground observations from Meteorological Aerodrome Report (METAR) for the weather phenomena of fog and clear sky, and the LUT from RTM. For this study, “dawn” was defined as within 2 h after sunrise (67° < SZA < 86°). The detailed information for the datasets is described below. For SZA, the average of LSF cases in 2018 was ~79° and that of clear-sky conditions was ~78°. In 2019, both SZA values were ~80°. Also, the meanings of acronyms used in this study are noted in Table A1.

2.1. Satellite Observations

The geostationary satellites of Himawari-8 and FY-4A are located at different longitudes over the equator, and they show a VZA difference of 40.4˚ near the nadir of Himawari-8 (Figure 1). The satellite information and application in this study with respect to the channel are described in Table 1. The primary payloads of FY–4A and Himawari-8 are the Advanced Geostationary Radiation Imager (AGRI; [38,39]) and AHI; [40], respectively. The instruments provided significant advancements in terms of the number of bands, spatial resolution, and temporal frequency. Their solar and thermal channels are separately presented in the table. There were two visible (VIS) channels of solar radiation centered at ~0.65 µm and ~1.6 µm. In the 0.65 µm channel, clouds appeared bright due to their high reflectance, while the signals from the land/sea surface were generally weak. The 1.6 µm channel was in the infrared (IR) region and allowed distinction between snow/ice-covered land and signals from other sources. This channel was sensitive to the cloud phase, and water clouds reflected more energy in this spectral range than the ice clouds.
Of the six channels used, the four non-solar channels were classified as thermal parts. The 3.7 µm lies in the spectrum of two outgoing energy sources (e.g., solar reflectance and thermally emitted radiation). This channel is not only sensitive to the cloud phase but also useful for the detection of water clouds at night. In the region around 8.5 µm, there is a little water vapor absorption, but it is still in the window region. The IR window channel around 11 µm is affected primarily by radiation from the Earth’s surface or cloud tops. The 13.3–13.5 μm band is used as cloud masks and cloud-top height.
Level 1B data as digital count values were used to convert the VIS and near-infrared (NIR) data to visible reflectance and the others (mid-IR, longwave-IR1, longwave-IR2, longwave-IR3) to brightness temperature. In order to observe the two satellite observations on the same-sized grid, collocation was performed by interpolating the Himawari-8 observations to FY-4A along with the ground observations shown in Table 2 and Figure 2. The analysis of this study is not so sensitive to the grid-size of either 4 km by 4 km or the original pixel (2 km by 2 km), but the data smoothing effect tends to increase in the former grid. To compare the dual-satellite observations by pixel, the spatial resolution at the channels noted above was set to 4 km by 4 km. Thus, the Himawari-8 data in a two-dimensional array of longitude and latitude near Japan (122.5-144E, 23.5-44N) were interpolated into the FY-4A array. This method keeps the FY-4A accuracy without degradation. The observations for seven variables (BTD13.5-8.5, R0.65, BTD3.7-11, NDSI, ΔR0.65, ΔNDSI, and ΔBTD3.7-11) were analyzed for LSF detection. Additionally, the difference value between Land Surface Temperature (LST) and 11 μm brightness temperature (BT11) was used to analyze cloud-top height as an acronym of LST-BT11. In this study, BTD13.5-8.5 was the brightness temperature difference between 13.5 μm and 8.5 μm. Based on the main weighting function peak in the Satellite Infrared Spectrometer (SIRS-B) channel at 750 cm−1 [42], the peak of AHI at 13.28 μm (753 cm−1) was located near the 850–900 hPa level. Meanwhile, the 8.5 μm AHI surface channel was used for the analysis of the cloud phase. Reflectance values (R0.6 and R1.6) were normalized by dividing by cos ( SZA ) and set to unity when the values were greater than 1 due to normalization (e.g., large SZA ≥ ~80°). The difference in values between the two collocated satellites is defined with Equations (1)–(3) as follows:
Δ R 0.65   =   R 0.65 FY - 4 A   R 0.65 Himawari text 8
Δ BTD 3.7 - 11   =   BTD 3.7 11 FY - 4 A   BTD 3.7 11 Himawari - 8
Δ NDSI   =   NDSI FY - 4 A   NDSI Himawari - 8
NDSI [43] was calculated from the visible reflectance values of 0.65 μm and 1.6 μm as shown in Equation (4).
NDSI   =   R 0.65     R 1.6 R 0.65   +   R 1.6
Four variables (BTD13.5-8.5, R0.65, BTD3.7-11, and NDSI) were obtained from a single satellite’s (i.e., Himawari-8 only) observations, while three variables (ΔR0.65, ΔNDSI, and ΔBTD3.7-11) were derived from both satellites. Also, LST data were given from ground observations.
The NDSI index has been widely used to distinguish between snow and clouds since snow has a lower reflectance at 1.6 μm [43]. The studies [14,44] presented the NDSI values for various surface conditions (e.g., water clouds, ice clouds, and vegetation,) to eliminate snow pixels, which have relatively large R0.65 values. Wu and Li [45] determined both the NDSI and R0.67 thresholds for daytime sea fog detection using MODIS data. In their study, the NDSI was used for LSF detection over inland China to divide vegetation and clouds during the warm season, because vegetation has larger R1.6 and smaller R0.65 than clouds. Ryu and Hong [46] proposed a regression-based method for sea fog detection from a combination of NDSI and reflectance in the AHI green band (0.51 μm).
The BTD13.5-8.5 value which is the difference in brightness temperature between 13.5 and 8.5 μm was used to detect higher ice clouds [47,48], because the 13.5 μm channel belongs to a CO2 absorption band that affects atmospheric transmission for the calculation of cloud altitudes. Thus, this value is sensitive to atmospheric altitude. This result indicates that the BTD13.5-8.5 value can be utilized as an indicator to classify LSF among the weather phenomena of clear sky, LSF, and high clouds.

2.2. Ground Observations

In this study, ground observation data from METAR were regarded as the ground truth for LSF and clear-sky cases from March to November in 2018 and 2019. The METAR is a weather status report designed for aviation based on information from an Automatic Weather Station (AWS), such as visibility, cloud amount, cloud base height present weather, ambient temperature, and relative humidity at hourly intervals. Figure 2 shows the locations of METAR stations in Japan including the “Ulleung” island station in South Korea. The information for a total of 31 stations over the study region is also described in Table 2. The LSF occurrences were selected only in cases where the current weather conditions were reported as “FG” (i.e., fog). In METAR data, fog appears in various forms depending on accompanying weather phenomena and full or partial spatial covering. Only cases of “FG” that covered the entire observation area without precipitation were selected in this study. The clear-sky cases as a contrast group of LSF cases were adopted under conditions with the abbreviation “CAVOK” or cloud amount with “FEW” (1 to 2 oktas). CAVOK stands for Ceiling (or clouds) And Visibility OK, indicating no clouds below 1500 m or the highest minimum sector altitude and no cumulonimbus or towering cumulus at any level, with a visibility of 10 km or more and no significant weather change [49]. The total observation numbers for LSF and clear-sky occurrences used in this study were 347 and 229 during 2018 and 2019, respectively.
Weather maps (i.e., SYNOP chart) provided by the Korea Meteorological Administration (KMA) were used to verify the presence of LSF and upper cloud information detected by satellites and to obtain extensive information on air pressure patterns and wet areas. However, since the weather maps were produced at three-hour intervals, a direct comparative analysis was not possible due to the ordinary rapid generation and dissipation of fog. Therefore, weather map analysis was performed only when the map time coincided with the ground or satellite observation time. If there was a time difference between fog occurrence time and weather map, the BT11 distribution, which is approximately related to cloud top temperature, was used instead.

2.3. Radiative Transfer Model

Optical characteristics of LSF can be estimated using the RTM simulation with satellite observations (e.g., [50]). In this study, the Santo Barbara DISORT (SBDART) of a plane-parallel radiative transfer model was used [51] and the LUT was calculated from RTM in advance for estimation. The cloud model in the SBDART yielded a radiance value relating to cloud optical properties (e.g., Fog Optical Thickness, FOT; Cloud Optical Thickness, COT; Effective Radius, ER) based on the Mie theory. This allowed us to understand satellite-observed spectral values in the LSF environment and their optical difference between the two satellites. In addition, the simulation provided insight into the effects of higher clouds over the LSF layer.
The Spectral Response Function (SRF) difference between the two satellites was less significant than that of the VZA [9]. Therefore, the values at the central wavelength of each channel (https://www.data.jma.go.jp accessed on 3 June 2019 for Himawari-8; http://satellite.nsmc.org.cn accessed on 5 June 2019 for FY-4A) were used instead of the SRF integrated over its entire wavelength range, of which the calculation process was time-consuming. Table A2 shows the input variables used for the RTM run in this study with their acronyms. The LUT was set up for five conditions similar to those in the study by [36] (their Figure A1). LSF1 at 0–1 km and 0–2 km without higher clouds above the LSF layer shown in Cases A and B. Cases C to E were with LSF2 accompanying the clouds. Case C and D were the clouds at 4–6 km above the 2 km fog layer, composed of either water particles or ice, respectively. Also, in Case E, clouds of ice particles exist 8–10 km above LSF1. We obtained simulated values from the LUT based on the input data of three sun or satellite orbit angles (SZA, VZA, and RAA) for each LSF case.

2.4. Formulation of Probability Index (PI)

Figure 3 shows the flow diagram for LSF detection near Japan at dawn during the period from March to November of 2018–2019, based on the near-simultaneous satellite observations of Himawari-8 and FY-4A. After collocating 31 ground station sites onto ~4 km × 4 km satellite observation grids, satellite data were collected when fog occurred at the ground stations. The optimal thresholds of eight satellite-derived variables were derived from iteration processes to discriminate between LSF and clear-sky cases based on a skill-score test. The variables were BTD13.5-8.5, ΔNDSI, R0.65, ΔBTD3.7-11, LST-BT11, ΔR0.65, BTD3.7-11, and NDSI. Using the top four skill-score variables, a total of 16 cases could be extracted from their statistical combinations (Table 3 and Table 4). In addition, the weights were given and merged into a total of five classes for the probability of LSF occurrence for the 16 cases. To calculate PI in real time, calculation of five PI classes and the AP (i.e., assigned probability) value in each class were conducted using the long-term data of the other training period (i.e., the 2018 control period in this current study). The introduction of the AP values is necessary to quantify the LSF occurrence by its probability. The magnitude order of AP was set based on the top four “POD minus FAR” values in the skill score test of the pre-processing. “POD minus FAR” could comprehensively explain both POD and FAR. This “PI formulation” is explained in detail in Section 3.3.
Utilizing these climatological values during the training period, a threshold test of the real-time dual satellite observations was carried out and provided the PI class (or no fog) in a spatial grid. The AP as a function of PI class (i.e., Classobs) was equivalent to the PIobs value, indicating the LSF probability from the real-time dual satellite observations. Fog detection methods from simultaneous observations of multiple satellites have been developed in two previous studies [9,36]. We analyzed various channel information obtained from advanced satellites and tried to present the probability index of fog occurrence more precisely in this study compared to previous studies. Table 5 shows a summary of the DSM differences between previous and the present studies.
The maximum number of PI class (i.e., 16) will be available based on the probability of the four variables if the long-term dataset of LSF and clear-sky occurrences is given in the future (Table 4).

3. Results

The processes for deriving the optimal thresholds and their verification will be discussed in Section 3.1 and Section 3.2, based on two-year observations and the RTM simulation. These sections explain the selection of optimal parameters for fog detection. The process of developing PI and its application to case studies to address the PI merits in LSF detection are presented in Section 3.3.

3.1. OptimalLSF Thresholds from the 2018 Control Data

To derive optimal LSF thresholds, the average values with standard deviation (±1σ) of satellite observations were investigated for LSF and clear-sky cases at dawn during the control (or training) period of 2018 (Figure 4). The seven satellite-observed values (BTD13.5-8.5, R0.65, ΔNDSI, ΔBTD3.7-11, NDSI, ΔR0.65, and BTD3.7-11) and the values of LST-BT11 with ground observations of land surface temperature were analyzed to determine the range of each variable in LSF and clear-sky conditions. With the separation of LSF from the clear sky, the LSF threshold can be relatively reliable for LSF detection. Separation without overlap of the weather phenomena was noted for R0.65, which was traditionally utilized for fog detection. In addition to the two variables (ΔNDSI and ΔBTD3.7-11) derived from the dual-satellites of Himawari-8 and FY-4A, LSF and clear sky were distinct in the BTD13.5-8.5 available from Himawari-8 only. This result suggests that the threshold values derived from dual satellites (ΔNDSI and ΔBTD3.7-11) can be worthy for LSF detection in addition to single satellite observations (BTD13.5-8.5 and R0.65). However, other thresholds (LST-BT11, BTD3.7-11, ΔR0.65, and NDSI) have uncertainties for LSF detection due to mixed information that cannot discriminate between LSF and clear sky.
Using the data for the control period of 2018, the optimal thresholds were derived through the iteration process. The lower boundary value (i.e., Thlower in Table 6) that distinguished LSF and clear sky was obtained through iteration, while the upper boundary value (i.e., Thupper) was adjusted to the +1 σ range of the average value to avoid middle or high cloud contamination. The results of statistical verification for each threshold variable are shown in the contingency table and the index definition (Table A3), based on the five skill score indices (HSS, CSI, POD, PC, and FAR) shown in Table 3. HSS can be an overall skill indicator of detection accuracy [52], although the difference between POD and FAR (i.e., POD minus FAR) was analyzed in this study as an important indicator of detection accuracy. The difference values were higher in magnitude order for BTD13.5-8.5, ΔNDSI, R0.65, and ΔBTD3.7-11. In these four variables, the skill scores of POD minus FAR (0.79–0.82) were roughly the same, within a 3% difference. Meanwhile, the POD minus FAR values for the other four variables did not exceed 0.47. Therefore, the threshold values of the top four variables were more useful for LSF detection than were those of the other four thresholds.
Figure 5 shows scatter plots of the top four variables (BTD13.5-8.5, ΔBTD3.7-11, ΔNDSI, and R0.65) with respect to LST-BT11 for the weather phenomena of LSF (red circle) and clear sky (blue triangle) near Japan at dawn during 2018. The value of LST-BT11 tended to be higher under middle/high clouds than in LSF/clear-sky conditions. The weather phenomena of LSF and clear sky were categorized well using the top four variables, with high verification scores of POD minus FAR.
To theoretically support these observational results, RTM simulation was performed considering various optical characteristics (e.g., altitude and optical thickness of the fog layer, altitude and optical thickness of upper clouds, effective radii of fog/cloud particles) in Table A2, and the results were shown in Figure 6. The figure shows the simulation results in the domains of a) BTD13.5-8.5 vs. LST-BT11, b) ΔBTD3.7-11 vs. LST-BT11, c) ΔNDSI vs. LST-BT11, and d) R0.65 vs. LST-BT11. The average values of the Himawari-8 angles (i.e., SZA = 79.5°, RAA = 74.5°, and VZA = 41.5°) for 176 LSF cases during 2018 were used as RTM input. The satellite-observed mean values are included in the figure for LSF (black asterisk) and clear sky (gray-cross). The clear-sky values from the observations and simulations (pink asterisk) are depicted using black circles to show the approximate agreement between them. Five colors (yellow, blue, red, navy, and green) denote the heights (FH or CH: 1, 2, 4–6, 8–10 km) of the fog and upper cloud layers with respect to effective radii (ER: 2, 4, 8, 16 μm) and water/ice phase. The RTM details are given in Table A2.
The DSM variables (i.e., ΔNDSI and ΔBTD3.7-11), as well as BTD13.5-8.5 and R0.65 from the single satellite observations (i.e., Himawari-8 only), had a significant tendency to separate LSF and clear sky at dawn (Figure 6a–d) due to the unique satellite orbit in space (i.e., VZA difference between the two satellites). The distributions of observed mean values for the two meteorological phenomena (LSF and clear sky) were consistent within the ordinate range of the simulation. The difference in “LST-BT11” between observations (5.49 K) and simulations (0.82 K) for clear-sky conditions could be explained in terms of the effect of water vapor and optically thin clouds during the warm season [53]. Considering that the LST-BT11 value is 19.7 K for LSF and 5.5 K for clear sky, the average LSF top height is 2–2.5 km, consistent with that of [11].
However, the simulated values of ΔBTD3.7-11 on the abscissa were somewhat out of range with the average of its observations (Figure 6b). According to the simulation, the upper clouds above the fog layer could cause errors in the LSF detection of satellite observations. The upper clouds generally show relatively large values of LST-BT11, compared to the lower clouds. Overall, the simulations were in good agreement with fog detection in DSM.

3.2. Verification of Satellite-Observed Thresholds during the Experimental Period of 2019

The eight LSF threshold values derived from the control period of 2018 for LSF detection at dawn were applied to the experimental data during 2019. Figure 7 shows statistical verification for the threshold values of eight satellite-observed variables in terms of a) POD minus FAR, b) POD, c) FAR, and d) CSI. The verification was performed based on the total ground-observations of 171 fog and 110 clear-sky occurrences at dawn near Japan. The thresholds of the top four variables (BTD13.5-8.5, R0.65, ΔNDSI, and ΔBTD3.7-11) in the verification skill-scores were substantially better than those of the other four variables. The “POD minus FAR” scores of the top four variables during the experimental period were lower by ~5% than those during the training period (Table 3). In particular, the “ΔNDSI” threshold was most sensitive to the interannual dependence of threshold accuracy. However, more long-term data seemed to be required to estimate interannual variability in view of the two years of data in this study.
Figure 8 presents the “±1σ” range of LSF and clear-sky cases for each variable in the same way as Figure 4 so that the biennial changes in satellite-observed variables under the weather phenomena were obtained from a comparison of the results shown in the two Figures. Since the threshold values were optimized for the training period data, statistical verification results decreased slightly (Figure 7 and Table 3). However, there were no significant differences in skill scores for the top four variables (BTD13.5-8.5, ΔBTD3.7-11, ΔNDSI, and R0.65) between the two years. In addition, the interannual change of NDSI among the entire eight variables was noticeable. Although the skill scores of NDSI itself were not good, ΔNDSI showed excellent results due to the DSM advantage. This result was due to not only the large VZA difference between the two satellites but also to the relatively long optical path of FY-4A, which was viewed from outside of the nadir. Therefore, the DSM was likely to cancel the noise error and amplify the signal for LSF detection. This tendency was similar to that of the ΔBTD3.7-11 case in a previous study [36].

3.3. PI development and Case Study

The probability index of LSF was developed to present fog occurrence as a probability concept rather than an alternative one (either no fog or fog), using the top four variables that showed excellent results (POD ~0.8; POD minus FAR ~0.7) during the experimental period. In Table 4, five LSF classes were derived from 16 probability sets of the variables based on the threshold test for detection of either LSF or clear-sky phenomenon. For instance, Case 1, which passed the test for LSF detection by all four variables, corresponded to Class 1 in terms of the set “P( A B C D ) ” with an assigned probability of 1.0. For Class 1, in which all four variables passed their own thresholds, the probability of LSF occurrence was classified as “very high”. In the same way, three threshold values among the four variables were achieved for those in Class 2. In addition, any combination of the three variables was classified as a “high” possibility of LSF occurrence. In this regard, the number of cases that passed the threshold test for the four satellite-observed variables (i.e., BTD13.5-8.5 in gray rectangle, ΔNDSI in orange, ΔBTD3.7-11 in blue, and R0.65 in green) was shown in the detection of a) LSF and b) clear sky in the so-called Venn diagram shown in Figure 9. In other words, the zones of red (Class 1 for the occurrence possibility of the two weather phenomena), yellow (Class 2), blue (Class 3), and white (Class 4) were defined and were also explained in the Table.
When the LSF class in Table 4 was determined for a total of 171 LSF cases during the experimental period of 2019, 99 passed all four thresholds (i.e., Class 1). In addition, 43 cases passed three thresholds (i.e., Class 2). These results are described in detail in Table 6. Satellite-derived PI distributions of the five classes in the Table were weighted at 0.25 probability intervals for detection of LSF and clear-sky cases in 2019, and detection was quantitatively presented in five steps rather than two (i.e., yes or no). The probability values were also shown in terms of “assigned probability” in Table 4. If the PI results were evaluated for LSF cases during 2019, the 171 fog occurrences represented a detection score of 140.25. In clear-sky conditions, the score was 93.75 of 110 (Table 6). This indicated excellent detection of realistic PI values in the spatial distribution, rather than either “fog” or “no fog” from the traditional method. In other words, when LSF detection was performed using a single threshold, the existence of LSF was generally classified with a probability of either 0 or 1. Meanwhile, the PI in this study showed five classified probability values based on the four satellite-observed variables in DSM, for useful and reliable LSF detection.
The PI in the case study was applied to two foggy scenes at dawn near Japan (128–145 E, 30–45 N) to compare the results with those of the conventional method (i.e., using a single threshold value of either R0.65 or BTD3.7-11) in the LSF spatial distribution. Figure 10 shows the fog probability distribution at dawn (06:00 LST) on July 9, 2019, near Japan using the detection methods of a) R0.65 threshold, b) BTD3.7-11 threshold, and c) LSF PI in this study. The BT11 distribution for approximate cloud-top temperature and the SYNOP map for station-observed fog occurrences were shown in Figure 10d,e, respectively. Fog occurrences (pink triangles in Figure 10a–d) were reported (station numbers 20, 26, and 31 in Table 2) at three stations by ASOS (automated surface observing system).
The R0.65 threshold successfully detected two of the three fog occurrences at a detection rate of 0.66 (Figure 10a); the rate was 0.33 with the BTD3.7-11 threshold (Figure 10b). The PI produced by combining the four variables using dual-satellite observations showed a probability of 0.66, higher than that in the BTD3.7-11 case but similar to that in the R0.65 case. Although PI accuracy was the same as that of the R0.65, it provided LSF spatial distribution in terms of probability in five classes. To examine the LSF-related information over the study area, the BT11 distribution and the weather map chart were analyzed simultaneously (Figure 10d,e). The BT11 was useful to estimate the approximate cloud top temperature. In addition, the map allowed the identification of the distributions of foggy and humid areas and lower or upper clouds. The BT11 values were either high (≥ ~290 K) or low (≤ ~230 K) in areas where PI indicated no or low probability of fog due to almost clear-sky or optically thin cloud conditions. In Figure 10c, the LSF PI values indicated that the LSF occurrences were “very high” or “high” in areas where fog actually occurred on the weather map (pink triangle). For areas at station numbers 20 and 31 where fog detection failed using R0.65 or BTD3.7-11, PI indicated a relatively low probability of fog but was the only method to generally succeed in the detection (PI = 0.75).
Figure 11 is the same as Figure 10 except for 05:30 LST on 4 July 2019, and without a SYNOP map (not available at the time). According to the ASOS report, there were three fog occurrences at the stations (in Table 2; 6, 16, and 27). The time difference between the fog occurrence on the report and that on the weather map was greater than 30 minutes, so the low cloud and fog areas were analyzed with the distribution of BT11 without the map. The PI and R0.65 successfully detected all three fog occurrences, showing a detection probability of one. However, the BTD3.7-11 threshold detected only one of the three, a probability of 0.33, only successful at station 6. These results show that PI was able to detect fog more accurately than BTD3.7-11 and was more realistic (or natural) and excellent than R0.65 in the LSF spatial distribution. In particular, detailed information on the LSF PI in five classes was useful for aviation and navigation. So far, it has been difficult to identify the occurrence of sea fog in the open ocean despite its frequent occurrence due to the limited number of ground observations. The spatial PI distribution of LSF over a large area could be useful to monitor fog phenomena over an extensive region including the ocean in almost real-time.

4. Discussion

The DSM used in this study is applicable to LSF (or fog) detection based on advanced satellite observations, including GEO-KOMPSAT-2A (GK2A; [54]) in the future, for real-time monitoring. Further improvement in nowcasting and forecasting (e.g., advective fog) can be expected. The Advanced Meteorological Imager (AMI) data of GK2A will be useful for the LSF analysis over the Korean Peninsula because of its nadir point. The DSM from FY-4A and Himawari-8 were applied to the LSF detection near Japan at dawn due to the Himawari-8 nadir. Similarly, the method can be utilized for the detection over the inland region of China under the FY-4A nadir at dusk, although LSF generally occurs at dusk less than at dawn due to daytime solar heating.
Various applications of the DSM in the meteorological and environmental fields are possible to observe fog targets in a three-dimensional structure from a stereoscopic view, similar to observations with two eyes [37]. Thus, this method can be applied to any worldwide region where simultaneous observations by two geostationary-orbit satellites are available without additional expense. However, the difference in VZA between the two satellites is required to be between 40° and 50° to achieve optimal LSF detection. For this reason, the detection accuracy of ΔR0.65, which had been excellent in the COMS and FY-2D (VZA difference ~46.5°) [9], was significantly reduced in this study. It seems that ΔR0.65 is more sensitive to both VZA and seasonality than to the top-four variables in this study. Here the variables of ΔNDSI and ΔBTD3.7-11 were newly utilized for LSF detection in terms of PI, despite the decrease in VZA, with the BTD13.5-8.5 from the spectral difference of Himawari-8 only.
The RTM simulation of this study may have some limitations in the diverse weather conditions of multilayer, broken, or mixed-phase clouds [9]. Also, errors in the LSF detection due to the higher (or upper) clouds above the LSF layer have been discussed in detail, based on the simulation [9,36]. However, the adverse issues may not be so critical for the LSF (under ~2.5km height) in relatively stable atmospheric conditions, considering the ratio (13–20%) of satellite-observed multi-layer cloudy pixels to all cloudy pixels [55] and the vertical properties of ice-over-water clouds [56]. In addition, for a comparison of the two GEO data in this study, the coarse-grid averaging of Himawari-8 observations from 2 km-by-2 km to 4 km-by-4 km may miss the detection of optically thin LSF phenomena in a small horizontal scale less than 4 km, particularly in the tests of BTD13.5-8.5 and R0.65.

5. Conclusions

In this study, the dual-satellite method was applied to advanced satellites (Himawari-8 and FY-4A), and newly developed thresholds with additional channels and LSF spatial distribution in terms of PI were provided to improve its detection in Japan at dawn. In addition, the optical characteristics of LSF were estimated using the LUT from the RTM simulation under various conditions and were consistent with satellite-derived observations. The PI of LSF in the DSM was derived by combining the major threshold components (i.e., BTD13.5-8.5, ΔNDSI, ΔBTD3.7-11, and R0.65) based on the top-four “POD minus FAR” scores (≥ 79%) among the eight variables.
The case study was carried out to address the PI merits of the LSF spatial distribution. The PI values in this study were similar to R0.65 simply based on the success rate of fog detection due to the use of R0.65 as one of the four equally-weighted variables for constructing the PI. Thus, the fog detection rate calculated from PI was expected to be similar to that from R0.65. However, the spatial distributions revealed the advantages of PI. Since multi-stage fog detection was possible in PI over a wide area, the LSF distribution in the DSM could be presented in more detail than that of the conventional methods of R0.65 and/or BTD3.7-11. More specifically, R0.65 only reported two classes of either presence or absence of fog, whereas PI presented fog occurrence probability as five classes. The classes can be increased to a maximum of 16 if a long-term dataset is available in the future, as shown in the statistical combination (Table 4). In this case, different probability weights would be assigned with respect to the variables. Thus, the LSF PI can be used as more practical and useful information for navigation and aviation. In particular, we believe that LSF PI is valuable for open oceanic areas where ground observations are limited in view of both accuracy and spatial distribution.

Author Contributions

Conceptualization, J.-H.Y., J.-M.Y.; methodology, J.-H.Y., J.-M.Y.; software, J.-H.Y.; validation, J.-H.Y. and J.-M.Y.; formal analysis, J.-H.Y. and J.-M.Y.; investigation, J.-H.Y. and J.-M.Y.; data curation, J.-H.Y.; writing—original draft preparation, J.-H.Y. and J.-M.Y.; writing—review and editing, J.-H.Y., J.-M.Y. and Y.-S.C.; visualization, J.-H.Y.; supervision, J.-M.Y.; project administration, J.-M.Y.; funding acquisition, J.-M.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. 2019R1F1A1059984).

Acknowledgments

We thank the employees at the CMA (for FY-4A) and the NMSC/KMA (for Himawari-8) for providing the satellite data. We are grateful to the anonymous reviewers for the constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. List of acronyms used in this study.
Table A1. List of acronyms used in this study.
AcronymsOriginal Words (or Details)AcronymsOriginal Words (or Details)
AGRIAdvanced Geostationary Radiation ImagerGK2AGEO-KOMPSAT-2A
AHIAdvanced Himawari ImagerGOCIGeostationary Ocean Color Imager
APassigned probabilityHRhit rate
ASOSautomated surface observing systemHSSPI estimated from satellite-observed scene
AVHRRAdvanced Very High-Resolution RadiometerIRinfrared
AWSautomatic weather stationKMAKorea Meteorological Administration
BRDFbidirectional reflectance functionLSFlow stratus and fog
BT11brightness temperature at ~11 µmLUTlook-up table
BT3.7brightness temperature at ~3.7 µmMETARMETeorological Aerodrome Report
BTD13.5-8.5difference between BT13.5 and BT8.5MODISModerate Resolution Imaging Spectroradiometer
BTD3.7-11difference between BT3.7 and BT11NDSInormalized difference snow index
CAVOKCeiling (or clouds) And Visibility OKPCpercentage correct
CHcloud heightPIprobability index
COMSKorean Communication, Ocean, and Meteorological SatellitePODprobability of detection
COTcloud optical thicknessR0.65reflectance at ~0.65 µm
CSIglobal telecommunications systemRAArelative azimuth angle
DSMdual satellite methodRTMradiative transfer model
EReffective radiusSIRS-BSatellite Infrared Spectrometer
FARfalse alarm ratioSRFspectral response function
FGfogSYNOPsurface synoptic observations
FHfog heightSZAsolar zenith angle
FOTfog optical thicknessTHLOWERlower threshold
FY-2DChinese FengYun-2DTHUPPERupper threshold
FY-4AChinese FengYun-4AVISvisible
GEOgeostationary-orbit satelliteVZAsatellite viewing zenith angle
Table A2. Input variables of SBDART for the LUT product.
Table A2. Input variables of SBDART for the LUT product.
Input variableContents
atmospheric profileMid-latitude summer, US62
wavelength (λ):three channels of VIS, SWIR, & IR1 for COMS & FY-2D0.55–0.90, 3.5–4.0, 10.3–11.3 μ m
solar zenith angle (SZA)0   SZA 80 ° at 10 ° intervals, and 85 °
surface typeOcean, Vegetation
fog height (FH)Water fog at 0–1 km or 0–2 km
upper cloud height (CH) above the fog layerWater/ice cloud (4–6 km), Ice cloud (8–10 km)
fog optical thickness (FOT) 0, 0.5, 1, 2, 4, 8, 16, 32, 64
cloud optical thickness (COT)0, 4, 8, 16, 32
effective radius of fog (FER)4, 8, 16, 32 μ m
effective radius of cloud (CER)2, 4, 8, 16 μ m
flux computation stream32
vertical resolution1 km
viewing zenith angle (VZA) 0   VZA 90 ° at 10 ° intervals
relative azimuth angle (RAA)0 RAA 180 ° at 30 ° intervals
boundary layer aerosol typeUrban
vertical optical depth of boundary layer aerosolsnominally at 0.55 μm0.2
Table A3. Contingency Table and definitions for the statistical skill test.
Table A3. Contingency Table and definitions for the statistical skill test.
SYNOP
LSFClear sky
Satellite observationLSFab
Clear skycd
POD = a a + c         CSI = a a + b + c         FAR = b a + b         HSS = 2 ( a d b c ) ( a + c ) ( c + d ) + ( a + b ) ( b + d )         PC = a + d a + b + c + d
POD-FAR = a a + c b a + b

References

  1. Kim, W.; Yum, S.S.; Kim, C.K. Numerical simulation of sea fog over the Yellow Sea: Comparison between UM + PAFOG and WRF + PAFOG coupled systems. Asia -Pac. J. Atmos. Sci. 2020, 56, 89–105. [Google Scholar] [CrossRef]
  2. Izett, J.G.; van de Wiel, B.J.H. Why does fog deepen? An analytical perspective. Atmosphere 2020, 11, 865. [Google Scholar] [CrossRef]
  3. Niu, S.; Lu, C.; Yu, H. Fog research in China: An overview. Adv. Atmos. Sci. 2010, 27, 639–662. [Google Scholar] [CrossRef]
  4. Zhang, G.; Bian, L.; Wang, J.; Yang, Y.; Yao, W.; Xu, W. Boundary layer structure of fog occurred in Beijing and its surrounding areas. Sci. China (D) 2005, 35 (Suppl. I), 73–83. (In Chinese) [Google Scholar]
  5. Wan, J.H.; Jiang, L.; Xiao, Y.F.; Sheng, H. Sea fog detection based on dynamic threshold algorithm at dawn and dusk time. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2019, XLII-3/W9, 159–163. [Google Scholar] [CrossRef] [Green Version]
  6. Liu, P.; Tang, M.; Yu, H.; Zhang, Y. Influence of arctic oscillation on frequency of wintertime fog days in eastern China. Atmosphere 2020, 11, 162. [Google Scholar] [CrossRef] [Green Version]
  7. Yu, H.; Li, T.; Liu, P. Influence of ENSO on frequency of wintertime fog days in eastern China. Clim. Dyn. 2019, 52, 5099–5113. [Google Scholar] [CrossRef]
  8. Hu, Z.Z.; Kumar, A.; Jha, B. How much of monthly mean precipitation variability over global land is associated with SST anomalies? Clim. Dyn. 2020, 54, 701–712. [Google Scholar] [CrossRef]
  9. Yoo, J.M.; Choo, G.H.; Lee, K.H.; Wu, D.L.; Yang, J.H.; Park, J.D.; Choi, Y.S.; Shin, D.B.; Jeong, J.H.; Yoo, J.M. Improved detection of low stratus and fog at dawn from dual geostationary (COMS and FY-2D) satellites. Remote Sen. Environ. 2018, 211, 292–306. [Google Scholar] [CrossRef]
  10. Egli, S.; Thies, B.; Bendix, J. A hybrid approach for fog retrieval based on a combination of satellite and ground truth data. Remote Sens. 2018, 10, 628. [Google Scholar] [CrossRef] [Green Version]
  11. Cermak, J.; Bendix, J. A novel approach to fog/low stratus detection using Meteosat 8 data. Atmos. Res. 2008, 87, 279–292. [Google Scholar] [CrossRef]
  12. Zhang, S.; Yi, L. A comprehensive dynamic threshold algorithm for daytime sea fog retrieval over the Chinese adjacent seas. Pure Appl. Geophys. 2013, 170, 1931–1944. [Google Scholar] [CrossRef]
  13. Bendix, J.; Thies, B.; Cermak, J.; Nauss, T. Ground fog detection from space based on MODIS daytime data-a feasibility study. Weather Forecast. 2005, 20, 989–1005. [Google Scholar] [CrossRef]
  14. Cermak, J. SOFOS-A New Satellite-Based Operational Fog Observation Scheme. Ph.D. Thesis, Philipps-University, Marburg, Germany, 2006. [Google Scholar]
  15. Heo, K.Y.; Park, S.; Ha, J.J.; Shim, J.S. Algorithm for sea fog monitoring with the use of information technologies. Meteorol. Appl. 2012. [Google Scholar] [CrossRef]
  16. Müller, M.D.; Masbou, M.; Bott, A.; Janjic, Z. Fog Prediction in a 3D Model with Parametrized Microphysics, Proc. Symp. on Nowcasting and Very Short Range Forecasting; World Weather Research Programme: Toulouse, France; Available online: http://www.meteo.fr/cic/wsn05/DVD/index.html (accessed on 26 June 2005).
  17. Maronga, B.; Bosveld, F.C. Key parameters for the life cycle of nocturnal radiation fog: A comprehensive large-eddy simulation study. Q. J. R. Meteorol. Soc. 2017, 143, 2463–2480. [Google Scholar] [CrossRef]
  18. Bendix, J.A. satellite-based climatology of fog and low-level stratus in Germany and adjacent areas. Tmos. Res. 2002, 64, 3–18. [Google Scholar] [CrossRef]
  19. Jedlovec, G. Automated Detection of Clouds in Satellite Imagery. In Advances in Geoscience and Remote Sensing; Jedlovec, G., Ed.; Intech: Rijeka, Croatia, 2009. [Google Scholar]
  20. Eyre, J.R.; Brownscombe, J.L.; Allam, R.J. Detection of fog at night using Advanced Very High Resolution Radiometer (AVHRR) imagery. Meteorol. Mag. 1984, 113, 266–271. [Google Scholar]
  21. Ellord, G.P. Advances in the detection and analysis of fog at night using GOES multispectral infrared imagery. Weather Forecast. 1995, 10, 606–619. [Google Scholar] [CrossRef] [Green Version]
  22. Cermak, J.; Bendix, J. Dynamical nighttime fog/low stratus detection based on Meteosat SEVIRI data: A feasibility study. Pure Appl. Geophys. 2007, 164, 1179–1192. [Google Scholar] [CrossRef]
  23. Gao, S.; Wu, W.; Zhu, L.; Fu, G.; Huang, B. Detection of nighttime sea fog/stratus over the Huang-Hai sea using MTSAT-1R IR data. Acta Oceanol. Sin. 2009, 28, 23–35. [Google Scholar]
  24. Chaurasia, S.; Gohil, B.S. Detection of day time fog over India using INSAT-3D data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 4524–4530. [Google Scholar] [CrossRef]
  25. Kim, S.H.; Suh, M.S.; Han, J.H. Development of fog detection algorithm during nighttime using Himawari-8/AHI satellite and ground observation data. Asia -Pac. J. Atmos. Sci. 2018. [Google Scholar] [CrossRef]
  26. Shin, D.; Kim, J.H. A new application of unsupervised learning to nighttime sea fog detection. Asia -Pac. J. Atmos. Sci. 2018, 54, 527–544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Han, J.-H.; Suh, M.-S.; Yu, H.-Y.; Roh, N.-Y. Development of fog detection algorithm using GK2A/AMI and ground data. Remote Sens. 2020, 12, 3181. [Google Scholar] [CrossRef]
  28. Bendix, J.; Thies, B.; Nauß, T.; Cermak, J. A feasibility study of daytime fog and low stratus detection with TERRA/AQUA-MODIS over land. Meteorol. Appl. 2006, 13, 111–125. [Google Scholar] [CrossRef]
  29. Kim, D.; Park, M.S.; Park, Y.J.; Kim, W. Geostationary Ocean Color Imager (GOCI) marine fog detection in combination with Himawari-8 based on the decision tree. Remote Sens. 2020, 12, 149. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, J.L.; Li, S.M.; Liu, X.L.; Wu, X.J. An analysis of the fog distribution in Beijing for the 2001–2005 period using NOAA and FY data. Atmos. Res. 2010, 96, 575–589. [Google Scholar] [CrossRef]
  31. Drönner, J.; Korfhage, N.; Egli, S.; Mühling, M.; Thies, B.; Bendix, J.; Freisleben, B.; Seeger, B. Fast cloud segmentation using convolutional neural networks. Remote Sens. 2018, 10, 1782. [Google Scholar] [CrossRef] [Green Version]
  32. Andersen, H.; Cermak, J. First fully diurnal fog and low cloud satellite detection reveals life cycle in the Namib. Atmos. Meas. Tech. 2018, 11, 5461–5470. [Google Scholar] [CrossRef] [Green Version]
  33. Schreiner, A.J.; Ackerman, S.A.; Baum, B.A.; Heidinger, A.K. A multispectral technique for detecting low-level cloudiness near sunrise. J. Atmos. Ocean. Technol. 2007, 24, 1800–1810. [Google Scholar] [CrossRef]
  34. Lee, J.; Chung, C.; Ou, M. Fog detection using geostationary satellite data: Temporally continuous algorithm. Asia-Pac. J. Atmos. Sci. 2011, 47, 113–122. [Google Scholar] [CrossRef]
  35. Schillings, C.; Mannstein, H.; Meyer, R. Operational method for deriving high resolution direct normal irradiance from satellite data. Solar Energy. 2004, 76, 475–484. [Google Scholar] [CrossRef]
  36. Yang, J.H.; Yoo, J.M.; Choi, Y.S.; Wu, D.; Jeong, J.H. Probability index of low stratus and fog at dawn using dual geostationary satellite observations from COMS and FY-2D near the Korean Peninsula. Remote Sens. 2019, 11, 1283. [Google Scholar] [CrossRef] [Green Version]
  37. Jensen, J.R. Remote Sensing of the Environment: An Earth Resource Perspective, 2nd ed.; Prentice Hall, Brigham Young University: Upper Saddle River, NJ, USA, 2000; pp. 150–151. [Google Scholar]
  38. Yang, J.; Zhang, Z.; Wei, C.; Lu, F.; Guo, Q. Introducing the new generation of Chinese geostationary weather satellites, Fengyun-4. Bull. Am. Meteorol. Soc. 2017, 98, 1637–1658. [Google Scholar] [CrossRef]
  39. Zhu, J.; Shu, J.; Guo, W. Biases characteristics assessment of the Advanced Geosynchronous Radiation Imager (AGRI) measurement on board Fengyun–4A Geostationary Satellite. Remote Sens. 2020, 12, 2871. [Google Scholar] [CrossRef]
  40. Kazumori, M. Assimilation of Himawari-8 clear sky radiance data in JMA’s global and mesoscale NWP systems. J. Meteorol. Soc. Jpn. 2018, 96B, 173–192. [Google Scholar] [CrossRef]
  41. Akimoto, Y.; Kusaka, H. A climatological study of fog in Japan based on event data. Atmos. Res. 2015. [Google Scholar] [CrossRef]
  42. Smith, W.L.; Woolf, H.M.; Fleming, H.E. Retrieval of atmospheric temperature profiles from satellite measurements for dynamical forecasting. J. Appl. Meteorol. Climatol. 1972, 11, 113–122. Available online: https://journals.ametsoc.org/view/journals/apme/11/1/1520-0450_1972_011_0113_roatpf_2_0_co_2.xml (accessed on 12 September 2020). [CrossRef]
  43. Dozier, J. Spectral signature of alpine snow cover from the Landsat thematic mapper. Remote Sens. Environ. 1989, 28, 9–22. [Google Scholar] [CrossRef]
  44. Kulkarni, A.V.; Singh, S.K.; Mathur, P.; Mishra, V.D. Algorithm to monitor snow cover using AWiFS data of RESOURCESAT-1 for the Himalayan region. Int. J. Remote Sens. 2006, 27, 2449–2457. [Google Scholar] [CrossRef]
  45. Wu, W.; Li, S. Automatic sea fog detection over Chinese adjacent oceans using Terra/MODIS data. Int. J. Remote Sens. 2014, 35, 7430–7457. [Google Scholar] [CrossRef]
  46. Ryu, H.S.; Hong, S. Sea Fog Detection Based on Normalized Difference Snow Index Using Advanced Himawari Imager Observations. Remote Sens. 2020, 12, 1521. [Google Scholar] [CrossRef]
  47. Wylie, D.P.; Menzel, W.P. Eight years of global high cloud statistics using HIRS. J. Clim. 1999, 12, 170–184. [Google Scholar] [CrossRef]
  48. Li, J.; Schmidt, C.C.; Nelson, J.P., III; Schmit, T.J.; Menzel, W.P. Estimation of total atmospheric ozone from GOES sounder radiances with high temporal resolution. J. Atmos. Ocean. Technol. 2001, 18, 157–168. [Google Scholar] [CrossRef]
  49. UK Met Office. Aerodrome Weather Report–World Meteorological Organization Archived 2012-02-24 at the Wayback Machine. 2012. Available online: https://web.archive.org/web/20120518165852/http://www.metoffice.gov.uk/media/pdf/4/1/Get_Met_2012_amends_tag_Artwork.pdf (accessed on 20 June 2019).
  50. Wetzel, M.A.; Borys, R.D.; Xu, L.E. Satellite microphysical retrievals for land-based fog with validation by balloon profiling. J. Appl. Meteor. 1996, 35, 810–829. [Google Scholar] [CrossRef] [Green Version]
  51. Ricchiazzi, P.; Yang, S.; Gautier, C.; Sowle, D. SBDART: A research and teaching software tool for plane-parallel radiative transfer in the Earth’s atmosphere. Bull. Am. Meteorol. Soc. 1998, 79, 2101–2114. [Google Scholar] [CrossRef] [Green Version]
  52. Egli, S.; Thies, B.; Dronner, J.; Cermak, J.; Bendix, J. A 10 year fog and low stratus climatology for Europe based on Meteosat Second Generation data. Q. J. R. Meteorol. Soc. 2016, 143, 530–541. [Google Scholar] [CrossRef]
  53. Prabhakara, C.; Kratz, D.P.; Yoo, J.M.; Dalu, G.; Verneka, A. Optically thin cirrus clouds: Radiative impact on the warm pool. Quant. Spectrosc. Radiat. Transfer. 1993, 49, 467–483. [Google Scholar] [CrossRef]
  54. Chung, S.R.; Ahn, M.H.; Han, K.S.; Lee, K.T.; Shin, D.B. Meteorological products of Geo-KOMPSAT 2A (GK2A) satellite. Asia-Pac. J. Atmos. Sci. 2020, 56, 185. [Google Scholar] [CrossRef] [Green Version]
  55. Marchant, B.; Platnick, S.; Meyer, K.; Wind, G. Evaluation of the MODIS collection 6 multilayer cloud detection algorithm through comparisons with CloudSat Cloud Profiling Radar and CALIPSO CALIOP products. Atmos. Meas. Tech. 2020, 13, 3263–3275. [Google Scholar] [CrossRef]
  56. Teng, S.; Liu, C.; Zhang, Z.; Wang, Y.; Sohn, B.J.; Yung, Y.L. Retrieval of ice-over-water cloud microphysical and optical properties using passive radiometers. Geophys. Res. Lett. 2020, 47, e2020GL088941. [Google Scholar] [CrossRef]
Figure 1. Viewing zenith angles (VZAs) of two geostationary satellites (Himawari-8 and FY-4A) available for near-simultaneous observations of the low stratus and fog (LSF) at dawn near Japan. The VZA difference between the satellites is 40.4° at the Himawari-8 nadir point.
Figure 1. Viewing zenith angles (VZAs) of two geostationary satellites (Himawari-8 and FY-4A) available for near-simultaneous observations of the low stratus and fog (LSF) at dawn near Japan. The VZA difference between the satellites is 40.4° at the Himawari-8 nadir point.
Remotesensing 13 01042 g001
Figure 2. Locations of 31 meteorological stations near Japan used for LSF verification during 2018–2019. The stations are located in Japan except for the 31st Ulleung station in South Korea.
Figure 2. Locations of 31 meteorological stations near Japan used for LSF verification during 2018–2019. The stations are located in Japan except for the 31st Ulleung station in South Korea.
Remotesensing 13 01042 g002
Figure 3. Flow diagram for LSF detection in terms of probability index (PI) near Japan at dawn during the period of March to November of 2018–2019, based on near-simultaneous satellite observations of Himawari-8 and FY-4A. The meaning of the acronyms in the diagram is explained in Table A1.
Figure 3. Flow diagram for LSF detection in terms of probability index (PI) near Japan at dawn during the period of March to November of 2018–2019, based on near-simultaneous satellite observations of Himawari-8 and FY-4A. The meaning of the acronyms in the diagram is explained in Table A1.
Remotesensing 13 01042 g003
Figure 4. Average values of eight satellite-observed variables with their standard deviation (±1σ) for two weather conditions (black for LSF and red for clear sky) at dawn during the control period of 2018 near Japan. (a) Brightness temperatures (BTD13.5-8.5, ΔBTD3.7-11, BTD3.7-11, and LST-BT11). (b) Reflectances (or NDSI) of ΔNDSI, NDSI, R0.65, and ΔR0.65.
Figure 4. Average values of eight satellite-observed variables with their standard deviation (±1σ) for two weather conditions (black for LSF and red for clear sky) at dawn during the control period of 2018 near Japan. (a) Brightness temperatures (BTD13.5-8.5, ΔBTD3.7-11, BTD3.7-11, and LST-BT11). (b) Reflectances (or NDSI) of ΔNDSI, NDSI, R0.65, and ΔR0.65.
Remotesensing 13 01042 g004
Figure 5. Scatter diagrams of (a) BTD13.5-8.5, (b) ΔBTD3.7-11, (c) ΔNDSI, and (d) R0.65 with respect to LST-BT11 for the weather phenomena of LSF (red circle) and clear sky (blue triangle) at dawn during 2018 near Japan. The values in the ordinate of the gray-shaded bands indicate the LSF threshold ranges for each variable. The satellite-observed value of one positive standard deviation (i.e., +1σ; upper boundary) for LSF detection is used to remove the middle or high-level clouds without accompanying LSF, except for ΔNDSI. The negative value of -1σ (lower boundary) is needed for ΔNDSI of Figure 5c based on the results shown in Figure 4b.
Figure 5. Scatter diagrams of (a) BTD13.5-8.5, (b) ΔBTD3.7-11, (c) ΔNDSI, and (d) R0.65 with respect to LST-BT11 for the weather phenomena of LSF (red circle) and clear sky (blue triangle) at dawn during 2018 near Japan. The values in the ordinate of the gray-shaded bands indicate the LSF threshold ranges for each variable. The satellite-observed value of one positive standard deviation (i.e., +1σ; upper boundary) for LSF detection is used to remove the middle or high-level clouds without accompanying LSF, except for ΔNDSI. The negative value of -1σ (lower boundary) is needed for ΔNDSI of Figure 5c based on the results shown in Figure 4b.
Remotesensing 13 01042 g005
Figure 6. RTM LSF simulation in the domain of (a) BTD13.5-8.5 vs. LST-BT11, (b) ΔBTD3.7-11 vs. LST-BT11, (c) ΔNDSI vs. LST-BT11, and (d) R0.67 vs. LST-BT11 at dawn during 2018 near Japan under the different LSF conditions. The yellow and blue symbols in the figures denote the fog layer at 0–1 km and 0–2 km without higher clouds. Middle clouds of water/ice at 4–6 km height are shown in red and navy, respectively. Green asterisk means higher cloud (8–10 km) above the 0–2km fog layer. The average values of the Himawari-8 angles (i.e., SZA, RAA, and VZA) for 176 LSF cases were used as RTM input. The RTM conditions are described in Table A2. The satellite-observed mean values are shown for LSF (black asterisk) and clear sky (gray-cross). The clear-sky values from the observations and simulations (pink asterisk) are depicted in black circles. Some simulated data out of range in the domain were not shown for a comparison with the observations of Figure 5.
Figure 6. RTM LSF simulation in the domain of (a) BTD13.5-8.5 vs. LST-BT11, (b) ΔBTD3.7-11 vs. LST-BT11, (c) ΔNDSI vs. LST-BT11, and (d) R0.67 vs. LST-BT11 at dawn during 2018 near Japan under the different LSF conditions. The yellow and blue symbols in the figures denote the fog layer at 0–1 km and 0–2 km without higher clouds. Middle clouds of water/ice at 4–6 km height are shown in red and navy, respectively. Green asterisk means higher cloud (8–10 km) above the 0–2km fog layer. The average values of the Himawari-8 angles (i.e., SZA, RAA, and VZA) for 176 LSF cases were used as RTM input. The RTM conditions are described in Table A2. The satellite-observed mean values are shown for LSF (black asterisk) and clear sky (gray-cross). The clear-sky values from the observations and simulations (pink asterisk) are depicted in black circles. Some simulated data out of range in the domain were not shown for a comparison with the observations of Figure 5.
Remotesensing 13 01042 g006aRemotesensing 13 01042 g006b
Figure 7. Statistical verification of (a) POD minus FAR, (b) POD, (c) FAR, and (d) CSI with respect to the LSF threshold values of eight satellite-observed variables during the experimental period of 2019. The verification was performed based on the total ground-observations of 171 fog and 110 clear-sky occurrences at dawn near Japan.
Figure 7. Statistical verification of (a) POD minus FAR, (b) POD, (c) FAR, and (d) CSI with respect to the LSF threshold values of eight satellite-observed variables during the experimental period of 2019. The verification was performed based on the total ground-observations of 171 fog and 110 clear-sky occurrences at dawn near Japan.
Remotesensing 13 01042 g007
Figure 8. (a) Same as Figure 4a except for the period of 2019. (b) Same as Figure 4b except for the period of 2019.
Figure 8. (a) Same as Figure 4a except for the period of 2019. (b) Same as Figure 4b except for the period of 2019.
Remotesensing 13 01042 g008
Figure 9. The number of cases that have passed the threshold test of the four satellite-observed variables (i.e., BTD13.5-8.5 in gray rectangle, ΔNDSI in orange, ΔBTD3.7-11 in blue, and R0.65 in green) for detection of (a) LSF and (b) clear sky at dawn during 2019 near Japan. The zones of red (Class 1 for the possibility of the weather phenomena), yellow (Class 2), blue (Class 3), and white (Class 4) are defined in Table 4.
Figure 9. The number of cases that have passed the threshold test of the four satellite-observed variables (i.e., BTD13.5-8.5 in gray rectangle, ΔNDSI in orange, ΔBTD3.7-11 in blue, and R0.65 in green) for detection of (a) LSF and (b) clear sky at dawn during 2019 near Japan. The zones of red (Class 1 for the possibility of the weather phenomena), yellow (Class 2), blue (Class 3), and white (Class 4) are defined in Table 4.
Remotesensing 13 01042 g009
Figure 10. Spatial distributions of fog probability at dawn (06:00 LST) on 9 July 2019, near Japan from detection methods of (a) R0.65 threshold, (b) BTD3.7-11 threshold, and (c) LSF PI in this study. (d) BT11 distribution for cloud-top temperatures, and (e) SYNOP map for station-observed fog occurrences. Fog occurrences at the stations (pink triangles in Figure 9a,b) were reported by ASOS.
Figure 10. Spatial distributions of fog probability at dawn (06:00 LST) on 9 July 2019, near Japan from detection methods of (a) R0.65 threshold, (b) BTD3.7-11 threshold, and (c) LSF PI in this study. (d) BT11 distribution for cloud-top temperatures, and (e) SYNOP map for station-observed fog occurrences. Fog occurrences at the stations (pink triangles in Figure 9a,b) were reported by ASOS.
Remotesensing 13 01042 g010
Figure 11. (a) Same as Figure 10a except for the date and time (05:30 LST on 4 July 2019), (b) Same as Figure 10b except for the date and time. (c) Same as Figure 10c except for the date and time. (d) Same as Figure 10d except for the date and time, and no SYNOP map (not available at the time).
Figure 11. (a) Same as Figure 10a except for the date and time (05:30 LST on 4 July 2019), (b) Same as Figure 10b except for the date and time. (c) Same as Figure 10c except for the date and time. (d) Same as Figure 10d except for the date and time, and no SYNOP map (not available at the time).
Remotesensing 13 01042 g011
Table 1. Information from the dual geostationary satellites of Himawari-8 and FY-4A near Japan, available for near-simultaneous observations of low stratus and fog (LSF) at dawn during 2018–2019. The meaning of acronyms is explained in Table A1.
Table 1. Information from the dual geostationary satellites of Himawari-8 and FY-4A near Japan, available for near-simultaneous observations of low stratus and fog (LSF) at dawn during 2018–2019. The meaning of acronyms is explained in Table A1.
Satellite (Nation, Lon at nadir)
Himawari-8 (Japan, 140.7 °E)FY-4A (China, 104.7 °E)
Channel
(Abbreviation)
Application of this StudyWavelength (μm)Central Wavelength (μm)Resolution at Nadir (km)Wavelength (μm)Central Wavelength (μm)Resolution at Nadir (km)
VIS (R0.65)Fog, cloud0.63–0.660.640.50.55–0.750.650.5–1
NIR (R1.6)Cloud, snow1.60–1.621.6121.58–1.641.612
SWIR (BT3.7)Fog, cloud3.74–3.963.8923.50–4.003.752
MIR (BT8.5)Cloud-top phase, SST, LST8.44–8.768.5928.00–9.008.504
LWIR1 (BT11)SST, LST, cloud-top temperature11.1–11.311.24210.3–11.310.804
LWIR3 (BT13.5)Cloud, air temperature13.2–13.413.28213.2–13.813.504
Table 2. Thirty-one meteorological stations of the Automatic Synoptic Observing System (ASOS) near Japan used for the LSF analysis in this study.
Table 2. Thirty-one meteorological stations of the Automatic Synoptic Observing System (ASOS) near Japan used for the LSF analysis in this study.
Station
Number
Station
in Japan
Lat (°N)Lon (°E)Height
(m)
Station
Number
Station
in Japan
Lat (°N)Lon (°E)Height
(m)
1AKENO (JASDF)34.53136.67917KUMAMOTO (CIV/JAS)32.84130.86196
2AKITA AIRPORT39.62140.229618KUSHIRO AIRPORT43.04144.1998
3AOMORI AIRPORT40.73140.6920119MATSUYAMA AIRPORT33.83132.707
4ASAHIKAWA AIRPORT43.67142.4521120MIYAZAKI AIRPORT31.88131.459
5CHITOSE42.77141.693021NEW HIROSHIMA34.44132.926
6FUKUI AIRPORT36.14136.22822NEW TOKYO INTL AIRPORT35.77140.3944
7FUKUOKA/
ITAZUKE
33.58130.451223OKAYAMA AIRPORT34.76133.86242
8FUKUSHIMA ARPT37.23140.4337524OMINATO (JASDF)41.23141.1310
9HACHINOHE (JMSDF)40.55141.474925SAGA AIRPORT33.15130.305
10HAKODATE AIRPORT41.77140.823626SENDAI AIRPORT38.14140.925
11HYAKURI (JASDF)36.18140.413527SHIMOFUSA (JMSDF35.80140.0133
12IZUMO AIRPORT35.41132.89528TAKAMATSU AIRPORT34.22134.02188
13KAGOSHIMA AIRPORT31.80130.7227529TOKUSHIMA(JMSDF)34.13134.6111
15KANSAI INTL34.43135.23830Yonaguni24.47122.9819
15KISARAZU (JGSDF)35.40139.916 Island station
in South Korea
16KOCHI AIRPORT33.55133.671031Ulleung37.48130.90223
Table 3. Statistical verification of eight threshold values for LSF detection at dawn near Japan from the contingency table during the period of March to November of 2018, using the Himawari-8 and FY-4A satellites. POD: Probability of Detection, CSI: Critical Success Index, HSS: Heidke Skill Scores, PC: Percentage Correct, FAR: False Alarm Ratio.
Table 3. Statistical verification of eight threshold values for LSF detection at dawn near Japan from the contingency table during the period of March to November of 2018, using the Himawari-8 and FY-4A satellites. POD: Probability of Detection, CSI: Critical Success Index, HSS: Heidke Skill Scores, PC: Percentage Correct, FAR: False Alarm Ratio.
Satellite-Derived ThresholdSkill Score
PODCSIHSSPCFARPOD-FAR
−24 K < BTD13.5-8.5 < −10 K0.8810.8330.7850.8950.0610.820
−0.1 < ΔNDSI < 0.30.8810.8290.7780.8920.0660.815
0.19 < R0.65 < 0.520.8470.8190.7740.8880.0390.808
7 K < ΔBTD3.7-11 < 19 K0.8470.8050.7530.8780.0570.790
4.5 K < LST-BT11 < 37.5 K0.6990.5750.3720.6920.2360.463
0.29 < ΔR0.65 < 0.550.4600.4360.3370.6440.1100.350
3 K < BTD3.7-11 < 15 K0.4660.3910.1670.5660.2930.173
−0.18 < NDSI < 0.200.3300.3090.2000.5590.1710.159
Table 4. Five LSF classes derived from 16 cases of the four variables (A: BTD13.5-8.5, B: ΔNDSI, C: R0.65, D: ΔBTD3.7-11) based on the threshold test for detection of either LSF or clear-sky phenomena. For instance, Case 1, which passed the test for LSF detection by all four variables, corresponds to the “very high” probability of Class 1 in terms of the probability set “P ( A B C D ) ”. The superscript “c” in A C , which consists of the elements that are not in A in set theory, means the complement of a set A .
Table 4. Five LSF classes derived from 16 cases of the four variables (A: BTD13.5-8.5, B: ΔNDSI, C: R0.65, D: ΔBTD3.7-11) based on the threshold test for detection of either LSF or clear-sky phenomena. For instance, Case 1, which passed the test for LSF detection by all four variables, corresponds to the “very high” probability of Class 1 in terms of the probability set “P ( A B C D ) ”. The superscript “c” in A C , which consists of the elements that are not in A in set theory, means the complement of a set A .
Assigned
Probability for 16 Cases
LSF Class (Possibility of Occurrence)
1 (Very High)2 (High)3 (Medium)4 (Low)5 (None)
1P ( A B C D )
0.75 P ( A B C D C )
0.75 P ( A B D C C )
0.75 P ( A C D B C )
0.75 P ( B C D A C )
0.5 P ( A B C C D C )
0.5 P ( A C B C D C )
0.5 P ( A D B C C C )
0.5 P ( B C A C D C )
0.5 P ( B D A C C C )
0.5 P ( C D A C B C )
0.25 P ( A B C C C D C )
0.25 P ( B A C C C D C )
0.25 P ( C A C B C D C )
0.25 P ( D A C B C C C )
0 P ( A C B C C C D C )
Table 5. Comparison of the results of this study with those of two previous studies of the Dual Satellite Method (DSM) for LSF detection at dawn [9,36].
Table 5. Comparison of the results of this study with those of two previous studies of the Dual Satellite Method (DSM) for LSF detection at dawn [9,36].
DetailsDSM Study
Yoo et al. (2018) [9]Yang et al. (2019) [36]This Study
Geostationary satellites
(Number of channels)
COMS (6) and
FY-2D (6)
COMS (6) and
FY-2D (6)
Himawari-8 (16) and
FY-4A (14)
Main research areaSouth KoreaSouth KoreaJapan
Spatial grid resolution~4 km~4 km~4 km
Satellite viewing zenith angle (VZA) difference at nadir46.5°46.5°40.4°
Period2013–20162013–20162018–2019
SeasonJune to AugustApril to AugustMarch to November
DSM variableΔR0.67ΔR0.67 and ΔBTD3.7-11ΔBTD3.7-11 and ΔNDSI
Additional variableSuggests R0.67R0.67BTD13.5-8.5 and R0.65
Probability index of LSFNoYesYes
Number of variables for PI calculationNoneThreeFour
Number of LSF occurrence classes3 or 485 extensible to 16 *
* The maximum class will be available based on the probability combination of the four variables (Table 4) if the long-term dataset of the LSF and clear-sky occurrences is given in the future.
Table 6. Satellite-derived PI distributions of five classes weighted at 0.25 intervals for detection of LSF and clear-sky cases in 2019.
Table 6. Satellite-derived PI distributions of five classes weighted at 0.25 intervals for detection of LSF and clear-sky cases in 2019.
LSF
Class
Assigned
Probability
Thlower < LSF < ThupperEither Clear Sky < Thlower or Clear Sky > Thupper
Number of
Actual
Occurrences
Sum of Detection ProbabilityNumber of Actual
Occurrences
Sum of Detection
Probability
119999.007979.00
20.754332.2586.00
30.52110.50126.00
40.2520.50112.75
5060.0000.00
Total171140.2511093.75
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, J.-H.; Yoo, J.-M.; Choi, Y.-S. Advanced Dual-Satellite Method for Detection of Low Stratus and Fog near Japan at Dawn from FY-4A and Himawari-8. Remote Sens. 2021, 13, 1042. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13051042

AMA Style

Yang J-H, Yoo J-M, Choi Y-S. Advanced Dual-Satellite Method for Detection of Low Stratus and Fog near Japan at Dawn from FY-4A and Himawari-8. Remote Sensing. 2021; 13(5):1042. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13051042

Chicago/Turabian Style

Yang, Jung-Hyun, Jung-Moon Yoo, and Yong-Sang Choi. 2021. "Advanced Dual-Satellite Method for Detection of Low Stratus and Fog near Japan at Dawn from FY-4A and Himawari-8" Remote Sensing 13, no. 5: 1042. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13051042

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop