Engineering PapersSearch

Engineering topics

Leonid Petrov

Publications and source records attributed to Leonid Petrov.

Investigating Biases in VLBI Clock Functions and Position Solutions

A systematic bias causing a statistically significant drift between group delay and phase delay measurements in geodetic Very Long Baseline Interferometry (VLBI)has long been present. The cause of this drift remains unknown, although it is readily reproduced in VLBI experiments including a short baseline such as those conducted with the twin 13-meter telescopes at the Wettzell observatory, given the station names WETTZ13N and WETTZELL. This statistical incompatibility has been an obstacle in adopting phase delays in routine VLBI processing. The aim of this study is to provide additional evidence to aid in identifying the source of this systematic bias. A series of four regular VLBI experiments are processed with both group delays and phase delays, and post fit residuals, estimated positions, and clock functions are presented with and without phase calibration applied. These results are also computed for individual intermediate frequency bands.

Geodesy

Automated Detection of Spurious Signals in VLBI Phase Calibration Data

In this memorandum, a set of processing strategies for automatic masking of phase calibration tones is outlined as implemented in the software package PIMA in the task Generate Phase Calibration Mask, or GEPM. The task relies on a robust procedure of cleaning phase calibration data before employing several mathematical strategies designed to selectively identify spurious signals from phase calibration tones. These strategies were derived as more rigorous implementations of heuristics traditionally used by analysts in manually identifying problematic phase calibration data. The task is intended to automate the process of generating a phase calibration mask and in so doing increase the speed and regularity of VLBI analysis. At the outset of the project, a series of goals were identified to evaluate the success of this mask generation. This included the development of an algorithm to identify and mask short-term (defined as less than 10 seconds in length) spurious signals affecting phase calibration data, the development of an algorithm to identify and mask phase calibration tones affected by constant radio-frequency interference, the implementation of a detection scheme for identifying large jumps in phase calibration phase caused by clock breaks, and finally a method of reporting phase calibration health metrics to the user. Processing of over 10 experiments has demonstrated that the final form of the task GEPM as detailed in this technical memorandum satisfies each of these conditions and satisfactorily performs the task of automatic phase calibration data masking, although the wide variety in quality and characteristics of phase calibration data makes a single solution to the problem quite difficult. A series of user inputs have therefore been defined to assist analysts in tailoring automatic masking to specific stations and phase calibration generators. The vast majority of applicable code was written in FORTRAN to increase execution speed and ease of interfacing with the existing code base in PIMA, but a wrapper function in Python was also written to allow for a simpler method of interacting with GEPM and inputting relevant parameters. In addition to this report, thorough documentation was added to the already existing repositories associated with PIMA as a whole.

Geodesy

VGOS VLBI Intensives Between MACGO12M and WETTZ13S for the Rapid Determination of UT1-UTC

In this work, we present a status update and preliminary results of the designated research and development VLBI Intensive program VGOS-INT-S, observed between MACGO12M and WETTZ13S for the rapid determination of the Earth’s phase of rotation, expressed via UT1-UTC. Since 2021, 27 Intensive sessions have been observed successfully utilizing a special observation strategy alternating between high- and low-elevation scans for improved determination of delays caused by the neutral atmosphere. Between the end of January and mid of March 2022, VGOS-INT-S was among the most accurate Intensive programs. During this time, eight sessions were observed with an average formal error of 3.1 µs and a bias w.r.t. IERS C04 of 1.1 µs. Later, the session performance decreased due to multiple technical difficulties.

VLBI

A Concept of Precise VLBI/GNSS Ties with Micro-VLBI

We present here a concept of measuring local ties between collocated GNSS and VLBI stations using the microwave technique that effectively transforms a GNSS receiver to an element of a VLBI network. This is achieved by modifying the signal chain that allows to transfer voltage of the GNSS antenna to a digitizer via a coaxial cable. We discuss the application of this technique to local tie measurement. We have performed observations with a GNSS antenna and FD-VLBA radiotelescope and detected a strong interferometric signal from both radiogalaxies and GNSS satellites.

Geodesy

Implementing a VLBI Time Delay Model for Earth-orbiting Satellites: Partial derivatives and Verification

This document describes the partial derivatives for commonly estimated parameters including antenna positions, satellite position, and satellite velocity from a near-field VLBI delay model for Earth-orbiting satellites. This model was presented in the Journal of Geodesy by Jaron & Nothnagel (2019). In this context, we present a streamlined version of the near-field VLBI delay. From this simplified model, we deduce a delay rate expression and calculate partial derivatives, maintaining only the terms with significant impact on the computed derivatives’ magnitude. To verify the accuracy of the simplified model and the partial derivatives computed from it, we have created a simple simulation in Matlab of an Earth-orbiting satellite at the altitude of a typical Global Navigation Satellite Systems (GNSS) satellite. From this simulation, we compare the simplified and original delays, and we verify the magnitude and direction of the partial derivatives against the numerically computed derivatives from the original delay model. The partial derivatives and simplified VLBI delay model detailed here are implemented in Fortran in the open-source library VTD.

Geodesy

Single-Band VLBI Absolute Astrometry

The ionospheric path delay impacts single-band very long baseline interferometry (VLBI) group delays, which limits their applicability for absolute astrometry. I consider two important cases: when observations are made simultaneously in two bands, but delays in only one band are available for a subset of observations and when observations are made at one band design. I developed optimal procedures of data analysis for both cases using Global Navigation Satellite System (GNSS) ionosphere maps, provided a stochastic model that describes ionospheric errors, and evaluated their impact on source position estimates. I demonstrate that the stochastic model is accurate at a level of 15%. I found that using GNSS ionospheric maps as is introduces serious biases in estimates of declination and I developed a procedure that almost eliminates them. I found serendipitously that GNSS ionospheric maps have multiplicative errors and have to be scaled by 0.85 in order to mitigate the declination bias. A similar scale factor was found in comparison of the vertical total electron content from satellite altimetry against GNSS ionospheric maps. I favor interpretation of this scaling factor as a manifestation of the inadequacy of the thin shell model of the ionosphere. I showed that we are able to model the ionospheric path delay to the extent that no noticeable systematic errors emerge and we are able to assess adequately the contribution of the ionosphere-driven random errors on source positions. This makes single-band absolute astrometry a viable option that can be used for source position determination.

Radio astrometry

The Event Horizon Explorer Mission Concept

The Event Horizon Explorer (EHE) is a mission concept to extend the Event Horizon Telescope via anadditional space-based node. We provide highlights and overview of a concept study to explore the feasibility ofsuch a mission. We present science goals and objectives, which include studying the immediate environment aroundsupermassive black holes, and focus on critical enabling technologies and engineering challenges. We provide anassessment of their technological readiness and overall suitability for a NASA Medium Explorer (MIDEX) class mission.

Peter Kurczynski

Powerful quasars with young jets in multi-epoch radio surveys

Energetic feedback driven by the large-scale (100’s of kpc) lobes of classical radio galaxies is known to play an important role in shaping galaxy evolution. However, the prevalence of young and compact jets – and their impact on the interstellar medium – remains an open question. Multi-epoch radio surveys with cadences of years to decades offer a promising means of identifying even faint (mJy-level) jets that are compact and potentially young on the basis of variability. Recently, a comparison of images from the Very Large Array Sky Survey (VLASS) and the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey has revealed a population of distant (0.2 < z < 3.2) quasars that have brightened dramatically in the past 1–2 decades. These quasars appear to have transitioned from “radio-quiet” nondetections in FIRST to “radio-loud” detections in VLASS. Extensive multiband follow-up observations with the VLA from 1 to 18GHz have revealed compact (sub-kpc) radio sources that are consistent with young jets that were recently triggered. Here, we summarize the status of our on-going study of quasars with newborn jets identified in the radio time domain.

galaxies

The Use of Astronomy VLBA Campaign MOJAVE for Geodesy

We investigated the suitability of the astronomical 15 GHz VLBA observing program MOJAVE-5 for estimation of geodetic parameters, such as station coordinates and Earth orientation parameters. We processed contemporary geodetic dual-band RV and CN experiments observed at 2.3 GHz and 8.6 GHz starting on September 2016 through July 2020 as reference dataset. We showed that the baseline length repeatability from MOJAVE-5 experiments is only a factor of 1.5 greater than from the dedicated geodetic dataset and still below 1 ppb. The wrms of the difference of estimated EOP with respect to the reference IERS C04 time series are a factor of 1.3 to 1.8 worse. We isolated three major differences between the datasets in terms their possible impact on the geodetic results, i.e. the scheduling approach, treatment of the ionospheric delay, and selection of target radio sources. We showed that the major factor causing discrepancies in the estimated geodetic parameters is the different scheduling approach of the datasets. We conclude that systematic errors in MOJAVE-5 dataset are low enough for these data to be used as an excellent testbed for further investigations on the radio source structure effects in geodesy and astrometry.

Hana Krasna

Quasars That Have Transitioned from Radio-quiet to Radio-loud on Decadal Timescales Revealed by VLASS and FIRST

We have performed a search over 3440 deg^2 of Epoch 1(2017–2019)of the Very Large Array Sky Survey to identify unobscured quasars in the optical(0.2 2500%)but roughly steady fluxes over a few months at 3 GHz are inconsistent with extrinsic variability due to propagation effects, thus favoring an intrinsic origin. We conclude that our sources are powerful quasars hosting compact/young jets. This challenges the generally accepted idea that “radio-loudness” is a property of the quasar/AGN population that remains fixed on human timescales. Our study suggests that frequent episodes of short-lived AGN jets that do not necessarily grow to large scales may be common at high redshift. We speculate that intermittent but powerful jets on subgalactic scales could interact with the interstellar medium, possibly driving feedback capable of influencing galaxy evolution

Kristina Nyland

The Wide-field VLBA Calibrator Survey: WFCS

This paper presents the results of the largest very long baseline interferometry (VLBI) absolute astrometry campaign to date of 13,645 radio source observations with the Very Long Baseline Array. Of these, 7220 have been detected, including 6755 target sources that have never been observed with VLBI before. This makes the present VLBI catalog the largest ever published. The positions of the target sources have been determined with the median uncertainty of 1.7 mas, and 15,542 images of 7171 sources have been generated. Unlike previous absolute radio astrometry campaigns, observations were made at 4.3 and 7.6 GHz simultaneously using a single wide-band receiver. Because of the fine spectral and time resolutions, the field of view was 4'–8'—much greater than the 10''–20'' in previous surveys. This made possible the use of input catalogs with low position accuracy and the detection of a compact component in extended sources. Unlike previous absolute astrometry campaigns, both steep- and flat-spectrum sources were observed. The observations were scheduled in the so-called filler mode to fill the gaps between other high-priority programs. This was achieved by the development of the totally automatic scheduling procedure.

Leonid Petrov

Single-Band Absolute Astrometry

The ionospheric contribution affects path delay. When simultaneous dual-band observations are used for absolute astrometry, the residual contribution of the ionosphere is at a level of several picoseconds. But there are two cases when we need to process single band group delays: a) some observations provided usable data only for one band; b) an entire experiment used only one band. I have developed a novel approach to utilize GNSS TEC maps to get the best solutions for these two cases and provide realistic estimates of residual errors. Applications of this technology is discussed. In particular, the impact of residual ionospheric path delays at K-band on estimates of source positions is quantitatively evaluated and conclusions about advantages and disadvantages of high frequency absolute astrometry are made.

Leonid Petrov

Advances in Estimation of the Earth Orientation Parameters From Observation

Historically, a two-step procedure for estimation of Universal Time and polar motion from observations was adopted: night averages from individual observatories formed a raw non-equidistant time series that was edited for outliers, smoothed, and formed the final data product. This approach was by inertia used for processing modern space geodesy observations. I will present an alternative one-step approach and its application for processing a dataset of VLBI observations since 1990.0 through present. The essence of this approach is to represent the Earth orientation parameters in a form of an expansion into basic functions and estimate of the expansion coefficients directly from observations. A further extension of this approach is to assimilate not only space geodesy data, but atmospheric and ocean angular momentum from assimilation numerical weather models and their forecasts. That makes a smooth bridge between Earth rotation estimation and prediction. I will talk about my experience of running data analysis in such a way on a semi-operational basis, present results, and discuss advantages and disadvantages of this approach.

Leonid Petrov

From Formal Errors Towards Realistic Uncertainties

Evaluation of uncertainties of geodetic parameter estimates is the problem that is not yet solved in a satisfactory way. A direct evaluation of the uncertainties derived from a least square solution is labeled "formal" and is usually biased, sometimes up to an order of magnitude. Customary, the use of formal errors for scientific analysis is discouraged. We claim that the root of the problem is neglecting off-diagonal elements in the variance-covariance matrix of the noise in the data. A careful reconstruction of the full variance-covariance matrix, including the off-diagonal terms greatly improves realism of uncertainty estimates derived from least squares. We processed the dataset of VLBI group delays and built a prior variance-covariance of the atmosphere-driven noise based on analysis of the output of NASA high-resolution numerical weather models. We found that the uncertainties of parameter estimates derived from this least square solution that uses such variance-covariance matrices become much closer to realistic errors. We consider approaches for implementation of this method in routine data analysis of space geodesy data.

Leonid Petrov

Improvements in VLBI Data Analysis for Advanced in Geodesy

I will provide an overview of advances in VLBI data analysis implemented at NASA GSFC. That include the use of B-splines for modeling non-linear station position evolution, the use of all non-tidal mass loadings in data reduction, the use of slant a priori path delays from the output of numerical weather models, the use of B-splines for evaluation of the Earth orientation parameters as a continuous function, the use of elevation-dependent weighting, processing both 24 hr and 1 hr experiments in a single least square solutions and other innovations. I will show how this changes in data analysis affect space geodesy results, such as estimation of station positions, Earth orientation parameters, and source coordinates.

Leonid Petrov