Atmospheric Retrievals and Assessment for Microwave Observations from Chinese FY-3C Satellite during Hurricane Matthew

: The evolution process of hurricane Matthew (NO. 8, 2016) was simulated using the mesoscale Weather Research and Forecasting (WRF) model at temporal resolution of 5 min and spatial resolution of 15 km. The atmospheric temperature and humidity proﬁles were retrieved accordingly for diagnostic analysis of the short-term heavy rainstorm. The satellite-based microwave observations from Microwave Humidity and Temperature Sounder (MWHTS) instrument on board the FY-3C polar-orbiting satellite were matched with the WRF grid points. In particular, the in-orbit calibration and data quality control are detailed, and an innovative method combining artiﬁcial neural network (ANN) and 1-D variational approach is presented to derive the high-performance retrieval proﬁles. Results show that the root-mean-square errors of the retrieved temperature and water vapor density proﬁles are 0.75 K and 0.41 g / m 3 , respectively. In addition, this study used both the retrievals and radiance from MWHTS as input to the WRF Data Assimilation (WRFDA) model to forecast the track and intensity of hurricane Matthew. The forecast results were cross-compared with the best track to verify the radiance quality and performance of the retrievals, especially for the 118 GHz channel, which was ﬁrstly used in meteorological satellite. study on the retrieved atmospheric temperature and water vapor using satellite observations provided the reference value for assessing the detection ability of 118.75 GHz oxygen line. provided by large-scale analysis from operational centers are often ill-defined, too weak, and sometimes misplaced. The methodology presented in this study is a new attempt to provide more confidence on the satellite-based atmospheric profile applications in hurricane domains. In addition to the atmospheric temperature and humidity profiles, precipitation, as one of the MWHTS L2 products, will also be an available choice to improve the hurricane storm track density. For the study case, using microwave radiance of MWHTS is recommended, which can improve the forecasting track by 50 km combining atmospheric temperature and humidity information.


Introduction
Several satellites in operation are equipped with high-frequency microwave sensors, including the Advanced Microwave Sounding Unit B (AMSU-B) [1] or Microwave Humidity Sounder (MHS) on the National Oceanic and Atmospheric Administration (NOAA) and Metop satellite series [2], Advanced Technology Microwave Sounder (ATMS) on the Suomi National Polar-orbiting Partnership (SNPP) and Joint Polar Satellite System (JPSS)-1 satellites [3], and MicroWave Humidity and temperature Sounder (MWHTS) on the Chinese FengYun-3 C (FY-3C) satellite [4]. Satellite-based observations provide expanded opportunities for atmospheric profile retrievals and hydro-meteorological monitoring through measuring brightness temperature over land and ocean at global scale.
Accurate profile and hydrometeor retrievals are crucial for applications in extreme weather events such as hurricanes. Microwave observations from radiometers can be used to produce physically reasonable retrievals through estimating the direct interaction between the radiation and atmospheric parameters in the atmospheric columns. A physical-based retrieval algorithm was introduced by Petty [5], aiming to invert multichannel microwave radiances to determine the physical information on temperature and hydrometeors. The Goddard Profiling Algorithm (GPROF)

Description of Instrument
The FY-3C satellite is polar orbiting with an equatorial crossing time of 10:00 (descending) and carries the microwave humidity and temperature sounders, MWHTS, which is similar to the ATMS temperature and humidity sounder onboard the Suomi-NPP satellite, except that MWHTS has additional channels at 118 GHz and lacks 23-60 GHz. The MWHTS has been in operation onboard the FY-3C satellite since 30 September 2013 [15]. It is a cross-track scanning radiometer that has 98 Fields of View (FOV) with a swath width of 2650 km, wider than that for MHS (2310 km) and ATMS (2580 km). The absolute visiting period is 5.5 days. The limitation threshold of matching domain (i.e., accuracy in domain selection for comparison) is 0.05 degree for both latitude and longitude in the spatial dimension, and 30 min in the temporal dimension. The MWHTS instrument is a selfcalibrating total-power passive microwave system. It can operate in three different modes: pointing mode, constant speed scanning mode and variable speed scanning mode. It is a cross-track scanner with variable speed and the antenna reflectors execute one complete revolution every 8/3 s. As illustrated in Figure 1, the scan view angle of earth field is within ±53.35° from the nadir direction. The cold space view angle is 106.1° and the zenith direction where the blackbody target is located is 180°. Figure 2 shows the scanning pixels of MWHTS, which has an instantaneous field-of-view of 1.1 or 2 degree and scans from nadir with 49 different viewing angles on each side, providing a nominal spatial resolution of 15 km at its nadir and 41X27 km for first and last pixel.

Description of Instrument
The FY-3C satellite is polar orbiting with an equatorial crossing time of 10:00 (descending) and carries the microwave humidity and temperature sounders, MWHTS, which is similar to the ATMS temperature and humidity sounder onboard the Suomi-NPP satellite, except that MWHTS has additional channels at 118 GHz and lacks 23-60 GHz. The MWHTS has been in operation onboard the FY-3C satellite since 30 September 2013 [15]. It is a cross-track scanning radiometer that has 98 Fields of View (FOV) with a swath width of 2650 km, wider than that for MHS (2310 km) and ATMS (2580 km). The absolute visiting period is 5.5 days. The limitation threshold of matching domain (i.e., accuracy in domain selection for comparison) is 0.05 degree for both latitude and longitude in the spatial dimension, and 30 min in the temporal dimension. The MWHTS instrument is a selfcalibrating total-power passive microwave system. It can operate in three different modes: pointing mode, constant speed scanning mode and variable speed scanning mode. It is a cross-track scanner with variable speed and the antenna reflectors execute one complete revolution every 8/3 s. As illustrated in Figure 1, the scan view angle of earth field is within ±53.35° from the nadir direction. The cold space view angle is 106.1° and the zenith direction where the blackbody target is located is 180°. Figure 2 shows the scanning pixels of MWHTS, which has an instantaneous field-of-view of 1.1 or 2 degree and scans from nadir with 49 different viewing angles on each side, providing a nominal spatial resolution of 15 km at its nadir and 41X27 km for first and last pixel.  The MWHTS has 15 channels and is mainly designed to provide information on atmospheric temperature and humidity profiles, rainfall and typhoon monitoring [16]. The instrument measures thermal radiation at microwave frequencies ranging 89-190.0 GHz. Atmospheric temperature profiles Remote Sens. 2019, 11, 896 4 of 18 are obtained primarily based on the measurements at channels near 118.75 GHz, which is an oxygen absorption band. In particular, the MWHTS sounding channels 2-9 respond to the thermal radiation at various altitudes as according to their weighting function distributions. As window channels, channels 1 and 10 mainly detect radiation background, clouds, precipitation, and surface properties. Since the satellite provides a nominal spatial resolution of 15 km at its nadir, the temperature perturbations from synoptic to mesoscale can be reasonably depicted. The instrument also contains five channels (11)(12)(13)(14)(15) that provide atmospheric humidity profiles, as shown in Table 1. The fifteen channels from MWHTS are designed for sounding atmospheric temperature and water vapor in all weather conditions except for heavy precipitation [17,18]. Since these instruments are designed as cross-track scanners, the instrument viewing angle is different from pixel to pixel along a scan line, as shown in Figure 1. A combination of several imaging channels at 89 and 150 GHz frequencies is utilized to determine cloud liquid and ice water content since they directly respond to the emission from liquid droplets and the scattering from ice particles. The MWHTS channels are non-polarized at nadir and mix-polarized off nadir. That is, the observed brightness temperature is only affected by one polarization state at nadir, but it is affected by both horizontal and vertical polarization off nadir. The polarization mixture reduces the angular dependence of the surface emissivity.

Simulation and Validation
In this study, we used the National Centers for Environmental Prediction (NCEP) 6-h reanalysis data [19] as input to the Weather Research and Forecasting (WRF) model. Specifically, the input data included the vertical profiles of temperature, specific humidity and pressure, the surface parameters of surface skin temperature, and the 2-m wind speed and directions. Along with the atmospheric and surface conditions, the sensor's (MWHTS) zenith and scan angles were used as additional input to the Atmospheric Radiative Transfer Simulator (ARTS) [20]. The NCEP/WRF data field were shaved at a horizontal resolution 0.05 • × 0.05 • and 29 vertical levels, with the highest vertical level located near 0.1 hPa.
In ARTS, the atmosphere is gridded in latitude, longitude, and pressure coordinates. Because radiative transfer calculations with scattering are generally much more expensive than clear sky radiative transfer, scattering calculations are confined to a sub-domain of the modeled atmosphere called "cloud box". In this research, we considered spherical and cylindrical ice particles, with either completely random orientation or horizontally aligned with random azimuthal orientation. The scattering properties of individual particles were derived using the random and fixed matrix code described by Mishchenko [21]. Because of the computational expense of the T-matrix calculations, single-scattering properties were pre-calculated, imported in XML format, and interpolated as required.
Calibration error analysis and quality control of MWHTS observations are critical for both profile retrievals and radiance assimilation. Before the launch of FY-3C, the intensive thermal vacuum (T/V) tests for MWHTS were carried out in a 2-m T/V chamber, and the basic parameters were obtained, such as the receiver nonlinearity and sensitivity of each channel of MWHTS [22]. Furthermore, system nonlinear error correction and bias correction of warm and cold targets were derived after the T/V data analysis and were used for calibration processing, which played an important role in producing Level 1 and Level 2 data. The calibration accuracy consisted mainly of contributions from four aspects: uncertainty of cold target, warm target, nonlinear term and the random noise of receiver. Table 2 shows the specifications of calibration accuracy for different channels. The interested reader is also referred to the work of He et al. [4] for more details about the calibration procedure. Table 2. Primary calibration results of MWHTS: ∆T W and ∆T C are uncertainty of warm and cold target, respectively; ∆T NL is uncertainty caused by nonlinearity; ∆T SYS is uncertainty from system noise; and T CAL is calibration uncertainty with system noise considered.

Channel
∆T W (K)

Clustering Method
The atmospheric states during hurricane development are not unified or stationary, and transition from one to another state appears frequently in long periods within the physical balance. As such, classification of atmospheric state is a key work before subsequent retrieving and assimilation analysis. We classified the atmospheric states in inner domain using a K-means clustering method [23] based on Integrated Water Path (IWP) integration estimated by the trapezoidal rule, and peak of height and width of the cloud mixing ratio. In particular, five classes were identified and labeled according to the hydrometer distribution and density. For illustration purposes, Figure 3 presents the maps of peak of Remote Sens. 2019, 11, 896 6 of 18 height and width of cloud mixing ratio and IWP (g/m 2 ) at 1200 UTC on 3 October 2016. Figure 4 shows the corresponding maps of pressure, rain density distribution, and k-means clustering results based on these three observations. Using the same criterion and algorithm, the four-day period of hurricane Matthew was classified into five categories at each 6-h interval for further retrieving and assimilation.
Remote Sens. 2019, 11, x FOR PEER REVIEW 6 of 18 purposes, Figure 3 presents the maps of peak of height and width of cloud mixing ratio and IWP (g/m 2 ) at 1200 UTC on 3 October 2016. Figure 4 shows the corresponding maps of pressure, rain density distribution, and k-means clustering results based on these three observations. Using the same criterion and algorithm, the four-day period of hurricane Matthew was classified into five categories at each 6-h interval for further retrieving and assimilation.

Retrieval Methodology
In satellite remote sensing, both statistical and physical retrieval methodologies are applied to deal with the regression problems and retrieve the atmospheric temperature profiles [12]. From the statistical point of view, the Neural Network (NN) algorithm is a typical approach used in the community [24], while the one-dimensional variational (1D-var) method is commonly used as a representative physical approach [25]. Both can obtain the nonlinear relationship between the satellite observation and the atmospheric profiles. In particular, the NN, which generally includes a few key steps, namely initialization, training, optimization and simulation of neural networks, is a powerful tool to reflect the nonlinear relation between the input and target spaces. Using the Back Propagation (BP) NN based on the error (BP) learning algorithm, one can learn and establish functions when the analytical relationships between inputs and outputs are complicated. For the 1D-var algorithm, it mainly includes two parts under the general term of physical concept: the radiative transfer model for brightness temperature simulations and the minimization of the cost function used for weighting the relative contribution of background information and satellite observations (see the right part of Equation (1)).
where x b is the background state variable; B is the background covariance matrix; R is the sum of the sensor noise and the covariance error in the radiance simulations; T represents the matrix transpose; H is the forward operator that simulates the satellite observations at the atmospheric state variable x; and O stands for the observations. Assuming the errors in the observations and the background information are neither biased nor correlated Gaussian distributions and there is a local linearity around x, the minimization of the cost function can be solved by: This results in the following solution to x, where H(x) is the tangent linear function of H at point x; n is the iteration index. From Equation (3), the final solution x n+1 is affected by the a priori information, including the background covariance matrix B, the background state variable x b , the initial guess, the bias (O − H(X)) between observation and simulation, and the observation covariance matrix R.
In this study, a combined methodology based on the above two methods was proposed and validated through the case study for hurricane Matthew. For background error covariance matrix B, we calculated it using the atmospheric temperature and humidity profiles within 20-40N and 60-80W from 28 September to 10 October 2016, where and when hurricane Matthew was formed.
The formula of covariance is given by: where σ 2 ij is arbitrary element in the covariance matrix B at row i and column j, respectively; and N is the number of atmospheric profiles used in this study. We took the mean of the atmospheric temperature and humidity profiles used in Equation (4) to calculate the background error covariance Bx b in the retrieval system. The initial guess of x n when n = 1 (i.e., x 1 ) in Equation (3) was obtained from NWP model outputs.
The brightness temperatures simulated using ARTS model were taken as the input vector O L . The length L of the input vector is 15 for all the MWHTS sounding channels. The atmospheric temperature and humidity profiles extracted from NWP model were taken as the output vector X M , and the length M of the output vector equalled to 32 pressure levels for profiles. The pairs of input/output vector were then used to train the BP NNs, and the steepest descent method was adopted in the training process. In particular, 90% of the random pairs (training data) were used to determine the weights and biases in the NN model, and the other 10% (testing data) of pairs were used to validate and evaluate the trained model. With detailed comparison and analysis, the number of hidden nodes in the hidden layer was determined as 28 for MWHTS in this study.
Before and after the launch of MWHTS, a series of tests were carried out, such as sensitivity, nonlinear correction factor, calibration accuracy and Allan deviation between observations and simulation model for each channel. In this study, all of these factors were considered when computing the observation error covariance matrix R. With the assumption that the measurements in one channel are uncorrelated to those from others, the diagonal elements of R can be expressed as: where r is the square root of the diagonal elements of R; f is the instrument channel noise (including all factors that cause calibration inaccuracy); and e is the forward model error.
In an actual case, the errors in the observations and the background information are correlated with bias, thus it no longer works for Equation (2). The convergence criterion adopted in this study is: When the relative difference of the cost function within two iterations was less than 0.01, the iteration process was terminated with a reasonable result. Otherwise, if the maximum iteration number reached 100, we assumed a mathematical solution failure occurred for the case.

Study Domain and WRF Configuration
As shown in Figure 5, the study area is located within 20-40N and 50-80W over the North Atlantic Ocean. For demonstration purpose, we selected the hurricane Matthew event for case study. Hurricane Matthew was a powerful and devastating tropical cyclone, which became the first Category 5 Atlantic hurricane since Hurricane Felix in 2007. This event evolved from 28 September to 10 October 2016, wrought widespread destruction and catastrophic loss of life during its journey across the Western Atlantic, including parts of Haiti, Cuba, Dominican Republic, the Lucayan Archipelago, the southeastern United States and Canadian Maritimes.
Numerous studies focus on various forecasting methods, including numerical model of dynamic prediction, objective statistical forecast, and subjective experience forecast. As mentioned in Section 3.1, this study used NCEP 6-h reanalysis data as input to WRF. The WRF model configuration consisted of a 15-km resolution parent domain and a 5-km nested grid with 32 vertical levels from the surface to 100 millibars. The model was configured with a Lambert Conformal (LCC) projection, centered at 29.25N, −65.878W. Figure 5 shows the modeling domains for the WRF simulations with the inner nested domains focused over the Atlantic Ocean. The Kain-Fritsch cumulus parameterization scheme was used on the 15-km domain only, while convection was explicitly simulated on the 5-km nested domains. The WSM6 microphysics scheme, the Yonsei University (YSU) PBL scheme, and the Noah land surface model were implemented for both domains. The NCEP 1-degree final analysis (FNL) data were used for the initial and lateral boundary conditions along with the NCEP Real Time Gridded 1/12-degree sea surface temperature data. Three-dimensional data assimilation was applied in all simulations using the FNL gridded analysis data. By using the WRF model, the evolution process and thermodynamic and hydrometric features of tropical cyclone with the highest resolution of 5 km were obtained and integrated to 72 h (three days).

Data Matching
The satellite observations used for retrieval in this study were collected from 28 September to 10 October 2016. The Level 1 data of FY-3C meteorological satellite are recorded once every 102 min for each orbit. The data can be obtained from http://satellite.nsmc.org.cn/portalsite/Data/FileShow.aspx. The required WRF inputs were obtained from the NCEP archival database. These inputs were: (1) 6h reanalysis data for atmospheric profiles; (2) 1-degree final analysis data for initial and lateral boundary conditions; and (3) real-time gridded data for sea surface temperature [19]. Each WRF forecast started with the initial and boundary conditions provided by NCEP reanalyses, creating 24 forecasts (with time interval of 15 min) at 5 km spatial grids for each initial atmospheric state. Hence, each pixel of MWHTS observations could be matched in the chosen domain within 30 min.
Considering that the invalid pixels should be deleted in the retrieval process, prior to retrieving, the first step Was to carry out quality control and screening. In particular, the footprint of MWHTS was flagged as invalid if brightness temperature for that pixel is less than 50 K or greater than 400 K. For data matching, the observations were screened if latitude is greater than 40N or smaller than 20N, or if the longitude is smaller than 50W or larger than 80W. The data from FY-3C satellite and WRF runs were matched in both spatial and temporal perspectives. The matched datasets were then used for further retrievals. Figure 6 illustrates sample matching data for this particular study case.
In addition, since the Level 1 datasets of FY-3C MWHTS are archived with surface type (ocean and land) and the Level 2 dataset includes rain detection product with same view in time and space, we considerrf both the surface type and rain and no-rain conditions in data assimilation for the following retrieving. In total, after matching the FY-3C MWHTS and WRF output data, 502,400 matched data pairs were obtained.

Data Matching
The satellite observations used for retrieval in this study were collected from 28 September to 10 October 2016. The Level 1 data of FY-3C meteorological satellite are recorded once every 102 min for each orbit. The data can be obtained from http://satellite.nsmc.org.cn/portalsite/Data/FileShow.aspx. The required WRF inputs were obtained from the NCEP archival database. These inputs were: (1) 6-h reanalysis data for atmospheric profiles; (2) 1-degree final analysis data for initial and lateral boundary conditions; and (3) real-time gridded data for sea surface temperature [19]. Each WRF forecast started with the initial and boundary conditions provided by NCEP reanalyses, creating 24 forecasts (with time interval of 15 min) at 5 km spatial grids for each initial atmospheric state. Hence, each pixel of MWHTS observations could be matched in the chosen domain within 30 min.
Considering that the invalid pixels should be deleted in the retrieval process, prior to retrieving, the first step Was to carry out quality control and screening. In particular, the footprint of MWHTS was flagged as invalid if brightness temperature for that pixel is less than 50 K or greater than 400 K. For data matching, the observations were screened if latitude is greater than 40N or smaller than 20N, or if the longitude is smaller than 50W or larger than 80W. The data from FY-3C satellite and WRF runs were matched in both spatial and temporal perspectives. The matched datasets were then used for further retrievals. Figure 6 illustrates sample matching data for this particular study case.
In addition, since the Level 1 datasets of FY-3C MWHTS are archived with surface type (ocean and land) and the Level 2 dataset includes rain detection product with same view in time and space, we considerrf both the surface type and rain and no-rain conditions in data assimilation for the following retrieving. In total, after matching the FY-3C MWHTS and WRF output data, 502,400 matched data pairs were obtained. Remote Sens. 2019, 11, x FOR PEER REVIEW 11 of 18 Figure 6. Schematic of matching between FY-3C MWHTS and WRF output (for 150 GHz channel). The matching threshold is 0.25° for latitude and longitude and 30 min for time. The matched pixels are based on MWHTS since the WRF output is gridded with high temporal and spatial resolution as the model configured.

Retrievals and Analysis
In this study, the brightness temperatures were simulated by ARTS model and validated by corrected brightness temperatures of MWHTS with observation error covariance added to matrix R. The atmospheric state vectors of temperature and humidity profiles were retrieved by the combined algorithm with a prior background error covariance matrix B estimated by ANN neural network based on historical datasets in the same domain.

Retrievals and Analysis
In this study, the brightness temperatures were simulated by ARTS model and validated by corrected brightness temperatures of MWHTS with observation error covariance added to matrix R. The atmospheric state vectors of temperature and humidity profiles were retrieved by the combined algorithm with a prior background error covariance matrix B estimated by ANN neural network based on historical datasets in the same domain.

Retrievals and Analysis
In this study, the brightness temperatures were simulated by ARTS model and validated by corrected brightness temperatures of MWHTS with observation error covariance added to matrix R. The atmospheric state vectors of temperature and humidity profiles were retrieved by the combined algorithm with a prior background error covariance matrix B estimated by ANN neural network based on historical datasets in the same domain.   Simulation datasets using ARTS were calculated, as shown in Figure 8, where the input data included three-dimensional atmospheric hydrometeor variables (e.g., cloud liquid water, cloud ice, graupel, rain, and snow mixing ratios) and thermodynamic state (e.g., temperature and water vapor density) variables at 460 mb from WRF model runs with NCEP FNL data.
Biases were calculated to evaluate the retrieval results and quantify the associated errors with respect to the WRF outputs and NWP model reanalysis, which were assumed to be truths. According to the channel specification of the MWHTS, the temperature and humidity retrievals were for 32 pressure levels from 1000 to 10 hPa, and totally nine sampled levels were chosen, as illustrated in Figure 9. The corresponding distribution of retrievals errors including root mean square error (RMSE) are shown in Figure 10.
For calculation of the background error covariance, we used NNs method and validated the retrievals with reanalysis from NWP model. Combining the NNs and the 1D-var retrieval algorithm, the temperature and water vapor density retrieval validation results are shown in Figure 11. Simulation datasets using ARTS were calculated, as shown in Figure 8, where the input data included three-dimensional atmospheric hydrometeor variables (e.g., cloud liquid water, cloud ice, graupel, rain, and snow mixing ratios) and thermodynamic state (e.g., temperature and water vapor density) variables at 460 mb from WRF model runs with NCEP FNL data.
Biases were calculated to evaluate the retrieval results and quantify the associated errors with respect to the WRF outputs and NWP model reanalysis, which were assumed to be truths. According to the channel specification of the MWHTS, the temperature and humidity retrievals were for 32 pressure levels from 1000 to 10 hPa, and totally nine sampled levels were chosen, as illustrated in Figure 9. The corresponding distribution of retrievals errors including root mean square error (RMSE) are shown in Figure 10.
For calculation of the background error covariance, we used NNs method and validated the retrievals with reanalysis from NWP model. Combining the NNs and the 1D-var retrieval algorithm, the temperature and water vapor density retrieval validation results are shown in Figure 11. According to the retrieval results during hurricane Matthew, it is concluded that MWHTS can be used to retrieve temperature and humidity at levels from surface to 100 hPa due to the weighting function distribution of multi-channel observation. This also offers a possibility for improving the retrieval accuracy of temperature and humidity profiles in the atmospheric layers using the MWHTS observations.
The retrieved temperature profiles from the combined algorithm show that the retrieval accuracy of MWHTS at levels from surface to 100 hPa was evenly distributed with a temperature bias less than 0.75 K, which is consistent with the peaks of 118 GHz channels. For the water vapor profile, the accuracy of surface and 750-550 was worse than other levels, and, due to the distribution of water vapor in troposphere, the upper levels had relatively smaller root mean square errors than the lower layers. The comparison study on the retrieved atmospheric temperature and water vapor using satellite observations provided the reference value for assessing the detection ability of 118.75 GHz oxygen line. According to the retrieval results during hurricane Matthew, it is concluded that MWHTS can be used to retrieve temperature and humidity at levels from surface to 100 hPa due to the weighting function distribution of multi-channel observation. This also offers a possibility for improving the retrieval accuracy of temperature and humidity profiles in the atmospheric layers using the MWHTS observations.
The retrieved temperature profiles from the combined algorithm show that the retrieval accuracy of MWHTS at levels from surface to 100 hPa was evenly distributed with a temperature bias less than 0.75 K, which is consistent with the peaks of 118 GHz channels. For the water vapor profile, the accuracy of surface and 750-550 was worse than other levels, and, due to the distribution of water vapor in troposphere, the upper levels had relatively smaller root mean square errors than the lower layers. The comparison study on the retrieved atmospheric temperature and water vapor using satellite observations provided the reference value for assessing the detection ability of 118.75 GHz oxygen line. Remote Sens. 2019, 11, x FOR PEER REVIEW 13 of 18

Data Assimilation with Retrievals and Radiance
We also conducted the following experiments: (1) Control trial (CTRL): Basic assimilation trial in which no radiance data from polar-orbiting instruments were assimilated. Assimilated observations were therefore only conventional observations. (2) Radiance trial (DA_MWHTS_118GHz and DA_ MWHTS_118GHz): Same as CTRL but with MWHTS 118 GHz and 183 GHz observations from FY-3C satellite added separately. (3) Retrievals trial (DA_ temperature and humidity profiles): Same as CTRL but with MWHTS retrievals added.
To explore the application of microwave sounding data in numerical prediction of typhoons and improve typhoon forecast, we assimilated data directly for numerical forecasting of the track and intensity for the case study of hurricane Matthew based on the WRF-3DVar system. The forecasted tracks were compared with other satellite observations. Results show that the initial fields of the numerical model were improved much more due to direct assimilation of FY-3C MWHTS observations than the assimilation of conventional observations alone. The model initial fields were

Data Assimilation with Retrievals and Radiance
We also conducted the following experiments: (1) Control trial (CTRL): Basic assimilation trial in which no radiance data from polar-orbiting instruments were assimilated. Assimilated observations were therefore only conventional observations. (2) Radiance trial (DA_MWHTS_118GHz and DA_ MWHTS_118GHz): Same as CTRL but with MWHTS 118 GHz and 183 GHz observations from FY-3C satellite added separately. (3) Retrievals trial (DA_ temperature and humidity profiles): Same as CTRL but with MWHTS retrievals added.
To explore the application of microwave sounding data in numerical prediction of typhoons and improve typhoon forecast, we assimilated data directly for numerical forecasting of the track and intensity for the case study of hurricane Matthew based on the WRF-3DVar system. The forecasted tracks were compared with other satellite observations. Results show that the initial fields of the numerical model were improved much more due to direct assimilation of FY-3C MWHTS observations than the assimilation of conventional observations alone. The model initial fields were more reasonable in reflecting the initial situation of typhoon circulation as well as temperature and humidity conditions, and typhoon central position was also adjusted. Therefore, the FY-3C microwave data could efficiently improve the numerical prediction of typhoons, and radiance assimilation showed more improvement than retrievals. Combining the simulation and observation of satellite-based microwave payloads, the prediction of hurricane path compared with the actual path from real monitoring, and the results are shown in Figure 12.
Remote Sens. 2019, 11, x FOR PEER REVIEW 16 of 18 more reasonable in reflecting the initial situation of typhoon circulation as well as temperature and humidity conditions, and typhoon central position was also adjusted. Therefore, the FY-3C microwave data could efficiently improve the numerical prediction of typhoons, and radiance assimilation showed more improvement than retrievals. Combining the simulation and observation of satellite-based microwave payloads, the prediction of hurricane path compared with the actual path from real monitoring, and the results are shown in Figure 12.

Conclusion and future work
Theoretically, satellite radiance assimilation is superior to the retrieval assimilation because the uncertainties in the radiance observations are much smaller and more justified in direct radiance assimilation than in retrieval assimilation. This study demonstrated that a 118-GHz aircraft or satellite observations could provide useful information on the vertical thermal structure of the atmosphere. In addition, combined with the 183 GHz and window channels, water vapor profiles could be accurately retrieved and used in the data assimilation model, which can validate the positive impacts on forecast of extreme weather such as the hurricanes. This study also provided incremental gains in the analysis and forecast accuracy and brought new capacity to the 118 GHz channel through the deliverable of new cloud ice information in data assimilation model over only assimilating the 183 GHz radiance.
Nevertheless, accurate prediction of hurricane track and intensity change is still a challenging problem. Due to the lack of in situ data for hurricane initialization, initial vortices provided by largescale analysis from operational centers are often ill-defined, too weak, and sometimes misplaced. The methodology presented in this study is a new attempt to provide more confidence on the satellitebased atmospheric profile applications in hurricane domains. In addition to the atmospheric temperature and humidity profiles, precipitation, as one of the MWHTS L2 products, will also be an available choice to improve the hurricane storm track density. For the study case, using microwave radiance of MWHTS is recommended, which can improve the forecasting track by 50 km combining atmospheric temperature and humidity information.

Acknowledgments:
The authors would like to express their gratitude to the anonymous reviewers for their comments that improved this manuscript.

Conflicts of Interest:
The authors declare no conflicts of interest.

Conclusions and Future Work
Theoretically, satellite radiance assimilation is superior to the retrieval assimilation because the uncertainties in the radiance observations are much smaller and more justified in direct radiance assimilation than in retrieval assimilation. This study demonstrated that a 118-GHz aircraft or satellite observations could provide useful information on the vertical thermal structure of the atmosphere. In addition, combined with the 183 GHz and window channels, water vapor profiles could be accurately retrieved and used in the data assimilation model, which can validate the positive impacts on forecast of extreme weather such as the hurricanes. This study also provided incremental gains in the analysis and forecast accuracy and brought new capacity to the 118 GHz channel through the deliverable of new cloud ice information in data assimilation model over only assimilating the 183 GHz radiance.
Nevertheless, accurate prediction of hurricane track and intensity change is still a challenging problem. Due to the lack of in situ data for hurricane initialization, initial vortices provided by large-scale analysis from operational centers are often ill-defined, too weak, and sometimes misplaced. The methodology presented in this study is a new attempt to provide more confidence on the satellite-based atmospheric profile applications in hurricane domains. In addition to the atmospheric temperature and humidity profiles, precipitation, as one of the MWHTS L2 products, will also be an available choice to improve the hurricane storm track density. For the study case, using microwave radiance of MWHTS is recommended, which can improve the forecasting track by 50 km combining atmospheric temperature and humidity information.