Profile and Precipitation Retrievals and Validation Based on Geostationary SubMillimeter Atmospheric Sounder

The key factors of temporal and spatial resolution for typhoon observation played an important role in the design of radiometer used for observing typhoon. The NCEP (National Centres for Environmental Prediction) operational global analysis data prepared operationally every six hours were used as the initial field for mesoscale weather research and forecasting model (WRF) and drove the model to output atmospheric parameters such as hydrometeor content, temperature and humidity profiles at different time, which were inputs for the Atmospheric Radiative Transfer Simulator (ARTS) to calculate brightness temperature observed from geostationary earth orbit at oxygen absorption and water absorption band. The atmospheric humidity and temperature profiles of typhoon domain were retrieved from geostationary sub-millimetre atmospheric sounder. The results show that the profile retrievals using BP-NN algorithm have a best agreement with those from radiosonde, which is less than 20% and 1 K of root mean square error, respectively. For precipitation rate retrievals, much better agreement with rain gauge and ECMWF datasets, the RMS is between 0.84 to 32.4 mm/h for sea surface 0.89 and 36.13 mm/h for land surface according to the classification by precipitation type.


Introduction
The current polar-orbiting meteorological satellite observation system ensures that the observation period is six How to cite this paper: He, J.Y., Zhang, S.W., Liu hours, which is unable to monitor typhoon and other fast changing situations.Only instruments on geostationary or comparable platforms can view regional disaster at the 15-min interval that is necessary to monitor rapidly evolving typhoon or cyclone events.This paper discusses the abilities of passive microwave/ sub-millimeter sensor which is being developed by our institute to retrieve atmospheric temperature and humidity profiles.Furthermore, we also carry out the work of retrieving surface precipitation rates and hydrometeor water paths.
Geostationary atmospheric sounder based on interferometric technology is the newest field of microwave remote sensing.Equipped with a sub-millimeter atmospheric remote sensing instruments on the geostationary orbit, meteorological satellite platform will increase the frequency of observations, at the same time, improve the cloud detection capability [1]- [3].

Instrument Description
According to the atmospheric sounding theory and gasous absorption coefficients from 0 -1000 GHz, the final frequencies are chosen to derive atmospheric temperature and humidity profiles for geostationary sub-millimeter sounder, which are list in Figure 1 and Table 1, Table 2 shows characteristics of ultimately selected frequencies and their applications.
Since atmospheric absorbing characteristics of water vapor and oxygen, using satellite-borne microwave radiometer to derive atmospheric temperature and humidity profiles is possible.According to the principle above,   when the radiometer operated at frequencies 50 -1000 GHz, the contribution of surface background noise can be reduced to negligible magnitude, because the transmittance of the atmosphere at these frequencies is approximately equal to 0. Atmospheric opacity thickness with pressure integration can be expressed as: ( ) ( ) where ( ) z α is atmospheric at altitude z, z' is height from surface to satellite.
For oxygen absorbing channels, ( ) z τ is basically due to the contribution of oxygen.Equation (1) can be expressed as the weighted integral of temperature: where, ( ) Furthermore, water vapor weighing function can be expressed as: Weighting function is the weight of atmospheric radiance at height z from surface to the height of the satellite.The plane parallel atmosphere was divided into N layers, the absorbing coefficient of atmospheric parameters was assumed uniform in each layer, and then attenuation contributions of entire atmosphere including the surface layer can be accumulated.Atmospheric absorbing coefficients in each channel can be calculated combining MPM93 [4] and PWR04 [5] model.Therefore, we can choose the central frequency and bandwidth for each channel of geostationary sub-millimeter atmospheric sounder according to the weighting functions, which is displayed in Figure 2.

Simulation of Full-Disk Model of Earth
Brightness temperatures are simulated by the validated global reference physical model, NCEP/WRF/ARTS, composed of the US National Center for Environment Prediction (NCEP) analyses, the new generation National Center for Atmospheric Research/Penn State Mesoscale Model (WRF) and the Atmospheric Radiative Transfer Simulator, ARTS, which is a software for performing simulations of atmospheric radiative transfer [6]- [9].The output of WRF model consists of temperature and humidity profiles, pressure profiles, surface mask, and surface parameters and so on.3 shows.Therefore the results show that: 1) the simulated total number of Typhoons is close to that from Best-track dataset; 2) The simulated intensity of typhoons is comparable to that from Best-track data set, but the time limit needs further exploration.

Observation from MWHTS
Until now, no observation data from geostationary satellite can be used to realize the profile and precipitation rates retrievals.To solve this problem, this paper uses the observing data from microwave humidity and temperature sounder onboard FY-3C (FY-3C MWHTS) which plays an important role in monitoring extreme climate, especially for typhoon since September 30 th , 2013, such as typhoon "kujira", "Linfa", "chan-hom" and "Nangka", monitoring their procedure of generating, evolution, strengthen and die out (as showing in Figure 4 and

Retrieval Algorithm
ANN is essentially a nonlinear statistical regression between a set of predictors (in this case the observation vectors X) and a set of predictands (in this case profiles of atmospheric temperature Z) [13].The structure of the ANN is shown in Figure 6.In this paper, we construct a three layer ANN model.The layers 1, 2, and 3 represent the input layer, the hidden layer, and the output layer, respectively.The neurons of the input layer are represented by vector X i (X 1 , X 2 , X 3 , •••, X L ), where L is the number of the input neurons.The neurons of the middle layer are represented by vector , where M is the number of the hidden neurons.The neurons of the output layer are represented by vector where N is the number of the output neurons.
The node in the hidden layer can be expressed as: ( ) where, ij w is the weighting of the connection between the j th hidden neuron and the i th input neuron and j b denotes the bias in the j th neuron of the hidden layer.The linear function is applied between the output layer and the hidden layer.Where, S denotes the sigmoid function: ( ) ( ) According to the characteristics of sigmoid function, the values of both input and output layers should be transformed to the range [0, 1].Outputs can be expressed as: where jk w is the weight of the connection between the j th hidden neuron and the k th output neuron; k b is the bias in the k th neuron of the output layer.
For geostationary sub-millimeter atmospheric sounder, several kinds of ANN are used according to the types of surface and sky, like land and sea, clear-sky, cloudy sky, rainy sky, typhoon-sky and so on.The schematic J. Y.He et al.

Retrieval Results and Analysis
Figure 8 shows the profile retrievals of temperature and humidity (T&H) using simulated brightness temperatures from microwave radiometer onboard gestational platform, which also gives the root mean square error distribution, respectively.For the temperature profile, when there is a thin inversion layer, large deviations occur from the neural network retrieval model.It is mainly because of its own shortcoming for nonlinear neural network.When the inversion layer is thick enough, the neural network inverse model can be well reflected for the details of atmospheric temperature changes.Water vapor varies significantly from time to time and from space to space.Certainly the fact is that water vapor is the source of all clouds and precipitation which would be enough to explain its retrieving difficulties.Therefore compared to temperature, water vapor and relative humidity profiles are retrieved with relatively larger difficulty with challenge.The retrievals show that the RMS of atmospheric temperature profile is less than 2.5 K, RMS of atmospheric relative humidity profile is better than 20%, which can be converted to atmospheric absolute humidity, and the RMS is less than 0.4 g/m 3 as shown in Figure 8.
For precipitation the primary radiometric signal at frequencies around 183 GHz from precipitating scenes results from the scattering by ice hydrometeors.This scattering can result in significant brightness temperature  depressions (several 10's K) relative to non-precipitating surroundings, and is therefore a sensitive proxy for the presence of precipitation at the surface.
According to the global difference distributions of brightness temperatures on Jan 20, 2015, a) channel 10 plus channel 15, b) for channel 10 plus channel 9, and c) for channel 9 plus channel 15 for FY-3C MWHTS, using neural network method, the precipitation and rain detection can be derived, which is shown in Figure 9 and Table 3.

Summary and Conclusions
According to the retrievals and analysis, geostationary sub-millimeter atmospheric sounder will play an important role in studying global climate and is the main remote sensing instrument for meteorology and disaster.It works in all weather and all day providing the observation of brightness temperature which can be used to retrieve temperature and humidity profiles and precipitation rate.
There is a well agreement in temperature and humidity profiles between radiosonde and retrievals.Compared to radiosonde, the retrievals show that the RMS of atmospheric temperature profile is less than 1 K, RMS of atmospheric relative humidity profile is less than 20%, which can be converted to atmospheric absolute humidity,  and the RMS is less than 0.4 g/m 3 .For the recent work, the surface is classified the surface mark as land, sea and coastal, using the observing data to detect weather it is rain or not and then test and validate the accuracy of the rain rate.Therefore, the work provides algorithms and data analysis of temperature and humidity profiles and precipitation distribution.So, the work will play an important role in the design and development of following meteorological satellites.
Because it is not enough, the authors are doing further improvement, like considering the surface covered by snow, ice and rainforest.Also, the rain rate is excepted to be retrieved more accurately.The radar data are also needed to validate the accuracy of rain detection.All of above will be described in future paper.
J. Y.He et al.

Figure 4 .
Figure 4. Brightness temperatures distribution at different pressure-level at 118 GHz.

Figure 5 )
Figure 5), especially locating the typhoon eye area clearly and intuitively with the resolution of 15 kilometers, and also predicting the heavy rainfall caused by typhoon for South China [10]-[12].

Figure 9 .
Figure 9. Rain detection according to the brightness temperature analysis.

Table 2 .
Characteristics of ultimately selected frequencies.

Table 3 .
Rain detection of typhoon domain by geostationary sub-millimeter atmospheric sounder between Jan. 1 to Oct. 31, 2014.