International terrestrial reference frame (ITRF) input data, generated by Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR), Very Long Baseline Interferometry (VLBI), and Doppler Orbitography and Radiopositioning integrated by satellite (DORIS) combination centers (CCs), are considered to be relatively high-quality and accurate solutions. Every few years, these input data are submitted to the three ITRS combination centers, namely Institut Géographique National (IGN), Deutsches Geodätisches Forschungsinstitut at the Technische Universität München (DGFI-TUM), and Jet Propulsion Laboratory (JPL), to establish a multi-technique combined terrestrial reference frame (TRF). Generally, these solutions have undergone three rounds of outlier removal: the first at the technique analysis centers during solution generations and the second during the technique-specific combination by the CCs; ITRS CCs then perform a third round of outlier removal and preprocessing during the multi-technique combination of TRFs. However, since the primary objective of CCs is to release the final TRF product, they do not emphasize the publication of analytical preprocessing results, such as the outlier rejection rate. In this paper, our specific focus is on assessing the precision improvement of ITRF input data from 2014 to 2020, which includes evaluating the accuracy of coordinates, the datum accuracy, and the precision of the polar motions, for all four techniques. To achieve the above-mentioned objectives, we independently propose a TRF stacking approach to establish single technical reference frameworks, using software developed by us that is different from the ITRF generation. As a result, roughly 0.5% or less of the SLR observations are identified as outliers, while the ratio of DORIS, GNSS, and VLBI observations are below 1%, around 2%, and ranging from 1% to 1.2%, respectively. It is shown that the consistency between the SLR scale and ITRF has improved, increasing from around −5 mm in ITRF2014 datasets to approximately −1 mm in ITRF2020 datasets. The scale velocity derived from fitting the VLBI scale parameter series with all epochs in ITRF2020 datasets differs by approximately 0.21 mm/year from the velocity obtained by fitting the data up to 2013.75 because of the scale drift of VLBI around 2013. The decreasing standard deviations of the polar motion parameter (XPO, YPO) offsets between Stacking TRFs and 14C04 (20C04) indicate an improvement in the precision of polar motion observations for all four techniques. From the perspective of the weighted root mean square (WRMS) in station coordinates, since the inception of the technique, the station coordinate WRMS of DORIS decreased from 30 mm to 5 mm for X and Y components, and 25 mm to 5 mm for the Z component; SLR WRMS decreased from 20 mm to better than 10 mm (X, Y and Z); GNSS WRMS decreased from 4 mm to 1.5 mm (X and Y) and 5 mm to 2 mm (Z); while VLBI showed no significant change.

ITRF input data which are integrated by GNSS, SLR, VLBI, DORIS combination centers are considered to be relatively high-quality and accurate solutions. However, when utilizing these inputs, one still needs to identify outliers, rescale inaccurate covariance matrix and evaluate the precision of the observed datum information. To achieve the above-mentioned objectives, we propose a terrestrial reference frame (TRF) stacking approach to establish the single technical reference frameworks for the ITRF2014 and ITRF2020 datasets of all four technologies. As a result, roughly 0.5% or less of the SLR observations are identified as outliers, while the ratio of DORIS, GNSS, and VLBI observations are below 1%, around 2%, and ranging from 1% to 1.2%, respectively. The post-rescaling covariance scale factors are 25.07, 27.25, 18.84, 6.98 for GNSS, SLR, VLBI, and DORIS in ITRF2014 datasets, and 8.95, 14.9, 16.8, 7.78 in ITRF2020 datasets, respectively. It is shown that the consistency between the SLR scale and ITRF has improved, increasing from around -5mm in ITRF2014 datasets to approximately -1mm in ITRF2020 datasets. The scale velocity derived from fitting the VLBI scale parameter series with all epochs in ITRF2020 datasets differs by approximately 0.21mm/year from the velocity obtained by fitting the data up to 2013.75 because of the scale drift of VLBI at around 2013. The decreasing standard deviations of the Polar motion parameter (XPO, YPO) offsets between Stacking TRFs and 14C04 (20C04) indicating an improvement in the precision of polar motion observations for all of the four techniques. From the perspective of the Weighted Root Mean Square in station coordinates, the measurement precision of GNSS, SLR, and DORIS techniques has improved, while VLBI shows no significant change.