Assessment of the Improvement in Observation Precision of GNSS, SLR, VLBI, and DORIS Inputs from ITRF2014 to ITRF2020 Using TRF Stacking Methods

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.