Impact of FY-3D MWRI Radiance Assimilation in GRAPES 4DVar on Forecasts of Typhoon Shanshan

+ Author Affiliations + Find other works by these authors
  • Corresponding author: Wei HAN, hanwei@cma.gov.cn
  • Funds:

    Supported by the National Natural Science Foundation of China (41675108), National Key Research and Development Program (2018YFC1506700), and Second Tibetan Plateau Scientific Expedition and Research Program (2019QZKK0105)

  • doi: 10.1007/s13351-020-9122-x

PDF

  • In this study, Fengyun-3D (FY-3D) MicroWave Radiation Imager (MWRI) radiance data were directly assimilated into the Global/Regional Assimilation and PrEdiction System (GRAPES) four-dimensional variational (4DVar) system. Quality control procedures were developed for MWRI applications by using algorithms from similar microwave instruments. Compared with the FY-3C MWRI, the bias of FY-3D MWRI observations did not show a clear node-dependent difference from the numerical weather prediction background simulation. A conventional bias correction approach can therefore be used to remove systematic biases before the assimilation of data. After assimilating the MWRI radiance data into GRAPES, the geopotential height and humidity analysis fields were improved relative to the control experiment. There was a positive impact on the location of the subtropical high, which led to improvements in forecasts of the track of Typhoon Shanshan.
  • 加载中
  • Fig. 1.  Mean (O–B) of each channel for the ascending (curve with solid square) and descending (curve with empty circle) data for the FY-3D MWRI. Statistics were calculated for the 2-day data after applying the quality control procedures described in Section 3.2.

    Fig. 2.  Flow chart showing the experiment setup used in this study. The experiment contains a 1-day cycle assimilation module and a 5-day forecast module. The abbreviation “AN” represents an analysis field.

    Fig. 3.  OMB probability density function (PDF) profiles of the (a) 19-V, (b) 24-V, and (c) 37-V channels. The shadow/black bars correspond to the PDF of the OMB before/after the quality control (QC) procedures.

    Fig. 4.  Results of the quality control experiments for the 37-V channel. (a) FY-3D MERSI image of Typhoon Shanshan on 5 August 2018. (b) The pixels contaminated by oceanic clouds, which are detected by the cloud detection scheme in Section 3.2, and the pixels with LWP > 0.1 mm are plotted. (c) Image of the brightness temperature (Tb) in the MWRI 37-V channel at 0600 UTC 5 August within a ±3-h time window from the observations without quality control. (d) The remaining pixels after quality control of the data in part (c). (e) The background Tb (K) corresponding to part (c). (f) The OMB (K) corresponding to part (d).

    Fig. 5.  Difference of Tb between the 19-V and 24-V channels around North America on 4 August 2018.

    Fig. 6.  OMB PDF profiles of the (a) 19-V, (b) 24-V, and (c) 37-V channels. The shadow/black bars correspond to the PDF of the OMB before/after the bias correction (BC) procedures. The solid and dashed lines are the fitted normal distribution profiles calculated by the black and shadow bars, respectively.

    Fig. 7.  Specific humidity analysis increments in the (a, b) 850-hPa and (c, d) 23°N sections from the (a, c) EXP3 and (b, d) CTL3 experiments.

    Fig. 8.  The 120-h typhoon track forecasts given by (a) CTL1 and EXP1, (b) CTL2 and EXP2, and (c) CTL3 and EXP3. The forecasts given by the control and sensitivity experiments are plotted as squares and triangles, respectively. The observed typhoon track is represented by circles. The time is marked along the track for each 24 h.

    Fig. 9.  The subtropical high represented by the contour line of 588 dagpm at 500 hPa given by (a) observational results in MICAPS (Meteorological Information Comprehensive Analysis and Process System) and (b) the forecast results in CTL3 and EXP3 at 72 h after forecast initialization. In (a), the Asian low-pressure trough is marked by the red arrowheads. In (b), the contour line of 587.1 dagpm at the same height and same time is also plotted because of some unavoidable deviation; the results of CTL3 and EXP3 are shown by dashed green line and solid blue line, respectively.

    Table 1.  FY-3D MWRI channel characteristics (see http://satellite.nsmc.org.cn/PortalSite/StaticContent/DeviceIntro_FY-3_MWRI.aspx)

    Channel characteristicFrequency (GHz)
    10.6518.723.836.589.0
    PolarizationVHVHVHVHVH
    Abbreviation10-V10-H19-V19-H24-V24-H37-V37-H89-V89-H
    Bandwidth (MHz)1802004004003000
    NEΔT (K)0.50.50.50.50.8
    IFOV (km2)51 × 8530 × 5027 × 4518 × 309 × 15
    Pixel (km2)40 × 11.240 × 11.220 × 11.220 × 11.210 × 11.2
    Range (K)3–340
    Number of scan positions266
    Scanning techniqueConical
    Swath width (km)1400
    Zenith angle (°)45 ± 0.1
    Incident angle (°)53.1
    Scan rate (s)1.8 ± 0.1
    Note: NEΔT: noise-equivalent brightness temperature (Tb).
    Download: Download as CSV

    Table 2.  Coefficients ${a_0}$, ${a_1}$, and ${a_2}$ in Eq. (6) for the MWRI

    LWPchana0a1a2
    10-V−3.87 4.470.09
    19-V−1.94 3.030.37
    37-V−0.97 2.740.39
    89-H−0.37−3.081.68
    Download: Download as CSV

    Table 3.  Threshold of the LWP corresponding to the three channels assimilated in this study

    Channel10-V19-V37-V
    LWP (mm)0.30.250.1
    Download: Download as CSV
  • [1]

    Alishouse, J. C., S. A. Snyder, J. Vongsathorn, et al., 1990: Determination of oceanic total precipitable water from the SSM/I. IEEE Trans. Geosci. Remote Sens., 28, 811–816. doi: 10.1109/36.58967.
    [2]

    Arakawa, A., and W. H. Schubert, 1974: Interaction of a cumulus cloud ensemble with the large-scale environment, Part I. J. Atmos. Sci., 31, 674–701. doi: 10.1175/1520-0469(1974)031<0674:IOACCE>2.0.CO;2.
    [3]

    Bao, Y. S., F. Mao, J. Z. Min, et al., 2014: Retrieval of bare soil moisture from FY-3B/MWRI data. Remote Sens. Land Resour., 26, 131–137. (in Chinese) doi: 10.6046/gtzyyg.2014.04.21.
    [4]

    Bettenhausen, M. H., C. K. Smith, R. M. Bevilacqua, et al., 2006: A nonlinear optimization algorithm for WindSat wind vector retrievals. IEEE Trans. Geosci. Remote Sens., 44, 597–610. doi: 10.1109/TGRS.2005.862504.
    [5]

    Bouttier, F., and P. Courtier, 2002: Data assimilation concepts and methods March 1999. Proceedings of ECMWF Meteorological Training Course Lecture Series, ECMWF, Bracknell, 1–58.
    [6]

    Chen, D. H., J. S. Xue, X. S. Yang, et al., 2008: New generation of multi-scale NWP system (GRAPES): General scientific design. Chinese Sci. Bull., 53, 3433–3445. doi: 10.1007/s11434-008-0494-z.
    [7]

    Chen, H., and Y. Q. Jin, 2012: In-orbit intercalibration of FY-3B/MWRI and applications for monitoring drought and flooding. J. Remote Sens., 16, 1024–1034. (in Chinese) doi: 10.11834/jrs.20121299.
    [8]

    Chen, L. S., 1979: On the causal analysis of typhoon tracks which turning direction westward suddenly over the sea area near the eastern China. Chinese J. Atmos. Sci., 3, 289–298. (in Chinese) doi: 10.3878/j.issn.1006-9895.1979.03.11.
    [9]

    Chen, X. M., Q. J. Liu, and J. C. Zhang, 2007: A numerical simulation study on microphysical structure and cloud seeding in cloud system of Qilian Mountain Region. Meteor. Mon., 33, 33–43. (in Chinese) doi: 10.3969/j.issn.1000-0526.2007.07.004.
    [10]

    Connor, L. N., and P. S. Chang, 2000: Ocean surface wind retrievals using the TRMM microwave imager. IEEE Trans. Geosci. Remote Sens., 38, 2009–2016. doi: 10.1109/36.851782.
    [11]

    Dai, Y. J., X. B. Zeng, R. E. Dickinson, et al., 2003: The common land model. Bull. Amer. Meteor. Soc., 84, 1013–1024. doi: 10.1175/BAMS-84-8-1013.
    [12]

    Dou, F. L., D. W. An, and J. R. Li, 2014: Sea surface wind speed retrieval based on FY-3B Microwave Imager. Remote Sens. Technol. Appl., 29, 984–992. (in Chinese).
    [13]

    Feng, C. C., and H. Zhao, 2015: Identification of radio-frequency interference signal from FY-3B microwave radiation imager over ocean. J. Remote Sens., 19, 465–475. (in Chinese) doi: 10.11834/jrs.20154056.
    [14]

    Ferraro, R. R., F. Z. Weng, N. C. Grody, et al., 1996: An eight-year (1987–1994) time series of rainfall, clouds, water vapor, snow cover, and sea ice derived from SSM/I measurements. Bull. Amer. Meteor. Soc., 77, 891–906. doi: 10.1175/1520-0477(1996)077<0891:AEYTSO>2.0.CO;2.
    [15]

    Gaiser, P. W., K. M. St Germain, E. M. Twarog, et al., 2004: The WindSat spaceborne Polarimetric microwave radiometer: Sensor description and early orbit performance. IEEE Trans. Geosci. Remote Sens., 42, 2347–2361. doi: 10.1109/TGRS.2004.836867.
    [16]

    Geer, A. J., K. Lonitz, P. Weston, et al., 2018: All-sky satellite data assimilation at operational weather forecasting centres. Quart. J. Roy. Meteor. Soc., 144, 1191–1217. doi: 10.1002/qj.3202.
    [17]

    Giorgi, F., Y. Huang, K. Nishizawa, et al., 1999: A seasonal cycle simulation over eastern Asia and its sensitivity to radiative transfer and surface processes. J. Geophys. Res. Atmos., 104, 6403–6423. doi: 10.1029/1998JD200052.
    [18]

    Grody, N. C., and R. R. Ferraro, 1992: A comparison of passive microwave rainfall retrieval methods. Proceeding of the 6th Conference on Meteorology and Oceanography, American Meteorological Society, Atlanta, 60–65.
    [19]

    Guo, L., H. Sheng, J. Wang, et al., 2017: Retrieving near sea surface air temperature by AMSR2 radiometer. Adv. Mar. Sci., 35, 124–130. (in Chinese) doi: 10.3969/j.issn.1671-6647.2017.01.013.
    [20]

    Han, W., and N. Bormann, 2016: Constrained adaptive bias correction for satellite radiance assimilation in the ECMWF 4D-Var system. Technical Memorandum No. 783, ECMWF, Shinfield Park, Reading, 26 pp.
    [21]

    Hargens, U., C. Simmer, and E. Ruprecht, 1992: Remote sensing of cloud liquid water during ICE’89. Proceedings of Specialist Meeting on Microwave Radiometry and Remote Sensing Applications, IEEE, Boulder, Colorado, 27–31.
    [22]

    Harris, B. A., and G. Kelly, 2001: A satellite radiance-bias correction scheme for data assimilation. Quart. J. Roy. Meteor. Soc., 127, 1453–1468. doi: 10.1002/qj.49712757418.
    [23]

    Hollinger, J. P., J. L. Peirce, and G. A. Poe, 1990: SSM/I instrument evaluation. IEEE Trans. Geosci. Remote Sens., 28, 781–790. doi: 10.1109/36.58964.
    [24]

    Hong, S. Y., and H. L. Pan, 1996: Nonlocal boundary layer vertical diffusion in a medium-range forecast model. Mon. Wea. Rev., 124, 2322–2339. doi: 10.1175/1520-0493(1996)124<2322:NBLVDI>2.0.CO;2.
    [25]

    Huang, W., Y. L. Hao, J. Wang, et al., 2013: Brightness temperature data comparison and evaluation of FY-3B microwave radiation imager with AMSR-E. Period. Ocean Univ. China, 43, 99–111. (in Chinese) doi: 10.16441/j.cnki.hdxb.2013.11.015.
    [26]

    JAXA, 2013: GCOM-W1 “SHIZUKU” Data Users Handbook, Japan Aerospace Exploration Agency. Tsukuba, Japan, 125 pp. Available online at https://gportal.jaxa.jp/gpr/assets/mng_upload/GCOM-W/GCOM-W1_SHIZUKU_Data_Users_Handbook_EN.pdf. Accessed on 19 August 2020.
    [27]

    Kawanishi, T., T. Sezai, Y. Ito, et al., 2003: The advanced microwave scanning radiometer for the Earth observing system (AMSR-E), NASDA’S contribution to the EOS for global energy and water cycle studies. IEEE Trans. Geosci. Remote Sens., 41, 184–194. doi: 10.1109/TGRS.2002.808331.
    [28]

    Kazumori, M., Q. H. Liu, R. Treadon, et al., 2008: Impact study of AMSR-E radiances in the NCEP global data assimilation system. Mon. Wea. Rev., 136, 541–559. doi: 10.1175/2007MWR2147.1.
    [29]

    Kazumori, M., A. J. Geer, and S. J. English, 2014: Effects of all-sky assimilation of GCOM-W1/AMSR2 radiances in the ECMWF system. Technical Memo 732, ECMWF, Reading, 1–34.
    [30]

    Krasnopolsky, V. M., L. C. Breaker, and W. H. Gemmill, 1995: A neural network as a nonlinear transfer function model for retrieving surface wind speeds from the special sensor microwave imager. J. Geophys. Res. Oceans, 100, 11,033–11,045. doi: 10.1029/95JC00857.
    [31]

    Kummerow, C., J. Simpson, O. Thiele, et al., 2000: The status of the Tropical Rainfall Measuring Mission (TRMM) after two years in orbit. J. Appl. Meteor., 39, 1965–1982. doi: 10.1175/1520-0450(2001)040<1965:TSOTTR>2.0.CO;2.
    [32]