1. Introduction
During the last decades, tremendous progress has been made in developing consistent long-term grid datasets for our understanding of climate variability and climate change when reanalyses of past meteorological observations using modern data assimilation systems developed for numerical weather prediction (NWP) have been used across a wide range of applications, and have become a mainstay of many types of atmospheric research (Schubert et al. 1993; Kalnay et al. 1996; Kistler et al. 2001; Kanamitsu et al. 2002; Gibson et al. 1997; Simmons and Gibson 2000; Uppala et al. 2005). These products provide basic meteorological variables, such as tropospheric pressure heights, humidity and winds, as well as 2-m surface temperature, precipitation, and radiation fluxes. Among those elements, the 2-m surface temperature plays the most important role in investigating the global or regional climate change and variability (Simmons et al. 2004; Schär et al. 2004; Frauenfeld et al. 2005). The reanalysis, however, had problems that made them suboptimal or even unusable for some applications. Perhaps the most serious problem for climate applications was that, while the assimilation system remained unchanged, changes in the observing systems did produce spurious changes in the perceived climate (Bengtsson et al. 2004a, b; Fiorino 1999; WCRP 2000; Trenberth et al. 2001). Therefore, the evaluation and validation of the reanalyzed products, whenever possible, using independent observations, are critical for the latter’s proper application (Smith et al. 2001).
China is a big country with a significant portion of land territory of the world characterized by remarkable topographical gradients and complexity. Recently, China’s daily surface air temperature data with a dense observation network and consistent observation practices have become available to the scientific community. In this paper, we attempt to study the effects of station elevation in calibrating and evaluating reanalyzed surface temperature in mainland China, which are important for evaluating the reliability of reanalyzed products and the performance of numerical model simulation over this region.
2. Data and method
The reanalyzed monthly mean 2-m surface temperature for the period 1979–2001 is obtained from the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR; Kalnay et al. 1996; Kistler et al. 2001) reanalysis at a resolution of approximately 1.875° × 1.875° (192 × 94 Gaussian grid points) and the 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40; Simmons and Gibson 2000; Uppala et al. 2005) on a 2.5° × 2.5° latitude–longitude grid. A daily temperature dataset of 597 stations in mainland China over the same period is utilized for comparison and calibration. A quality control procedure by the National Meteorological Center of the China Meteorological Administration (CMA) was applied to Chinese daily datasets to guarantee the reliability and consistency of the observations (Zhai et al. 1999).
To investigate the topographic effects on reanalysis surface temperature quantitatively, China’s mainland territory is divided into three regions: eastern coastal plain–hill regions, Inner Mongolia Plateau–Loess Plateau–Yunnan Guizhou Plateau, and Tibetan Plateau with an average elevation of less than 500 m, 500–3000 m, and over 3000 m, respectively. Figure 1 is the division of the three regions and the location of the meteorological stations. The density of the stations is lower in the sparsely populated high mountainous and desert areas of west and northwest China.
Since the resolutions of the two reanalyzed grid boxes do not match each other, the 2-m surface temperature from the reanalyzed grid was transformed to a finer 0.5° × 0.5° latitude–longitude grid using bilinear interpolation at first. The elevation of the reanalyzed grid was processed in the same way. To evaluate the reanalysis temperature based on consistent spatial resolution, the station values of the observed temperature and elevation were also converted onto the same regular grid adopting the kriging interpolation technique (Phillips et al. 1992; Goovaerts 1997; Hunter and Meentemeyer 2005).
3. Results
Shown in Fig. 2 are the differences of elevation and surface temperature between the ERA-40 reanalysis and the observations (reanalysis minus observations) before and after topographic correction in mainland China. It is evident that the differences of, not only the elevation, but also the temperature are generally larger in the region with a higher altitude, and vice versa. To investigate the origin of the temperature errors, the differences of elevation between the reanalysis and the observations are presented in Fig. 2a, and the differences of the temperature in summer [June–August (JJA)] and winter [December–February (DJF)] between ERA-40 and the observations are shown in Figs. 2b,c. It is meaningful that the differences of the temperature are nearly inversely proportional to that of the elevation, especially in summer, which means that the interpolated temperature/elevation from native reanalysis or station values is highly dependent on the altitudes of each grid point and topographic complexity. This result may be more obvious in region III where the largest topographic gradient exists.


As a tentative step, “topographical correction” was conducted for the interpolated temperature of reanalyzed products and observations at each grid. Though such a correction may be very simple, the errors introduced by elevation bias can be greatly removed and a dramatic improvement has been achieved for the interpolated temperature (Figs. 2c,e). For example, the Std of the temperature differences between ERA-40 and the observations in summer is reduced by an average of 43.6%. In winter, the reduction is 47.1% in region III, but Std is actually increased in regions I and II (Table 1). Such results indicate that the effects of the topographic correction have obvious regional dependency (i.e., the larger the region’s topographic gradient, the better the corrected effects are expected to be for the interpolated surface temperature).
Figure 3 is the same as Fig. 2 except for the differences between the NCEP–NCAR reanalysis and the observations. The correlation of the elevation differences to the surface temperature is somewhat decreased in some areas like regions I and II, which is obvious especially in winter. However, some interesting results can also be obtained similar to the results mentioned above. The correlation coefficients of the temperature bias with the elevation differences are negative and are over the 99.9% confidence level for all regions and both seasons. Table 1 suggests that the accuracy of the surface temperature in ERA-40 is more dependent on model topography than is the case for the NCEP–NCAR reanalysis. The Std of surface temperature differences is reduced by an average of 38.3% for the three regions in summer, but the reduction is only 11.7% in winter.
Dramatic improvements have been achieved through a topographical correction for the interpolated surface temperature, especially in the region with higher altitude and complex terrain, which indicates that the difference of elevation between the interpolated reanalyzed grid cell and the station mainly accounts for the interpolated surface temperature bias. However, the effects of topographic correction also show a seasonal and regional dependency to a certain extent. For example, in northwest Xinjiang (a province of China) and southwest of the Tibetan Plateau, the differences of interpolated surface temperature increase after topographic correction rather than decrease in DJF, and the effects of topographic correction become worse inversely. Besides the overall elevation and the complexity of the terrain, the station density together with the interpolation method can introduce biases when irregular observations are converted into gridded data. Hence, there are many uncertainties in the interpolated observations in the sparsely populated high mountainous and desert areas of west and northwest China due to the lack of the observational coverage. In addition, inconsistencies in both the ERA-40 and NCEP–NCAR assimilated system are the other important error sources for the interpolated reanalysis. That is, the biases of interpolated surface temperature are results of both the observations and the reanalysis. Therefore, to evaluate the reanalysis dataset in west and northwest China in terms of the biases mentioned above is not convincing enough.
According to the above analysis, the ERA-40 calibrated surface temperature is closer to the observations than the NCEP–NCAR counterparts in most areas of China, especially in west and northwest China. While the NCEP–NCAR reanalysis does not assimilate surface observations such as temperature, moisture, and wind (Kalnay and Cai 2003), ERA-40 does incorporate these variables when available. These differences of treatment for surface observations between two assimilated systems could produce many inconsistencies in their products certainly. Therefore, how to improve evaluation for the different reanalysis datasets is very crucial for their proper application to weather and climate studies. The topographical correction may be one of the important ways to evaluate the uncertainties of reanalysis surface temperature.
4. Conclusions
Based on China’s daily temperature and two reanalyzed datasets from 1979 to 2001, the topographical effects on interpolated surface temperature errors and their correction are discussed in this paper. The major conclusions are summarized as follows:
Obvious biases of interpolated temperature, either from reanalyzed products or from the station observations, are found in mainland China, and ERA-40 is closer to the observations than the NCEP–NCAR reanalysis due to the different treatment of surface observations.
The interpolated temperature errors could be mainly (but not solely) attributed to the biases of interpolated elevation because the former varies very well with local elevation and topographical complexity, which indicates the necessity of calibrating reanalysis surface temperature before evaluating their accuracy and utilizing them in climate change studies.
Outstanding improvements of interpolated surface temperature have been achieved through a “topographical correction,” especially in the boundary areas with a larger altitude difference. The effects of topographic correction for the interpolated surface temperature have a seasonal dependence. The effects of the topographic correction are better in summer than in winter, especially in the sparsely populated high mountainous and desert areas of west and northwest China.
Acknowledgments
Special thanks are given to the ECMWF public Web server (http://data.ecmwf.int/data/) and CIRES Climate Diagnostics Center (http://www.cdc.noaa.gov/) for providing the free reanalysis datasets. This research was supported by National Key Basic Research Program of China (2006CB400500) and the National Natural Science Foundation of China (Grants 40405014 and 40645028).
REFERENCES
Bengtsson, L., K. I. Hodges, and S. Hagemann, 2004a: Sensitivity of the ERA-40 reanalysis to the observing system: Determination of the global atmospheric circulation from reduced observations. Tellus, 56A , 456–471.
Bengtsson, L., S. Hagemann, and K. I. Hodges, 2004b: Can climate trends be calculated from reanalysis data? J. Geophys. Res., 109 .D11111, doi:10.1029/2004JD004536.
Cosgrove, B. A., and Coauthors, 2003: Real-time and retrospective forcing in the North American Land Data Assimilation Systems (NLDAS) project. J. Geophys. Res., 108 .8842, doi:10.1029/2002JD003118.
Fiorino, M., 1999: The impact of observing system changes on the climate-scale variability and temperature in the ECMWF and NCEP reanalysis. Second Int. Conf. of Reanalysis, Reading, England, WMO, 65–68.
Frauenfeld, O. W., T. Zhang, and M. C. Serreze, 2005: Climate change and variability using European Centre for Medium-Range Weather Forecasts reanalysis (ERA-40) temperatures on the Tibetan Plateau. J. Geophys. Res., 110 .D02101, doi:10.1029/2004JD005230.
Gibson, J. K., P. Kallberg, S. Uppala, A. Noumura, A. Hernandez, and E. Serrano, 1997: ERA description. ECMWF Re-Analysis Project Rep. Series 1, 77 pp.
Goovaerts, P., 1997: Geostatistics for Natural Resources Evaluation. Oxford University Press, 467 pp.
Hunter, R. D., and R. K. Meentemeyer, 2005: Climatologically aided mapping of daily precipitation and temperature. J. Appl. Meteor., 44 , 1501–1510.
Kalnay, E., and M. Cai, 2003: Impact of urbanization and land–use change on climate. Nature, 423 , 6939. 528–531.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77 , 437–471.
Kanamitsu, M., W. Ebisuzaki, J. Woollen, S-K. Yang, J. J. Hnilo, M. Fiorino, and G. L. Potter, 2002: NCEP–DOE AMIP-II Reanalysis (R-2). Bull. Amer. Meteor. Soc., 83 , 1631–1643.
Kistler, R., and Coauthors, 2001: The NCEP–NCAR 50-Year Reanalysis: Monthly means CD-ROM and documentation. Bull. Amer. Meteor. Soc., 82 , 247–267.
Phillips, D. L., J. Dolph, and D. Marks, 1992: A comparison of geostatistical procedures for spatial analysis of precipitation in mountainous terrain. Agric. For. Meteor., 58 , 119–141.
Schär, C., P. L. Vidale, D. Luthi, C. Frei, C. Häberli, M. A. Liniger, and C. Appenzeller, 2004: The role of increasing temperature variability in European summer heatwaves. Nature, 427 , 332–336.
Schubert, S. D., R. B. Rood, and J. Pfaendtner, 1993: An assimilated dataset for earth science applications. Bull. Amer. Meteor. Soc., 74 , 2331–2342.
Simmons, A. J., and J. K. Gibson, 2000: The ERA-40 Project plan. ERA-40 Project Rep. Series 1, ECMWF, 63 pp.
Simmons, A. J., and Coauthors, 2004: Comparison of trends and low-frequency variability in CRU, ERA-40, and NCEP/NCAR analyses of surface air temperature. J. Geophys. Res., 109 .D24115, doi:10.1029/2004JD005306.
Smith, S. R., D. M. Legler, and K. V. Verzone, 2001: Quantifying uncertainties in NCEP reanalysis using high-quality research vessel observations. J. Climate, 14 , 4062–4072.
Trenberth, K. E., D. P. Stepaniak, J. W. Hurrell, and M. Fiorino, 2001: Quality of reanalysis in the Tropics. J. Climate, 14 , 1499–1510.
Uppala, S. M., and Coauthors, 2005: The ERA-40 re-analysis. Quart. J. Roy. Meteor. Soc., 131 , 2961–3012.
WCRP, 2000: Proceedings of the second WCRP International Conference on Reanalyses. WMO/Tech. Doc. 959, 465 pp.
Zhai, P., A. Sun, F. Ren, X. Liu, B. Gao, and Q. Zhang, 1999: Changes of climate extremes in China. Climatic Change, 42 , 1. 203–218.
Division of the three subregions in mainland China with topography (m) and station distribution.
Citation: Journal of Climate 21, 6; 10.1175/2007JCLI1463.1
The differences in elevation and in temperature between ERA-40 and observed values (ERA-40 − Obs) interpolated from the meteorological station using kriging interpolation over China: (a) The difference for elevation (m) and the differences for (b),(d) uncorrected and (c),(e) corrected temperature (°C).
Citation: Journal of Climate 21, 6; 10.1175/2007JCLI1463.1
As in Fig. 2 but between the NCEP–NCAR reanalysis and observed values.
Citation: Journal of Climate 21, 6; 10.1175/2007JCLI1463.1
Correlation and Std of elevation and temperature difference between interpolated reanalysis value and observations in three regions.