1. Introduction
The use of satellite information has improved numerical weather prediction results and increased forecasts’ reliability [
1,
2]. As a result, the contribution of satellite data to numerical forecasting systems has steadily grown. Statistics from the European Centre for Medium-Range Weather Forecasts (ECMWF) show that 91.41% of data input to the assimilation system come from satellite observations [
3]. Having analyzed and assessed the impact of the initial decade of use of the TIROS Operational Vertical Sounder (TOVS) on weather forecasting, the World Meteorological Organization (WMO) concluded that significant improvements in weather forecasting capabilities could exclusively improve the precision of global atmospheric vertical temperature and humidity profiles to equal that of radio soundings [
4]. Among the various types of satellite remote sensing data, high-spectral infrared data have emerged as crucial for reducing forecast errors and enhancing assimilation effects in all satellite observations [
5].
Spaceborne infrared hyperspectral sounders, which boast a high detection capability, a high spectral resolution, and spectral coverage of CO
2 and H
2O absorption bands, predominantly provide information on atmospheric temperature and humidity [
6]. Their high spectral resolution provides a high vertical resolution, allowing for fine-scale detection of meteorological elements (especially temperature and humidity) throughout the atmospheric troposphere. They are the most important satellite data source, reducing forecast errors and improving assimilation effectiveness [
5]. In particular, an infrared hyperspectral sounder mounted on a geostationary mereological satellite platform facilitates continuous observations of the covered area with a high temporal resolution. As small- and medium-scale weather systems evolve relatively quickly, these high temporal resolution instruments are important for detecting high-impact weather events that are strongly associated with atmospheric thermodynamic instability, as well as convective systems that develop over a short period of time [
7]. Therefore, observations from geostationary hyperspectral infrared sounders are crucial for weather monitoring and numerical prediction models [
8].
FY-4A is a Chinese second-generation geostationary meteorological satellite, positioned at 104.7°E over the equator. It has a more than 20-fold increased observation efficiency and can perform 24 h Earth observations from high altitudes [
9]. The Geostationary Interferometric Infrared Sounder (GIIRS) is the world’s first precision infrared hyperspectral instrument specifically designed to measure the vertical structure of the atmosphere in geostationary orbit via infrared interferometric spectroscopy. The GIIRS covers the mid-wave infrared (1650–2250 cm
−1) and longwave infrared (700–1130 cm
−1) bands at a spectral resolution of 0.625 cm
−1. It comprises 1650 observation channels, with 689 channels in the longwave and 961 channels in the mid-wave band.
Atmospheric clouds exhibit a high forecast sensitivity [
10], and their role in the radiative balance of the Earth’s atmospheric system is determined by their optical properties [
11]. Clouds consist of water droplets and ice crystals, both of which strongly absorb infrared radiation, and can be approximated as black bodies [
12]. If clouds are present in the field of view, the infrared radiation emitted by the atmosphere below clouds cannot be distinguished from the radiance from clouds [
13]. Cloud detection, which determines if there is any cloud cover in each satellite instrument’s instantaneous field of view, is a fundamental and essential step for atmospheric parameter retrieval and cloud microphysical property retrieval. Cloud detection also plays a critical role in clear-sky and cloud-area satellite data assimilation. Therefore, in practical applications of infrared data, it is often necessary to first eliminate cloud-contaminated scenes [
14].
The common methods for cloud detection can be grouped into four types: finding ‘clear FOVs’ [
15]; ‘clear channels’ [
16]; ‘cloud clearing’ [
17]; and ‘cloud product matching between different instruments’ [
18]. The ‘clear FOV’ approach aims to find completely cloud-free fields of view (FOVs) by excluding observations from all channels within cloudy FOVs. For example, Wang et al. selected optimal channels for the FY-4A GIIRS’s mid-wave infrared band and used the minimum residual method to calculate the effective cloud cover for cloud classification [
19]. The ‘clear channel’ method does not completely eliminate clouds from the satellite’s FOV. Instead, it assigns cloud-free channel heights to different channels based on their varying sensitivities to clouds. This method, proposed by McNally et al. (2003), is employed by the ECMWF for cloud detection using an Atmospheric Infrared Sounder (AIRS) and an Infrared Atmospheric Sounding Interferometer (IASI) [
20]. The ‘cloud clearing’ method combines moderate-resolution spectroradiometer imaging data with AIRS FOV radiance data equivalent to a clear sky and was proposed by Li et al. Finally, as an example of cloud product matching, Guan et al. (2007). performed spatial and temporal matching of high-resolution Moderate-resolution Imaging Spectroradiometer (MODIS) imager cloud detection products to AIRS’s FOVs for successful cloud detection [
21]. In addition, based on the clear channel method, Li et al. added Advanced Himawari Imager (AHI) cloud mask products to select clear sky areas and remove high-level clouds. This two-step cloud detection method serves as a complementary method for clear sky areas [
22].
These cloud detection algorithms are mainly used on polar-orbiting satellites, which have limitations in capturing rapidly developing weather systems due to their poor temporal resolution. In contrast, geostationary satellites offer high spectral and temporal resolutions. They can continuously monitor rapid and intense changes in water vapor levels and temperatures during severe weather events when compared to infrared hyperspectral detection instruments on polar-orbiting satellites. At present, the GIIRS Level 2 (L2) operational cloud mask (CLM) product relies on the Advanced Geostationary Radiation Imager (AGRI) mounted on the same satellite [
23]. The high-spatial-resolution imager cloud mask product is spatially and temporally matched to the low-spatial-resolution GIIRS’s FOV. Cloud detection in the GIIRS’s FOV is based on the proportion of AGRI clear/cloud pixels within it. This method only processes clear-sky FOVs and eliminates scenes contaminated by clouds, which wastes a large amount of observation data. Thus, it is desirable to develop a cloud detection algorithm that is based solely on GIIRS observations and does not rely on other instruments.
A cloud detection algorithm based solely on FY-4A GIIRS observations (advanced infrared hyperspectral cloud detection, hereafter referred to as AIRCD) is proposed in this paper according to the basic principles of cloud detection proposed by the Cross-track Infrared and Microwave Sounding Suite Environmental Data Records (CrIMSS EDR) algorithm [
24]. The algorithm employs the following four steps to extract cloud information and determine cloud coverage in the GIIRS’s FOV: (1) clear FOV identification, (2) cloud amount estimation, (3) thermal contrasting, and (4) cloud classification.
3. Cloud Detection Algorithm Based on GIIRS Observations
As mentioned in
Section 2.1, for the cloud detection algorithm, adjacent 2 × 2 FOVs were selected as samples according to the arrangement of the FY-4A GIIRS’s scanning detectors. Each FOR contains 16 × 2 cloud detection clusters, as shown by the black circles in
Figure 1.
The physical basis of cloud detection lies in the fact that when detecting the wavelength in a specific infrared channel, the detected thermal energy of an object is mainly related to its own temperature [
26]. What sets clouds apart from other underlying surfaces is that they have high reflectance and low brightness temperatures. During cloud detection tasks, brightness temperature differences can be used to distinguish between clouds and a clear sky. The principle is that a single scattering by clouds results in changes in spectral information. Additionally, there is a nonlinear variation in the Planck function during this process; in other words, it expresses the temperature difference between clouds and different surface types. The absorption and emission of other factors are not considered in this process. Changes in spectral information are often influenced by the microphysical properties of the clouds. Therefore, infrared brightness temperature differences can distinguish between clouds and clear skies [
27].
The most significant infrared absorption bands for atmospheric CO2 are the 4.3 μm and 15 μm bands. When clouds are present in the FOV, the attenuation of brightness temperature is significantly greater in the shortwave infrared channel than in the longwave infrared channel; this discrepancy arises from the strong scattering effect of ice particles and can be effectively used for cloud detection. Therefore, we based the AIRCD algorithm on the dual gas absorption bands of CO2 in the GIIRS, which are the longwave infrared band at 709.5–746.0 cm−1 and the shortwave infrared band at 2190–2250 cm−1.
A flow chart of the cloud detection algorithm is shown in
Figure 2. The process is divided into four steps: clear FOV identification, cloud amount estimation, thermal contrasting, and cloud classification. The classification results within each cloud detection cluster are divided into the following categories: clear, partly cloudy, or overcast.
3.1. Clear FOV Identification
We observed no significant differences in the overall observed radiance in each FOV when the cloud detection cluster was in either overcast or clear sky conditions (i.e., the scene was relatively homogeneous) [
28]. Therefore, to determine whether a cloud detection cluster consisting of 2 × 2 FOVs was clear, the number of clear FOVs within the cluster was first calculated. The clear radiance of each FOV was simulated using the CRTM. The root-mean-square error of radiances (
) between the observed and simulated clear radiations for all channels in the longwave infrared band of 709.5–746.0 cm
−1 was calculated as follows:
where
and
denote the radiances in FOV
and channel
for the observed and simulated clear FOVs, respectively, and ‘
’ is the 59 channels in the 709.5–746.0 cm
−1 band.
Ideally, if the FOV is clear, the radiance observed by the satellite should be the same as the simulated clear radiance. However, the difference between the observed and simulated radiance can be affected by instrument noise and the simulation error of the radiative transfer model. Thus, the conditions used to determine a clear FOV were as follows:
where
denotes the noise-equivalent radiance in FOV
and channel
and comes from the GIIRS Level 1 dataset.
was obtained by averaging the
of all channels within the 709.5–746.0 cm
−1 band for each FOV
j. If
, the FOV is clear; otherwise, it is overcast [
24]. The number of clear FOVs within the cloud detection cluster is denoted as
in Equations (1) and (2).
3.2. Cloud Amount Estimation
In this step, a principal component analysis is used to determine the number of cloudy FOVs within a cloud detection cluster and to initially distinguish between clear, overcast, and partly cloudy conditions. For the four FOVs within the cloud detection cluster, the GIIRS-observed radiance of all channels in the longwave infrared band ranges from 709.5 cm
−1 to 746.0 cm
−1, forming a 4 × 59 matrix. A principal component analysis was applied to this matrix to obtain the four orthogonal principal components and their corresponding eigenvalues. Susskind et al., (2003) indicated that principal components with larger eigenvalues are informative, whereas others are typically associated with the observed noise [
29]. The cloud amount is approximately equal to the number of significant principal components within the cloud detection cluster that corresponds to the cloud features.
First, we calculated the residual standard deviation (RSD) of the last few eigenvalues and estimated the observed noise,
:
where
is the eigenvalue and
,
, where n takes values from one in sequence. When
, the number of significant principal components is
and the cloud amount
.
Second, the chi-square,
was calculated in turn from the observed radiance and radiance for the first n principal component reconstructions,
:
The reconstructed radiance, was computed using the mat function with four principal components and their corresponding eigenvalues. As the number of principal components, n, used for radiance reconstruction increases, the minimum n that satisfies corresponds to the significant principal components related to the cloud. At this point, the reconstructed radiance, was already approximately representative of the original radiance, and the cloud amount . After comparing the values of and obtained in the above two steps, we used the larger value as the cloud amount, , of the cloud detection cluster.
3.3. Thermal Contrasting
The difference in observed radiance between FOVs is more significant when partial cloud contamination exists within the cloud detection cluster. In contrast, the observed radiance between FOVs is similar under overcast or clear sky conditions. Thus, to further distinguish partly cloudy conditions, a thermal contrast test is required.
The average observed radiances for the 709.5–746.0 cm
−1 and 2190–2250 cm
−1 channels of the four FOVs within the cloud detection cluster were sorted in descending order. The FOV with the highest average radiance was labeled as the ‘warmest’ FOV, whereas that with the lowest average radiance was labeled as the ‘coldest’ FOV. The absolute difference in observations between the ‘warmest’ and ‘coldest’ FOVs was then calculated and compared with the GIIRS channel noise.
In Equation (6),
is the absolute difference between the observed radiance,
in channel
of the ‘warmest’ FOV and the observed radiance,
of the ‘coldest’ FOV in the regions of 709.5–746.0 cm
−1 and 2190–2250 cm
−1. In Equation (7),
is the noise variance in channel
of the ‘warmest’ FOV, and the coefficient 4.246 is from reference [
24]. We defined the number of channel thermal contrasts,
, within each cloud detection cluster. The initial value was zero; when
,
increased by one. By iterating through all channels of the GIIRS in the spectral regions of 709.5–746.0 cm
−1 and 2190–2250 cm
−1, the total number of channels with a radiance difference greater than the instrument noise was obtained as
.
3.4. Cloud Classification
Following the previous three steps, we used
(the number of clear FOVs within each cloud detection cluster),
(the number of cloudy FOVs), and
(the number of channel thermal contrasts) to categorize the cloud detection clusters as clear, partly cloudy, or overcast according to the process shown in
Figure 3. If
≤ 1, this indicates that there are very few clouds within this cluster, and then we further assess the value of
. If
> 2, this means that most of the FOVs are clear, so the cluster is identified as clear sky; otherwise, it is overcast. If
, this indicates the presence of clouds in this cluster, and then
and
are used to continue the analysis. If
and
, the cluster is overcast; otherwise, it is partly cloudy.
4. Results and Discussion
The AIRCD algorithm was applied to GIIRS radiance observations from 1200 UTC to 1340 UTC on 13 May 2022, as an example.
Figure 4 illustrates the crucial parameters involved in the algorithm for the cloud detection cluster located at [43.09°E, 94.12°N].
Figure 4a shows the eigenvalues obtained after the principal component analysis. The eigenvalues decreased rapidly as n increased, with larger eigenvalues indicating a stronger correlation between the principal components and cloud features.
Figure 4b shows the calculated RSD values when the first n eigenvalues were considered. When n = 3, the RSD value satisfies the condition of being less than the observation noise; consequently,
.
Figure 4c shows the calculated
values when considering the first n eigenvalues, with the diagonal representing the set threshold. When n = 3, the threshold was satisfied; thus,
. Taken together, these results indicate that
; thus, the cloud detection cluster was classified as partly cloudy according to the threshold classification in
Figure 2.
Figure 5a shows the GIIRS brightness temperature in the longwave infrared window channel at 900 cm
−1 observed between 1200 UTC and 1340 UTC on 13 May 2022. The observed radiance in the infrared window channel is related to the target’s temperature; colder-toned areas correspond to lower brightness temperatures, indicating higher cloud heights and lower cloud-topped temperatures, whereas warmer-toned areas represent higher brightness temperatures, typically associated with clear sky areas or low-level clouds.
Figure 5b shows the results obtained using the AIRCD algorithm. The green, red, and blue areas indicate clear, partly cloudy, and overcast FOVs, respectively. The results show that the identified cloudy FOVs (both partly cloudy and overcast) were consistent with the cold-colored areas with lower brightness temperatures in
Figure 5a.
To further evaluate the performance of the AIRCD algorithm, we employed the closest AGRI’s real-time cloud detection product to the GIIRS observation time (
Figure 5c). The AGRI is an imager on the same satellite platform as the GIIRS that possesses highly reliable CLM products with a high spatial resolution of 4 km, as well as visible, near-infrared, and longwave infrared channels. The CLM products were categorized as clouds (blue), clear (green), probably clear (yellow), and probably clouds (red).
Figure 5d shows the GIIRS’s L2 operational CLM product, with green areas representing clear skies and blue areas indicating cloudy skies. A comparison of
Figure 5b–d shows that the GIIRS’s operational product identified fewer clear FOVs. The AIRCD algorithm identified clear sky areas that were consistent with the AGRI’s cloud mask product and the warm brightness temperature images in the window channel at 900 cm
−1; however, the GIIRS’s operational product identified more cloudy FOVs. The areas where the algorithmic products differed significantly from the AGRI’s real-time cloud detection products were primarily along the clear sky/cloud boundary and in convective cloud areas with a relatively low cloud coverage. This was attributed to the fact that the GIIRS’s cloud detection cluster has a lower spatial resolution (approximately 32 km) than the AGRI’s 4 km FOV over the same area. These differences were more evident when there was a lack of homogeneity within a large FOV. For example, in the case of fine-scale cumulus clouds, many cloud-free areas within the FOV were classified as partly cloudy at a coarse resolution. Consequently, the AIRCD algorithm identified slightly more cloudy areas than the AGRI’s product with a higher spatial resolution.
In the quantitative assessment process, considering its inconsistency between the AGRI and the GIIRS, AGRI’s CLM product should not be used directly as a reference value to quantify the AIRCD algorithm’s performance. Therefore, spatial matching between the AGRI and the GIIRS is necessary before evaluating the AIRCD’s identification results. The AGRI’s pixels were considered to fall within the GIIRS’s FOV if the distance between the central latitude and longitude of the GIIRS and AGRI was less than a certain threshold (Equation (8)), indicating spatial matching [
30]. In Equation (8), R is the Earth’s radius (6371 km) and ‘arccos’ is the inverse cosine transformation. If FOV deformations due to changes in the satellite’s observation angle are not considered, the threshold should be the radius of the GIIRS’s FOV, which is 8 km. However, as the satellite scanning angle increases, the FOV gradually becomes elliptical and may approach an egg shape, making it difficult to define using mathematical equations. Therefore, a threshold of 9 km was used in this study [
31].
The spatial collocations of the GIIRS and AGRI data are shown in
Figure 6. The large circle represents the GIIRS’s FOV that must be matched, whereas the small circles represent the number of AGRI pixels within it. The spatial resolution of the GIIRS’s FOV (16 km) was coarser than that of the AGRI (4 km). Approximately 4 × 4 AGRI pixels fall within one GIIRS’s FOV, as illustrated in
Figure 6.
The AGRI’s cloud mask is categorized as clear, probably clear, probably cloud, and cloud. The GIIRS’s cloud labels were categorized into three groups according to the proportion of clear or cloud pixels within the GIIRS’s FOV: (1) clear, if over 80% of the matched AGRI pixels were clear or probably clear; (2) overcast, if matched AGRI cloud pixels made up 87.5% or more, or combined cloud and probably cloud pixels totaled 100% (with at least 75% cloud pixels); or (3) partly cloudy, when neither of the above two situations applied.
According to the matching steps described above, by collocating the GIIRS’s FOVs with the AGRI pixels, 49,536 pairs of samples were obtained over a two-hour observation of China and its surrounding areas. An analysis of the monthly atmospheric circulation and weather reports published by the National Satellite Meteorological Centre reveals that severe convective weather primarily occurs in China from April to October every year [
32]. This study was tested using one month of GIIRS observations from May 2022. Observations were conducted ten times daily, covering China and its surrounding areas, with no data available from 1600 UTC to 1900 UTC [
33].
Given that cloud detection is a binary classification task, there are four potential outcomes: True Positive (TP), which represents predictions that are true and match actual positive samples; True Negative (TN), which means predictions that are true and match actual negative samples; False Positive (FP), indicating predictions that are false but match positive samples; False Negative (FN), meaning predictions that are false and match actual negative samples. These four metrics have been used in this study to assess the cloud classification results. Clear sky instances are considered positive samples, while partly cloudy and overcast conditions are considered negative samples.
Table 1 details the specific application of the classification metrics in the AIRCD algorithm.
Equations (9) and (10) are used to assess the performance of the cloud detection results using the false alarm rate (FAR) [
34] and hit rate (HR) [
35].
In this study, clear sky instances are considered positive samples. The FAR, often referred to as specificity, represents the proportion of clear-sky FOVs classified as partly cloudy or overcast by the AIRCD algorithm out of all the negative samples. HR, also known as sensitivity, corresponds to the proportion of FOVs consistently classified as clear sky by the AIRCD algorithm among all positive samples. It is generally accepted that an algorithm is effective when the HR approaches one and the FAR approaches zero.
In the infrared spectral region, clouds have strong absorption and scattering characteristics, impacting the detection performance of atmospheric profiles and the underlying surfaces [
36]. Therefore, the classification statistics are based on whether the solar zenith angle of FY-4A/GIIRS exceeds 75 degrees, dividing the data into day, night, and all-day periods and incorporating eight different surface types as provided by GIIRS L2 [
37]. The quantitative statistical results, presented in
Table 2, evaluate the classification performance of the AIRCD algorithm. To facilitate comparisons, both the FAR and HR results were multiplied by 100.
Regarding all-day observations, the AIRCD algorithm consistently exceeded a 68% hit rate for clear sky classifications over various surface types, including the deep ocean, shallow ocean, land, and ephemeral water. However, the lowest hit rate was 56% over deep inland water. The false alarm rates for clear sky classifications such as shallow inland water, ephemeral (intermittent) water, and deep inland water exceed 15%. The hit rates for partly cloudy classifications over ephemeral (intermittent) water and deep inland water are notably lower than for other surface types, falling below 30%. The false alarm rates for partly cloudy classifications ranged between 15 and 25%, peaking at 36.7% for deep inland water. Despite exhibiting a high identification performance in overcast areas, the AIRCD algorithm’s hit rate dipped below 70% over shallow inland water, ephemeral (intermittent) water, and deep inland water, with the lowest rate at 36.7% for deep inland water. The false alarm rates for overcast classifications remained relatively stable, concentrated between 6 and 12%. Notably, the performance of all cloud classifications was most significantly impacted by deep inland water for the all-day observations.
When distinguishing between daytime and nighttime periods, the clear sky hit rate over the ocean and land was higher during the night. In contrast, for surface types characterized by inland water, the algorithm exhibited more coincident clear sky classification during the daytime compared to the night. Notably, the hit rate for clear sky observations was similar over land for both day and night. Regarding ephemeral (intermittent) water, the false alarm rate for clear sky classifications increased at night. The algorithm consistently had a higher partly cloudy hit rate during the day for all surface types, a trend that is paradoxically accompanied by an increased false alarm rate. The hit rate for overcast classifications was more consistent during the day, with the exception of shallow inland water and ephemeral (intermittent) water. As for the overcast false alarm rate, excluding inland water, it increases during the day for the ocean and land.
The above results indicate that the AIRCD algorithm has more consistency in clear sky classifications over ocean and land, with an increased hit rate and false alarm rate noted during the nighttime. For partly cloudy classifications, both the hit rate and false alarm rate rise during the daytime. In contrast, the overcast classifications show a higher hit rate during the day, coupled with an increased false alarm rate over the ocean and a decreased rate over land. However, there is a decline in the classification hit rate when the surface is inland water. This decrease could be attributed to the FOVs over such surfaces, including shallow, ephemeral (intermittent), and deep inland water, comprising merely 0.86% of the total samples, which highlights the unreliability of the statistical results. In addition, the radiation of the inland water body is affected by various factors such as chlorophyll concentrations, suspended matter, and colored substances in the water [
38]. For the infrared spectrum, pigments and suspended matter can affect the temperature of the inland water by absorbing solar energy, possibly causing a slight increase in temperature, which indirectly affects the observation of infrared radiation. Another reason is the uneven underlying surface within the GIIRS observation. Due to the relatively smaller area of the inland water body, there may also be land surface besides water within one FOV. These factors collectively increase the complexity of cloud detection and classification over inland water. Consequently, it is necessary to perform separate cloud detection calculations for surface types categorized as inland water and discuss them separately to enhance the algorithm’s identification and classification performance under these specific circumstances.
The AIRCD algorithm’s classification performance was evaluated for the GIIRS’s L2 CLM operational products through cross-comparison. By utilizing the spatial matching method previously mentioned, the AGRI’s CLM pixels were matched to each FOV of the GIIRS; cloud labels served as the reference value for comparison. To assess the cloud detection results in May 2022, four evaluation methods were employed using all the FOVs from ten complete observations each day. The daily statistics accounted for the following conditions: FOVs where both the AIRCD algorithm and GIIRS L2 CLM had inconsistent classifications (termed as ‘both missed’) and consistent classifications (termed as ‘both hit’) as a percentage of the total samples. Furthermore, FOVs where the AIRCD algorithm performed consistent classifications but the GIIRS L2 CLM did not (termed as ‘Only AIRCD hit’) and vice versa (termed as ‘Only L2 CLM hit’) were also considered as a percentage of the total samples.
Figure 7 illustrates that the AIRCD algorithm and L2 CLM have a 40–60% consistency, with both missed rates hovering around 16–18%. Notably, the peak and valley of ‘both hit’ occurred on May 9th and May 6th, respectively. The cloud coverage (percentage of cloud FOVs to total samples after AGRI spatial matching to GIIRS) was 63.3% and 50.9%, respectively. These values correspond with the above-mentioned example, where the classification performance of GIIRS L2 was higher when the cloud coverage was higher. This indicates an over-identification of clouds in the L2 operational products. The AIRCD algorithm demonstrates a recognition performance approximately 10% higher than that of L2 CLM. The increased recognition performance highlights the consistency of identification under clear sky conditions, which enhances the utilization of clear-sky data. Furthermore, it also points to the consistency of identification under partly cloudy conditions. The algorithm’s consistency with the partly cloudy labels accounted for 40.6% of all partly cloudy labels, which subsequently, through the use of the cloud cleaning algorithm, can be used to calculate the equivalent clear sky radiation in partly cloudy areas to extract clear sky areas.