Abstract

FY-3A is the second Chinese Polar Orbital Meteorological Satellite with global, three-dimensional, quantitative, and multispectral capabilities. Its missions include monitoring global disasters and environment changes. This study describes some basic parameters and major technical indicators of the FY-3A and evaluates data quality and drought monitoring capability of the Medium-Resolution Imager (MERSI) onboard the FY-3A. Data obtained with the MERSI was compared with that of the MODerate-resolution Imaging Spectroradiometer (MODIS), imaged at the same time period and geographic zone. In addition, the Temperature/Vegetation Drought Index (TVDI), a highly accurate and stable monitoring model, was used to monitor drought condition with MERSI and MODIS sensors. It is found in the study that the relative accuracy of data, obtained with these two devices, was consistent with the acceptable overall accuracy of 93.8. Furthermore, spatial resolution of MERSI is superior as compared to that of MODIS. Therefore, FY-3A MERSI can serve a reliable and new data source for drought monitoring.

1. Introduction

1.1. Drought Monitoring

Drought is a normal, recurrent feature of climate. It occurs almost everywhere, although its features vary from region to region. In the most general sense, drought originates from a deficiency of precipitation over an extended period of time, resulting in a water shortage for some activity, group, or environmental sector (NDMC define). Consequently, it is required to demonstrate the distribution and degree of drought condition timely, which is crucial for drought warning and resisting effectively [13].

At present, remote sensing methods for drought monitoring are mainly classified into four categories: Vegetation Index-based [4, 5], Temperature-based [6], Vegetation and Temperature-based [7], and Cloud-based [8]. The representative indicies include Vegetation Supply Water Index (VSWI) [9],Temperature/Vegetation Dryness Index (TVDI) [10, 11], Apparent Thermal Inertia Index (ATI) [12, 13], and Cloud Parameters Index (CPI) [14]. As FY-3A satellite is a polar-orbital, its temporal resolution is hard to meet the requirements of cloud-based method, so CPI is not utilized in this study. ATI requires the index of temperature variation between day and night. However, polar-orbital satellites cannot provide real-time temperature; consequently ATI method is not suitable for FY-3A to monitor drought. VSWI and TVDI need only two parameters—Normal Difference Vegetation Index (NDVI) and Land Surface Temperature (LST), which can be calculated from MERSI data. Meanwhile, taking the superiority of TVDI into account, it is be used in this study.

1.2. FY-3A MERSI

The FY-3A satellite is the second Chinese polar orbital meteorological satellite which provides global, three-dimensional, high-quantified, and multispectral remote sensed data. The satellite weighs 2400 kg, and has in-orbit dimensions of meters. Its orbit altitude is 836.4 km, with an obliquity of 98.753 degrees and an orbit time of 101.496 minutes. Its designed life-span is over 2 years. FY-3A has refined the sensor technology, which is multisensor, microwave, hundred-meter resolution, and global data reception rather than single-sensor, optical only, kilo-meter resolution, and domestic data reception in the field of meteorological satellite. FY-3A missions include monitoring natural disasters and environmental changes all over the world and offering meteorological information for agriculture. Since, drought monitoring has naturally become one of the important applications of FY-3A. Therefore, drought monitoring capability of FY-3A deserves a comprehensive evaluation. Using MODIS results obtained at the same time period as reference, data quality, and drought monitoring capability of MERSI is analyzed in the study (CMA).

With multispectral and high-resolution features, MERSI can detect atmosphere, land, and ocean through the reflection of visible channels and thermal infrared radiation. There are three channels located in the water vapor absorption window (0.9050.980 um) which increase detection capability for atmospheric water vapor, in particular water vapor at low levels. Visible bands with spatial resolution of 250 meters can be utilized to derive NDVI and thermal bands for LST. Hence, MERSI has an enormous potential for drought monitoring. Details about some spectral properties of MERSI are listed in Table 1.

It can be found in Table 1 that FY-3A MERSI has red, near-infrared, and thermal infrared channels with spatial resolution of 250 meters, so it can provide drought monitoring information of 250 meters. This should be a great improvement of meteorological satellite data in spatial resolution for drought monitoring application.

Both FY-3A and TERRA are polar-orbiting satellites. Channels setting of MERSI, including the center wavelength and the wave width, is basically consistent with that of MODIS, in particular the bands used for monitoring surface vegetation, atmospheric water vapor, and surface temperature. MODIS data is calibrated on orbit and it uses the complicated recorrecting technology to locate when it scans. Because of high-quality and effective monitoring, MODIS has become a widely used data source in drought monitoring. Therefore, MODIS is optimal reference data in analyzing and evaluating radiometric calibration, relative geometric location, and drought monitoring results of MERSI.

In order to enlarge the remote sensing data source, the objective of the present study is to evaluate the data quality and the drought monitoring capability of FY-3A MERSI data with TERRA MODIS data as a reference.

2. Materials and Methods

2.1. Study Area

The case study area is near Bohai Bay in North China (area inside the red frame, Figure 1) with latitude range from to and longitude range from to respectively. There are a variety of land covers in the area, including vegetation, bare area, and waters.

2.2. MERSI and MODIS Data

The data of the study were extracted from MERSI on October 2008, at 03:00 AM (GMT) (Figure 2(a)), and from TERRA-MODIS of the same area 15 minutes earlier (Figure 2(b)). The image of MERSI was 2000×2084 pixels, of TERRA-MODIS 1354×2030 pixels. To make it comparable with the result of MODIS data, FY-3A MERSI data was resampled to 1000-meter spatial resolution. These two kinds of data were provided by China Meteorological Administration (CMA) and MODIS Data Receiving Station in Wuhan University, respectively.

The orange lines in Figure 2 denote orbit central lines. Their spatial resolution has achieved the anticipated design requirement and the texture information is plentiful. The study area is in the lower right corner of MERSI data (Figure 2(a)) and at the upper left side of MODIS data (Figure 2(b)).

2.3. Methods
2.3.1. NDVI

The Normalized Difference Vegetation Index is a satellite-derived global vegetation indicator obtained from the red, near-infrared (NIR) ratio of vegetation reflectance in the electromagnetic spectrum [15, 16]

The parameters, and , are the reflectivity of near-infrared and red channels, respectively.

NDVI provides information on vegetation productivity and phenology over large temporal and spatial scales and has been widely used in the recent ecological studies as a proxy for vegetation productivity and phenology [17].

2.3.2. LST

The Land Surface Temperature can be calculated from the Brightness Temperature (BT) of the thermal infrared channels based on the thermal transmission equation. Because FY-3A MERSI contains only one thermal infrared channel, LST is replaced by BT of the thermal infrared channel in this paper [18].

2.3.3. TVDI

Sandholt and others extracted water stress indices (i.e., the index of temperature vegetation drought) on the basis of the simplified NDVI-LST feature space. (minimum LST) is a straight-line paralleling with NDVI axis, and (maximum LST) has a linear relation with NDVI. The TVDI can be obtained with formula (2) on the basis of NDVI-LST feature space [19, 20]:

In fact, in the NDVI-LST feature space varies under different vegetation coverage conditions [2123]. Formula (2) is obtained by calculating the four Vertex Coordinates of trapezium theoretically based on the assumption that NDVI-LST feature space is trapezoidal. Therefore, and can be obtained through Least-squares linear fitting when NDVI-LST feature space is simplified into a triangle. In traditional studies, the function between and NDVI is defined as wet side equation while that of is accordingly called dry side equation. The formulas are as follows:

where and are the coefficients of the dry side equation, and and are the coefficients of the wet side equation, respectively.

In the experiment, a serial points (e.g., take ), () is extracted, where is the serial number of these point. Thus, the error matrix equation can be established

where is the error matrix, and is the number of points.

The normal equation can be deduced from the error matrix equation based on the Least-squares algorithm

where is the transposed operation of matrix.

The coefficients can be calculated based on the matrix operation:

where –1 is the inversion operation of matrix.

Last, the relation square (-square, ) demonstrates that this statistic measures how successful the fit is in explaining the variation of the data. A value closer to 1 indicates a better fit

where is the average value of all and is the average of all NDVI.

3. Result and Discussion

After calibration, MERSI and MODIS data are projected with the latitude and longitude information recorded in scan. The relative accuracy of geometric location and radiometric calibration of the two kinds of data are also analyzed. NDVI is obtained through red and near infrared channels, and LST in the clear sky or cloud top temperature can be retrieved through thermal infrared and vapor channels. Furthermore, it is necessary to compare and analyze the distribution character of the dry and wet sides of the two images in the NDVI-LST feature space. User’s accuracy, producer’s accuracy, and the overall accuracy are evaluated for TVDI with MERSI and MODIS data.

3.1. Data Preprocessing
3.1.1. Calibration

The visible channels of MERIS are directly calibrated with slope-intercept form, and revised solar altitude and Earth-Sun Distance at the imaging moment. Digital number (DN) of infrared channels is already radiation value, so the brightness temperature can be gained with Planck Formula.

The visible data of MODIS is calibrated with slope-intercept form, while the infrared channel is calibration with Planck Formula to get brightness temperature directly [2426]. The in-lab calibration accuracy of MERSI visible and infrared channels is 7% and 1 K, respectively, but the radiation calibration accuracy of MODIS data is much higher.

Both record the latitude and longitude values in imaging; they are translated into geographic coordinates in the corresponding projection through strict transformation calculation method [27, 28]. At the same time, gray value of a pixel is extracted to a corresponding location in the new image. The new projection is Lambert Equal-angle Intercross Taper Projection (Lambert Conformal Conic). Its specific parameters include Spheroid of Krasovsky, Datum of Krasovsky, the first standard parallel latitude of the second standard parallel latitude of the central meridian latitude of the central meridian longitude of and the origin coordinates of (5000000.0, 5000000.0) (measured in meters). The results of radiometric calibration and projection of the two data are shown in Figure 3.

It is obvious in Figure 3 that the two images have the similar color effect and the color of MERSI is slightly lighter than that of MODIS. Moreover, there is small difference in atmospheric conditions, for example, the sizes of cloud cluster are a little different at the Bohai Bay region. For region B in Figures 3(a) and 3(b), it can be easy noticed that the resolution of MERSI data, located at the right side of original images, is coarser than that of MODIS which locates at the center of the original scene. Furthermore, the partial textures information of region B in MERSI data is poorer than that of MODIS. For region A in Figures 3(a) and 3(b), it is vice versa. This is mainly due to the different positions of the study area in the original images. Some statistic parameters can be seen in Table 2, and the gray histograms of selected channels are shown in Figure 4.

The population characteristic value different of red channel reflectivity between MERSI and MODIS data is 0.05 and the population variance is 0.02. In other words, red channels of two sensors have same sensitivity to different types of ground objects. Of course, those asynchronous images and the differences of central wavelength and wave width will also have some effect on the scaling result, and this is beyond the scope of this study. The similar situation also happens to the reflectivity of near infrared channel. Its population characteristic value difference between the MERSI and MODIS data is 0.06 and the population variance is 0.02. As can be seen from the reflectivity histogram of near infrared channel of MERSI and MODIS data, the histogram has two-peak values. According to the reflection spectral characteristics of ground objects, the lower peak value of reflectivity mainly results from the lower reflectivity of the surface waters, as most of the study area is covered by waters. Furthermore, LST in the clear sky or cloud top temperatures retrieved from two images are of little difference, and this indicates that MERSI has a strong capability for LST retrieving.

3.1.2. Relative Geometric Location Analysis

MODIS data is of high orbital location accuracy, so it can serve as reference for analyzing geometric location accuracy of MERSI data. Fifteen control points (balanced distributed in study area) are selected manually in the study. The geometric location differences of these points are used to analyze geometry location accuracy of MERSI data. Figure 5 shows that the spatial distribution of the points and Table 3 shows the geometric location difference of each point between two images.

It is shown in Table 3 that MERSI data has quite good location accuracy when compared with MODIS data. Its overall difference is 0.4666 pixels, with direction difference of 0.415 pixels and direction difference of 0.2131 pixels. The maximum difference point for is the point no. 12 with a small difference of mainly because of manual point-choosing. Actual spatial resolution of the points at the side is coarser than those near nadir, so the points at the left side of MODIS (no. 1, no. 2, no. 8, no. 12, and no. 14) have poor matching with their corresponding points near nadir of MERSI.

Based on the analysis of calibration results, MERSI can retrieve similar reflectivity, LST, and cloud top temperature as MODIS does. The geometric location difference between MERSI and MODIS is less than half pixel, while difference of some areas reaches one pixel, but it is still within the tolerance. In other words, the geolocation of MERSI is acceptable. Therefore, the MERSI data has potential for drought monitoring on basis of calibration and relative geometric location analysis.

3.2. Cloud Detection

NDVI and LST, the two key parameters for TVDI method, cannot be retrieved when remote sensing data is covered by clouds, so TVDI is inapplicable for the cloud-covered areas. Consequently, to identify whether an image is qualified for drought monitoring, cloud detection should be performed before deducing parameters [29]. In addition, areas covered by water need to be identified and ignored in case of error in establishing NDVI-LST feature space.

Cloud in visible band shows high reflectance, and remote sensing image with 0.66 um is ideal for distinguishing borders between land and cloud [30]. The spectrum feature of cloud in NIR band (near 0.936 um) is mainly relevant to water vapor content in atmosphere, so it can help to indicate the moisture feature, namely, absorption window. Because 0.66 um and 0.936 um clouds form sharp contrast with spectrum features of all kinds of ground objects, cloud information can be obtained by normalization, which not only gives prominence to cloud information and partially eliminate ill effects caused by solar altitude, satellite sweep angle and atmosphere route radiation. Normalized Cloud Detection Index (CDI) [12, 21] is expressed as

As cloud refection spectrum has reflection characteristics at 0.66 um and its absorption window at 0.936 um is influenced by water vapor, the CDI is positive (CDI > 0), because soil reflection spectra has little difference in reflection properties between 0.66 um and 0.936 um, the CDI is close to 0; vegetation reflection spectra has low reflection at 0.66 um and high reflection at 0.936 um, the CDI is negative (CDI < 0).

As the reflectance of water is higher than vegetation near the red band and the reflectance of vegetation is obviously higher than water round the NIR band (0.841 um–0.876 um)   [31, 32], a typical bimodal distribution appears in the histogram in Figures 4(c) and 4(d). Thus, NDVI is adopted to enhance land-water contrast.

Based on the cloud detection and water extraction algorithms discussed above, it can be concluded that MERSI and MODIS data can have cloud and water detection with the same algorithm and the results is shown in Figure 6.

Figure 6 shows that the results of MERSI and MODIS are unanimous in most areas. However, as a result of difference in imaging time, there is a slight difference in the cloud detection at the top left corner of images due to small atmospheric condition change.

3.3. Extraction of Drought Feature Parameters

For MERSI and MODIS, NDVI is usually calculated by red band and NIR band. Meanwhile, MERSI has only one TIR band (1012.5 um with 250 m resolution in nadir point), which can be used to extract LST in the clear sky or cloud top temperature. Qin Zhihao deduced a simple, feasible LST algorithm with high accuracy according to the equation of land surface heat conduction [33, 34] and this method can be quoted to process MERSI data.

From previous studies, it is known that the distribution of the wet and dry sides is significant in the two-dimensional feature space composed by NDVI and LST from MERSI. It is necessary to analyze the characteristics of the NDVI-LST feature space and determine the wet and dry sides.

LST is calculated with the NDVI datasets and single-window algorithm, of different NDVI conditions is extracted on the basis of 0.01 NDVI step (NDVI is among 0.3 to 0.7 for medium vegetation coverage), and then NDVI-LST feature space is obtained (Figure 7).

It is shown in Figure 7 that and have better correlation with NDVI in MERSI than in MODIS, but their linear correlations are quite consistent. This may be due to the difference of cloud detection, that is to say, the cloud top temperature may be considered as the minimum LST in the clear sky in the process of selecting MODIS LST samples. From Figure 7, MODIS has smaller and higher for the same NDVI value. This demonstrates that MODIS data is doing better in retrieving LST in the clear sky than MERSI. Based on the principle of least squares solver fitting, the relationship functions are regressed in the NDVI-LST feature space. Because wet and dry side fitting of the low and high vegetation coverage (with NDVI ranging from 0 to 0.30 and from 0.70 to 1, resp.) is quite instable, only middle NDVI of 0.30 to 0.70 is used in regression.

The fitting equations of wet and dry sides based on MERSI are as follows.

Dry side equation:

Wet side equation:

The above formulas demonstrate that if the slope of dry side is less than zero, will decrease with the increase of vegetation coverage. Oppositely, the minimum LST is vice versa.

The fitting -squares of MERSI wet and dry sides are 0.8307 and 0.9541, respectively. It shows that changes of and are quite consistent with the changes of NDVI, and the mechanism of drought is well demonstrated. The fitting -squares of MODIS are 0.6852 and 0.6904, respectively.

3.4. Drought Monitoring

It is found through the above analysis that the wet and dry sides possess high-fitting relevance in the NDVI-LST feature space of MERSI and it means that MERSI is suitable for drought monitoring. TVDI is used to analyze the drought monitoring capability of FY-3A MERSI in the study. MODIS has been proved to be one of the most successful data for drought monitoring and is naturally utilized as reference data. The monitoring results and comparison of MERSI and MODIS are as shown in Figure 8.

Table 4 is comparison of TVDI confusion matrices between FY-3A MERSI and TERRA MODIS. User’s accuracy, producer’s accuracy, and overall accuracy of MERSI are analyzed with MODIS TVDI as a reference. The overall accuracy contains the correct and acceptable accuracy, where the correct accuracy indicates that the monitoring results of MERSI TVDI and MODIS TVDI are unanimous. The definition of the acceptable accuracy is that the difference of grade results is lower than one level. For instance, it is acceptable if one pixel is shown as severe drought in MERSI TVDI and its corresponding point shown as medium drought in MODIS TVDI. It is clear in Table 4 that MERSI has almost the same drought monitoring result with MODIS. In other words, absolutely correct accuracy is 72.1274 and the acceptable accuracy is as high as 93.8306. Thus, MERSI is as strong as MODIS in drought monitoring.

4. Conclusion

The following conclusions can be reached on the basis of above comparison and analysis.

(1) FY-3A MERSI data enjoys a high quality. The relative difference of calibration results between MERSI and MODIS is negligible, and temperature retrieving capability of MERSI is as good as MODIS. In addition, it is also proved that their cloud detection results have strong consistency. Furthermore, the overall difference of geo-location between MERSI and MODIS is about 0.4666 pixels, in other words, MERSI data possess wonderful geo-location capability.(2)The decreases with the increase of vegetation coverage, vice versa for The fitting R-squares of wet and dry sides extracted from MERSI are as high as 0.8307 and 0.9541, respectively. It means that changes of and have high consistency with those of NDVI, and the mechanism of drought is well presented. While, the fitting -squares of MODIS data (for wet and dry side) are, respectively, 0.6852 and 0.6904 which are lower than that of MERSI. From the NDVI-LST feature space of MERSI, the scatter points’ distribution of wet side is highly correlated with that of dry side, and this is a clear demonstration of drought information.(3)The TVDI confusion matrices indicate that FY-3A MERSI possesses the similar capability to MODIS, with absolute accuracy of 72.1274 and acceptable accuracy of 93.8306. Therefore, it is possible for us to obtain the drought products in operation with 250 m spatial resolution and global scale with MERSI as a new kind of data.

Acknowledgments

The authors are grateful to the members of the MODIS Data Receiving Station, Wuhan University, China. This paper is granted by FY-3A Satellite Application in Drought Monitoring and Early Warning in Northwest China which is subtask of FY-3A Satellite Application Projection from China Meteorological Administration (CMA). Meanwhile, this paper is also granted the AMD University Plan.