{'fraction_non_alphanumeric': 0.04030801730893727, 'fraction_numerical': 0.033766056765136804, 'mean_word_length': 4.955762987012987, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Impact of the Observation Data Assimilation on Atmospheric Reanalyses over Tibetan Plateau and Western Yunnan-Guizhou Plateau

: Three modern atmospheric reanalyses with different input observation data (NOAA– CIRES 20th Century Reanalysis (20CR), Japanese 55-year Reanalysis (JRA-55), and JRA-55C) were compared against the independent radiosonde observations over the Tibetan Plateau (TP) and the western Yunnan–Guizhou Plateau (YGP) from the China-Japan Meteorological Disaster Reduction Cooperation (JICA/Tibet) Center Project in the summer of 2018 to investigate the effects of the assimilation of the observation data on the quality and accuracy of the reanalyses in the troposphere. The results indicate that the mean biases and mean root-mean-square errors of horizontal wind, temperature, and speciﬁc humidity signiﬁcantly decreased when comparing the 20CR reanalysis method (assimilating only surface pressure) to the JRA-55C method (assimilating conventional surface and upper-air observations) and the JRA-55 (assimilating conventional and satellite observations), and the bias spreads of these aboveground variables in JRA-55C and JRA-55 were cut to almost half of those observed in 20CR. However, the mean biases and uncertainties varied little from JRA-55C to JRA-55. This means that the assimilation of conventional observation data plays a vital role in the quality of reanalyses for the troposphere over these data-sparse plateaus. It was also found that the temperature and speciﬁc humidity near the ground over TP showed larger mean biases and bias spans than those over YGP, likely due to the sparser surface observation over TP. Author Contributions: Data curation, Y.Z. and X.B.; Formal analysis, X.B.; Investigation, X.B.; Methodology, F.Z. and X.B.; Validation, Y.C.; Writing—original draft, X.B.; Writing—review & editing,


Introduction
The Tibetan Plateau (TP), the largest and highest plateau on earth, plays an essential role in the regional and even the global climate [1][2][3][4][5][6]. The Yunnan-Guizhou Plateau (YGP), located southeast of TP, also affects the climate and weather (e.g., Indian monsoon, precipitation) over Central Asia [7,8]. The use of observations alone is very limiting in climate and weather studies over these plateaus due to the spatial and temporal sparsity of the data. As a better choice, assimilation reanalyses, which combine the benefits of both forecast models and observations, have been frequently used in atmospheric science studies for TP regions [4,[9][10][11][12].
The quality and performance of reanalysis datasets are related to many factors, such as the model structure, assumptions, and data assimilation methods, as well as the input observation data sources, etc. [13][14][15][16]. The quantity and quality of the observation data assimilated as part of the reanalyses and consisting of conventional and satellite data are major factors influencing the performance of reanalyses [16]. Reanalyses, especially in the stratosphere, are sensitive to changes in the number and type of input satellite observations [17][18][19][20]. The performance of various atmospheric reanalyses of TP regions has been assessed in previous studies [13,14,[21][22][23]. It was found that updating the data assimilation systems, improving the quality, and increasing the amount of input observations have a positive effect on the quality of reanalysis data across these regions. The quality of almost all the reanalyses was improved following the updating of the satellite observing system at the end of 1998 [24,25]. Bao and Zhang [14] found that the quality of reanalysis over TP had improved significantly from the early period (1998) to the later period (2008 and 2015-2016) mainly due to the huge increase in the amount and quality of input satellite observations. They also note that the Japanese 55-year Reanalysis (JRA-55) [26], a global atmospheric reanalysis with a~55 km resolution led by the Japan Meteorological Agency (JMA), showed relatively better quality than other reanalyses in 1998 because the soundings collected in the field campaign were assimilated into the JRA reanalysis. With its use of conventional and satellite observation data assimilation, JRA-55 is one of the widely used modern atmospheric reanalysis products. It has been verified to have good quality over TP [14]. However, few studies have documented to what extent the assimilation of observation data affects the reanalysis products. The quantitative differences between reanalyses assimilating conventional data with and without satellite data over TP and YGP are still unknown.
The NOAA-CIRES 20th Century Reanalysis (20CR) [27], spanning 1851 to 2012, used an ensemble Kalman filter (EnKF) approach to assimilate only surface pressure. JRA-55 and 20CR essentially represent the reanalysis product performance with and without the assimilation of observation data, respectively. Additionally, JRA-55C [28], which uses the same data assimilation system as JRA-55 to assimilate only conventional surface and upper-air observations, can provide insights on the quality of reanalysis in the absence of satellite observations. The China-Japan Meteorological Disaster Reduction Cooperation Research Center Project (JICA/Tibet Project [29]) deployed numerous radiosonde observations over TP and western YGP in the summer of 2008 ( Figure 1). None of these sounding observations were introduced to any of the reanalyses. This dataset can thus be used to evaluate the sensitivity of reanalyses over TP and YGP according to different input observation sets. For this study, these different sets consist of (a) only surface pressures with the 20CR reanalysis, (b) only conventional observations with the JRA-55 reanalysis, and (c) both conventional and satellite data with JRA-55C. Moreover, the quality of reanalysis products over YGP is seldom discussed. It is also not very clear whether reanalysis products over YGP show better or worse performance compared with those over TP. The spatial density of conventional observation data also affects the quality of reanalyses. Reanalyses show different quality levels for regions with different data densities; for example, the root-meansquare errors (RMSEs) are lower for horizontal winds and relative humidity over data-rich Southern China, compared to those from the data-sparse TP [30]. This study is an extension of our earlier studies [13,14].

Data and Methodology
Three reanalysis datasets-consisting of 20CR, JRA-55, and JRA-55C-are intercompared in this study. 20CR [26] is a reanalysis project jointly conducted by the National Oceanic and Atmospheric Administration's (NOAA) Physical Sciences Laboratory and the Cooperative Institute for Research in Environmental Sciences (CIRES) at the University of Colorado, which applies an ensemble Kalman filter scheme and assimilates only surface pressure observations into NOAA's Global Forecast System. The version used in this study is V2c, with a grid resolution of 75 km, covering 1851 to 2012. JRA-55 [27] is a global atmospheric reanalysis led by the Japan Meteorological Agency (JMA). It uses a 4D-Var data assimilation method to assimilate both conventional and satellite observations. JRA-55 spans 1958 to the present and has ~55km of spatial resolution and 60 vertical levels. JRA-55C [28] stems from the same system as JRA-55, with the assimilation of only conventional surface and upper-air observations. It produces a reanalysis dataset, which is more homogeneous and not influenced by variations of the satellite observing systems throughout history.
To evaluate the quality of these three reanalyses, the enhanced independent radiosonde sounding observations recorded during the intensive observing period (IOP) of the JICA/Tibet project from 20 June to 19 July 2018 were used in this study. Sounding observations were gathered every 6 h (0000, 0600, 1200, 1800 UTC) at 15 locations covering the Tibetan Plateau and the western Yunnan-Guizhou Plateau ( Figure 1). There are nine sparsely and highly distributed stations (>3000 m) located on TP, and six relatively closer and lower stations (1300-2400 m) located on YGP. The 0000 and 1200 UTC observations

Data and Methodology
Three reanalysis datasets-consisting of 20CR, JRA-55, and JRA-55C-are intercompared in this study. 20CR [26] is a reanalysis project jointly conducted by the National Oceanic and Atmospheric Administration's (NOAA) Physical Sciences Laboratory and the Cooperative Institute for Research in Environmental Sciences (CIRES) at the University of Colorado, which applies an ensemble Kalman filter scheme and assimilates only surface pressure observations into NOAA's Global Forecast System. The version used in this study is V2c, with a grid resolution of 75 km, covering 1851 to 2012. JRA-55 [27] is a global atmospheric reanalysis led by the Japan Meteorological Agency (JMA). It uses a 4D-Var data assimilation method to assimilate both conventional and satellite observations. JRA-55 spans 1958 to the present and has~55km of spatial resolution and 60 vertical levels. JRA-55C [28] stems from the same system as JRA-55, with the assimilation of only conventional surface and upper-air observations. It produces a reanalysis dataset, which is more homogeneous and not influenced by variations of the satellite observing systems throughout history.
To evaluate the quality of these three reanalyses, the enhanced independent radiosonde sounding observations recorded during the intensive observing period (IOP) of the JICA/Tibet project from 20 June to 19 July 2018 were used in this study. Sounding observations were gathered every 6 h (0000, 0600, 1200, 1800 UTC) at 15 locations covering the Tibetan Plateau and the western Yunnan-Guizhou Plateau ( Figure 1). There are nine sparsely and highly distributed stations (>3000 m) located on TP, and six relatively closer and lower stations (1300-2400 m) located on YGP. The 0000 and 1200 UTC observations at Naqu, Lasa, Changdu, Ganzi, Xichang, Tengchong, Kunming, and Mengzi, which belong to the standard observing network and were assimilated in JRA-55 and JRA-55C, were excluded for completely independent verifications. Four types of radiosondes, these being the GTS1 (Tengchong, Kunming, Hongyuan, Ganzi, Dingri, Xichang, Naqu, Lasa), 59-701 (Mengzi, Linzhi, Changdu), Vaisala RS92 (Gaize, Litang, Dali), and TD2-A digital sonde (Lijiang), were applied for the enhanced sounding observations during the 2008 JICA/Tibet Project. In the World Meteorological Organization (WMO) intercomparison of high-quality radiosonde systems implemented at Yangjiang, China in 2010, GTS1 and Vaisala RS92 had the higher scores of 4.5-5 for temperature (<16 km), humidity (T > −40 • C), wind (troposphere) in the intercomparison of 13 radiosonde systems from nine countries. The score for evaluating operational performance was set to range one to six, but no radiosonde system scored more than five [31]. Bao and Zhang [14] have verified three of them (GTS1, Vaisala RS92, and 59-701) as having high accuracy and good quality and stability; the interinstrument difference between these radiosondes was small in wind speed (<± 0.2 m/s) and temperature (<± 0.2 • C) in the mid-upper troposphere and slightly large in relative humidity (~10%) above −40 • C. The TD2-A radiosondes were also frequently used during that period.
The enhanced sounding observation data provided by the JICA/Tibet team include four quality-controlled conventional above ground variables (wind direction, wind speed, temperature, and dewpoint depression) interpolated to eight standard vertical levels (600, 500, 400, 300, 250, 200, 150, and 100 hPa) over TP and nine levels (700, 600, 500, 400, 300, 250, 200, 150, and 100 hPa) over YGP. As in Bao and Zhang [14], we used the wind direction and wind speed to derive the horizontal meridional and zonal winds (U and V). The specific humidities (Q) were calculated from the temperatures and dewpoint temperatures in accordance with the method used in the European Centre for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim) forecast model [32,33].
We used the same methods as in Bao and Zhang [13,14] to make a direct comparison between the gridded reanalyses and discrete soundings. Firstly, we interpolated the pressure coordinate reanalysis data to the standard pressure levels at each of the sounding locations at the same synoptic times. Secondly, we rejected the abnormal values under a subjective criterion of the differences between observation and JRA-55 (i.e., if ∆U > 30 m s −1 , and/or ∆V > 30 m s −1 , as well as ∆T > 20 • C, Q will be rejected when the temperature T is eliminated). Thirdly, we compared the sounding observations and the three reanalyses in terms of mean, bias, root-mean-square error (RMSE), and spread of the reanalysis-minusobservation (R-O).  (Figure 2b,d,f,h) show the averaged vertical profiles of U, V, T, and Q from the independent soundings at TP (YGP) and the three reanalysis products. At a glance, 20CR provided significantly different mean profiles of these above ground variables, with the exception of temperature, as compared to the other reanalyses and soundings. As mentioned in Bao and Zhang (2019), the mean observed zonal winds were predominantly westerlies with a maximum of~10 m/s around 200 hPa over TP (Figure 2a). Meanwhile, the averaged observed zonal winds over YGP were weak westerlies near the ground and easterlies with a maximum of~14 m/s near the tropopause (Figure 2b). The mean observed meridional winds over TP and YGP showed a similar vertical profile, both of them with northerly winds throughout almost all levels except the very weak southerly flows near the ground (Figure 2c,d). JRA-55 and JRA-55C closely followed well those averaged horizontal winds over TP and YGP. The mean horizontal winds of 20CR differed notably from the verifying observed sounding mean and the other reanalyses. 20CR significantly undervalues the zonal winds throughout the vertical profile, overestimates (undervalues) the meridional winds near the tropopause over TP and YGP (in the mid-troposphere (250-700 hPa) over YGP) (Figure 2a-d). The mean temperature profiles and the mean specific humidity profiles of each reanalysis and soundings all show expected decreases with increasing height (Figure 2e-h). The mean temperature profiles from the sounding and the three reanalysis products were hardly distinguishable from each other. Therefore, these reanalyses represent very well the mean dry, static stability (Figure 2e,f), as also found in Bao and Zhang [13,14]. The mean specific humidity profiles of JRA-55 and JRA-55C nearly overlapped with the mean soundings over TP and YGP, but the mean specific humidity from 20CR as wetter (larger) than the soundings observed over TP and YGP (Figure 2g,h).

Mean Profiles
Atmosphere 2021, 12, x FOR PEER REVIEW 5 of 13 dervalues) the meridional winds near the tropopause over TP and YGP (in the mid-troposphere (250-700 hPa) over YGP) (Figure 2a-d). The mean temperature profiles and the mean specific humidity profiles of each reanalysis and soundings all show expected decreases with increasing height (Figure 2e-h). The mean temperature profiles from the sounding and the three reanalysis products were hardly distinguishable from each other. Therefore, these reanalyses represent very well the mean dry, static stability (Figure 2e,f), as also found in Bao and Zhang [13,14]. The mean specific humidity profiles of JRA-55 and JRA-55C nearly overlapped with the mean soundings over TP and YGP, but the mean specific humidity from 20CR as wetter (larger) than the soundings observed over TP and YGP (Figure 2g,h).

Mean Biases and RMSEs
In order to quantitatively assess the quality and utility of these reanalysis products, the mean biases and RMSEs between each reanalysis and the observed soundings (R-O) were calculated (Figures 3 and 4).

Mean Biases and RMSEs
In order to quantitatively assess the quality and utility of these reanalysis produc the mean biases and RMSEs between each reanalysis and the observed soundings (Rwere calculated (Figures 3 and 4).

Horizontal Wind
Not surprisingly, 20CR had the largest mean meridional wind biases. It underestimated the westerlies at almost all levels, with a maximum bias of~−4 m/s between 200-400 hPa over TP and a relatively smaller maximum bias of~3 m/s between 200-300 hPa over YGP. The mean meridional wind biases for JRA-55C and JRA-55 were rather small and mostly within 1 m/s over TP, and 2 m/s over YGP throughout the vertical column (Figure 3a,b). The zonal wind biases were also the largest among the three reanalysis products over both regions, with the bias over YGP being larger. The  Figure 3d). Furthermore, the difference in V biases between JRA-55 and JRA-55C over YG was larger than that over TP. The zonal wind biases for JRA-55 and JRA-55C were similar over TP below 300 hPa, but negative and positive, respectively, above that level (Figure 3c). JRA-55 and JRA-55C had negative (positive) zonal wind biases below (above) 300 hPa, whereas the zonal wind biases from JRA-55 was slightly larger (smaller) than JRA-55C, below (above) 300 hPa over YGP (Figure 3d). The obvious differences in the meridional and zonal wind between 20CR and two other reanalyses were also shown through RMSE (R-O) (Figure 4a-d). 20CR had nearly twice the RMSE for meridional wind compared to JRA-55 and JRA-55C. The meridional wind RMSE of 20CR peaked at 8 m/s at 150 hPa over TP and about 7 m/s at 150-200 hPa over YGP, whereas the meridional wind RMSEs of JRA-55 and JRA-55C were~3 m/s, with little variation over TP and YGP. The zonal wind RMSEs showed similarities with the meridional wind RMSEs. JRA-55 had slightly smaller meridional and zonal wind RMSEs than JRA-55C, with the relatively obvious difference between 150-400 hPa. On the whole, 20CR had the largest RMSEs for horizontal winds, whereas JRA-55 had the smallest RMSEs, with tiny differences compared to JRA-55C

Temperature
The largest temperature biases over TP and YGP were also seen for 20CR (Figure 3e,f). 20CR has a positive (warm) bias compared to sounding observations below 150 hPa, where the positive (warm) was around 2 • C (1.5 • C) over TP (YGP), and a negative (cold) bias of around −2 • C (−1.5 • C) near 100 hPa. JRA-55C had the smallest negative (cold) bias over TP, which was about −2 • C near the surface and which generally decreased with height, and was almost negligible around 400 hPa and above 200 hPa. JRA-55 had a similar bias with JRA-55C, but was slightly larger than JRA-55C above 400 hPa. The T biases over YGP showed similar features to those over TP, but all three reanalyses had smaller bias below 500 hPa, except near the surface where the bias was around zero (Figure 3e,f).
The temperature RMSEs of 20CR and the two other reanalyses showed great differences over TP and YGP (Figure 4e,f). Clearly, the largest was for 20CR, from a maximum of~3.2 • C at 600 hPa to a valley value of~2 • C around 300 hPa, and then rise to~3 • C at 100 hPa. The smaller overall RMSEs were for JRA-55 and JRA-55C, both of which were very close to each other. The RMSEs were~4 • C near the surface and around 1.2 • C above 500 hPa over TP. 20CR over YGP also had the largest RMSE, from~1.8 • C at 700 hPa, increasing to 2.6 • C at the tropopause. The temperature RMSEs of JRA-55 and JRA-55C were only one-third to one-half of that of 20CR.
Overall, these temperature biases and RMSEs largely diminished from 20CR (assimilating surface pressure only) to JRA-55C (assimilating conventional data). Although assimilating both conventional and satellite data, JRA-55 did not show a notable reduction, but it showed slightly larger biases above 400 hPa than JRA-55C. Comparing radiosonde temperatures in China with the temperatures from eight reanalysis datasets, Guo et al. [34] also found that 20CR had obviously higher temperatures than observations at almost all vertical levels except 200 hPa, and significant differences with other reanalyses assimilating conventional and satellite data.

Specific Humidity
The mean biases in specific humidity (Q) also showed huge differences between 20CR and the two other reanalyses, whereas JRA-55 and JRA-55C agree closely with each other (Figure 3g,h). 20CR had the overall largest bias in moisture content over TP and YGP. The bias in 20CR had a positive maximum of~2.2 g/kg at 600 hPa (700 hPa) and decreased with height to almost zero above 200 hPa over TP (YGP). The mean Q biases in JRA-55 and JRA-55C over TP were moist (positive) overall, with the maximum value of~0.9 g/kg near the surface and near-zero above 500 hPa. The mean Q biases in JRA-55 and JRA-55C over YGP varied between −0.2 g/kg and 0.2 g/kg. In contrast with the maximum bias over TP, the biases near the surface in JRA-55 and JRA-55C was almost negligible near the surface, indicating that JRA-55 and JRA-55C had larger biases of specific humidity over the higher TP.
Similar to wind and temperature, specific humidity over TP and YGP also showed greatly different RMSE profiles between 20CR and the two others, whereas the RMSE profiles were similar in JRA-55 and JRA-55C (Figure 4g,h). 20CR had the largest RMSE of specific humidity, and its maximum value was~3.2 (3) g/kg near the surface, which was nearly 1.5 times that of JRA-55 and JRA-55C over TP (YGP). Note that the positive maxima in JRA-55 and JRA-55C near the surface over TP (2.1 g/kg) were larger than that over YGP (1.5 g/kg), which means that these analyses were not consistent over the different regions with different heights and different observational data densities.

Variations of Uncertainty and Bias
The reanalysis-to-reanalysis variations in bias spreads for TP and YGP are summarized in Figure 5. The results further support the evidence of the prominent decrease in bias spread from 20CR to JRA-55 and JRA-55C. The zonal and meridional wind biases for 20CR over TP and YGP were significantly greater than those for JRA-55C and JRA-55. Such differences were more significant in the upper troposphere, where the span of wind bias for JRA-55C and JRA-55 was cut to almost half of that for 20CR (Figure 5a,b,e,f). For temperature, the bias spreads presented a relatively small reduction in the trend from 20CR to JRA-55 and JRA-55C (Figure 5c,g). The temperature biases of 20CR were mostly positive between 150-500 hPa (150-600 hPa) over TP (YGP). Almost 80% of the temperature biases of JRA-55C and JRA-55 were negative at almost all altitude ranges over TP, whereas these temperature biases showed a more normal distribution over YGP. For specific humidity, more than 85% of the humidity biases of 20CR were more than zero over TP and YGP. By comparison, the specific humidity biases of JRA-55 and JRA-55C had a smaller span and were near to a normal distribution. Moreover, the bias spreads of temperature and specific humidity near the surface were the maximum throughout the vertical column over TP, which was significantly greater than that over YGP (Figure 5d,h).
lating upper-air observations. Moreover, 20CR exhibited a similar performance with the older generation reanalyses (prior to 2010) of the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) and the Japan Meteorological Agency Climate Data Assimilation System (JCDAS) in the Antarctic Peninsula region [37]. However, compared with the other new-generation reanalyses, 20CR has obviously larger biases and larger RMSE for temperature, specific humidity, and wind in the Antarctic Peninsula region [38].  Evidently, the 20CR reanalysis (without conventional and satellite data assimilation) had the largest analysis uncertainty, according to the bias spreads of these four aboveground variables, then the two other reanalysis methods. The differences between JRA-55C and JRA-55 were rather small for the quality of these variables, although the satellite data was also used in JRA-55. He et al. [35] have developed a method of regional reanalysis over TP based on an ensemble-based data assimilation system and found prototype reanalysis assimilating non-radiance observations had better accuracy in analyzing atmospheric humidity than the leading modern atmospheric reanalyses-the fifth generation ECMWF atmospheric reanalysis of the global climate (ERA5) and the ERA-Interim analysis.
The significant effect of conventional observations on a multi-year reanalysis (the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2)) was also noted in some recent studies which employed the forecast sensitivity-based observation impact (FSOI) technique [36,37]. Diniz and Todling [36] indicated that conventional observations played a major role in reducing forecast errors throughout the nearly 40-year reanalysis period. Diniz et al. [36] applied the same FSOI technique to MERRA-2 over the Amazon basin and found that more than half of the forecast error reduction was related to conventional observations after 1999. Excepting the differences in input observations, differences in the forecast models and data assimilation methods also affect the quality of reanalyses. These differences were much smaller than those observed between 20CR and JRA-55 in 2008. This is an example showing that the differences in forecast models and data assimilation methods are minor in impact on the performance of reanalyses as compared to differences of assimilating only surface observations and of also assimilating upper-air observations. Moreover, 20CR exhibited a similar performance with the older generation reanalyses (prior to 2010) of the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) and the Japan Meteorological Agency Climate Data Assimilation System (JCDAS) in the Antarctic Peninsula region [37]. However, compared with the other new-generation reanalyses, 20CR has obviously larger biases and larger RMSE for temperature, specific humidity, and wind in the Antarctic Peninsula region [38].

Conclusions
In this study, we evaluated and compared three reanalysis products that assimilated different observations, namely, 20CR (using only surface pressure), JRA-55C (using only conventional observations), and JRA-55 (using conventional and satellite observations), by validating them against independent sounding observations from the 2008 JICA/Tibet project. Our purpose was to determine the influence of the assimilation of observed data (especially conventional data) on the quality of reanalysis products in the troposphere over the Tibetan Plateau and the western Yunnan-Guizhou Plateau. Secondly, in this study, we explored the degrees to which these reanalyses differently impacted the atmospheric variables over the high Tibetan Plateau (>3000 m) and the relatively low Yunnan Guizhou Plateau (about 1300-2400 m).
It was found that 20CR, with no use of conventional and satellite observation, had the largest mean biases and uncertainties for the horizontal wind, temperature, and specific humidity over TP and YGP. 20CR had significantly overvalued temperatures and specific humidities throughout almost all levels, along with undervalued westerlies and upperlevel easterlies. For JRA-55C, assimilating conventional surface and upper-air observations, and JRA-55, assimilating conventional and satellite data, the mean biases and RMSEs of these aboveground variables were generally cut down to 50-70% of those found in 20CR. The bias spreads had a smaller span and neared a normal distribution compared to 20CR. However, although the RMSEs of horizontal winds in JRA-55 were slightly smaller than those in JRA-55C, a notable decreasing trend for the mean biases and RMSEs of these variables was not present from JRA-55C to JRA-55, although satellite data was also included in JRA-55. The biases between JRA-55 and JRA-55C had no significant differences, especially for temperature and humidity. This shows that the assimilation of conventional observation data has a significant effect on quality improvements for reanalysis products in the troposphere over these data-sparse plateaus. Adding satellite data assimilation to conventional data assimilation has a weak positive impact on the upper-level horizontal winds but has a negligible impact on temperature and specific humidity. However, this does not mean that the effect of satellite data assimilation by itself is small; the result will be more reliable if the reanalysis data assimilated only satellite observations added into the comparison. In addition, the reanalyses showed similar overall performances over TP and YGP. In general, JRA-55 and JRA-55C over both TP and YGP were able to reproduce well the vertical profiles of the mean horizontal wind, temperature, and specific humidity against the verifying radiosonde observations with small mean biases, although 20CR had up to twice the biases of JRA-55 and JRA-55C for these variables. However, compared to YGP, the temperature and specific humidity near the surface in JRA-55 and JRA-55C showed greater uncertainties by having higher RMSEs and larger spans of bias spread. This is probably due to the fact that the near-surface observational data for TP were more sparse than those for YGP.