An Algorithm to Retrieve Total Precipitable Water Vapor in the Atmosphere from FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) Data

: The atmosphere has substantial e ﬀ ects on optical remote sensing imagery of the Earth’s surface from space. These e ﬀ ects come through the functioning of atmospheric particles on the radiometric transfer from the Earth’s surface through the atmosphere to the sensor in space. Precipitable water vapor (PWV), CO 2 , ozone, and aerosol in the atmosphere are very important among the particles through their functioning. This study presented an algorithm to retrieve total PWV from the Chinese second-generation polar-orbiting meteorological satellite FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) data, which have three near-infrared (NIR) water vapor absorbing channels, i.e., channel 16, 17, and 18. The algorithm was improved from the radiance ratio technique initially developed for Moderate-Resolution Imaging Spectroradiometer (MODIS) data. MODTRAN 5 was used to simulate the process of radiant transfer from the ground surfaces to the sensor at various atmospheric conditions for estimation of the coe ﬃ cients of ratio technique, which was achieved through statistical regression analysis between the simulated radiance and transmittance values for FY-3D MERSI-2 NIR channels. The algorithm was then constructed as a linear combination of the three-water vapor absorbing channels of FY-3D MERSI-2. Measurements from two ground-based reference datasets were used to validate the algorithm: the sun photometer measurements of Aerosol Robotic Network (AERONET) and the microwave radiometer measurements of Energy’s Atmospheric Radiation Measurement Program (ARMP). The validation results showed that the algorithm performs very well when compared with the ground-based reference datasets. The estimated PWV values come with root mean square error (RMSE) of 0.28 g / cm 2 for the ARMP and 0.26 g / cm 2 for the AERONET datasets, with bias of 0.072 g / cm 2 and 0.096 g / cm 2 for the two reference datasets, respectively. The accuracy of the proposed algorithm revealed a better consistency with ground-based reference datasets. Thus, the proposed algorithm could be used as an alternative to retrieve PWV from FY-3D MERSI-2 data for various remote sensing applications such as agricultural monitoring, climate change, hydrologic cycle, and so on at various regional and global scales.


Introduction
The total precipitable water vapor (PWV) in the atmosphere is an important variable to determine the dynamics of the atmospheric movements. PWV is also required for quantitative remote sensing, such as land surface temperature (LST) retrieval from thermal infrared data [1,2]. This is because the radiance transferring from the ground to the sensor in space is strongly affected by the atmosphere, in which PWV plays a critical role in the atmosphere even though it only contains a small proportion in comparison to the total atmospheric mass. Generally, there are two ways through which PWV affects the satellite remote sensing of the Earth's ground surface: First, water vapor is an important absorber to the radiance at several spectral ranges for remote sensing. Second, the emittance of the atmosphere in the thermal infrared range is closely related to the content of PWV in that the atmosphere has at the specific moment of the satellite overpass. Therefore, retrieval of PWV has been an important study area for quantitative remote sensing [3][4][5][6][7]. Accurate information of PWV in the atmosphere can help to improve the understanding of weather patterns, climate change, hydrological cycles, atmospheric dynamics and the chemical composition of aerosols [8,9].
There are several ways to estimate PWV from ground-based stations such as radiosonde observations, upward-looking microwave radiometers, sun photometer, and so on. Lower spatial resolution (i.e., point features) of ground-based PWV measurements limits its application over regional and global scales. To overcome this inconvenience, several remote sensing satellites such as the Moderate-Resolution Imaging Spectroradiometer (MODIS) [10,11], Advanced Very High-Resolution Radiometer (AVHRR) [12] and FengYun (FY-3C and FY-3D) have been launched in recent decades to monitor Earth's processes. Using the data from the above-mentioned remote sensing satellites, we can estimate the PWV at both regional and global scales for various remote sensing applications. Several studies have been published on PWV retrieval from MODIS data over land, cloud, and oceans [13][14][15].
Since launch, the Aqua and Terra instruments carried by MODIS have successfully been in operation in an orbit for more than 18 and 21 years, respectively, which were originally designed for 6 years. Even so, both instruments of MODIS and their onboard calibration devices have been functional up-to-now and can be expected to continue operation for the near future. However, the remote sensing community always look forward to alternative data resources to enhance the earth observations at better temporal and spatial scales. Chinese second-generation polar-orbiting meteorological satellite FY-3D onboard with MERSI-2 can be an alternative data source for global observations since it is already in-service to provide data since November 2018. FY-3D MERSI-2 is equipped with three near-infrared (NIR) water vapor absorption channels (channel 16, 17, and 18) and one neighboring non-absorption channel 4, which confirms the basic requirements for PWV estimation. The spatial resolution of these channels is 250 m for channel 4 and 1000 m for channel 16, 17, and 18. However, very few studies are published on PWV estimation from FY-3D MERSI-2 data.
There are several methods to calculate PWV from satellite data, such as ratio technique, regression slope, split-window difference of the thermal bands, look-up table derived from atmospheric radiative transfer models and the HITRAN2000 spectroscopic database [1,11,13,14,16]. As stated in Sobrino et al. [14], the methods for PWV estimation from remote sensing data should meet the following three essential requirements: first, it can minimize the sensibility among various types of the ground surfaces; second, it can minimize the sensibility to the noise resulted from the variations of different channels; third, it can minimize the sensibility caused by variations in other components of the atmosphere [14]. Considering these requirements, this study utilized the ratio technique to develop an algorithm for PWV estimation from FY-3D MERSI-2 data. The differential absorption method assumes that transmission in a single water vapor absorption channel can be calculated considering the radiance values of both absorbing and non-absorbing channels.
Similarly, for the development of this algorithm, the ratio of three NIR channels and the neighboring non-absorbing channel were used to estimate atmospheric PWV. The work of Sobrino et al. [14] successfully showed the application of the algorithm utilizing such ratio technique with MODTRAN code and the reflectance data of different ground surfaces for the MODIS satellite and archived accuracy with root mean square error (RMSE) of 0.48 g/cm 2 [17]. The detailed description of the ratio technique and its applications have been validated and described in several studies [17][18][19].
Based on abovementioned foundation, the main objective of this study was to develop an algorithm to retrieve PWV from FY-3D MERSI-2 data. Since FY-3D MERSI-2 has slightly similar NIR water vapor absorbing bands and a neighboring non-absorbing band to those of MODIS for atmospheric water vapor monitoring, we intended to propose an algorithm similar to that for MODIS to retrieve PWV from FY-3D MERSI-2 data. In order to have a logical presentation of the study, we structured the remaining part of the paper as follows: Section 2 outlines the data and methodology; Section 3 describes the results and analysis, especially the validation with two ground-based reference datasets; discussion and conclusion are given in Sections 4 and 5, respectively.

The FY-3D MERSI-2 Data
FY-3D is the fourth satellite in the FY-3 series, which was launched in November 2017 and began to provide data in November 2018. The satellite was designed to operate as a primary afternoon orbiting over an altitude of 836 km above the Earth's surface. Its orbital inclination angle is 98.75 • to the equator. FY-3D completes 14 orbital observations of the Earth's surface at global scales twice a day. The MERSI-2 instrument onboard with FY-3D has been greatly improved from MERSI-1, with a higher accuracy of onboard and lunar-calibration capabilities. MERSI-2 is equipped with 25 channels: 6 visible, 3 shortwave IR, 10 visible/NIR, and 6 thermal IR channels with the scanning range of +55.4 • .
The spectral range of MERSI-2 is between 0.412 µm and 12.0 µm, and its spatial resolution at nadir is 250 m for channels 1-4 and 24-25, and 1000 m for the remaining channels. The capability of daily global observation makes it a very high potential satellite with many applications at both regional and global scales, such as climate change, numerical weather prediction, space weather prediction, and ecosystem monitoring. An algorithm to accurately and quickly retrieve PWV from MERSI-2 data is beneficial to the above-mentioned applications.   Table 1 that the spectral range and center wavelengths of both sensors are slightly similar, which implies that it may be possible to develop an algorithm similar to MODIS for FY-3D MERSI-2 by using its spectral information in MODTRAN. The Level-1B radiance data of FY-3D MERSI-2 was used in this study, which were downloaded via the Fengyun Satellite Data Archiving and Order Portal (http://satellite.nsmc.org.cn) over the period of August-December 2019.

Theoretical Basis for PWV Algorithm Development
In order to obtain PWV, the solar radiation reflected by the surface, the absorption and scattering characteristics of the atmosphere and the surface around 1 µm region of wavelength must be taken into account, because absorption and scattering vary with wavelength and viewing angle. The total radiance at a downward looking sensor geometry can be written in a simplified form, as Fraser and Kaufman [20] and Hansen and Travis [21]: where λ is the wavenumber/wavelength, L Sensor (λ) is total radiance reached at sensor, L Path (λ) is the path scattered radiance, L Sun (λ) is above atmosphere solar radiance, ρ(λ) is bi-directional surface reflectance and τ(λ) is combined/total transmittance, considered as the total atmospheric transmittance resulted from sun to surface and surface to sensor [11].

The Ratio Technique for PWV Estimation
The ratio technique is simply detecting the attenuation of the reflected solar radiation which is affected by the water vapor absorption by passing through the atmosphere. A number of atmospheric elements, i.e., carbon monoxide, carbon dioxide, ozone, nitrogen oxide, methane and other gases, are relatively un-changeable in the atmosphere; such atmospheric elements are controlled by default standard atmospheric profiles of MODTRAN 5. However, PWV is the most important element in the atmosphere. Hence, the changing of atmospheric transmittance directly related to changing of PWV in the profile. So, Equation (1) can be written as: In Equation (2), the L path (λ) is directly reflected solar radiation which is usually a few percent because the Rayleigh scattering is negligible near the 1µm region. So, Equation (2) can be further simplified as: where L sensor (λ) is a Level-1B sensor radiance of MERSI-2 and ρ(λ) is a surface reflectance. The total transmittance τ(λ) can be calculated if the L sun (λ) is known. Gao and Kaufman [22], Kaufman and Gao [11] and Sobrino et al. [14] used the non-absorption/window channel to replace L sun (λ) ρ(λ), which makes the ratio technique the only possible way to get the total atmospheric transmittance τ(λ) of absorption bands. In this way the PWV can be computed by building a linear relationship between simulated radiance ratios and PWV with the help of radiative transfer code MODTRAN 5. It is clearly observed from Figure 2 that the MERSI-2 channel 4 at the wavelength range of 0.865 µm is a window channel, and the channels at wavelength range of 0.935, 0.94, and 0.905 µm are water absorption channels.

Algorithm Development for PWV Estimation
In the ratio technique it was supposed that the surface reflectance differs almost linearly with wavelength for various ground surface types around the water absorption spectra. Therefore, the reflectance values of 10 surface types such as forest, snow cover, desert, farm, old grass, cloud deck, ocean, maple leaf, burnt grass, and decayed grass were used as the input parameters in MODTRAN 5. The range of PWV from 0.3 to 3.5 g/cm 2 with a step size of 0.3 g/cm 2 was set as the input for the simulation, the water vapor value 0.3 g/cm 2 reveals drier atmospheres and 3.5 g/cm 2 is for wet atmospheres [14]. The radiances of NIR channels are mainly affected by volume of PWV, so the remaining parameters are set as the default in six standard atmospheres of Earth, including mid-latitude summer and winter, sub-Arctic summer and winter, tropical, and 1976 U.S. standard in MODTRAN. The R 16 , R 17 and R 18 are the ratios of channel 16, 17 and 18 as follows: where L i are the simulated radiances for MERSI-2 channel 4, 16, 17, and 18 using six standard atmospheres and 10 ground surface types. After the 960 simulations of MODTRAN 5, a look-up table was generated, and then the polynomial expressions ( Figure 3) were developed with the help of statistical regression analysis as follows:  As the PWV absorption is a function of wavelength, it can vary between each of the MERSI-2 NIR channels. Considering this issue, the sensitivity of each NIR channel under similar atmospheric conditions could be different. For example, the channel 17 is most sensitive under dry conditions due to its strong absorption, while channel 16 has weak absorption and it is most sensitive under humid conditions [11]. Hence, under the same atmospheric condition, the amount of PWV from the three NIR channels cannot be the same. Therefore, the mean PWV value was obtained according to the following equation: where f 16 , f 17 , and f 18 are weighting functions defined according to: where ∆τ i and ∆W are a difference (maximum and minimum) of transmissivities and water vapor in channel (i) from the six standard atmospheres of Earth [11].
The values obtained for f i are f 16 = 0.208, f 17 = 0.433 and f 18 = 0.359, respectively. Finally, Equation (10) can be written as: Equation (13) is a final equation that we developed to estimate the total PWV from MERSI-2 images.

Framework and Technical Procedures
The framework of the study is described in Figure 4 with the detailed steps of PWV estimation from MERSI-2 data in the development of the algorithm. The primary approach of the study was based on following steps: first, the MODTRAN 5 [23] radiative transfer model to simulate sensor reaching total radiance in absorption and non-absorptions bands for different surface types and atmospheric conditions. The important input parameters were configured to drive MODTRAN 5 shown in Table 2, while other less important parameters, which were not shown in Table 2, were set as default. Second, determination of the coefficients for polynomial expressions using the output information of the MODTRAN 5 by statistical regression analysis was carried out; then, the weight of each band was estimated using the maximum and minimum values of water vapor and atmospheric transmittance. Finally, we performed a validation analysis for the retrieved water vapor to evaluate the applicability and suitability of the algorithm.

Validation of the Algorithm
In order to assure the accuracy of the proposed algorithm, we compared the retrieved PWV amounts with ground-based datasets, using the following assessment criteria: root mean square error (RMSE), bias, and error percentage (EP), generally computed as follows: Two ground-based datasets were used for validation of the algorithm: the upward looking microwave radiometer data provided by the U.S. Department of Energy's Atmospheric Radiation Measurement Program (ARMP) and the sun-photometer data of Aerosol Robotic Network (AERONET). The ARMP's microwave radiometer measurements were conducted at the Southern Great Plains (SGP) in Oklahoma, USA [24]. The microwave radiometer instrument measured the radiation emitted by the atmospheric water vapor and liquid water at frequencies of 23.8 and 31.4 GHz [25]. The microwave radiometer measurement is one of the most precise methods to determine the water vapor from the ground. The uncertainty of the measured water vapor from microwave radiometer data is expected to be in the range of 0.3 mm [26].
The sun photometer measurements of AERONET were used for validation, which were downloaded from http://aeronet.gsfc.nasa.gov. The retrieval produced water vapor values with a consistent dry bias of approximately 5-6% and an estimated uncertainty of 12-15% [27]. The AERONET program is a federation of ground-based remote sensing aerosol networks established by NASA and PHOTONS [28]. The details of sixteen sun photometer sites shown in Table 3 were used for validation, which were randomly selected from the different regions of North America over different types of land covers.

Sensitivity Analysis
In this section, the sensitivity analysis of the MERSI-2 algorithm was performed on the basis of different surface types and atmospheric conditions applied to Equation (13) as follows: where the change in water vapor (∆W i ) of each channel is equal to:  (20) The change in the radiance ratio (∆R i ) was obtained according to Kaufman and Gao [11] as, where σ[R i ] is the standard deviation of radiance ratios (Ri) for all surface covers which are considered in simulation; the denominator Ri (3.5 g/cm 2 ) and Ri (0.3 g/cm 2 ) are the mean value of radiance ratios for wet and dry atmospheres from surface covers considered. To this end, the sensitivity analysis applied to Equation (13) gives a standard deviation of 0.356 g/cm 2 for wet atmospheres to 0.11 g/cm 2 for dry atmospheres.

Precipitable Water from FY-3D NIR Algorithm
The newly developed algorithm was practically applied to retrieve PWV from FY-3D MERSI-2 data. Figure 5 explains an example of PWV retrievals from the MERSI-2 NIR bands. These subsets of images were randomly selected from different lower and higher elevated regions of North America, which include different kinds of land covers. The left panel of Figure 5a,c shows a set of the true color composited images of MERSI-2 bands centered at 0.645 µm (red), 0.86 µm (green), and 0.47 µm (blue). The right panel of Figure 5b,d shows a set of "pixel-based" PWV images processed by the newly developed algorithm. It is easy to demonstrate the spatial distribution of PWV; the most remarkable movement of water vapor can be observed, especially over coastal and mountainous regions as compared to the lower elevated areas. The spatial patterns of water vapor distributions for every location were drastically different.

Validation and Error Analysis
For the better validation of the newly developed algorithm, the FY-3D images between August and December 2019 were used, which mostly cover the sixteen AERONET and single microwave radiometer (MWR) sites referred to in Section 2.6. The programming script was developed to read the data from images by using latitude/longitude as a control condition, matching with the observation data. If the value of PWV extracted by the proposed algorithm is comparable to the PWV obtained from the ground station, we set that pixel as a center pixel, and get the values for positioning pixel 3 × 3 matrix. It is worth mentioning that in some cases, the comparison with ground-based observational data points far from the estimated PWV values. However, this situation was found near cloudy areas, and for this reason, each image was deeply studied because there is not any available cloud mask product to detect the cloudy pixels for better validation. The algorithm proposed in this study is only applicable over clear land surface pixels, except the clouds and oceans with sun-glint. Far from the clouds, 1206 data points for sun-photometer and 103 for microwave radiometer were obtained from 410 satellite images. The validation results are shown in Figures 6 and 7. In each figure, the frequency distribution of error, RMSE, offset and slope of the linear regression, bias, correlation coefficient, and the sample size as the number of data points are shown.  In the left panel of Figure 6, the comparison of the retrieved PWV from the proposed algorithm and the ground-based microwave radiometer measurements is shown; only 100% valid pixels of FY-3D MERSI-2 were considered, to exclude the influence of clouds. It was observed that the microwave radiometer dataset shows almost perfect agreement, even the sample size is relatively low due to cloud contamination. All 103 matched samples were presented in this section based on scatterplots and histogram, which can clearly depict the difference in the distribution of the PWV values. We find that the scatter fitted line is closer to the 1:1 line when the water vapor values are between 1.5 and 2.5 g/cm 2 . On the other hand, the deviations of the PWV become very large when the values of PWV are higher than 2 g/cm 2 . The validation results with microwave radiometer measurements ( Figure 6) produce an RMSE of 0.28 g/cm 2 (12.6%), a bias of 0.072 g/cm 2 , and a correlation coefficient of 0.95.
The comparison between retrieved PWV and sun photometer observational values are shown in Figure 7. Here, again, only 100% valid pixels of FY-3D MERSI-2 were considered. Moreover, there was a large number of available samples which enhanced the validation analysis. Compared to the results of the microwave radiometer in Figure 6, the scatters in Figure 7 are much more focused around the 1:1 line, and the frequency exhibits a much narrower and higher distribution. Furthermore, the deviations in retrieved PWV become very large when the observational values are higher than 1.2 g/cm 2 . The validation results produced a very good consistency with sun-photometer measurements with RMSE of 0.26 g/cm 2 , a bias of 0.096 g/cm 2 , and a correlation coefficient of 0.93.

Comparison with MOD05 Water Vapor Product of MODIS
In order to validate the accuracy of the proposed algorithm, the pixel-based water vapor comparison was performed with the most widely used water vapor product (MOD05) of MODIS satellite. Figure 8 shows the inter-comparison results of PWV between the MOD05 NIR product and the newly developed algorithm of FY-3D MERSI-2. The data acquired for comparison were investigated over North America on 21 August 2019, at UTC 20:45 for FY-3D MERSI-2 and UTC 18:25 for MODIS. Both sensors have almost similar spatial resolution, but their image acquisition time was slightly different. In the upper panel of Figure 8, the retrieved PW data of FY-3D MERSI-2 are plotted against the MODIS product, which reveals significant differences. The comparison shows a larger scatter and wider histogram in the data against the MOD05 PWV product, as evidenced by the lower correlation coefficients 0.58 and increased RMSE up to 0.52 g/cm 2 . A bias of −0.112 g/cm 2 was found, which produces large RMSEs over the North American conditions.

Comparison with Other Combinations
The FY-3D MERSI-2 satellite equipped with three NIR water vapor channels is discussed in Section 2.1. Equation (13) was made with a combination of three NIR channels, and the retrieved water vapors were compared with two kinds of in-situ ground datasets to see how well the algorithm would perform. In this section, we further investigate the accuracy of three other combinations of NIR channels, such as:  (24) Similar AERONET ground-based stations were used in this section, which were initially investigated for the validation of the three-channel NIR algorithm. The comparisons of three different combinations of NIR channels such as band 16 with 17, band 17 with 18, and band 16 with 18 were shown in Figure 9. The RMSE of 0.38, 0.28 and 0.34 g/cm 2 was found for band 16 with 17, band 17 with 18, and band 16 with 18, respectively, which are equated to be in the range between 28 and 64% of error. The line representing a 1:1 correspondence is included in each plot to illustrate slopes or offsets which may exist between the data in the comparisons, and the histogram shows how data are frequently distributed.

Discussion
This study presents an applicable algorithm to estimate PWV from FY-3D MERSI-2 data over clear sky land pixels. The suitability of the proposed algorithm in terms of the error percentage, compared to different ground-based stations and MODIS water vapor product, was investigated. The results showed that the agreement between MERSI-2 and in-situ ARMP and AERONET measurements was acceptable, with an RMSE of 0.28 g/cm 2 (12.6%) and 0.26 g/cm 2 (23.7%), respectively. These findings were comparable to the accuracy obtained for the MODIS institutional product, MOD05, reported by Gao and Kaufman [22] with a 5-10% error range. Similarly, Sobrino et al. [14] reported~21% accuracy from a validation study conducted for the MODIS NIR algorithm using the ratio technique. A bias of 0.09 g/cm 2 , confirmed by Gao and Kaufman [22] and Sobrino et al. [14] for MODIS PWV validations, was also comparable to our proposed method's biases (0.072 to 0.097 g/cm 2 ). During the pixel-based comparison between the MOD05 water vapor product and FY-3D MERSI-2 algorithm, the RMSE reached up to 0.52 g/cm 2 , which was relatively higher than other comparisons and far from the acceptable limits. The number of possible reasons can explain this larger RMSE: First, the three-channel NIR algorithm used in this study and the MOD05 NIR algorithm differ slightly. The MODIS water vapor product was developed by ratios of apparent reflectance in the NIR channels, rather than ratios of measured radiance as utilized in this study. Second, the MOD05 NIR algorithm was developed by consideration of two non-absorbing channels (0.865 and 1.240) rather than MERSI-2, which utilized only a single non-absorbing channel [17]. Third, both instruments, MERSI-2 and MODIS, look at the surface from different angles, and their spectral response functions are slightly different. These differences lead to slightly different results from each of the respective algorithms which do not provide the best inter-comparison results between the MERSI-2 and MODIS. It was observed that the three-channel algorithm is the best way to calculate PWV with acceptable limits from FY-3D MERIS-2 data. This is also because under a similar atmosphere, every NIR band has different water vapor sensitivities. Another finding of this study suggested that the strong absorption band 17 is most sensitive in a moisture deficient environment, while the weak absorption band 16 is most sensitive in humid environmental conditions [11]. Furthermore, there were numerous sources of errors for PWV retrieval from NIR bands that were reported in the literature by Kaufman and Gao [11]; Gao et al. [10]; Bouffies et al. [29]; Gao et al. [30]. Along with the uncertainties in the spectral reflectance of the surface, these errors include spectral calibrations and sensor radiometric, pixel registration and so on. Through the sensitivity analysis of FY-3D data, it was revealed that the amount of deviations depends on the surface targets in different atmospheres.

Conclusions
FY-3D is a second generation polar-orbiting meteorological satellite of China. FY-3D has provided various products to the remote sensing community. However, estimations of water vapor from the FY-3D MERSI-2 data attained were less focused. This study provides an algorithm to retrieve PWV from FY-3D MERSI-2 data over clear sky land surface pixels. During the development of the ratio technique, the simulated radiance values of three water vapor absorbing channels and one non-absorbing channel were used. The simulation of radiance from the ground surface to the sensor in space was investigated through the MODTRAN 5 radiative transfer model. From the output information of MODTRAN, a look-up table was generated and then polynomial expression were developed with the help of regression analysis. Several regression equations were constructed based on various combinations of the ratios against the PWV estimation, and the best one was then selected as the possible algorithm as W = 0.208W 16 + 0.433W 17 + 0.359W 18 , where W 16 , W 17 , and W 18 are polynomial expressions of PWV estimation for FY-3D MERSI-2 bands 16, 17, and 18, respectively.
Additionally, two ground-based datasets, ARMP and AERONET, were used in the study to validate the proposed algorithm. The results of the validation showed that the proposed algorithm has acceptable accuracy. The estimated PWV values come with an RMSE of 0.28 g/cm 2 (12.6%) for ARMP measurements and 0.26 g/cm 2 (23.7%) for AERONET measurements. Moreover, the proposed algorithm is very easy to apply and reproducible. It only requires the Level 1B radiance data of NIR bands from FY-3D MERSI-2 for PWV retrieval. In this sense, the proposed algorithm can be used as an alternative to retrieve PWV from FY-3D MERSI-2 data at global scales for various remote sensing applications. The future line of work is to develop an advanced method to retrieve PWV over land, oceans and clouds with more detailed validation analysis throughout the Earth.