PSR J1012+5307, a millisecond pulsar in orbit with a helium white dwarf (WD), has been timed with high precision for about 25 years. One of the main objectives of this long-term timing is to use the large asymmetry in gravitational binding energy between the neutron star and the WD to test gravitational theories. Such tests, however, will be eventually limited by the accuracy of the distance to the pulsar. Here, we present VLBI (very long baseline interferometry) astrometry results spanning approximately 2.5 years for PSR J1012+5307, obtained with the Very Long Baseline Array as part of the MSPSRPI project. These provide the first proper motion and absolute position for PSR J1012+5307 measured in a quasi-inertial reference frame. From the VLBI results, we measure a distance of $0.83^{+0.06}_{-0.02}$kpc (all the estimates presented in the abstract are at 68% confidence) for PSR J1012+5307, which is the most precise obtained to date. Using the new distance, we improve the uncertainty of measurements of the unmodeled contributions to orbital period decay, which, combined with three other pulsars, places new constraints on the coupling constant for dipole gravitational radiation $\kappa_D=(-1.7\pm1.7)\times 10^{-4}$ and the fractional time derivative of Newton's gravitational constant $\dot{G}/G = -1.8^{\,+5.6}_{\,-4.7}\times 10^{-13}\,{\rm yr^{-1}}$ in the local universe. As the uncertainties of the observed decays of orbital period for the four leading pulsar-WD systems become negligible in $\approx10$ years, the uncertainties for $\dot{G}/G$ and $\kappa_D$ will be improved to $\leq1.5\times10^{-13}\,{\rm yr^{-1}}$ and $\leq1.0\times10^{-4}$, respectively, predominantly limited by the distance uncertainties.

Newton's gravitational constant $G$ may vary with time at an extremely low level. The time variability of $G$ will affect the orbital motion of a millisecond pulsar in a binary system and cause a tiny difference between the orbital period-dependent measurement of the kinematic distance and the direct measurement of the annual parallax distance. PSR J0437$-$4715 is the nearest millisecond pulsar and the brightest at radio. To explore the feasibility of achieving a parallax distance accuracy of one light-year, comparable to the recent timing result, with the technique of differential astrometry, we searched for compact radio sources quite close to PSR J0437$-$4715. Using existing data from the Very Large Array and the Australia Telescope Compact Array, we detected two sources with flat spectra, relatively stable flux densities of 0.9 and 1.0 mJy at 8.4 GHz and separations of 13 and 45 arcsec. With a network consisting of the Long Baseline Array and the Kunming 40-m radio telescope, we found that both sources have a point-like structure and a brightness temperature of $\geq$10$^7$ K. According to these radio inputs and the absence of counterparts in the other bands, we argue that they are most likely the compact radio cores of extragalactic active galactic nuclei rather than Galactic radio stars. The finding of these two radio active galactic nuclei will enable us to achieve a sub-pc distance accuracy with the in-beam phase-referencing very-long-baseline interferometric observations and provide one of the most stringent constraints on the time variability of $G$ in the near future.

The massive photometric data collected from multiple large-scale sky surveys offer significant opportunities for measuring distances of celestial objects by photometric redshifts. However, catastrophic failure is still an unsolved problem for a long time and exists in the current photometric redshift estimation approaches (such as $k$-nearest-neighbor). In this paper, we propose a novel two-stage approach by integration of $k$-nearest-neighbor (KNN) and support vector machine (SVM) methods together. In the first stage, we apply KNN algorithm on photometric data and estimate their corresponding z$_{\rm phot}$. By analysis, we find two dense regions with catastrophic failure, one in the range of z$_{\rm phot}\in[0.3,1.2]$, the other in the range of z$_{\rm phot}\in [1.2,2.1]$. In the second stage, we map the photometric input pattern of points falling into the two ranges from original attribute space into a high dimensional feature space by Gaussian kernel function in SVM. In the high dimensional feature space, many outlier points resulting from catastrophic failure by simple Euclidean distance computation in KNN can be identified by a classification hyperplane of SVM and further be corrected. Experimental results based on the SDSS (the Sloan Digital Sky Survey) quasar data show that the two-stage fusion approach can significantly mitigate catastrophic failure and improve the estimation accuracy of photometric redshifts of quasars. The percents in different |$\Delta$z| ranges and rms (root mean square) error by the integrated method are $83.47\%$, $89.83\%$, $90.90\%$ and 0.192, respectively, compared to the results by KNN ($71.96\%$, $83.78\%$, $89.73\%$ and 0.204).