Next Article in Journal
Drifting Effects of NOAA Satellites on Long-Term Active Fire Records of Europe
Next Article in Special Issue
Impact of Cyclonic Ocean Eddies on Upper Ocean Thermodynamic Response to Typhoon Soudelor
Previous Article in Journal
Detection and Validation of Tropical Peatland Flaming and Smouldering Using Landsat-8 SWIR and TIRS Bands
Previous Article in Special Issue
Variability of the Great Whirl and Its Impacts on Atmospheric Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Surface Air Specific Humidity and Air–Sea Latent Heat Flux Using FY-3C Microwave Observations

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Key Laboratory of Earth Observation, Hainan Province, Sanya 572029, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(4), 466; https://doi.org/10.3390/rs11040466
Submission received: 21 January 2019 / Revised: 6 February 2019 / Accepted: 20 February 2019 / Published: 24 February 2019

Abstract

:
Latent heat flux (LHF) plays an important role in the global hydrological cycle and is therefore necessary to understand global climate variability. It has been reported that the near-surface specific humidity is a major source of error for satellite-derived LHF. Here, a new empirical model relating multichannel brightness temperatures ( T B ) obtained from the Fengyun-3 (FY-3C) microwave radiometer and sea surface air specific humidity ( Q a ) is proposed. It is based on the relationship between T B , Q a , sea surface temperature (SST), and water vapor scale height. Compared with in situ data, the new satellite-derived Q a and LHF both exhibit better statistical results than previous estimates. For Q a , the bias, root mean square difference (RMSD), and the correlation coefficient (R2) between satellite and buoy in the mid-latitude region are 0.08 g/kg, 1.76 g/kg, and 0.92, respectively. For LHF, the bias, RMSD, and R2 are 2.40 W/m2, 34.24 W/m2, and 0.87, respectively. The satellite-derived Q a are also compared with National Oceanic and Atmospheric Administration (NOAA) Cooperative Institute for Research in Environmental Sciences (CIRES) humidity datasets, with a bias, RMSD, and R2 of 0.02 g/kg, 1.02 g/kg, and 0.98, respectively. The proposed method can also be extended in the future to observations from other space-borne microwave radiometers.

Graphical Abstract

1. Introduction

Air–sea latent heat flux (LHF), combined with momentum flux and sensible heat flux, represents an important aspect of the atmosphere–ocean interaction, energy budget, water cycle variability, and global climate systems [1]. In addition, inland lakes influence greatly on surface energy budgets and have shown considerable changes worldwide [2,3]. This is particularly true in high-latitude regions, where anthropogenic climate change is predicted to be exceptionally rapid [4]. LHF can be estimated by the bulk aerodynamic formula that includes sea surface temperature (SST), surface winds, air temperature ( T a ), and near-surface air specific humidity ( Q a ). These flux-related variables can be obtained from three major sources: In situ measurements, satellite observations, and numerical weather prediction models. Several previous studies have used in situ measurements from ships or buoys to estimate LHF [5,6]. Since the temporal and spatial resolution of ship measurements is limited, large uncertainties are still present in the estimation of LHF. The atmospheric reanalysis products can be strongly influenced by variations in the type and amount of data assimilated [7].
Depending on the data fusion techniques, blended global surface heat flux datasets have been derived from multi-satellite information. These mainly include the Hamburg Ocean–Atmosphere Parameters and Fluxes from Satellite Data (HOAPS-3; [8,9]), Japanese Ocean Flux Data Sets with Use of Remote Sensing Observations (J-OFURO; [10]), Goddard Satellite-Based Surface Turbulent Flux (GSSTF-3; [11]), and French Research Institute for Exploitation of the Sea (IFREMER) turbulent flux products [12]. Compared with in situ and reanalysis flux products, satellite observations provide much higher spatial resolution and a homogeneous time series of atmospheric state variables for surface flux computation. However, a large portion of the errors in satellite-derived LHF is found to be associated with uncertainties in the estimation of Q a [13].
Several authors have investigated the estimation of Q a from satellite measurements. Using the column water vapor (W) as an independent parameter is a common way to determine the Q a from satellite measurements [14,15]. This Q a −W relationship is widely used in LHF-retrieval algorithms. However, errors in Q a originating from the relationship might result in large errors in estimating LHF, especially on hourly and daily scales. Schulz et al. (1993) developed a method to estimate the bottom-layer (500 m height) precipitable water from the brightness temperature [16]. Thereafter, several similar studies based on Schulz’s model have been proposed. For example, Schlüssel et al. (1995) (hereafter SC95) directly regressed the remotely sensed brightness temperature and Q a to avoid error propagation [17]. Later, this approach was updated using improved training data [12,18]. The difference in the frequencies used for Q a retrieval in each algorithm is considered to be a cause of statistical difference between each product. High-frequency channels, 85–89 GHz, had not been used for Q a retrieval in previous studies because the polarized radiation was considered to have no significant information about the boundary layer [16]. Iwasaki et al. (2012) (hereafter IW12) concluded that the use of the brightness temperature determined by the 85 GHz polarized radiation can considerably reduce the Q a bias, particularly in the wet regions [19]. Information on the vertical profile of humidity is also useful for improving the estimation accuracy of Q a [20,21]. Moreover, methods using artificial intelligence such as neural networks and genetic algorithm have also been proposed [22,23,24]. More recently, Tomita et al. (2018) (hereafter TO18) developed a method to estimate the Q a using vertical water profile information and W [25]. However, the TO18 algorithm ignored the effect of SST, which is highly correlated with W, on Q a .
In this study, we developed an improved method to estimate Q a based on the TO18 algorithm. Vertical water profiles and other meteorological parameters (SST and W) were used to improve retrieval accuracy. The remaining parts of this paper are organized as follows: Datasets used in this study are introduced in Section 2. Analysis procedures and retrieval updates of Q a are described in Section 3. The results and validation are presented in Section 4. The discussion and conclusions are included in Section 5 and Section 6, respectively.

2. Data Descriptions

2.1. Satellite Observations

The brightness temperature ( T B ) was obtained by the Microwave Radiation Imager (MWRI) onboard the new generation polar-orbiting meteorological satellite of China (Fengyun-3C, FY-3C), launched in September 2013. The MWRI instrument is a microwave conical-scanning imager following on the heritage of the Special Sensor Microwave Imager/Sounder (SSMI/S), Advanced Microwave Scanning Radiometer for EOS (AMSR-E), and the AMSR-2 instruments. MWRI has vertically and horizontally polarized radiation for 10.65, 18.7, 23.8, 36.5, and 89 GHz [26] (hereafter referred as 10 V/H, 19 V/H, 23 V/H, 37 V/H, and 89 V/H respectively; the suffixes V and H mean vertical and horizontal polarization, respectively). In this study, data from ten common channels were used; the correlation between multiple channels is shown in Table A1. After data preprocessing, the instantaneous T B data from January to December 2014 were used for algorithm development. The main characteristics of the MWRI are shown in Table 1.

2.2. In Situ Measurements

Two kinds of in situ measurements were used as training data to derive the Q a regression formula. One was the International Comprehensive Ocean–Atmosphere Data Set Release 3.0 (ICOADS 3.0) [27]. The other was from mooring buoy networks: The National Data Buoy Center (NDBC) buoys off the U.S. Atlantic, Pacific, and Gulf coasts maintained by the National Oceanic and Atmospheric Administration (NOAA); the European offshore data acquisition system (ODAS) buoys in the eastern Atlantic, maintained by the UK Met Office and Météo-France (MFUK); and the Tropical Atmosphere Ocean (TAO) buoys located in the tropical Pacific Ocean, maintained by NOAA’s Pacific Marine Environment Laboratory (PMEL).
Data from January to December 2014 that included sea surface temperature, relative humidity (dew point temperature), air temperature, and pressure were used in this study. The Q a was computed with the saturation vapor by using a pressure formula [28]. We removed the outlier values of Q a , which were above the upper quartile plus 1.5 × inter-quartile range (IQR) or below the lower quartile minus 1.5 × IQR [29]. After carrying out quality control for Q a , the Coupled Ocean-Atmosphere Response Experiment version 3.0 (COARE3.0) model was used to adjust the Q a data to standard height (10 m) using the metadata on the type of ship and the height of the meteorological sensors [30]. To match up in situ data with the satellite observations, we set criteria that the temporal and spatial differences had to be less than 30 min and 25 km, respectively. Finally, in total 23,200 pairs of ship/satellite collocated data were divided randomly into two equal groups (Samples 1 and 2). Data in Sample 1 were used for the development of the new algorithm, while Sample 2 was used for an additional evaluation of the retrieval algorithms as independent data.

2.3. NOAA CIRES

The instantaneous FY-3C Q a values are compared with NOAA CIRES Multi-Satellite Humidity. This product combines Special Sensor Microwave Imager/Special Sensor Microwave Imager Sounder (SSMI/SSMIS) and Advanced Microwave Sounding Unit (AMSU) satellite observations to create a new satellite humidity dataset. The retrievals were developed using research vessel data from the Shipboard Automated Meteorological and Oceanographic System (SAMOS) initiative and validated using the ICOADS dataset [31]. The NOAA data are available through the Center for Ocean-Atmospheric Prediction Studies (COAPS) Marine Data Center (https://mdc.coaps.fsu.edu/data).

2.4. ERA-Interim Reanalysis Data

The ERA-Interim is a global atmospheric reanalysis produced by the European Centre for Medium-Range Weather Forecasts (ECMWF) [32]. ERA-Interim products cover the period from 1979, and are continuously updated in real time. In this study, humidity, SST, and W were obtained from the ERA-Interim daily and gridded products, which have been used for theoretical research and validation. This data can be download from the website: http://apps.ecmwf.int/datasets/. The time series is from January 2014 to December 2014.

3. Methodology

There are different satellite remote sensing algorithms for air-–sea specific humidity and LHF, but these algorithms do not fully consider the effects of SST, W, and vertical moisture structure on Q a . In this study, a methodological framework was proposed to estimate Q a and LHF from FY-3C observations and meteorological auxiliary data (see Figure 1).

3.1. Relationship between SST and W

Kanemaru and Masunaga (2013) illustrated that the water vapor scale height ( H v ) is a good indicator of the vertical moisture gradient between the boundary layer and the free troposphere,
H v = W ρ a q v ,
where W is the column water vapor (kg/m2), ρ a is the density of dry air (1.2 kg/m3), and q v is the surface water vapor mixing ratio (g/kg) [33]. Here, values of W and q v were derived from ERA-Interim [32].
Tomita et al. (2018) presented the value of H v influencing the relationship between Q a and the brightness temperature. They proposed a method using additional data for H v as a proxy for the vertical moisture gradient to retrieve Q a beyond previous linear multi-regressions. The influence of W on Q a is considered in their method.
Kanemaru and Masunaga (2013) also found that the relationship between SST and W changes depending on the values of SST [33]. W shows contrasting sensitivities to SST depending on the SST. The SST–W relation resembles the C–C (Clausius–Clapeyron) relation [34,35] when SST is lower than 27 °C. For higher SST, the SST–W relation does not follow the Clausius–Clapeyron relation because the influence of W on H v is larger than the SST contribution. Hence, the SST should be considered when we study the relationship between W and Q a .
Figure 2a–e shows the relationship between the obtained SST and W with the change in H v values. The color depth represents the value of q v . Five groups of H v (≤1300 m, 1300–1800 m, 1800–2300 m, 2300–2800 m, and >2800 m) were temporally averaged for the year 2014, using the data from 60° S to 60° N. Lower H v corresponded to relatively lower W values. For example, when H v ≤ 1300 m, the maximum value of W was less than 35 kg/m2, while the maximum value reached 85 kg/m2 if H v > 2800 m. This may be because lower H v corresponds to atmospheric columns containing moisture in relatively lower layers, and higher H v corresponds to columns containing moisture up to relatively higher layers [25].
Furthermore, compared with higher H v values, the increase ratio of d W / d S S T was smaller at lower H v values (see Figure 2). The ratio changed with SST, the thick line shows a cubic regression curve calculated from SST and W. Figure 2f shows the correlation between the increase ratio and the SST at each H v range. When the SST was lower than 27 °C, the increase ratio had less change with the increase of SST at lower H v ranges. When the SST was near 27 °C, the ratio of d W / d S S T was significantly increased (especially at higher H v ranges), and then the increase ratio was reduced. The influence of the SST–W relation on Q a will be further discussed in Section 3.2.

3.2. Channel Sensitivity

Figure 3 and Figure 4 show the relationship between Q a and the brightness temperature ( T B , 23 V and T B , 89 H ) at five H v ranges. The brightness temperature dataset is from the FY-3C MWRI. The color depth indicates the density of data and the black line was calculated using all H v data. There was a basic positive linear relationship between brightness temperature and Q a , and the H v values influenced this relationship. When the H v was small (Figure 3a–c, under 2300 m; Figure 4a,b, under 1800 m), the relationship appeared linear. However, the relationship between Q a and brightness temperature appeared to curve at higher H v ranges (Figure 3d and e, over 2300 m; Figure 4c–e, over 1800 m), which is similar to the results of Tomita et al. (2018) [25].
The basic relationship between Q a and the brightness temperature can be understood from the semi-linear relationship between W and the brightness temperature. T B , 23 and T B , 89 are well related to columnar water vapor [19,36]. When the value of W is large, the atmospheric column contains more water vapor [15]. This results in a higher brightness temperature due to higher microwave radiation from the increased amount of water vapor, so the basic relationship appears linear.
To further explore the effect of SST on the above phenomena, we divided the data shown in Figure 3 into four groups by SST. Figure 5 shows the relationship between brightness temperature ( T B , 23 V ), Q a , and SST at each H v range. Red lines were calculated for each SST interval, and black lines were calculated using all data of each H v interval. Figure 5 demonstrates that the distribution of points is influenced by SST. In the SST ≤ 10 °C range, the slope of the scatter distribution was lower than the slope of all scatter data, and their difference was the largest in this interval. As the temperature increased, the difference between them decreased. Except for H v ≤ 1300 m, the range of 20–27 °C showed minimal difference, and then the difference increased again.
This phenomenon may be due to two reasons: (1) Difference in the amount of data between different temperature ranges has an impact on the results. For example, the ratio of the number of scatters in the range of 20–27 °C to the total number of scatters was the largest, which had a greater contribution to the slope of the fitting. In contrast, for the range of 0–10 °C and >27 °C the ratio was small and corresponded to a small contribution to the overall slope. (2) The linear relationship between Q a and T B was different at each SST range, which may be affected by the SST–W relation mentioned above. The difference in the effect of the SST increment on W further caused a change in the relationship between the Q a and T B .
In the high frequency (89 H GHz) channel, we did the same analysis as shown in Figure 6. Compared with 23 V GHz, the scatters exhibited a more curvilinear relationship in the higher H v ranges. There was a significant difference in the relationship between Q a and T B in different SST ranges. We further used the observation sensitivity ( d T B / d Q a ) to analyze this phenomenon quantitatively, and the results are shown in Figure 7.
Figure 7 shows that the observation sensitivity varies with SST on the same H v range, and the relationship is affected by H v values. For the 23 V GHz channel, when H v ≤ 2300 m, the maximum value of observation sensitivity appeared in the mid-SST range, and there was minimal difference between different SST ranges. When H v > 2300 m, observation sensitivity gradually decreased with the increase of SST, which resulted in a curve distribution of scatters. In the 89 H GHz channel, there was a similar trend. When H v ≤ 1300 m, the maximum value of observation sensitivity appeared at 10–20 °C. Observation sensitivity decreased as SST increased when H v was larger than 1300 m. In addition, Figure 7 shows that the observation sensitivity value of 23 V GHz is lower than 89 H GHz channel as a whole. Furthermore, the variation amplitude between different SST ranges was lower than 89 H GHz. This difference has an effect on the relationship between Q a and T B . Since the observation sensitivity is an inverse value of d Q a / d T B , the smaller variation of observation sensitivity in different SST ranges generates a closer relationship between Q a and T B . Hence, the relationship appears linear at 23 V GHz radiation but appears curved at 89 H GHz. Higher H v values exhibit scatter points of both channels in the form of a curve. A similar result was found for the 23 H GHz and 89 V GHz observation channels, but was unclear for 10, 19, and 37 GHz channels.
In general, it was found that observation sensitivity within the same H v range changed depending on the value of SST. This trend may contribute to the curved relationship between Q a and T B , particularly at higher H v values (Figure 5 and Figure 6). This result indicates that SST is an important parameter for Q a estimates and cannot be ignored. Previous algorithms based only on direct regression between T B and Q a might produce estimation errors for some ranges of SST.

3.3. Developed Algorithm

SST has been considered as an important near-surface meteorological parameter in some artificial intelligence algorithms when retrieving Q a and LHF [23,37]. From the relationship discussed in Section 3.1 and Section 3.2, we found clear connections between in situ Q a and T B , and W was also useful because of its relationship with Q a . In this study, an improved model was proposed based on the TO18 method. We suggested using all the brightness temperature data in order to obtain a more accurate estimation of the Q a (for more details of channel correlation and its impact on the regression please see Section 5 and Table A1). Moreover, because the relationship between Q a and T B appeared curved at higher H v ranges, the quadratic terms for T B , 23 V / H and T B , 89 V / H were used for inversion. Ultimately, the new regression equation had the following form:
Q a = c 0 + c 1 T B , 10 V + c 2 T B , 10 H + c 3 T B , 19 V + c 4 T B , 19 H + c 5 T B , 23 V + c 6 T B , 23 V 2 + c 7 T B , 23 H + c 8 T B , 23 H 2 + c 9 T B , 37 V + c 10 T B , 37 H + c 11 T B , 89 V + c 12 T B , 89 V 2 + c 13 T B , 89 H + c 14 T B , 89 H 2 + c 15 W SST ,
where Q a denotes the specific air humidity at 10 m in g/kg, T B is the brightness temperature in Kelvin, SST is in °C, c 0 15 are regression coefficients, and the p-values for these retrieved coefficients are listed in Table 2.
The p-values for regression coefficients indicate the significance in estimating the Q a . A small p-value (<0.05) indicates that the regression is significant. Overall, we found that all input parameters were significant in the regression model (see Table 2). But there were some exceptions, which were related to different H v intervals. For example, when H v ≤ 1300 m, the p-values of c 3 and c 14 were larger than 0.05. In comparison, when 1300 m < H v ≤ 3300 m, all p-values were less than 0.05, representing very significant parameters in improved models. For c 14 , the quadratic coefficient of T B , 89 H in extreme H v values (≤1300 m or >3300 m) was larger than 0.05. These demonstrate that in extremely large or very small H v , the term was weakly correlated with retrieved results. In this study, we deleted the terms with p-values greater than 0.05 when obtaining the retrieved coefficients. Finally, the values of coefficients are listed in Table A3.
The LHF is widely estimated using the COARE3.0 bulk aerodynamic formula [30], which is suitable for both satellite and in situ data. The LHF can be expressed as follows:
L H F = ρ a C e U L e ( Q a Q s ) ,
where C e is the turbulent exchange coefficient for latent heat flux, L e is the latent heat of evaporation, and Q s is surface humidity. U and Q a are the wind speed and specific air humidity relative to the sea surface at the height of 10 m.
In this study, Q a was estimated using satellite data and U was derived from FY-3C MWRI orbit products. Q s was calculated from SST assuming saturation at the surface, where a multiplier factor of 0.98 was used when considering the reduction in vapor pressure caused by a typical salinity of 34 psu. C e was computed by the COARE3.0 algorithm, which was mainly influenced by conditions of light wind, stable stratification, sea spray, and other sea-state information. In addition, FY-3C MWRI orbit products of rain rate were used as input data to estimate LHF. Detailed information of these products can be obtained from the website: http://www.nsmc.org.cn/en/NSMC/Home/Index.html.

4. Experimental Results

4.1. Evaluation with In Situ Data

4.1.1. Air Specific Humidity

In this section, satellite-derived Q a fields were compared with in situ observed Q a . We also compared the proposed model with three other existing empirical algorithms (SC95, IW12, and TO18). Details of these four Q a retrieval algorithms are given in Table 3.
Figure 8 shows scatter plots of the comparison of in situ observed Q a with four satellite-derived Q a values (the Sample 2 data). The proposed model has good performance, with a root mean square difference (RMSD) of 1.53 g/kg, and a correlation coefficient (R2) of 0.94. TO18 has an RMSD and R2 of the order of 1.95 g/kg and 0.92, respectively. SC95 has an RMSD and R2 of 2.88 g/kg, and 0.86, respectively. IW12 performed better (RMSD, 2.44 g/kg; R2, 0.89) than SC95 due to the use of the high-frequency 89 GHz channel.
Figure 9 shows the zonal averages of the bias (satellite derived Q a minus in situ Q a ), RMSD, and standard deviation difference (SDD) for the Sample 2 data. For the proposed model, the RMSD are less than 2.0 g/kg and biases are less than 1 g/kg over most of the latitudes. The bias (Figure 9a) for all algorithms tends to show positive values, especially in low and mid-latitude regions. This means that all four algorithms tend to overestimate Q a . Compared with the other three algorithms, the proposed model’s Q a shows the smallest bias. The accuracy of satellite-derived Q a was significantly improved by using SST.

4.1.2. Latent Heat Flux

To assess the improvement resulting from the proposed Q a model on satellite-derived LHF, the satellite-derived LHF was compared with in situ LHF, which was calculated from in situ measurements (i.e., Sample 2 data). Table 4 provides some statistical results of these comparisons. The statistics were averaged in the high (45°–60° N, 45°–60° S), mid (15°–45° N, 15°–45° S), and low (15° S–15° N) latitudinal regions.
Compared with the other three existing algorithms, the proposed method had good performance in estimating LHF. There were significant statistical improvements of LHF due to these improvements of Q a . For example, in mid-latitude regions, the bias, RMSD, and R2 for SC95 method were 0.66 g/kg, 3.16 g/kg, and 0.81, respectively, while those of the proposed model were 0.08 g/kg, 1.76 g/kg and 0.92, respectively. Corresponding parameters of LHF (SC95) were 14.10 W/m2, 51.96 W/m2, and 0.51, while those of the proposed model were 2.40 W/m2, 34.24 W/m2, and 0.87, respectively. The differences between satellite and in situ data were dependent on zonal locations, which were mainly influenced by the accuracy of Q a .

4.2. Q a Comparison with NOAA CIRES Datasets

The satellite-derived Q a fields were also compared with NOAA CIRES Multi-Satellite Humidity. To match up two Q a datasets, we set criteria that the temporal and spatial differences had to be less than 180 min and 25 km, respectively. Figure 10 shows the two datasets covering the common area on 6 October 2014.
Figure 11 shows the results of the comparison between NOAA CIRES Multi-Satellite Humidity and FY-3C humidity (total 693,836 collocated data) on 6 October 2014. Compared with the Multi-Satellite Q a dataset, the Q a derived from FY-3C satellite had a bias, RMSD, and R2 of 0.02 g/kg, 1.02 g/kg, and 0.98, respectively.

5. Discussion

In this study, we investigated the observation sensitivity ( d T B / d Q a ) variation with the SST at different H v ranges (as described in Section 3.2). This phenomenon indicates that previous Q a algorithms based only on direct regression between brightness temperature and Q a might not be accurate enough for some SST and humidity conditions.
Hence, we introduced the SST as a key input parameter into Q a estimation (Equation (2)). Radiations from 10 microwave channels of the FY-3C satellite were all used in the proposed algorithm. It is well known that T B with the same frequency but different polarizations, has a strong correlation. The Pearson correlation coefficients also clearly show that there are strong dependencies between some channels used in Equation (2) (Table A1). To better show the importance of each indicator in the proposed model, a stepwise regression was used (Table A2). The results of the stepwise regression show that when the number of parameters used in the regression model was larger than six, the improvements of R2 and residual mean square of the regression model were very small. However, it was also shown that including more channels could improve the proposed linear model, even if the newly added channel was highly correlated to existing channels. Besides, we regressed the model coefficients according to different H v ranges. In order to obtain the best fitness of the regression model, we chose to include all channels in the proposed model. Moreover, we decided to use quadratic terms for T B , 23 V / H and T B , 89 V / H because the relationship between Q a and brightness temperature appears curved at higher H v ranges. These measures ensure a superior output of our model for different humidity conditions in the global oceans.
It is worth mentioning that the original SC95, IW12, and TO18 algorithms were derived from different training data and may not be well suited for the FY-3C MWRI instrument. Thus, we derived new regression coefficients for three algorithms before using them to compare with the proposed method. Figure 8 shows that the proposed method has the lowest RMSD (1.53 g/kg) and the highest R2 (0.94). In addition, satellite-derived Q a were also compared with the NOAA CIRES Multi-Satellite humidity dataset, their bias, RMSD, and R2 were 0.02 g/kg, 1.02 g/kg, and 0.98, respectively. These comparison results indicate that adding SST as an input parameter can improve Q a and LHF estimation. There are regional differences of the outputs of the proposed model. The RMSD of Q a in mid-latitude regions (15°–45° N, 15°–45° S) was 1.76 g/kg, and the RMSD of high (45°–60° N, 45°–60° S) and low (15° S–15° N) latitude regions were 1.11 and 1.51 g/kg, respectively. This difference may result from the SST–W relation. Kanemaru and Masunaga (2013) concluded that the W is closely related with SST and dependent on regional position [33].
We obtained SST data from the FY-3C MWRI orbit products and H v data from atmospheric reanalysis. Therefore, any potential error and dependency for either variable may affect the accuracy of Q a . Furthermore, the training and testing dataset of in situ measurements was collected from the global sea area between 60° S and 60° N. This means that the proposed Q a model is also only valid for these regions. Although we believe that the form of this algorithm is universal in calculating global specific humidity, the retrieved coefficients are only available for FY-3C MWRI. For other space-borne radiometers, these coefficients need to be recalculated.

6. Conclusions

Latent heat flux represents an important aspect of the atmosphere–ocean interaction. Additionally, satellite-derived specific humidity has been a major source of error for the estimation of LHF. In order to improve the estimation accuracy of specific humidity, we proposed a new Q a model based on FY-3C microwave brightness temperature observations. We found the relationship of observation sensitivity ( d T B / d Q a ) changed depending on the range of SST, especially at larger H v , so SST was added into the estimation model. This improved model gives better estimates of Q a and LHF when compared with three existing methods (SC95, IW12, and TO18). The satellite-derived Q a were also compared with in situ data, and their bias, RMSD, and R2 were 0.17 g/kg, 1.53 g/kg, and 0.94, respectively. The satellite Q a were also compared with NOAA CIRES Q a datasets, with a bias, RMSD, and R2 of 0.02 g/kg, 1.02 g/kg, and 0.98, respectively. The differences between satellite LHF and in situ LHF were dependent on zonal locations: In high latitude regions (45°–60° N, 45°–60° S), the bias, RMSD, and R2 were 12.89 W/m2, 23.44 W/m2, and 0.72, respectively; in mid-latitude regions (15°–45° N, 15°–45° S), the bias, RMSD, and R2 were 2.4 W/m2, 34.24 W/m2, and 0.87, respectively; in low latitude regions (15° S–15° N), the bias, RMSD, and R2 were 4.87 W/m2, 29.05 W/m2, and 0.82, respectively, which were mainly influenced by the accuracy of Q a (see Table 4).
As shown in the study, SST as an additional input parameter greatly improved the estimation accuracy of Q a and LHF for the FY-3C satellite observation. The innovative aspect of the algorithm was based on the finding of a relationship between SST, W, and water vapor profile information, which provided a theoretical basis for satellite-derived Q a and LHF. Furthermore, SST and W were directly used as input parameters for estimating Q a . Both variables may introduce unknown redundant information to model outputs. Future work should explore the essence of the SST–W relationship and improve the simulation method of the vertical moisture structure.

Author Contributions

All authors made great contributions to this article. Conceptualization, X.Y.; methodology, X.Y. and Q.G.; software, Q.G. and S.W.; validation, Q.G.; writing—original draft, Q.G. and X.Y.; and writing—review and editing, X.Y. and Q.G.

Funding

This study is supported by the National Key Research and Development Program of China under Grant 2017YFB0502800, in part by the National Natural Science Foundation of China under Grant 41871268, and in part by the Key Research and Development Program of Hainan Province under Grant ZDYF2017167.

Acknowledgments

The authors are grateful to the National Satellite Meteorological Center, FENGYUN Satellite Data Center of China. The authors would also like to thank ECMWF, NDBC, NOAA, Météo-France, UK Met Office, PMEL, and Woods Hole Oceanographic Institution for providing numerical and in situ data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Values of Pearson’s correlation coefficient for the parameters in Equation (2).
Table A1. Values of Pearson’s correlation coefficient for the parameters in Equation (2).
Q a B1B2B3B4B5B6B7B8B9B10W·SST
Q a 1.000
B1 *0.4591.000
B20.2200.9161.000
B30.7440.8050.6831.000
B40.6700.8220.7400.9611.000
B50.9010.6240.4370.9100.8801.000
B60.8710.6270.4580.9160.9010.9931.000
B70.6760.6670.5670.9120.9350.8920.9101.000
B80.5430.6290.5810.8550.9130.7970.8350.9711.000
B90.8690.4030.1850.6950.6250.8740.8500.6950.5671.000
B100.8180.4390.2630.7610.7340.9080.9080.8180.7420.9491.000
W·SST0.9470.4600.2420.7740.7180.9100.8900.7260.6120.8210.8051.000
* B1–B10 are 10 V/H, 19 V/H, 23 V/H, 37 V/H, and 89 V/H GHz, respectively.
Table A2. Stepwise regression analysis for predicting the specific air humidity (Qa).
Table A2. Stepwise regression analysis for predicting the specific air humidity (Qa).
ModelAdjusted R SquareResidual Mean SquareF-Test
Fdf1df2Sig. F Change
1 *0.8972.544367104,793.588111,5980.000
20.9221.933543708,22.499111,5970.000
30.9281.77498051,823.867111,5960.000
40.9301.63641742,380.563111,5950.000
50.9351.48837737,518.207111,5940.000
60.9401.42801332,669.545111,5930.000
70.9421.39930928,620.837111,5920.000
80.9431.39462625,140.234111,5910.000
90.9441.38229722,546.893111,5900.000
100.9441.37866020,621.421111,5890.000
110.9451.36134319,567.475111,5880.000
* 1: Constant, WSST; 2: Constant, WSST, B9; 3. Constant, WSST, B9, B10; 4: Constant, WSST, B9, B10, B5; 5: Constant, WSST, B9, B10, B5, B7; 6: Constant, WSST, B9, B10, B5, B7, B6; 7: Constant, WSST, B9, B10, B5, B7, B6, B8; 8: Constant, WSST, B9, B10, B5, B7, B6, B8, B4; 9: Constant, WSST, B9, B10, B5, B7, B6, B8, B4, B1; 10: Constant, WSST, B9, B10, B5, B7, B6, B8, B4, B1, B2; 11: Constant, WSST, B9, B10, B5, B7, B6, B8, B4, B1, B2, B3.
Table A3. Regression coefficients of proposed model for different water vapor scale heights H v   (m).
Table A3. Regression coefficients of proposed model for different water vapor scale heights H v   (m).
Coefficients H v     1300 1300   <   H v   1800 1800   <   H v   2300 2300   <   H v   2800 2800   <   H v   3300 H v >   3300
C0−101.7520−74.1441−56.4953−46.2155−61.2600−86.3314
C10.0252−0.0103−0.0149−0.0089−0.0725−0.0519
C2−0.01250.00930.00430.00210.03250.0183
C30.0000−0.0110−0.0404−0.0725−0.1106−0.0247
C4−0.0358−0.01380.01050.01630.04180.0000
C5−0.20150.16040.36490.67170.91741.1411
C60.0005−0.0003−0.0009−0.0013−0.0019−0.0024
C70.10120.05480.0359−0.1322−0.2903−0.3309
C8−0.0002−2.40 × 10−57.54 × 10−50.00040.00080.0010
C9−0.0902−0.0450−0.0123−0.00970.06840.0000
C100.0235−0.0264−0.0228−0.0029−0.0442−0.0223
C110.99190.56580.2115−0.1695−0.20980.0000
C12−0.0016−0.0010−0.00030.00050.00058.82 × 10−5
C13−0.0743−0.0805−0.02820.10240.1952−0.0173
C140.00000.00011.12 × 10−5−0.0003−0.00050.0000
C150.01800.01210.00970.00740.00710.0061

References

  1. Gulev, S.K.; Latif, M.; Keenlyside, N.; Park, W.; Koltermann, K.P. North atlantic ocean control on surface heat flux on multidecadal timescales. Nature 2013, 499, 464–467. [Google Scholar] [CrossRef] [PubMed]
  2. Dutra, E.; Stepanenko, V.M.; Balsamo, G.; Viterbo, P.; Miranda, P.M.A.; Mironov, D.; Schar, C. An offline study of the impact of lakes on the performance of the ecmwf surface scheme. Boreal Environ. Res. 2010, 15, 100–112. [Google Scholar]
  3. Woolway, R.I.; Verburg, P.; Lenters, J.D.; Merchant, C.J.; Hamilton, D.P.; Brookes, J.; de Eyto, E.; Kelly, S.; Healey, N.C.; Hook, S.; et al. Geographic and temporal variations in turbulent heat loss from lakes: A global analysis across 45 lakes. Limnol. Oceanogr. 2018, 63, 2436–2449. [Google Scholar] [CrossRef]
  4. Bourassa, M.A.; Gille, S.T.; Bitz, C.; Carlson, D.; Cerovecki, I.; Clayson, C.A.; Cronin, M.F.; Drennan, W.M.; Fairall, C.W.; Hoffman, R.N.; et al. High-latitude ocean and sea ice surface fluxes: Challenges for climate research. Bull. Am. Meteorol. Soc. 2013, 94, 403–423. [Google Scholar] [CrossRef]
  5. Bunker, A.F. Computations of surface-energy flux and annual air-sea interaction cycles of north-atlantic ocean. Mon. Weather Rev. 1976, 104, 1122–1140. [Google Scholar] [CrossRef]
  6. Hasse, L.; Smith, S.D. Local sea surface wind, wind stress, and sensible and latent heat fluxes. J. Clim. 1997, 10, 2711–2724. [Google Scholar] [CrossRef]
  7. Bentamy, A.; Piolle, J.F.; Grouazel, A.; Danielson, R.; Gulev, S.; Paul, F.; Azelmat, H.; Mathieu, P.P.; von Schuckmann, K.; Sathyendranath, S.; et al. Review and assessment of latent and sensible heat flux accuracy over the global oceans. Remote Sens. Environ. 2017, 201, 196–218. [Google Scholar] [CrossRef]
  8. Andersson, A.; Fennig, K.; Klepp, C.; Bakan, S.; Grassl, H.; Schulz, J. The hamburg ocean atmosphere parameters and fluxes from satellite data—Hoaps-3. Earth Syst. Sci. Data 2010, 3, 143–194. [Google Scholar] [CrossRef]
  9. Andersson, A.; Klepp, C.; Fennig, K.; Bakan, S.; Grassl, H.; Schulz, J. Evaluation of hoaps-3 ocean surface freshwater flux components. J. Appl. Meteorol. Climatol. 2011, 50, 379–398. [Google Scholar] [CrossRef]
  10. Kubota, M.; Iwasaka, N.; Kizu, S.; Kondo, M.; Kutsuwada, K. Japanese ocean flux data sets with use of remote sensing observations (j-ofuro). J. Oceanogr. 2002, 58, 213–225. [Google Scholar] [CrossRef]
  11. Shie, C.L.; Chiu, L.S.; Adler, R.; Nelkin, E.; Lin, I.I.; Xie, P.P.; Wang, F.C.; Chokngamwong, R.; Olson, W.; Chu, D.A. A note on reviving the goddard satellite-based surface turbulent fluxes (gsstf) dataset. Adv. Atmos. Sci. 2009, 26, 1071–1080. [Google Scholar] [CrossRef]
  12. Bentamy, A.; Grodsky, S.A.; Katsaros, K.; Mestas-Nunez, A.M.; Blanke, B.; Desbiolles, F. Improvement in airsea flux estimates derived from satellite observations. Int. J. Remote Sens. 2013, 34, 5243–5261. [Google Scholar] [CrossRef]
  13. Tomita, H.; Kubota, M.; Cronin, M.F.; Iwasaki, S.; Konda, M.; Ichikawa, H. An assessment of surface heat fluxes from j-ofuro2 at the keo and jkeo sites. J. Geophys. Res. Oceans 2010, 115. [Google Scholar] [CrossRef]
  14. Liu, W.T.; Niiler, P.P. Determination of monthly mean humidity in the atmospheric surface-layer over oceans from satellite data. J. Phys. Oceanogr. 1984, 14, 1451–1457. [Google Scholar] [CrossRef]
  15. Liu, W.T. Statistical relation between monthly mean precipitable water and surface-level humidity over global oceans. Mon. Weather Rev. 1986, 114, 1591–1602. [Google Scholar] [CrossRef]
  16. Schulz, J.; Schluessel, P.; Grassl, H. Water-vapor in the atmospheric boundary layer over oceans from ssm/i measurements. Int. J. Remote Sens. 1993, 14, 2773–2789. [Google Scholar] [CrossRef]
  17. Schlussel, P.; Schanz, L.; Englisch, G. Retrieval of latent-heat flux and longwave irradiance at the sea-surface from ssm/i and avhrr measurements. In Satellite Monitoring of the Earth’s Surface and Atmosphere; Arnault, S., Ed.; Pergamon Press Ltd: Oxford, UK, 1995; Volume 16, pp. 107–116. [Google Scholar]
  18. Bentamy, A.; Katsaros, K.B.; Mestas-Nunez, A.M.; Drennan, W.M.; Forde, E.B.; Roquet, H. Satellite estimates of wind speed and latent heat flux over the global oceans. J. Clim. 2003, 16, 637–656. [Google Scholar] [CrossRef]
  19. Iwasaki, S.; Kubota, M. Algorithms for estimation of air-specific humidity using tmi data. Int. J. Remote Sens. 2012, 33, 7413–7430. [Google Scholar] [CrossRef]
  20. Jackson, D.L.; Wick, G.A.; Bates, J.J. Near-surface retrieval of air temperature and specific humidity using multisensor microwave satellite observations. J. Geophys. Res. Atmos. 2006, 111. [Google Scholar] [CrossRef]
  21. Jackson, D.L.; Wick, G.A.; Robertson, F.R. Improved multisensor approach to satellite-retrieved near-surface specific humidity observations. J. Geophys. Res. Atmos. 2009, 114. [Google Scholar] [CrossRef]
  22. Jones, C.; Peterson, P.; Gautier, C. A new method for deriving ocean surface specific humidity and air temperature: An artificial neural network approach. J. Appl. Meteorol. 1999, 38, 1229–1245. [Google Scholar] [CrossRef]
  23. Singh, R.; Joshi, P.C.; Kishtawal, C.M. A new technique for estimation of surface latent heat fluxes using satellite-based observations. Mon. Weather Rev. 2005, 133, 2692–2710. [Google Scholar] [CrossRef]
  24. Roberts, J.B.; Clayson, C.A.; Robertson, F.R.; Jackson, D.L. Predicting near-surface atmospheric variables from special sensor microwave/imager using neural networks with a first-guess approach. J. Geophys. Res. Atmos. 2010, 115. [Google Scholar] [CrossRef]
  25. Tomita, H.; Hihara, T.; Kubota, M. Improved satellite estimation of near-surface humidity using vertical water vapor profile information. Geophys. Res. Lett. 2018, 45, 899–906. [Google Scholar] [CrossRef]
  26. Yang, J.; Zhang, P.; Lu, N.M.; Yang, Z.D.; Shi, J.M.; Dong, C.H. Improvements on global meteorological observations from the current fengyun 3 satellites and beyond. Int. J. Digit. Earth 2012, 5, 251–265. [Google Scholar] [CrossRef]
  27. Freeman, E.; Woodruff, S.D.; Worley, S.J.; Lubker, S.J.; Kent, E.C.; Angel, W.E.; Berry, D.I.; Brohan, P.; Eastman, R.; Gates, L.; et al. Icoads release 3.0: A major update to the historical marine climate record. Int. J. Climatol. 2017, 37, 2211–2232. [Google Scholar] [CrossRef]
  28. Buck, A.L. New equations for computing vapor-pressure and enhancement factor. J. Appl. Meteorol. 1981, 20, 1527–1532. [Google Scholar] [CrossRef]
  29. Wilks, D.S. Statistical Methods in the Atmospheric Sciences, 2nd ed.; Academic Press: San Diego, CA, USA, 2006; 627p. [Google Scholar]
  30. Fairall, C.W.; Bradley, E.F.; Hare, J.E.; Grachev, A.A.; Edson, J.B. Bulk parameterization of air-sea fluxes: Updates and verification for the coare algorithm. J. Clim. 2003, 16, 571–591. [Google Scholar] [CrossRef]
  31. Jackson, D.L.; Wick, G.A. Near-surface air temperature retrieval derived from amsu-a and sea surface temperature observations. J. Atmos. Ocean. Technol. 2010, 27, 1769–1776. [Google Scholar] [CrossRef]
  32. Dee, D.P.; Uppala, S.M.; Simmons, A.J.; Berrisford, P.; Poli, P.; Kobayashi, S.; Andrae, U.; Balmaseda, M.A.; Balsamo, G.; Bauer, P.; et al. The era-interim reanalysis: Configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 2011, 137, 553–597. [Google Scholar] [CrossRef]
  33. Kanemaru, K.; Masunaga, H. A satellite study of the relationship between sea surface temperature and column water vapor over tropical and subtropical oceans. J. Clim. 2013, 26, 4204–4218. [Google Scholar] [CrossRef]
  34. Stephens, G.L. On the relationship between water-vapor over the oceans and sea-surface temperature. J. Clim. 1990, 3, 634–645. [Google Scholar] [CrossRef]
  35. Jackson, D.L.; Stephens, G.L. A study of ssm/i-derived columnar water-vapor over the global oceans. J. Clim. 1995, 8, 2025–2038. [Google Scholar] [CrossRef]
  36. Grody, N.C. Remote-sensing of atmospheric water-content from satellites using microwave radiometry. IEEE Trans. Antennas Propag. 1976, 24, 155–162. [Google Scholar] [CrossRef]
  37. Singh, R.; Joshi, P.C.; Kishtawal, C.M.; Pal, P.K. A new method for estimation of near surface specific humidity over global oceans. Meteorol. Atmos. Phys. 2006, 94, 1–10. [Google Scholar] [CrossRef]
Figure 1. Methodological framework to estimate global air–sea Q a and latent heat flux (LHF) using FY-3C satellite observations.
Figure 1. Methodological framework to estimate global air–sea Q a and latent heat flux (LHF) using FY-3C satellite observations.
Remotesensing 11 00466 g001
Figure 2. (ae) The relationship between sea surface temperature (SST), column water vapor (W), and surface water vapor mixing ratio q v (g/kg, color) on different ranges of water vapor scale height H v (m). Each figure shows a cubic regression curve calculated from SST and W. (f) SST compared to increase ratio (W/SST, kg/m2·°C) over five H v   ranges, calculated from daily data and averaged from the range of 60° S to 60° N.
Figure 2. (ae) The relationship between sea surface temperature (SST), column water vapor (W), and surface water vapor mixing ratio q v (g/kg, color) on different ranges of water vapor scale height H v (m). Each figure shows a cubic regression curve calculated from SST and W. (f) SST compared to increase ratio (W/SST, kg/m2·°C) over five H v   ranges, calculated from daily data and averaged from the range of 60° S to 60° N.
Remotesensing 11 00466 g002
Figure 3. Scatter diagrams showing the relationship between   T B , 23 V , Q a , and H v , using the Sample 1 data. (a) For H v ≤ 1300 m; (b) for 1300 m < H v ≤ 1800 m; (c) for 1800 m < H v ≤ 2300 m; (d) for 2300 m < H v ≤ 2800 m; (e) for H v > 2800 m; and (f) for all H v   ranges. Each figure shows a linear regression line calculated using all H v range data. Color depth indicates the density of data.
Figure 3. Scatter diagrams showing the relationship between   T B , 23 V , Q a , and H v , using the Sample 1 data. (a) For H v ≤ 1300 m; (b) for 1300 m < H v ≤ 1800 m; (c) for 1800 m < H v ≤ 2300 m; (d) for 2300 m < H v ≤ 2800 m; (e) for H v > 2800 m; and (f) for all H v   ranges. Each figure shows a linear regression line calculated using all H v range data. Color depth indicates the density of data.
Remotesensing 11 00466 g003
Figure 4. Scatter diagrams showing the relationship between T B , 89 H , Q a , and H v , using the Sample 1 data. (a) For H v ≤ 1300 m; (b) for 1300 m <   H v ≤ 1800 m; (c) for 1800 m < H v ≤ 2300 m; (d) for 2300 m < H v ≤ 2800 m; (e) for H v > 2800 m; and (f) for all H v   ranges. Each figure shows a linear regression line calculated using all H v range data. Color depth indicates the density of data.
Figure 4. Scatter diagrams showing the relationship between T B , 89 H , Q a , and H v , using the Sample 1 data. (a) For H v ≤ 1300 m; (b) for 1300 m <   H v ≤ 1800 m; (c) for 1800 m < H v ≤ 2300 m; (d) for 2300 m < H v ≤ 2800 m; (e) for H v > 2800 m; and (f) for all H v   ranges. Each figure shows a linear regression line calculated using all H v range data. Color depth indicates the density of data.
Remotesensing 11 00466 g004
Figure 5. The relationship between T B , 23 V , Q a , and SST. (ae) For H v ≤ 1300 m; (fj) for 1300 m <   H v ≤ 1800 m; (ko) for 1800 m < H v ≤ 2300 m; (pt) for 2300 m < H v ≤ 2800 m; and (uy) for H v > 2800 m. Each subfigure shows a linear regression black line calculated using all data of the same H v range; the red line was calculated using the corresponding SST range data.
Figure 5. The relationship between T B , 23 V , Q a , and SST. (ae) For H v ≤ 1300 m; (fj) for 1300 m <   H v ≤ 1800 m; (ko) for 1800 m < H v ≤ 2300 m; (pt) for 2300 m < H v ≤ 2800 m; and (uy) for H v > 2800 m. Each subfigure shows a linear regression black line calculated using all data of the same H v range; the red line was calculated using the corresponding SST range data.
Remotesensing 11 00466 g005
Figure 6. The relationship between T B , 89 H , Q a , and SST. (ae) For H v ≤ 1300 m; (fj) for 1300 m <   H v ≤ 1800 m; (ko) for 1800 m <   H v ≤ 2300 m; (pt) for 2300 m <   H v ≤ 2800 m; and (uy) for H v > 2800 m. Each subfigure shows a linear regression black line calculated using all data of the same H v range; the red line was calculated using the corresponding SST range data.
Figure 6. The relationship between T B , 89 H , Q a , and SST. (ae) For H v ≤ 1300 m; (fj) for 1300 m <   H v ≤ 1800 m; (ko) for 1800 m <   H v ≤ 2300 m; (pt) for 2300 m <   H v ≤ 2800 m; and (uy) for H v > 2800 m. Each subfigure shows a linear regression black line calculated using all data of the same H v range; the red line was calculated using the corresponding SST range data.
Remotesensing 11 00466 g006
Figure 7. SST compared with observation sensitivity from a linear fitting for each SST range. (a) For 23 V GHz and (b) for 89 H GHz.
Figure 7. SST compared with observation sensitivity from a linear fitting for each SST range. (a) For 23 V GHz and (b) for 89 H GHz.
Remotesensing 11 00466 g007
Figure 8. Comparisons between in situ observed and four satellite-derived Q a values (g/kg). (a) For SC95; (b) for IW12; (c) for TO18; and (d) for the proposed model.
Figure 8. Comparisons between in situ observed and four satellite-derived Q a values (g/kg). (a) For SC95; (b) for IW12; (c) for TO18; and (d) for the proposed model.
Remotesensing 11 00466 g008
Figure 9. Zonal distributions of the comparative statistics between in situ observed and satellite-derived Q a values for four algorithms, including (a) averaged bias, (b) root mean square difference (RMSD), and standard deviation difference (SDD, plotted as red lines). Statistics were averaged for each 2° bias.
Figure 9. Zonal distributions of the comparative statistics between in situ observed and satellite-derived Q a values for four algorithms, including (a) averaged bias, (b) root mean square difference (RMSD), and standard deviation difference (SDD, plotted as red lines). Statistics were averaged for each 2° bias.
Remotesensing 11 00466 g009
Figure 10. The common coverage area of NOAA CIRES Multi-Satellite Humidity (dark blue strip) and FY-3C humidity (colored strip) on 6 October 2014. Color depth indicates the value of specific air humidity. (a) Ascending orbit and (b) descending orbit.
Figure 10. The common coverage area of NOAA CIRES Multi-Satellite Humidity (dark blue strip) and FY-3C humidity (colored strip) on 6 October 2014. Color depth indicates the value of specific air humidity. (a) Ascending orbit and (b) descending orbit.
Remotesensing 11 00466 g010
Figure 11. Results of the comparison between NOAA CIRES Multi-Satellite Humidity and FY-3C humidity on 6 October 2014. Color depth indicates the density of data.
Figure 11. Results of the comparison between NOAA CIRES Multi-Satellite Humidity and FY-3C humidity on 6 October 2014. Color depth indicates the density of data.
Remotesensing 11 00466 g011
Table 1. Main parameters for FY-3C Microwave Radiation Imager (MWRI).
Table 1. Main parameters for FY-3C Microwave Radiation Imager (MWRI).
Channel NameCentral Frequency (GHz)Bandwidth (MHz)IFOVResolution
10 V/H *10.6518051 × 85 km
19 V/H18.720030 × 50 km
23 V/H23.840027 × 45 km
37 V/H36.590018 × 30 km
89 V/H89.046009 × 15 km
* V: Vertical polarization and H: Horizontal polarization.
Table 2. p-Values of retrieved coefficients for different water vapor scale heights H v (m).
Table 2. p-Values of retrieved coefficients for different water vapor scale heights H v (m).
p-Values H v     1300 1300   <   H v 1800 1800   <   H v   2300 2300   <   H v   2800 2800   <   H v   3300 H v   >   3300
C01.51 × 10−630.0000 0.0000 0.0000 0.0000 6.60 × 10−77
C13.27 × 10154.71 × 10696.58 × 10945.12 × 101410.0000 3.33 × 10−76
C22.41 × 10068.39 × 102272.40 × 10594.94 × 101520.0000 3.77 × 10−37
C30.532 9.67 × 101250.0000 0.0000 0.0000 6.14 × 10−08
C44.02 × 10111.68 × 103630.0272 0.0010 8.81 × 10220.2000
C51.04 × 10060.0000 0.0000 0.0000 0.0000 5.22 × 10−132
C64.71 × 10080.0000 0.0000 0.0000 0.0000 3.25 × 10−124
C71.17 × 10130.0000 2.71 × 10050.0000 0.0000 3.60 × 10−107
C83.86 × 10050.0000 5.87 × 101560.0000 0.0000 6.83 × 10−145
C91.53 × 10423.38 × 10566.78 × 10158.98 × 10066.73 × 10170.0640
C101.10 × 10090.0000 0.0000 3.65 × 10454.90 × 10889.94 × 10−19
C116.15 × 10590.0000 0.0000 5.15 × 102238.84 × 10900.4254
C123.89 × 10370.0000 1.10 × 102260.0000 4.14 × 101180.0214
C133.09 × 10100.0000 0.0000 2.01 × 102151.02 × 10−1500.0070
C140.647 0.0000 2.93 × 101410.0000 1.19 × 10−1880.6800
C151.51 × 10630.0000 0.0000 0.0000 0.0000 6.62 × 10−77
Table 3. Channels and meteorological parameters for each algorithm.
Table 3. Channels and meteorological parameters for each algorithm.
Algorithm10 GHz 19 GHz23 GHz37 GHz89 GHzParameters
SC95 19 V/H23 V37 V/H
IW12 19 V/H23 V37 V/H89 V/H
TO1810 V/H19 V/H23 V/H37 V/H89 V/HW
Proposed model10 V/H19 V/H23 V/H37 V/H89 V/HW and SST
SC95: Schlüssel et al. (1995); IW12: Iwasaki et al. (2012); and TO18: Tomita et al. (2018).
Table 4. Comparisons between satellite and in situ data for the four algorithms.
Table 4. Comparisons between satellite and in situ data for the four algorithms.
StatisticsSpecific Air HumidityLatent Heat Flux
Bias
(g/kg)
RMSD
(g/kg)
Correlation
(R2)
Bias
(W/m2)
RMSD (W/m2)Correlation
(R2)
SC95
45°–60° N
45°–60° S
−1.342.460.85−26.5039.240.48
15°–45° N
15°– 45° S
0.663.160.8114.1051.960.51
15° S–15° N1.722.960.7034.9848.140.49
IW12
45°–60° N
45°– 60° S
−1.162.290.86−22.7636.540.49
15°–45° N
15°–45° S
0.532.980.8411.8948.840.55
15° S–15° N1.412.410.7428.7442.610.54
TO18
45°– 60° N
45°– 60° S
−0.481.570.93−16.1831.870.63
15°– 45° N
15°–45° S
0.152.200.9010.9844.140.75
15° S–15° N0.441.830.828.4535.500.70
Proposed model
45°–60° N
45°–60° S
−0.141.110.95−12.8923.440.72
15°–45° N
15°–45° S
0.081.760.922.4034.240.87
15° S–15° N0.241.510.864.8729.050.82

Share and Cite

MDPI and ACS Style

Gao, Q.; Wang, S.; Yang, X. Estimation of Surface Air Specific Humidity and Air–Sea Latent Heat Flux Using FY-3C Microwave Observations. Remote Sens. 2019, 11, 466. https://doi.org/10.3390/rs11040466

AMA Style

Gao Q, Wang S, Yang X. Estimation of Surface Air Specific Humidity and Air–Sea Latent Heat Flux Using FY-3C Microwave Observations. Remote Sensing. 2019; 11(4):466. https://doi.org/10.3390/rs11040466

Chicago/Turabian Style

Gao, Qidong, Sheng Wang, and Xiaofeng Yang. 2019. "Estimation of Surface Air Specific Humidity and Air–Sea Latent Heat Flux Using FY-3C Microwave Observations" Remote Sensing 11, no. 4: 466. https://doi.org/10.3390/rs11040466

APA Style

Gao, Q., Wang, S., & Yang, X. (2019). Estimation of Surface Air Specific Humidity and Air–Sea Latent Heat Flux Using FY-3C Microwave Observations. Remote Sensing, 11(4), 466. https://doi.org/10.3390/rs11040466

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