CRA-40/Atmosphere—The First-Generation Chinese Atmospheric Reanalysis (1979–2018): System Description and Performance Evaluation

中国气象局第一代全球大气再分析CRA-40(1979–2018):系统介绍和性能评估

+ Author Affiliations + Find other works by these authors
  • Corresponding author: Lipeng JIANG, jianglp@cma.gov.cn; Chunxiang SHI, shicx@cma.gov.cn
  • Funds:

    Supported by the China Meteorological Administration Special Public Welfare Research Fund (GYHY201506002) and National Innovation Project for Meteorological Science and Technology (CMAGGTD003-5)

  • doi: 10.1007/s13351-023-2086-x

PDF

  • Atmospheric reanalysis reproduces the past atmospheric conditions through assimilation of historical meteorologi-cal observations with fixed version of a numerical weather prediction (NWP) model and data assimilation (DA) system. It is widely used in weather, climate, and even business-related research and applications. This paper reports the development of CMA’s first-generation global atmospheric reanalysis (RA) covering 1979–2018 (CRA-40; CRA refers to CMA-RA). CRA-40 is produced by using the Global Spectral Model (GSM)/Gridpoint Statistical Interpolation (GSI) at a 6-h time interval and a TL574 spectral (34-km) resolution with the model top at 0.27 hPa. A large number of reprocessed satellite data and widely collected conventional observations were assimilated during the reanalyzing process, including the reprocessed atmospheric motion vector (AMV) products from FY-2C/D/E/G satellites, dense conventional observations (at about 120 radiosonde and 2400 synoptic stations) over China, as well as MWHS-2 and GNSS-RO observations from FY-3C. The reanalysis fitting to observations is improved over time, especially for surface pressure with root-mean-square error reduced from 1.05 hPa in 1979 to 0.8 hPa, and for upper air temperature from 1.65 K in 1979 to 1.35 K, in 2018. The patterns of global analysis increments for temperature, specific humidity, and zonal wind are consistent with the changes in the observing system. Near surface temperature from the model’s 6-h forecast reflects the global warming trend reasonably. The CRA-40 precipitation pattern matches well with those of GPCP and other reanalyses. CRA-40 also successfully captures the QBO and its vertical and temporal development, hemispherical atmospheric circulation change, and moisture transport by the East Asian summer monsoon. CRA is now operationally running in near real time as a climate data assimilation system in CMA.
    大气再分析是利用固定的、最先进的数值预报模式和数据同化系统,以及各种来源的观测资料,对大气状况进行历史回算,得到尽可能稳定、均一和精确的长序列历史分析数据产品。本文主要介绍中国气象局第一代40年(1979年–准实时)全球大气再分析(CRA-40;注意本文CRA等同于CMA-RA)系统的系统配置和产品评估结果。CRA-40产品时间分辨率逐6小时,水平分辨率约34 km,模式层顶0.27 hPa。与国际上其他全球再分析产品相比,CRA-40在同化过程采用了更多的中国特色观测资料,比如中国风云2号系列卫星大气运动矢量重处理产品,中国特有常规观测资料,FY-3C微波湿度计和掩星观测资料,等。CRA-40降水产品的空间分布形态与GPCP降水分析产品及国际其他再分析降水产品比较一致;CRA-40近地面气温场能够成功反映出全球变暖趋势。另外,CRA-40产品能成功再现赤道平流层低层风场准两年振荡(QBO)、全球大气环流特征及东亚夏季风水汽输送特征等。目前,CRA-40系统已作为中国气象局业务系统准实时运行。
  • 加载中
  • Fig. 1.  Monthly counts of conventional observations for (a) surface pressure, (b) temperature, (c) specific humidity, and (d) wind used in CRA-40 (red line) and CFSR (black line) from 1979 to 2014.

    Fig. 2.  Timeline of conventional data, AMVs (i.e., SATWND), ocean surface winds, ozone, and GPSRO data used in CRA-40.

    Fig. 3.  Timeline of satellite radiance data using in CRA-40, including those from various microwave (MSU, AMSU-A, AMSU-B, MHS, MWHS-2, ATMS, and SSMIS) and infrared sensors (HIRS2–4, AVHRR3, AIRS, IASI, and CrIS) onboard various polar-orbiting satellites (e.g., n* for NOAA series, f17 for DMSP 17 satellite). Infrared radiances from geostationary GOES (g*) sounders are indicated by sndr*, while those from Meteosat-10 (m10) are indicated by SEVIRI.

    Fig. 4.  The 30-day mean (thin curves) and 360-day mean (thick curves) of the root-mean-square of first-guess (blue) (OMB) and analysis (red) (OMA) departures from observations in CRA-40 for (a) upper-air zonal wind from radiosondes, dropsondes, and PILOTs; (b) upper-air temperature from radiosondes and dropsondes; (c) upper-air specific humidity from radiosondes and dropsondes; and (d) surface pressure from SYNOP, buoys, ships, and METAR.

    Fig. 5.  Pressure–time diagrams of monthly and globally averaged analysis increments for (a) temperature (K), (b) zonal wind (m s−1), and (c) specific humidity (mg kg−1) on a logarithmic pressure scale.

    Fig. 6.  Pressure–time diagrams of three reanalyses’ monthly and globally averaged root-mean-square errors (RMSEs) for (a) geopotential height (gpm), (b) zonal wind (m s−1), (c) temperature (K), and (d) specific humidity (g kg−1) against ERA5. The RMSEs of CRA-40, CFSR, and JRA-55 are shown from top to bottom in each subplot.

    Fig. 7.  Twelve-month running averages of global-mean near surface air temperature (T2m) anomalies (K) from 1979 to 2019 for CRA-40/Atmosphere. The T2m anomaly is computed relative to the mean from 1981 to 2010.

    Fig. 8.  Thirty-year (1981–2010) averaged precipitation rates (mm day−1) in (a) CRA-40, (b) ERA5, (c) JRA-55, and (d) GPCP-v2.3.

    Fig. 9.  Time series of 360-day running average of global mean precipitation rates (mm day−1) from CRA-40 (red solid line), ERA5 (blue line), JRA-55 (green line), and GPCP-v2.3 (gray line). The red dash line is the corresponding evaporation from CRA-40.

    Fig. 10.  Time–pressure cross-section of zonal-mean of monthly-mean zonal wind (m s−1) averaged between 5°S and 5°N for (a) ERA5, (b) CRA-40, (c) JRA-55, and (d) CFSR.

    Fig. 11.  Longitude–height cross-section of summer (June, July, and August) zonal–vertical atmospheric circulation climatology (1980–2010) between 15° and 50°N (a) for ERA5 and the difference (b) between CRA-40 and ERA5, (c) between JRA-55 and ERA5, and also (d) between CFSR and ERA5.

    Fig. 12.  Forty-year mean of JJA water vapor transport (kg m−2 s−1) at 850 hPa (a) from ERA5, and (b) the differences between CRA-40 and ERA5, (c) between JRA-55 and ERA5, and (d) between CFSR and ERA5. The colors indicate the magnitude of the water vapor flux vector.

    Table 1.  Global atmospheric reanalyses

    GenerationReanalysisProducerPeriodResolutionMethodReference
    1stNASANASA1980–19952° × 2.5°OISchubert et al. (1993)
    R1NCEP + NCAR1948– T623DVarKalnay et al. (1996)
    ERA-15ECMWF1979–1994T106OIGibson et al. (1999)
    2ndR2NCEP + DOE1979– T623DVarKanamitsu et al. (2002)
    ERA-40ECMWF1957–2002TL1593DVarUppala et al. (2005)
    JRA-25JMA-CRIEPI1979– T1063DVarOnogi et al. (2007)
    3rdERA-InterimECMWF1979–2019TL2554DVarDee et al. (2011)
    CFSRNCEP1979– T3823DVarSaha et al. (2010)
    MERRANASA1979–20161/2° × 2/3°3DVarRienecker et al. (2011)
    MERRA-2NASA1980– 1/2° × 5/8°3DVarGelaro et al. (2017)
    JRA-55JMA1957– TL3194DVarKobayashi et al. (2015)
    4thERA5ECMWF1950– TL639Hybrid 4DVarHersbach et al. (2020)
    Download: Download as CSV

    Table 2.  Data sources of conventional observations for CRA-40

    Data supplierData type and suppliers’ identifierPeriod
    CMAChinese operational observationsSYNOP, radiosonde, and pilot1979–2018
    Buoy and ship2013–2018
    Wind profiler2007–2018
    Aircraft2003–2018
    Operational GTS dataSYNOP, radiosonde, aircraft, buoy, and ship1979–2018
    Intense observations of TIPEX III (Zhao et al., 2018)Radiosonde2015–2018
    NCAR/RDAds099.0SYNOP, radiosonde, pilot, dropsonde, aircraft, buoy, and ship1979–2014
    Wind profiler over US1987–2014
    ds735.0SYNOP, radiosonde, pilot, dropsonde, aircraft, buoy, and ship2015–2018
    ds351.0Aircraft1999–2018
    NCEIISDSYNOP and METAR1979–2018
    IGRA V2.0Radiosonde1979–2018
    ICOADS 3.0Buoy and ship1979–2018
    DSI-6380Aircraft1979–1998
    CEDAAMDAR reports collected by the Met OfficeAircraft2009–2018
    JODCNEAR-GOOS Regional Delayed Mode Data BaseBuoy and ship1979–2017
    Download: Download as CSV

    Table 3.  Reprocessed satellite datasets assimilated in CRA-40

    Observation typeSatellite productionPeriod (yyyymm)Source
    AMVGOES-8/9/10/11/12/13/14/15 -reprocessed199501–201307CIMSS
    METEOSAT-2/3/4/5/6/7/8/9 -reprocessed198205–200101,
    200403–201212
    EUMETSAT
    GMS-1/3/4/5/GOES-9/MTSAT-1R -reprocessed197901–197911,
    198703–200909
    JMA MSC
    FY-2C/D/E/G -reprocessed200506–201706CMA NSMC
    NOAA-7/9/10/11/12/14/15/16/17/18 polar wind -reprocessed198201–201412CIMSS
    MetOp-A polar wind -reprocessed200703–201412EUMETSAT
    Ocean surface windERS-1/2 AMI 25km -reprocessed199203–200101KNMI
    QuikSCAT Sea Winds 25km -reprocessed199907–200911KNMI
    MetOp-A ASCAT 25km -reprocessed200701–201403KNMI
    Oceansat-2 OSCAT 25km -reprocessed200912–201402KNMI
    GPSROCHAMP -reprocessed200105–200810CDAAC/UCAR
    SAC-C -postprocessed200701–201108CDAAC/UCAR
    COSMIC-1/2/3/4/5/6 -reprocessed200701–201404CDAAC/UCAR
    COSMIC-1/2/3/4/5/6 -postprocessed201405–201512CDAAC/UCAR
    GRACE-A/B -postprocessed200702–201609CDAAC/UCAR
    MetOp-A -reprocessed200710–201512CDAAC/UCAR
    MetOp-B -reprocessed201302–201512
    TerraSAR-X occultation -postprocessed200802–201608CDAAC/UCAR
    Download: Download as CSV
  • [1]

    Bengtsson, L., M. Kanamitsu, P. Kallberg, et al., 1982: FGGE 4-dimensional data assimilation at ECMWF. Bull. Amer. Meteor. Soc., 63, 29–43. doi: 10.1175/1520-0477-63.1.29.
    [2]

    Cannon, D. J., D. J. Brayshaw, J. Methven, et al., 2015: Using reanalysis data to quantify extreme wind power generation statistics: A 33 year case study in Great Britain. Renewable Energy, 75, 767–778. doi: 10.1016/j.renene.2014.10.024.
    [3]

    Chen, Z., Z. J. Zhou, Z. Q. Liu, et al., 2021: Bias adjustment and analysis of Chinese daily historical radiosonde temperature data. J. Meteor. Res., 35, 17–31. doi: 10.1007/s13351-021-9162-x.
    [4]

    Dee, D. P., S. M. Uppala, A. J. Simmons, et al., 2011: The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597. doi: 10.1002/qj.828.
    [5]

    Durre, I., X. G. Yin, R. S. Vose, et al., 2018: Enhancing the data coverage in the Integrated Global Radiosonde Archive. J. Atmos. Oceanic Technol., 35, 1753–1770. doi: 10.1175/JTECH-D-17-0223.1.
    [6]

    Freeman, E., S. D. Woodruff, S. J. Worley, et al., 2017: ICOADS Release 3.0: a major update to the historical marine climate record. Int. J. Climatol., 37, 2211–2232. doi: 10.1002/joc.4775.
    [7]

    Geer, A. J., F. Baordo, N. Bormann, et al., 2017: The growing impact of satellite observations sensitive to humidity, cloud and precipitation. Quart. J. Roy. Meteor. Soc., 143, 3189–3206. doi: 10.1002/qj.3172.
    [8]

    Gelaro, R., W. McCarty, M. J. Suárez, et al., 2017: The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). J. Climate, 30, 5419–5454. doi: 10.1175/JCLI-D-16-0758.1.
    [9]

    Gibson, J. K., P. Kallberg, S. M. Uppala, et al., 1999: ECMWF Reanalysis Project Report 1, ERA-15 Description (version 2). ECMWF Technical Report, Reading, UK, 1–74.
    [10]

    Haimberger, L., 2007: Homogenization of radiosonde temperature time series using innovation statistics. J. Climate, 20, 1377–1403. doi: 10.1175/JCLI4050.1.
    [11]

    Hersbach, H., B. Bell, P. Berrisford, et al., 2020: The ERA5 global reanalysis. Quart. J. Roy. Meteor. Soc., 146, 1999–2049. doi: 10.1002/qj.3803.
    [12]

    Hu, M., G. Q. Ge, H. Shao, et al., 2017: Grid-Point Statistical Interpolation (GSI) User’s Guide Version 3.6. Developmental Testbed Center, 1–157. Available online at http://www.dtcenter.org/sites/default/files/GSIUserGuide_v3.6.pdf. Accessed on 11 February 2023.
    [13]

    Jiang, L. P., C. X. Shi, T. Zhang, et al., 2020: Evaluation of assimilating FY-3C MWHS-2 radiances using the GSI global analysis system. Remote Sens., 12, 2511. doi: 10.3390/rs12162511.
    [14]

    Kalnay, E., M. Kanamitsu, R. Kistler, et al., 1996: The NCEP/NCAR 40-year reanalysis project. Bull. Amer. Meteor. Soc., 77, 437–472. doi: 10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.
    [15]

    Kanamitsu, M., W. Ebisuzaki, J. Woollen, et al., 2002: NCEP–DOE AMIP-II reanalysis (R-2). Bull. Amer. Meteor. Soc., 83, 1631–1644. doi: 10.1175/BAMS-83-11-1631.
    [16]

    Kleist, D. T., D. F. Parrish, J. C. Derber, et al., 2009: Introduction of the GSI into the NCEP Global Data Assimilation System. Wea. Forecasting, 24, 1691–1705. doi: 10.1175/2009WAF2222201.1.
    [17]

    Kobayashi, S., Y. Ota, Y. Harada, et al., 2015: The JRA-55 reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 5–48. doi: 10.2151/jmsj.2015-001.
    [18]

    Li, C. X., T. B. Zhao, C. X. Shi, et al., 2020: Evaluation of daily precipitation product in China from the CMA Global Atmospheric Interim Reanalysis. J. Meteor. Res., 34, 117–136. doi: 10.1007/s13351-020-8196-9.
    [19]

    Li, C. Y., and Z. X. Long, 1992: QBO and its influence on the general atmospheric circulation and the climate in East Asia. Chinese J. Atmos. Sci., 16, 167–176. doi: 10.3878/j.issn.1006-9895.1992.02.05. (in Chinese)
    [20]

    Liang, X., L. P. Jiang, Y. Pan, et al., 2020: A 10-yr global land surface reanalysis interim dataset (CRA-Interim/Land): Implementation and preliminary evaluation. J. Meteor. Res., 34, 101–116. doi: 10.1007/s13351-020-9083-0.
    [21]

    Liao, J., H. Y. Wang, Z. J. Zhou, et al., 2021: Integration, quality assurance, and usage of global aircraft observations in CRA. J. Meteor. Res., 35, 1–16. doi: 10.1007/s13351-021-0093-3.
    [22]

    McCormack, J. P., S. D. Eckermann, D. E. Siskind, et al., 2006: CHEM2D-OPP: A new linearized gas-phase ozone photochemistry parameterization for high-altitude NWP and climate models. Atmos. Chem. Phys., 6, 4943–4972. doi: 10.5194/acp-6-4943-2006.
    [23]

    NCDC, 2002: National Climatic Data Center Data Document for Data Set 6380 (DSI-6380): Worldwide Aircraft Reports. NCDC, Asheville, USA, 14 pp.
    [24]

    Onogi, K., J. Tsutsui, H. Koide, et al., 2007: The JRA-25 reanalysis. J. Meteor. Soc. Japan, 85, 369–432. doi: 10.2151/jmsj.85.369.
    [25]

    Park, C.-K., and S. D. Schubert, 1997: On the nature of the 1994 East Asian summer drought. J. Climate, 10, 1056–1070. doi: 10.1175/1520-0442(1997)010<1056:OTNOTE>2.0.CO;2.
    [26]

    Pawson, S., and M. Fiorino, 1999: A comparison of reanalyses in the tropical stratosphere. Part 3: inclusion of the pre-satellite data era. Climate Dyn., 15, 241–250. doi: 10.1007/s003820050279.
    [27]

    Rienecker, M. M., M. J. Suarez, R. Gelaro, et al., 2011: MERRA: NASA’s Modern-Era Retrospective Analysis for Research and Applications. J. Climate, 24, 3624–3648. doi: 10.1175/JCLI-D-11-00015.1.
    [28]

    Saha, S., S. Moorthi, H.-L. Pan, et al., 2010: The NCEP Climate Forecast System Reanalysis. Bull. Amer. Meteor. Soc., 91, 1015–1058. doi: 10.1175/2010BAMS3001.1.
    [29]

    Schubert, S. D., R. B. Rood, and J. Pfaendtner, 1993: An assimilated dataset for Earth science applications. Bull. Amer. Meteor. Soc., 74, 2331–2342. doi: 10.1175/1520-0477(1993)074<2331:AADFES>2.0.CO;2.
    [30]

    Sela, J., 2010: The Derivation of Sigma Pressure Hybrid Coordinate Semi-Lagrangian Model Equations for the GFS. NCEP Office Note 462, Maryland, USA, 1–31.
    [31]

    Sela, J. G., 2009: The Implementation of the Sigma Pressure Hybrid Coordinate into the GFS. NCEP Office Note 461, Maryland, USA, 1–25.
    [32]

    Serreze, M. C., A. P. Barrett, and J. Stroeve, 2012: Recent changes in tropospheric water vapor over the Arctic as assessed from radiosondes and atmospheric reanalyses. J. Geophys. Res. Atmos., 117, D10104. doi: 10.1029/2011JD017421.
    [33]

    Shao, H., J. Derber, X.-Y. Huang, et al., 2016: Bridging research to operations transitions: Status and plans of community GSI. Bull. Amer. Meteor. Soc., 97, 1427–1440. doi: 10.1175/BAMS-D-13-00245.1.
    [34]

    Smith, A., N. Lott, and R. Vose, 2011: The Integrated Surface Database: Recent developments and partnerships. Bull. Amer. Meteor. Soc., 92, 704–708. doi: 10.1175/2011BAMS3015.1.
    [35]

    Toreti, A., A. Maiorano, G. De Sanctis, et al., 2019: Using reanalysis in crop monitoring and forecasting systems. Agric. Syst., 168, 144–153. doi: 10.1016/j.agsy.2018.07.001.
    [36]

    Trenberth, K. E., J. T. Fasullo, and J. Mackaro, 2011: Atmospheric moisture transports from ocean to land and global energy flows in reanalyses. J. Climate, 24, 4907–4924. doi: 10.1175/2011JCLI4171.1.
    [37]

    Uppala, S. M., P. W. Kållberg, A. J. Simmons, et al., 2005: The ERA-40 re-analysis. Quart. J. Roy. Meteor. Soc., 131, 2961–3012. doi: 10.1256/qj.04.176.
    [38]

    Wallace, J. M., 1973: General circulation of the tropical lower stratosphere. Rev. Geophys., 11, 191–222. doi: 10.1029/RG011i002p00191.
    [39]

    Wang, J. C., H. J. Lu, W. Han, et al., 2017: Improvements and performances of the operational GRAPES_GFS 3DVar system. J. Appl. Meteor. Sci., 28, 11–24. doi: 10.11898/1001-7313.20170102. (in Chinese)
    [40]

    Wu, W.-S., R. J. Purser, and D. F. Parrish, 2002: Three-dimensional variational analysis with spatially inhomogeneous covariances. Mon. Wea. Rev., 130, 2905–2916. doi: 10.1175/1520-0493(2002)130<2905:TDVAWS>2.0.CO;2.
    [41]

    Xue, Y. H., J. Li, Z. L. Li, et al., 2020: Assessment of upper tropospheric water vapor monthly variation in reanalyses with near-global homogenized 6.5-μm radiances from geostationary satellites. J. Geophys. Res. Atmos., 125, e2020JD032695. doi: 10.1029/2020JD032695.
    [42]

    Yang, J. X., M. T. Huang, and P. M. Zhai, 2021: Performance of the CRA-40/Land, CMFD, and ERA-Interim datasets in reflecting changes in surface air temperature over the Tibetan Plateau. J. Meteor. Res., 35, 663–672. doi: 10.1007/s13351-021-0196-x.
    [43]

    Yao, B., C. Liu, Y. Yin, et al., 2020a: Evaluation of cloud properties from reanalyses over East Asia with a radiance-based approach. Atmos. Meas. Tech., 13, 1033–1049. doi: 10.5194/amt-13-1033-2020.
    [44]

    Yao, B., S. W. Teng, R. Z. Lai, et al., 2020b: Can atmospheric reanalyses (CRA and ERA5) represent cloud spatiotemporal characteristics? Atmos. Res., 244, 105091. doi: 10.1016/j.atmosres.2020.105091.
    [45]

    Zeng, Z. L., X. Wang, Z. M. Wang, et al., 2022: A 35-year daily global solar radiation dataset reconstruction at the Great Wall Station, Antarctica: First results and comparison with ERA5, CRA40 reanalysis, and ICDR (AVHRR) satellite products. Front. Earth Sci., 10, 961799. doi: 10.3389/feart.2022.961799.
    [46]

    Zhang, L., A. Kumar, and W. Q. Wang, 2012: Influence of changes in observations on precipitation: A case study for the Climate Forecast System Reanalysis (CFSR). J. Geophys. Res. Atmos., 117, D08105. doi: 10.1029/2011JD017347.
    [47]

    Zhang, L., Y. Z. Liu, Y. Liu, et al., 2019: The operational global four-dimensional variational data assimilation system at the China Meteorological Administration. Quart. J. Roy. Meteor. Soc., 145, 1882–1896. doi: 10.1002/qj.3533.
    [48]

    Zhang, S.-Q., G.-Y. Ren, Y.-Y. Ren, et al., 2021: Comprehensive evaluation of surface air temperature reanalysis over China against urbanization-bias-adjusted observations. Adv. Climate Change Res., 12, 783–794. doi: 10.1016/j.accre.2021.09.010.
    [49]

    Zhao, B., B. Zhang, C. X. Shi, et al., 2019: Comparison of the global energy cycle between Chinese Reanalysis Interim and ECMWF Reanalysis. J. Meteor. Res., 33, 563–575. doi: 10.1007/s13351-019-8129-7.
    [50]

    Zhao, P., and L. X. Chen, 2001: Climatic features of atmospheric heat source/sink over the Qinghai-Xizang Plateau in 35 years and its relation to rainfall in China. Sci. China Ser. D Earth Sci., 44, 858–864. doi: 10.1007/BF02907098.
    [51]

    Zhao, P., R. H. Zhang, J. P. Liu, et al., 2007: Onset of southwesterly wind over eastern China and associated atmospheric circulation and rainfall. Climate Dyn., 28, 797–811. doi: 10.1007/s00382-006-0212-y.
    [52]

    Zhao, P., X. D. Xu, F. Chen, et al., 2018: The third atmospheric scientific experiment for understanding the Earth–atmosphere coupled system over the Tibetan Plateau and its effects. Bull. Amer. Meteor. Soc., 99, 757–776. doi: 10.1175/BAMS-D-16-0050.1.
    [53]

    Zhao, Q. Y., and F. H. Carr, 1997: A prognostic cloud scheme for operational NWP models. Mon. Wea. Rev., 125, 1931–1953. doi: 10.1175/1520-0493(1997)125<1931:APCSFO>2.0.CO;2.
    [54]

    Zhou, T.-J., and R.-C. Yu, 2005: Atmospheric water vapor transport associated with typical anomalous summer rainfall patterns in China. J. Geophys. Res. Atmos., 110, D08104. doi: 10.1029/2004JD005413.
    [55]

    Zhou, X. J., P. Zhao, J. M. Chen, et al., 2009: Impacts of thermodynamic processes over the Tibetan Plateau on the Northern Hemispheric climate. Sci. China Ser. D Earth Sci., 52, 1679–1693. doi: 10.1007/s11430-009-0194-9.
    [56]

    Zhu, Y. Q., J. Derber, A. Collard, et al., 2014: Enhanced radiance bias correction in the National Centers for Environmental Prediction’s Gridpoint Statistical Interpolation data assimilation system. Quart. J. Roy. Meteor. Soc., 140, 1479–1492. doi: 10.1002/qj.2233.
    [57]

    Zhu, Y. Q., J. C. Derber, R. J. Purser, et al., 2015: Variational correction of aircraft temperature bias in the NCEP’s GSI analysis system. Mon. Wea. Rev., 143, 3774–3803. doi: 10.1175/MWR-D-14-00235.1.
    [58]

    Zhu, Y. Q., E. Liu, R. Mahajan, et al., 2016: All-sky microwave radiance assimilation in NCEP’s GSI analysis system. Mon. Wea. Rev., 144, 4709–4735. doi: 10.1175/MWR-D-15-0445.1.
  • Chunxiang SHI and Zhiquan LIU.pdf

  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

CRA-40/Atmosphere—The First-Generation Chinese Atmospheric Reanalysis (1979–2018): System Description and Performance Evaluation

    Corresponding author: Lipeng JIANG, jianglp@cma.gov.cn
    Corresponding author: Chunxiang SHI, shicx@cma.gov.cn
  • 1. National Meteorological Information Centre, China Meteorological Administration (CMA), Beijing 100081, China
  • 2. National Center for Atmospheric Research, Boulder, CO 80301, USA
  • 3. CMA Earth System Modeling and Prediction Centre, China Meteorological Administration, Beijing 100081, China
Funds: Supported by the China Meteorological Administration Special Public Welfare Research Fund (GYHY201506002) and National Innovation Project for Meteorological Science and Technology (CMAGGTD003-5)

Abstract: Atmospheric reanalysis reproduces the past atmospheric conditions through assimilation of historical meteorologi-cal observations with fixed version of a numerical weather prediction (NWP) model and data assimilation (DA) system. It is widely used in weather, climate, and even business-related research and applications. This paper reports the development of CMA’s first-generation global atmospheric reanalysis (RA) covering 1979–2018 (CRA-40; CRA refers to CMA-RA). CRA-40 is produced by using the Global Spectral Model (GSM)/Gridpoint Statistical Interpolation (GSI) at a 6-h time interval and a TL574 spectral (34-km) resolution with the model top at 0.27 hPa. A large number of reprocessed satellite data and widely collected conventional observations were assimilated during the reanalyzing process, including the reprocessed atmospheric motion vector (AMV) products from FY-2C/D/E/G satellites, dense conventional observations (at about 120 radiosonde and 2400 synoptic stations) over China, as well as MWHS-2 and GNSS-RO observations from FY-3C. The reanalysis fitting to observations is improved over time, especially for surface pressure with root-mean-square error reduced from 1.05 hPa in 1979 to 0.8 hPa, and for upper air temperature from 1.65 K in 1979 to 1.35 K, in 2018. The patterns of global analysis increments for temperature, specific humidity, and zonal wind are consistent with the changes in the observing system. Near surface temperature from the model’s 6-h forecast reflects the global warming trend reasonably. The CRA-40 precipitation pattern matches well with those of GPCP and other reanalyses. CRA-40 also successfully captures the QBO and its vertical and temporal development, hemispherical atmospheric circulation change, and moisture transport by the East Asian summer monsoon. CRA is now operationally running in near real time as a climate data assimilation system in CMA.

中国气象局第一代全球大气再分析CRA-40(1979–2018):系统介绍和性能评估

大气再分析是利用固定的、最先进的数值预报模式和数据同化系统,以及各种来源的观测资料,对大气状况进行历史回算,得到尽可能稳定、均一和精确的长序列历史分析数据产品。本文主要介绍中国气象局第一代40年(1979年–准实时)全球大气再分析(CRA-40;注意本文CRA等同于CMA-RA)系统的系统配置和产品评估结果。CRA-40产品时间分辨率逐6小时,水平分辨率约34 km,模式层顶0.27 hPa。与国际上其他全球再分析产品相比,CRA-40在同化过程采用了更多的中国特色观测资料,比如中国风云2号系列卫星大气运动矢量重处理产品,中国特有常规观测资料,FY-3C微波湿度计和掩星观测资料,等。CRA-40降水产品的空间分布形态与GPCP降水分析产品及国际其他再分析降水产品比较一致;CRA-40近地面气温场能够成功反映出全球变暖趋势。另外,CRA-40产品能成功再现赤道平流层低层风场准两年振荡(QBO)、全球大气环流特征及东亚夏季风水汽输送特征等。目前,CRA-40系统已作为中国气象局业务系统准实时运行。
    • Atmospheric reanalysis reproduces the past atmos-pheric conditions through the assimilation of historical meteorological observations and the fixed version of a numerical weather prediction (NWP) model and a data assimilation (DA) system. Since the mid-1990s, several agencies (mainly NOAA, NASA, ECMWF, and JMA) have successively conducted four rounds of global atmos-pheric reanalysis as listed in Table 1. Note that ECMWF named their latest reanalysis as ERA5 (i.e., the 5th gene-ration) with consideration of the FGGE reanalysis (Beng-tsson et al., 1982). These reanalysis products have been widely used in weather, climate, and related research and applications (e.g., Trenberth et al., 2011; Serreze et al., 2012). They are also used increasingly for business applications such as in energy (e.g., Cannon et al., 2015) and agriculture (e.g., Toreti et al., 2019) sectors. Most of these reanalysis systems are continuously running as operational “climate data assimilation system (CDAS)” for climate monitoring. More recent reanalysis products usually have been produced with more advanced NWP models and DA systems at higher resolutions and with more optimal use of observations.

      GenerationReanalysisProducerPeriodResolutionMethodReference
      1stNASANASA1980–19952° × 2.5°OISchubert et al. (1993)
      R1NCEP + NCAR1948– T623DVarKalnay et al. (1996)
      ERA-15ECMWF1979–1994T106OIGibson et al. (1999)
      2ndR2NCEP + DOE1979– T623DVarKanamitsu et al. (2002)
      ERA-40ECMWF1957–2002TL1593DVarUppala et al. (2005)
      JRA-25JMA-CRIEPI1979– T1063DVarOnogi et al. (2007)
      3rdERA-InterimECMWF1979–2019TL2554DVarDee et al. (2011)
      CFSRNCEP1979– T3823DVarSaha et al. (2010)
      MERRANASA1979–20161/2° × 2/3°3DVarRienecker et al. (2011)
      MERRA-2NASA1980– 1/2° × 5/8°3DVarGelaro et al. (2017)
      JRA-55JMA1957– TL3194DVarKobayashi et al. (2015)
      4thERA5ECMWF1950– TL639Hybrid 4DVarHersbach et al. (2020)

      Table 1.  Global atmospheric reanalyses

      In November 2013, China Meteorological Administration (CMA) launched the global reanalysis project, which is led by the National Meteorological Information Centre (NMIC) with participation of several other national centers of CMA and also external research institutes and universities. The goal was to produce China’s first-generation global atmospheric and land reanalysis (RA) products for a 40-yr (1979–2018) period (so this reanalysis is named CRA-40; here CRA actually refers to CMA-RA) and then the reanalysis system will become operational as a CDAS to support climate monitoring activities at the National Climate Centre (NCC) of CMA.

      At the time of proposing CRA-40, the CMA-developed GRAPES-GFS was not mature and the CMA’s operational global NWP system was based on an old version of ECMWF’s spectral model configured at TL639 spectral resolution (~31 km) and NCEP’s Gridpoint Statisti-cal Interpolation (GSI) DA system. The original plan was to use operational T639 + GSI at that time for CRA-40 atmospheric reanalysis. However, during the early stage of project execution it was found that the T639 model exhibited significant temperature bias near the model top, which greatly limited the optimal use of satellite radiance data. GRAPES-GFS began operational on 1 June 2016 with a 3DVar DA scheme (Wang et al., 2017) and was then upgraded to 4DVar on 1 July 2018 (Zhang et al., 2019). However, GRAPES-GFS’s DA system cannot assimilate satellite radiance data from old satellite sensors and relatively low model top (36 km) of GRAPES-GFS also restricts the optimal assimilation of satellite radiances. Therefore, it was decided to replace the T639 model with NCEP’s publically available Global Spectral Model (GSM), while keeping GSI as atmospheric DA system for CRA-40.

      An interim version of GFS/GSI-based CRA (CRA-Interim) for a 10-yr period (2007–2016) at a 6-h time interval and a TL574 spectral resolution (34 km) was produced in February 2018. A 10-yr (2007–2016) CRA-Interim/Land dataset with the same time and spatial (34-km) resolutions as CRA-Interim was also produced by forcing Noah land surface model using observation blended CRA-Interim precipitation and near surface temperature, pressure, humidity, wind, and radiation directly from CRA-Interim (Liang et al., 2020). Evaluation of CRA-Interim shows some advantages in terms of clouds (Yao et al., 2020a), precipitation (Li et al., 2020), upper troposphere water vapor (Xue et al., 2020), and the ability of depicting global energy cycle (Zhao et al., 2019). In late 2019, CRA-40/Atmosphere and CRA-40/Land were completed. This article documents CRA-40/Atmosphere. Section 2 describes the model and DA system used in CRA-40 as well as the production streams. Section 3 provides some details about the observations assimilated in CRA-40. Section 4 presents the DA diagnostics of CRA-40, followed by the evaluation of some aspects of the products in Section 5. The paper is concluded in Section 6.

    2.   Model, data assimilation, and production streams
    • As mentioned earlier, CRA-40 adopted NCEP’s GSM and the community GSI DA system for atmospheric reana-lysis. Configurations of GSM model and GSI DA system will be briefly described in the following subsections.

    • The atmospheric prediction model used by CRA-40 is GSM version 14, operational at NCEP from July 2017 to June 2019. The GSM-v14 was made publicly available by the NCEP’s Central Operation (NCO). This is the last operational implementation of GSM before the transition to FV3-based GFS at NCEP in June 2019. GSM-v14 was configured with TL1534 spectral resolution (~13-km grid spacing at the equator) at the NCEP’s operation, but CRA-40 configured it at a reduced spectral resolution TL574 (~34-km grid spacing at the equator) to lower the computational cost. In the vertical, there are 64 hybrid sigma pressure layers (Sela, 2009) with the model top around 0.27 hPa (~55 km). The GSM dynamic core adopts a two time-level semi-implicit and semi-Lagrangian time discretization (Sela, 2010). The semi-Lagrangian advection calculation and physics processes are performed on a linear, reduced Gaussian grid in the horizontal domain.

      GSM’s prognostic variables consist of vorticity, divergence, logarithm of surface pressure, specific humidity, virtual temperature, mixing ratios of cloud liquid water and cloud ice (Zhao and Carr, 1997), and O3 (McCormack et al., 2006) (see https://dtcenter.ucar.edu/GMTB/v4.1.0/sci_doc/GFS_OZPHYS.html). GSM-v14 replaced spectral file output (sigma files) with new NEMSIO (https://www.emc.ncep.noaa.gov/emc/pages/infrastructure/nems.php) binary files on model native grid. More details about physical parameterization schemes of GSM and specific changes in GSM-v14 can be found at https://www.emc.ncep.noaa.gov/emc/pages/numerical_forecast_systems/gfs/documentation.php. SST and sea ice are important lower boundary conditions for any atmospheric reanalyses. CRA-40 uses daily SST and sea ice data (0.25° × 0.25°) from CFSR for the period from 1979 to 2014. As CFSR did not continue after December 2014, we use NCEP operational daily Real-Time Global SST (RTG SST, 8 km × 8 km) and sea ice analysis (0.5° × 0.5°) for more recent years from 2015 to 2018.

    • The community version 3.6 of GSI (Wu et al., 2002; Kleist et al., 2009; Shao et al., 2016) was adopted for the final production of CRA-40/Atmosphere (http://www.dtcenter.org/community-code/gridpoint-statistical-interpolation-gsi/community-gsi-version-3-6-enkf-version-1-2), which was released in September 2017 by the Developmental Testbed Center (DTC). While more advanced DA techniques such as hybrid-3DEnVar or hybrid-4DEnVar are available in GSI, CRA-40 used 3DVar for the sake of its simplicity and computational efficiency. This choice is also constrained by the requirement that the production needs to be completed within a relatively short time period (several months) for all 40-yr products.

      The background and analysis files of the GSI-3DVar contain the grid point fields (instead of spectral coefficients) at a 34-km horizontal resolution and in the NEMSIO format. The analysis was carried out at each 0000, 0600, 1200, and 1800 UTC with assimilation of conventional observations, satellite retrievals, and satellite radiance data (see more details in Section 3). GSI-3DVar was run in the FGAT mode, in which assimilation time window is 3-h centered at analysis times and the three background files at −3-, 0-, and 3-h come from a 9-h GSM forecast initialized from the previous cycle’s analysis. Variational bias correction (VarBC) was applied to clear-sky satellite radiances (Zhu et al., 2014) and aircraft temperature data (Zhu et al., 2015). However, cloudy AMSU-A radiance assimilation capability (Zhu et al., 2016) available in GSI version 3.6 was not utilized. For radiance DA, version 2.2.3 of the Community Radiative Transfer Model (CRTM) was used. Channel selection of satellite radiances for various sensors over the reanalysis period is based upon CFSR, which is also a GSI-based reanalysis for its atmospheric component.

      GSI observation interface was modified to allow the assimilation of atmospheric motion vectors (AMVs) from Chinese FY-2 geostationary satellites, and GNSS radio occultation bending angle data and MWHS-2 microwave radiances (Jiang et al., 2020) from Chinese FY-3C polar-orbiting satellite. Another important development within GSI is to interface the model level ERA-Interim reanalysis data so that GSI can serve as a pre-evaluation tool for observations by comparing ERA-Interim with observations over the 40-yr period. The addition of this functionality in GSI and the observation pre-evaluation step prior to the reanalysis production greatly helped to implement CRA-40’s own quality control and blacklisting strategy for conventional observations, which is approved to be crucial for the success of CRA-40.

      The analysis variables of GSI consist of stream function, unbalanced part of velocity potential, temperature and surface pressure, and normalized relative humidity as well as O3 with the background error covariance model following Wu et al. (2002). The GSI’s static background error covariance statistics file of CRA-40 is taken from the NCEP’s operational implementation of GFS-v14 and is not changed over the 40-yr reanalysis period. Note that the inner loop minimization of GSI hybrid-4DEnVar incremental analysis at the NCEP’s GFS-v14 operation was done at T574 resolution (so the resolution of static background error covariances) even though the GSM-v14’s forecast background is at T1534 resolution. Observation errors were configured with relevant files in the community GSI package (Hu et al., 2017). It is worth mentioning that temperature and humidity at 2 m and zonal and meridional winds at 10 m are model diagnostic variables and are not analyzed by GSI. The analysis of those screen-level variables will come from CRA-40/Land.

    • CRA-40 covers a 40-yr period from January 1979 to December 2018. The CRA-40 production is divided into 10 parallel streams with a 4.5-yr period for each stream. This enables the whole 40-yr production being completed within one year as required by the project. There is a 6-month overlap between two adjacent streams and the first six months’ data of each stream are discarded so that different streams’ data can be joined together without discontinuity. For instance, stream-10 is run from 1 July 2014 to 31 December 2018. However, only the data between 1 January 2015 and 31 December 2018 in stream-10 are kept for the final delivered products and the data between 1 July 2014 and 31 December 2014 in the final products are actually from stream-9, which is run from 1 July 2010 to 31 December 2014. This multiple stream production strategy is common practice for most reanalysis projects even though the length of each stream and overlap period can be different among different reanalyses. Due to computational and storage cost, no longer range forecasts are carried out. Model level analysis is postprocessed into pressure level products by using UPP (https://dtcenter.org/community-code/unified-post-processor-upp). A list of CRA-40/Atmosphere products is given in Appendix A.

    3.   Observations
    • Conventional observations assimilated in CRA-40 are based on NCEP’s Prepbufr data used in CFSR (1979–2014, https://rda.ucar.edu/datasets/ds099.0) (Saha et al., 2010) and GDAS (2015–2018, https://rda.ucar.edu/datasets/ds735.0), but supplemented with additional data collected by NMIC.

      For synoptic observations over land (SYNOP), CRA-40 combined data from CFSR/GDAS, Integrated Surface Database (ISD; Smith et al., 2011), and more than 2400 CMA’s surface synoptic stations. No METAR data are used in CFSR before June 1997. All the hourly METAR data used in CRA-40 come from ISD. METAR data from 1997 to 2009 in CFSR contain more observations than ISD, and therefore the counts of surface pressure (Ps) used in CRA-40 are slightly less than CFSR. Surface observations over the ocean (buoy and ships) combined data from CFSR/GDAS, ICOADS R3.0 (Freeman et al., 2017), NMIC archived global and regional datasets, and the North–East Asian ocean observations from Japan Oceanographic Data Center (JODC; https://near-goos1.jodc.go.jp). While multiple observed parameters are contained in these datasets, only surface pressure (Ps) observations are assimilated over land. Over the ocean, surface pressure, temperature, humidity, and wind observations are all assimilated in GSI.

      Radiosonde observations combined six data sources: (1) NCEP’s CFSR/GDAS; (2) NCEI’s IGRA V2.0 dataset (Durre et al., 2018); (3) about 120 radiosonde stations operated by CMA; (4) about 80 pilot balloon stations operated by CMA from 1980 to 1991; (5) NMIC archived global radiosonde data from GTS; and (6) some field campaign datasets collected and archived by NMIC. Radiosonde temperature observations are bias corrected following RAOBCORE 1.4 (Haimberger, 2007) for non-Chinese stations, and Chen et al. (2021) for Chinese stations.

      Aircraft data (AMDAR, ACARS, and AIREP) are based upon the integrated global aircraft dataset developed by NMIC (1973–present, Liao et al., 2021), which combined data from CFSR/GDAS, NCAR archived reports (from October 1999, https://rda.ucar.edu/datasets/ds351.0), NCEI’s worldwide aircraft reports (1973–1998, DSI-6380; NCDC, 2002), AMDAR reports from the Centre for Environmental Data Analysis (CEDA) archive of Met Office (https://catalogue.ceda.ac.uk/uuid/33f44351f9ceb09c495b8cef74860726), and NMIC archived Chinese AMDAR data (from May 2003). Variational bias correction is applied to aircraft temperature observations in GSI (Zhu et al., 2015).

      Global wind profiler and pilot/pibal data are mainly from CFSR/GDAS and mostly available over US. Chinese wind profiler data are supplemented from 2007 with more than 100 stations in 2018.

      Duplicated data from multiple data sources are removed according to the priority order predetermined for different data sources. Each type of conventional data is converted into the NCEP’s “BUFR DUMP” file format, which is then preprocessed into NCEP’s “Prepbufr” file format along with the GFS background file to perform in-fly quality control and is subsequently assimilated in GSI. The monthly blacklists of radiosonde, aircraft, and surface data are pre-generated by evaluating the data quality for each station, airplane, ship, and buoy at different times. The differences between observations and the ERA-Interim reanalysis are considered as “quality feedback information” to estimate the quality of conventional observations. The blacklists are determined depending on the quality feedback information and carefully manual diagnostic.

      Availability of different types of conventional observations used in CRA-40 is given in Table 2. Figure 1 shows the monthly counts for all conventional observations (Ps, t, q, u) used in CRA-40 and CFSR. The two reanalyses differ mainly in the additional observations over China, the METAR data usage, and the assimilation of aircraft humidity.

      Data supplierData type and suppliers’ identifierPeriod
      CMAChinese operational observationsSYNOP, radiosonde, and pilot1979–2018
      Buoy and ship2013–2018
      Wind profiler2007–2018
      Aircraft2003–2018
      Operational GTS dataSYNOP, radiosonde, aircraft, buoy, and ship1979–2018
      Intense observations of TIPEX III (Zhao et al., 2018)Radiosonde2015–2018
      NCAR/RDAds099.0SYNOP, radiosonde, pilot, dropsonde, aircraft, buoy, and ship1979–2014
      Wind profiler over US1987–2014
      ds735.0SYNOP, radiosonde, pilot, dropsonde, aircraft, buoy, and ship2015–2018
      ds351.0Aircraft1999–2018
      NCEIISDSYNOP and METAR1979–2018
      IGRA V2.0Radiosonde1979–2018
      ICOADS 3.0Buoy and ship1979–2018
      DSI-6380Aircraft1979–1998
      CEDAAMDAR reports collected by the Met OfficeAircraft2009–2018
      JODCNEAR-GOOS Regional Delayed Mode Data BaseBuoy and ship1979–2017

      Table 2.  Data sources of conventional observations for CRA-40

      Figure 1.  Monthly counts of conventional observations for (a) surface pressure, (b) temperature, (c) specific humidity, and (d) wind used in CRA-40 (red line) and CFSR (black line) from 1979 to 2014.

    • Satellite observations represent a majority of dataset assimilated in “satellite-era” reanalyses (typically referred to those reanalyses from 1979), including radiances from microwave and infrared sensors onboard polar-orbiting and geostationary platforms, GPS radio occultation bending angle data, retrieved AMVs, surface ocean winds, and O3 concentration, as listed in Table 3 and Figs. 2 and 3. While satellite radiances and retrieved O3 data are from CFSR/GDAS BUFR files, CRA-40 assimilated the latestly reprocessed GPSRO data by UCAR’s COSMIC, AMVs by CIMSS, EUMETSAT, JMA, CMA, and surface ocean winds by KNMI. Unlike conventional observations, quality control procedures of satellite data are mostly performed within GSI. For reprocessed retrieval products (AMVs, surface ocean winds), quality flags provided by the processing agencies are also kept when converting raw retrieval data into GSI’s BUFR format files and used in data assimilation.

      Observation typeSatellite productionPeriod (yyyymm)Source
      AMVGOES-8/9/10/11/12/13/14/15 -reprocessed199501–201307CIMSS
      METEOSAT-2/3/4/5/6/7/8/9 -reprocessed198205–200101,
      200403–201212
      EUMETSAT
      GMS-1/3/4/5/GOES-9/MTSAT-1R -reprocessed197901–197911,
      198703–200909
      JMA MSC
      FY-2C/D/E/G -reprocessed200506–201706CMA NSMC
      NOAA-7/9/10/11/12/14/15/16/17/18 polar wind -reprocessed198201–201412CIMSS
      MetOp-A polar wind -reprocessed200703–201412EUMETSAT
      Ocean surface windERS-1/2 AMI 25km -reprocessed199203–200101KNMI
      QuikSCAT Sea Winds 25km -reprocessed199907–200911KNMI
      MetOp-A ASCAT 25km -reprocessed200701–201403KNMI
      Oceansat-2 OSCAT 25km -reprocessed200912–201402KNMI
      GPSROCHAMP -reprocessed200105–200810CDAAC/UCAR
      SAC-C -postprocessed200701–201108CDAAC/UCAR
      COSMIC-1/2/3/4/5/6 -reprocessed200701–201404CDAAC/UCAR
      COSMIC-1/2/3/4/5/6 -postprocessed201405–201512CDAAC/UCAR
      GRACE-A/B -postprocessed200702–201609CDAAC/UCAR
      MetOp-A -reprocessed200710–201512CDAAC/UCAR
      MetOp-B -reprocessed201302–201512
      TerraSAR-X occultation -postprocessed200802–201608CDAAC/UCAR

      Table 3.  Reprocessed satellite datasets assimilated in CRA-40

      Figure 2.  Timeline of conventional data, AMVs (i.e., SATWND), ocean surface winds, ozone, and GPSRO data used in CRA-40.

      Figure 3.  Timeline of satellite radiance data using in CRA-40, including those from various microwave (MSU, AMSU-A, AMSU-B, MHS, MWHS-2, ATMS, and SSMIS) and infrared sensors (HIRS2–4, AVHRR3, AIRS, IASI, and CrIS) onboard various polar-orbiting satellites (e.g., n* for NOAA series, f17 for DMSP 17 satellite). Infrared radiances from geostationary GOES (g*) sounders are indicated by sndr*, while those from Meteosat-10 (m10) are indicated by SEVIRI.

    4.   Data assimilation diagnostics
    • The 40-yr evolution of the root-mean-square (RMS) of the first-guess departures (OMB) and analysis departures (OMA) from observations is shown in Fig. 4. The statistics are averaged over globe and over all vertical levels for upper air observations. Generally, the fitting to observations is improved over time, which is more pronounced from surface pressure and upper air temperature with the RMS of OMB (OMA) reduced from ~1.35 (1.05) hPa in 1979 to ~1.0 (0.8) hPa in 2018 for surface pressure and from ~2.1 (1.65) K in 1979 to ~1.6 (1.35) K in 2018 for upper air temperature. The improving fit over time is believed to be resulted from more and more satellite radiances assimilated. This can be partly verified with a larger improving fit to surface pressure over ocean (not shown), where satellite data are the dominant dataset.

      Figure 4.  The 30-day mean (thin curves) and 360-day mean (thick curves) of the root-mean-square of first-guess (blue) (OMB) and analysis (red) (OMA) departures from observations in CRA-40 for (a) upper-air zonal wind from radiosondes, dropsondes, and PILOTs; (b) upper-air temperature from radiosondes and dropsondes; (c) upper-air specific humidity from radiosondes and dropsondes; and (d) surface pressure from SYNOP, buoys, ships, and METAR.

      The change of fit to radiosonde specific humidity observations over time is relatively small, which might be an indication of relatively small impact of satellite radiances from humidity-sensitive channels. Geer et al. (2017) demonstrated the growing impact of microwave humidity sounders with all-sky assimilation approach. However, only clear-sky radiances were assimilated in CRA-40. For upper air zonal wind, there is an apparent trend of degraded analysis fit over time from ~1990 (red line in Fig. 4a) although the first-guess fit keeps improving or remains in a similar level. While this wind analysis fit degradation trend can be an aspect to be improved in the future version of reanalysis, it is also likely an indication of assimilation effect of more and more satellite radiances, which are directly sensitive to temperature and moisture and thus reduce the analysis fit to wind observations. Note that the first-guess fit to wind observations, which is usually a better indication of data assimilation performance than the analysis fit, does not exhibit an obvious trend of degradation.

    • Pressure–time diagrams for monthly and globally averaged mean analysis increments (i.e., analysis minus background) for specific humidity, zonal wind, and temperature are displayed in Fig. 5. With bias-corrected observations and unbiased model forecast, mean increments should be small. Systematic biases of increments typically indicate an average conflict between observations and the model background. Abrupt changes in mean analysis increments are usually related to changes in the observing system.

      Figure 5.  Pressure–time diagrams of monthly and globally averaged analysis increments for (a) temperature (K), (b) zonal wind (m s−1), and (c) specific humidity (mg kg−1) on a logarithmic pressure scale.

      The pattern of temperature increments (Fig. 5a) is very similar to that of ERA5 [see Fig. 16a of Hersbach et al. (2020)], e.g., opposite increments at the top of the stratosphere indicating a model bias with respect to anchoring satellite observations which peak at those heights, positive increment band around 200 hPa, and negative increment band below 200 hPa. Abrupt change of temperature increments at the upper stratosphere in late 1998 is most likely related to the introduction of NOAA-15 AMSU-A in September 1998 (Fig. 3). Mean zonal wind increments (Fig. 5b) are relatively small with a maximum of ~0.1 m s−1 between 1 and 3 hPa. For specific humidity (Fig. 5c), sharper change from 2000 is related to the introduction of NOAA-15 AMSU-B (Fig. 3).

    5.   Performance evaluation
    • Pressure–time diagrams for monthly and globally averaged root-mean-square errors (RMSEs) of CRA-40, CFSR, and JRA-55 against ERA5 for geopotential height, zonal wind, temperature, and specific humidity are shown in Fig. 6. The RMSE of geopotential height (Fig. 6a) at high-levels (10–100 hPa) for CRA-40 is larger than that of CFSR, except for the period of 1998–2010. JRA-55 has a sudden increased RMSE for geopotential height at the layer of 10–100 hPa from 2000 to 2007. For zonal wind, all of three reanalyses have larger and periodically changing RMSEs above 100 hPa, reflecting the zonal wind high uncertainty at upper levels likely due to the lack of wind observations aloft. Similar to the time variation of analysis increment shown in Fig. 5a, there are also abrupt changes of temperature RMSE between 10 and 30 hPa in 1998, which may also be related to the evolution of observing systems. For specific humidity, CRA-40 has the smallest RMSE among three reanalyses when taking ERA5 as a reference (Fig. 6d).

      Figure 6.  Pressure–time diagrams of three reanalyses’ monthly and globally averaged root-mean-square errors (RMSEs) for (a) geopotential height (gpm), (b) zonal wind (m s−1), (c) temperature (K), and (d) specific humidity (g kg−1) against ERA5. The RMSEs of CRA-40, CFSR, and JRA-55 are shown from top to bottom in each subplot.

    • The trend analysis of near surface air temperature at 2 m (t2m) is one important application of reanalysis data. Figure 7 shows 12-month running averages of global mean t2m anomalies (K) from 1979 to 2018 for CRA-40/Atmosphere. The global warming trend is well reflected in Fig. 7. While t2m is part of CRA-40/Atmosphere products, it is worth noting that the t2m is not an analy-sis variable of GSI and it is actually from the model’s 6-h forecast background. Therefore, we suggest that users employ the t2m (also for humidity at 2 m and wind at 10 m) product from CRA-40/Land, in which those near surface variables are produced by merging the corresponding model background from CRA-40/Atmosphere with global near surface observations using the EnOI method (Liang et al., 2020). Further details on CRA-40/Land will be presented in a separate article.

      Figure 7.  Twelve-month running averages of global-mean near surface air temperature (T2m) anomalies (K) from 1979 to 2019 for CRA-40/Atmosphere. The T2m anomaly is computed relative to the mean from 1981 to 2010.

    • Figure 8 shows the climatology of global precipitation distributions in CRA-40, ERA5, JRA-55, and GPCP-v2.3 (2.5° × 2.5°) as an observational dataset. The global precipitation patterns of three reanalysis datasets match well with that of GPCP, while the precipitation magnitudes of three reanalyses are larger than that of GPCP. Over the ocean (e.g., ITCZ over Pacific), the magnitude of CRA-40 agrees more closely with that of ERA5, while JRA-55 exhibits stronger rainfall than both CRA-40 and ERA5. In some regions over land (e.g., central Africa, Southeast Asia, and north of South America), CRA-40 generally overestimates the precipitation compared to ERA5 and JRA-55, which is mainly caused by the sudden increase of precipitation over land since 1998 (not shown). Zhang et al. (2012) investigated the similar abrupt increase in precipitation from CFSR around 1998, and pointed out that it is related to the spin down of the increased initial moisture introduced by the assimilation of ATOVS data. It is also worth noting that these excessive precipitation regions in CRA-40 apparently correspond to regions with low correlation between ERA5/ERA-Interim and observed precipitation [see Fig. 24 of Hersbach et al. (2020)].

      Figure 8.  Thirty-year (1981–2010) averaged precipitation rates (mm day−1) in (a) CRA-40, (b) ERA5, (c) JRA-55, and (d) GPCP-v2.3.

      Figure 9 displays the time series of 360-day running average of global mean precipitation (mm day−1) from CRA-40, ERA5, JRA-55, and GPCP-v2.3 over 1979–2018. The magnitudes of mean daily precipitation are consistent with those shown in Fig. 8. JRA-55 has the strongest global mean daily precipitation with a mean value of 3.28 mm day−1, while ERA5 with a mean value of 2.91 mm day−1 agrees most closely with GPCP (2.69 mm day−1). CRA-40 with a value of 3.12 mm day−1 is roughly in the middle of JRA-55 and ERA5. Compared to ERA5 and JRA-55, one noticeable difference in CRA-40 is the increased precipitation amount from 2011, which is apparently associated with the wetter troposphere moisture increments in Fig. 5. The cause of excessive water vapor and precipitation from 2011 remains unclear, but we suspect that it might be related to the deficiency of satellite radiance DA for some sensors and/or channels. The mean evaporation of CRA-40 is lower than precipitation at most of years except for 2007–2009 and similar feature of evaporation vs. precipitation is also present in ERA5 (Hersbach et al., 2020). Unlike precipitation, the evaporation does not exhibit abnormal increase from 2011.

      Figure 9.  Time series of 360-day running average of global mean precipitation rates (mm day−1) from CRA-40 (red solid line), ERA5 (blue line), JRA-55 (green line), and GPCP-v2.3 (gray line). The red dash line is the corresponding evaporation from CRA-40.

    • The QBO of the zonal mean wind is a main feature of the tropical lower stratosphere, which is mainly manifested in the cyclical changes of east–west wind, with an average period of about two years (Wallace, 1973). Studies have shown that QBO is related to the Walker circulation, ENSO, tropospheric circulation, and climate change (Li and Long, 1992; Pawson and Fiorino, 1999). Figure 10 shows the zonal mean of monthly-mean zonal wind averaged between 5°S and 5°N as a function of pressure and time for ERA5, CRA-40, JRA-55, and CFSR. The results show that the QBO period and the corresponding vertical and time evolution characteristics are similar among four reanalyses. The main issue with CRA-40 is that the amplitudes of the westerly phase near 10 hPa and the easterly phase near 20 hPa are too strong when compared to other reanalyses. The issue is apparently reduced after 2005, when more AMSU-A radiance data available to constrain high-level state. We suspect that this deficiency is mainly caused by the GFS model (e.g., relatively low model top). Saha et al. (2010) pointed out that CFSR (based upon an older version of the GFS model) had a similar issue in an earlier version of the production, which was alleviated by assimilating ERA-40 stratospheric wind profiles above 20 hPa from 20°S and 20°N as bogus observations until December 1998.

      Figure 10.  Time–pressure cross-section of zonal-mean of monthly-mean zonal wind (m s−1) averaged between 5°S and 5°N for (a) ERA5, (b) CRA-40, (c) JRA-55, and (d) CFSR.

    • The hemispherical atmospheric circulation is highly related to the variability of the Tibetan Plateau heat source, and consequently influences greatly on the Asian monsoon, water and energy cycle, and precipitation variation (Zhao and Chen, 2001; Zhou et al., 2009). Based on the ECMWF reanalysis data from 1958 to 2001, Zhou et al. (2009) proposed a basic physical model to explain the influence of the Tibetan Plateau heat source on Northern Hemispherical atmospheric circulation. Following Zhou et al. (2009), we reproduced the summer zonal–vertical atmospheric circulation climatology (averaged from 1980 to 2010) across 15°–50°N, using the global reanalysis datasets from ECMWF (ERA5), NCEP (CFSR), JMA (JRA-55), and CMA (CRA-40). Figure 11a shows the longitude–pressure cross-section of zonal–vertical (15°–50°N) mean averaged summer (June, July, and August; JJA) wind of ERA5 covering 1980–2010. It can be seen that the updraft in the Tibetan Plateau separated into two parts. One part is uplifted to the lower stratosphere, moving westward and sinking over Europe, forming an anti-clock zonal–vertical circulation in the upper troposphere and the lower stratosphere. The other part flows eastward in the troposphere and sinks near the eastern Pacific Ocean. It is worth mentioning that, some of the eastward flow merges into the weak North American updraft, and continues to move eastward and then sink over the eastern Atlantic. Similar zonal–vertical atmospheric circulation characters are also found in CRA-40, CFSR, and JRA-55 (figure omitted). To see the details more clearly, Fig. 11b shows the differences between CRA-40 and ERA5, while Figs. 11c and 11d for JRA-55 and CFSR, respectively. It can be seen that, CRA-40 is closer to ERA5 than JRA-55 and CFSR over the Tibetan Plateau.

      Figure 11.  Longitude–height cross-section of summer (June, July, and August) zonal–vertical atmospheric circulation climatology (1980–2010) between 15° and 50°N (a) for ERA5 and the difference (b) between CRA-40 and ERA5, (c) between JRA-55 and ERA5, and also (d) between CFSR and ERA5.

    • Atmospheric water vapor transport is one of the most important features in the East Asian summer monsoon (EASM) system, and it has significant influences on summer rainfall patterns in China (e.g., Park and Schubert, 1997; Zhou and Yu, 2005; Zhao et al., 2007). Figure 12a shows the 40-yr (1979–2018) mean of summer (JJA) water vapor transport (in kg m−2 s−1) over East Asia at 850 hPa from ERA5. Three major water vapor pathways can be seen clearly for normal climatological conditions in East Asia, as pointed out by Zhou and Yu (2005). The first one is the strong transport from the Arabian Sea and the Bay of Bengal to eastern China; the second comes from the tropical western Pacific Ocean; and the third is a cross-equator flow straddling 105°–150°E. Figures 12b–d are similar to Fig. 12a, but show differences of CRA-40, JRA-55, and CFSR relative to ERA5, respectively. CRA-40 has smaller differences from ERA5 than JRA-55 and CFSR for the first and second transport pathways. The spatial correlations of CRA-40, JRA-55, and CFSR with ERA5 are 0.98, 0.98, and 0.97, and the root-mean-square errors are 77.1, 84.7, and 84.3 kg m−2 s−1 respectively.

      Figure 12.  Forty-year mean of JJA water vapor transport (kg m−2 s−1) at 850 hPa (a) from ERA5, and (b) the differences between CRA-40 and ERA5, (c) between JRA-55 and ERA5, and (d) between CFSR and ERA5. The colors indicate the magnitude of the water vapor flux vector.

    • Xue et al. (2020) evaluated the monthly variation of upper tropospheric water vapor (UTWV) in several reanalysis products using near-global homogenized water vapor channel (6.5 μm) radiances from geostationary satellites. They pointed out that the UTWV from ERA5, JRA-55, and CRA-40 agreed more closely to observations when compared to that from ERA-Interim, CFSR, and MERRA-2. Using MODIS data as the reference, Yao et al. (2020b) showed that CRA-40 slightly underestimated cloud cover, but captured well the influences of El Niño–Southern Oscillation (ENSO) on cloud spatiotemporal distributions like ERA5. With a newly reconstructed 35-yr daily solar radiation dataset for the Great Wall and Antarctica stations, Zeng et al. (2022) demonstrated that CRA-40 followed a relatively consistent trend with observations with a correlation coefficient of 0.977, but overestimated the solar radiation with a mean bias of 6.90 MJ m−2 during the austral summer half-year. The assessment of surface air temperature over China against urbanization-bias-adjusted observations indicated that CRA-40/Land performed the best among eight global reanalysis products for the period 1979–2015 in reproducing various aspects of climatological features in China (Zhang et al., 2021). The evaluation of surface air temperature over the Tibetan Plateau by Yang et al. (2021) showed that CRA-40/Land is spatially consistent with 22 in situ observations, and has a comparable performance to ERA-Interim.

    6.   Concluding remarks
    • A new global atmospheric reanalysis CRA-40 produced by CMA for the first time is documented. The original 6-h model output covers a period from 1979 to 2018, with a horizontal resolution of TL574 (~34 km) and 64 levels spanning from surface to 0.27 hPa. A large number of reprocessed satellite datasets and widely collected conventional observations were assimilated during the production. It is worth mentioning that, the 13-yr reprocessed AMV products from FY-2C/D/E/G are used in global reanalysis for the first time. Dense conventional observations over China (~120 radiosonde and ~2400 synoptic stations) are also used.

      The original model level output is postprocessed into 47 pressure level products listing in Appendix A using UPP. All the variables are then interpolated to four horizontal resolutions in longitude–latitude projection, including 0.25° × 0.25°, 0.5° × 0.5°, 1.0° × 1.0°, and 2.5° × 2.5°. Daily and monthly products are also produced based on the 6-h products. All these products are publicly released in http://data.cma.cn/CRA. CRA-40/Atmosphere is now running operationally in near real time as a climate data assimilation system for the climate monitoring. In conjunction with CRA-40/Atmosphere, a land product is also produced through a single simulation using the Noah land surface model driven by observations adjusted near the surface. User is recommended to use t2m, q2m, u10, v10, and precipitation from the CRA-40/Land products, as they have merged SYNOP observations using EnOI method. Details on the production of CRA-40/Land will be presented in a separate article.

      The next version of CMA’s global reanalysis will be based on CMA Global Forecast System (CMA-GFS) with the 4DVar DA technique. More satellite radiance data in early years (e.g., those from SSU) will be included to improve high-level temperature analysis. The reprocessed radiances from FY-3 series will also be used in the next generation reanalysis of CMA.

      Acknowledgments. As with other reanalyses, CRA-40 is a collective work of many people from multiple organizations. We express our gratitude to all members involved in the funding projects from the National Meteorological Information Centre (NMIC), National Climate Centre (NCC), National Meteorological Centre (NMC), and National Satellite Meteorological Centre (NSMC) of CMA, and the Institute of Atmospheric Physics (IAP) of Chinese Academy of Sciences, as well as the Nanjing University of Information Science and Technology (NUIST).

      We wish to thank the Central Operations (NCO) of NCEP for making the GSM model publicly available, Developmental Testbed Center (DTC) of NCAR for providing the GSI community software, and also the Research Data Archive (RDA) of NCAR for providing the CFSR input observations.

Reference (58)

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return