Accurately accounting for spectral structure in spectrometer data induced by instrumental chromaticity on scales relevant for detection of the 21-cm signal is among the most significant challenges in global 21-cm signal analysis. In the publicly available EDGES low-band data set, this complicating structure is suppressed using beam-factor based chromaticity correction (BFCC), which works by dividing the data by a sky-map-weighted model of the spectral structure of the instrument beam. Several analyses of this data have employed models that start with the assumption that this correction is complete. However, while BFCC mitigates the impact of instrumental chromaticity on the data, given realistic assumptions regarding the spectral structure of the foregrounds, the correction is only partial. This complicates the interpretation of fits to the data with intrinsic sky models (models that assume no instrumental contribution to the spectral structure of the data). In this paper, we derive a BFCC data model from an analytic treatment of BFCC and demonstrate using simulated observations that, in contrast to using an intrinsic sky model for the data, the BFCC data model enables unbiased recovery of a simulated global 21-cm signal from beam-factor chromaticity corrected data in the limit that the data is corrected with an error-free beam-factor model.

This article focuses on the new generation of geodetic very long baseline interferometry (VLBI), the VLBI global observing system (VGOS), and measurements carried out during the CONT17 campaign. It uses broadband technology that increases both the number and precision of observations. These characteristics make VGOS a suitable tool for studying the atmosphere. This study focuses on the effects of the ionosphere on VGOS signals using a model that incorporates and extends ideas originally published in Hobiger et al. (2006, https://doi.org/10.1029/2005RS003297 ). Our investigation revealed that the differential total electron content (dTEC) data product calculated with the VGOS post‐processing software had a sign error that fortunately, does not change the final values of the phase and group delay. Therefore, this study was a way to identify this problem within the dTEC product. After diagnosing and solving this problem, the underlying model was modified such that instead of considering a single unknown for the latitude gradient of the ionosphere, a time series of latitude gradients were considered that enhanced the resulting vertical total electron content (VTEC) estimates. For evaluation purposes, time series of VTEC at each station during the CONT17 campaign were compared with VTEC obtained from the global navigation satellite system (GNSS). The final agreement between VGOS and GNSS was between 1.1 and 5.9 TEC units (TECU).

We measured the components of the 31-m-long vector between the two very-long-baseline interferometry (VLBI) antennas at the Kokee Park Geophysical Observatory (KPGO), Hawaii, with approximately 1 mm precision using phase delay observables from dedicated VLBI observations in 2016 and 2018. The two KPGO antennas are the 20 m legacy VLBI antenna and the 12 m VLBI Global Observing System (VGOS) antenna. Independent estimates of the vector between the two antennas were obtained by the National Geodetic Survey (NGS) using standard optical surveys in 2015 and 2018. The uncertainties of the latter survey were 0.3 and 0.7 mm in the horizontal and vertical components of the baseline, respectively. We applied corrections to the measured positions for the varying thermal deformation of the antennas on the different days of the VLBI and survey measurements, which can amount to 1 mm, bringing all results to a common reference temperature. The difference between the VLBI and survey results are 0.2 ± 0.4 mm, −1.3 ± 0.4 mm, and 0.8 ± 0.8 mm in the East, North, and Up topocentric components, respectively. We also estimate that the Up component of the baseline may suffer from systematic errors due to gravitational deformation and uncalibrated instrumental delay variations at the 20 m antenna that may reach ± 10 and −2 mm, respectively, resulting in an accuracy uncertainty on the order of 10 mm for the relative heights of the antennas. Furthermore, possible tilting of the 12 m antenna increases the uncertainties in the differences in the horizontal components to 1.0 mm. These results bring into focus the importance of (1) correcting to a common reference temperature the measurements of the reference points of all geodetic instruments within a site, (2) obtaining measurements of the gravitational deformation of all antennas, and (3) monitoring local motions of the geodetic instruments. These results have significant implications for the accuracy of global reference frames that require accurate local ties between geodetic instruments, such as the International Terrestrial Reference Frame (ITRF).

Earth and Space Science Open Archive This preprint has been submitted to and is under consideration at Radio Science. ESSOAr is a venue for early communication or feedback before peer review. Data may be preliminary.Learn more about preprints preprintOpen AccessYou are viewing the latest version by default [v1]Deriving ionospheric TEC by VGOS data analysis during CONT17AuthorsSanamMotlaghzadehiDM MahdiAlizadehRogerCappalloRobertHeinkelmannHaraldSchuhiDSee all authors Sanam MotlaghzadehiDCorresponding Author• Submitting AuthorK N Toosi University of TechnologyiDhttps://orcid.org/0000-0003-0895-2893view email addressThe email was not providedcopy email addressM Mahdi AlizadehK N Toosi University of TechnologyTechnical University of Berlinview email addressThe email was not providedcopy email addressRoger CappalloMIT Haystack Observatoryview email addressThe email was not providedcopy email addressRobert HeinkelmannGerman Research Centre for Geosciences (GFZ)view email addressThe email was not providedcopy email addressHarald SchuhiDTechnical University of BerlinGerman Research Centre for Geosciences (GFZ)iDhttps://orcid.org/0000-0001-5443-0370view email addressThe email was not providedcopy email address

3C 279 is an archetypal blazar with a prominent radio jet that show broadband flux density variability across the entire electromagnetic spectrum. We use an ultra-high angular resolution technique – global Very Long Baseline Interferometry (VLBI) at 1.3 mm (230 GHz) – to resolve the innermost jet of 3C 279 in order to study its fine-scale morphology close to the jet base where highly variableγ-ray emission is thought to originate, according to various models. The source was observed during four days in April 2017 with the Event Horizon Telescope at 230 GHz, including the phased Atacama Large Millimeter/submillimeter Array (ALMA), at an angular resolution of ∼20 μas (at a redshift ofz = 0.536 this corresponds to ∼0.13 pc ∼ 1700 Schwarzschild radii with a black hole massMBH = 8 × 108 M⊙). Imaging and model-fitting techniques were applied to the data to parameterize the fine-scale source structure and its variation. We find a multicomponent inner jet morphology with the northernmost component elongated perpendicular to the direction of the jet, as imaged at longer wavelengths. The elongated nuclear structure is consistent on all four observing days and across different imaging methods and model-fitting techniques, and therefore appears robust. Owing to its compactness and brightness, we associate the northern nuclear structure as the VLBI “core”. This morphology can be interpreted as either a broad resolved jet base or a spatially bent jet. We also find significant day-to-day variations in the closure phases, which appear most pronounced on the triangles with the longest baselines. Our analysis shows that this variation is related to a systematic change of the source structure. Two inner jet components move non-radially at apparent speeds of ∼15 cand ∼20 c(∼1.3 and ∼1.7 μas day−1, respectively), which more strongly supports the scenario of traveling shocks or instabilities in a bent, possibly rotating jet. The observed apparent speeds are also coincident with the 3C 279 large-scale jet kinematics observed at longer (cm) wavelengths, suggesting no significant jet acceleration between the 1.3 mm core and the outer jet. The intrinsic brightness temperature of the jet components are ≲1010K, a magnitude or more lower than typical values seen at ≥7 mm wavelengths. The low brightness temperature and morphological complexity suggest that the core region of 3C 279 becomes optically thin at short (mm) wavelengths.

We present the design and implementation of an automated data calibration and reduction pipeline for very long baseline interferometric (VLBI) observations taken at millimeter wavelengths. These short radio wavelengths provide the best imaging resolution available from ground-based VLBI networks such as the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA) but require specialized processing owing to the strong effects from atmospheric opacity and turbulence, as well as the heterogeneous nature of existing global arrays. The pipeline builds on a calibration suite (HOPS) originally designed for precision geodetic VLBI. To support the reduction of data for astronomical observations, we have developed an additional framework for global phase and amplitude calibration that provides output in a standard data format for astronomical imaging and analysis. The pipeline was successfully used toward the reduction of 1.3 mm observations from the EHT 2017 campaign, leading to the first image of a black hole “shadow” at the center of the radio galaxy M87. In this work, we analyze observations taken at 3.5 mm (86 GHz) by the GMVA, joined by the phased Atacama Large Millimeter/submillimeter Array in 2017 April, and demonstrate the benefits from the specialized processing of high-frequency VLBI data with respect to classical analysis techniques.

We present the calibration and reduction of Event Horizon Telescope (EHT) 1.3 mm radio wavelength observations of the supermassive black hole candidate at the center of the radio galaxy M87 and the quasar 3C 279, taken during the 2017 April 5-11 observing campaign. These global very long baseline interferometric observations include for the first time the highly sensitive Atacama Large Millimeter/submillimeter Array (ALMA); reaching an angular resolution of 25 mu as, with characteristic sensitivity limits of similar to 1 mJy on baselines to ALMA and similar to 10 mJy on other baselines. The observations present challenges for existing data processing tools, arising from the rapid atmospheric phase fluctuations, wide recording bandwidth, and highly heterogeneous array. In response, we developed three independent pipelines for phase calibration and fringe detection, each tailored to the specific needs of the EHT. The final data products include calibrated total intensity amplitude and phase information. They are validated through a series of quality assurance tests that show consistency across pipelines and set limits on baseline systematic errors of 2% in amplitude and 1 degrees in phase. The M87 data reveal the presence of two nulls in correlated flux density at similar to 3.4 and similar to 8.3 G lambda and temporal evolution in closure quantities, indicating intrinsic variability of compact structure on a. timescale of days, or several light-crossing times for a. few billion solar-mass black hole. These measurements provide the first opportunity to image horizon-scale structure in M87.

We report on a search for low-frequency radio variability in 944 bright (>4 Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources, we find very low levels of variability with typical modulation indices <5 per cent. We detect 15 candidate low-frequency variables that show significant long-term variability (>2.8 yr) with time-averaged modulation indices (M) over bar = 3.1-7.1 per cent. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7 per cent of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milliarcsec (which we support with multiwavelength data). At 154 MHz, we demonstrate that interstellar scintillation time-scales become long (similar to decades) and have low modulation indices, while synchrotron-driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work, we infer that the low-frequency extragalactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than 1 yr.

The Event Horizon Telescope (EHT) is a very long baseline interferometry (VLBI) array that comprises millimeter- and submillimeter-wavelength telescopes separated by distances comparable to the diameter of the Earth. At a nominal operating wavelength of 1.3 mm, EHT angular resolution (lambda/D) is 25 micro-as, which is sufficient to resolve nearby supermassive black hole candidates on spatial and temporal scales that correspond to their event horizons. With this capability, the EHT scientific goals are to probe general relativistic effects in the strong-field regime and to study accretion and relativistic jet formation near the black hole boundary. In this Letter we describe the system design of the EHT, detail the technology and instrumentation that enable observations, and provide measures of its performance. Meeting the EHT science objectives has required several key developments that have facilitated the robust extension of the VLBI technique to EHT observing wavelengths and the production of instrumentation that can be deployed on a heterogeneous array of existing telescopes and facilities. To meet sensitivity requirements, high-bandwidth digital systems were developed that process data at rates of 64 gigabit/s, exceeding those of currently operating cm-wavelength VLBI arrays by more than an order of magnitude. Associated improvements include the development of phasing systems at array facilities, new receiver installation at several sites, and the deployment of hydrogen maser frequency standards to ensure coherent data capture across the array. These efforts led to the coordination and execution of the first Global EHT observations in 2017 April, and to event-horizon-scale imaging of the supermassive black hole candidate in M87.

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional uv coverage in the N–S direction, and leads to a spatial resolution of ∼30 μ as (∼3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of ∼4%–13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of ∼3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

Type III solar radio bursts are the Sun’s most intense and frequent nonthermal radio emissions. They involve two critical problems in astrophysics, plasma physics, and space physics: how collective processes produce nonthermal radiation and how magnetic reconnection occurs and changes magnetic energy into kinetic energy. Here magnetic reconnection events are identified definitively in Solar Dynamics Observatory UV-EUV data, with strong upward and downward pairs of jets, current sheets, and cusp-like geometries on top of time-varying magnetic loops, and strong outflows along pairs of open magnetic field lines. Type III bursts imaged by the Murchison Widefield Array and detected by the Learmonth radiospectrograph and STEREO B spacecraft are demonstrated to be in very good temporal and spatial coincidence with specific reconnection events and with bursts of X-rays detected by the RHESSI spacecraft. The reconnection sites are low, near heights of 5–10 Mm. These images and event timings provide the long-desired direct evidence that semi-relativistic electrons energized in magnetic reconnection regions produce type III radio bursts. Not all the observed reconnection events produce X-ray events or coronal or interplanetary type III bursts; thus different special conditions exist for electrons leaving reconnection regions to produce observable radio, EUV, UV, and X-ray bursts.

We present the first low-frequency Murchison Widefield Array (MWA) radio continuum maps of the Magellanic Clouds (MCs), using mosaics from the GaLactic Extragalactic All-Sky MWA (GLEAM) survey. In this paper, we discuss the overall radio continuum morphology between 76 and 227 MHz and compare them with neutral hydrogen maps, 1.4 GHz continuum maps and optical images. Variation of diffuse emission is noticeable across the Large Magellanic Cloud (LMC) but absent across the bar of the Small Magellanic Cloud (SMC). We also measure the integrated flux densities and derive the spectral indices for the MCs. A double power-law model with fixed alpha(1) = -0.1 fit between 19.7 MHz and 8.55 GHz yields alpha(0) = -0.66 +/- 0.08 for the LMC. Apower-law model yields alpha(85.5MHz)(8.55GHz) = -0.82 +/- 0.03 for the SMC. The radio spectral index maps reveal distinctive flat and steep spectral indices for the H IT regions and supernova remnants, respectively. We find strong correlation between H II regions and H, emission. Using a new 150 MHz-H-alpha relation as a star formation rate indicator, we estimate global star formation rates of 0.068-0.161 M-circle dot yr(-1) and 0.021-0.050 M-circle dot yr(-1) for the LMC and SMC, respectively. Images in 20 frequency bands, and wideband averages are made available via the GLEAM virtual observatory server.

A prototype broadband geodetic very long baseline interferometry system has been implemented, and measurements of the baseline length over approximately two years, between December 2014 and January 2017, have been made in the process of exercising the system, developing operational procedures, and assessing geodetic precision for the new broadband observing concept. In addition to developing a broadband signal chain and installing the instrumentation on both a new 12‐m antenna at the Goddard Geophysical and Astrophysical Observatory and the 18‐m Westford antenna at the Massachusetts Institute of Technology Haystack Observatory, it was necessary to develop new correlation and analysis procedures to process the four‐band, dual‐linear‐polarization data. A geodetic analysis of the data from 19 sessions that were observed during this period yielded a weighted root‐mean‐square deviation of the baseline length residuals about the weighted mean of 1.6 mm. These results validate several of the expectations set forth for the vision of the next‐generation geodetic very long baseline interferometry system.

Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilizing the first year of observations. The catalogue covers 24 831 square degrees, over declinations south of + 30 degrees and Galactic latitudes outside 10 degrees of the Galactic plane, excluding some areas such as theMagellanic Clouds. It contains 307 455 radio sources with 20 separate flux density measurements across 72-231 MHz, selected from a time-and frequency-integrated image centred at 200 MHz, with a resolution of approximate to 2 arcmin. Over the catalogued region, we estimate that the catalogue is 90 per cent complete at 170 mJy, and 50 per cent complete at 55 mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97 per cent above the detection threshold of 5 sigma, which itself is typically 50 mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low-frequency flux density scale of the southern sky to better than 10 per cent. This paper presents details of the flagging, imaging, mosaicking and source extraction/characterization, as well as estimates of the completeness and reliability. All source measurements and images are available online. (1) This is the first in a series of publications describing the GLEAM survey results.

Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak short-lived narrowband emission features, even during moderately quiet solar conditions. These nonthermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans, and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of −2.23 in the 12–155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1–2 s and possess bandwidths of about 4–5 MHz. Their occurrence rate remains fairly flat in the 140–210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts.

We present new radio continuum observations of NGC 253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC 253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of a central starburst and extended emission. The central component, corresponding to the inner 500 pc of the starburst region of the galaxy, is best modeled as an internally free–free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the spectrum of NGC 253 is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the southeast halo, and may be indicative of synchrotron self-absorption of shock-reaccelerated electrons or an intrinsic low-energy cutoff of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC 253 in our radio images. At 154–231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ∼8 kpc in the z -direction (from the major axis).

Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (< 300 MHz) remains relatively unexplored. Blind surveys with new wide-field radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble of faint, unresolved sources, one can in principle detect transients below the classical confusion limit to the extent that the classical confusion noise is independent of time. We develop a technique for detecting radio transients that is based on temporal matched filters applied directly to time series of images, rather than relying on source-finding algorithms applied to individual images. This technique has well-defined statistical properties and is applicable to variable and transient searches for both confusion-limited and non-confusion-limited instruments. Using the Murchison Widefield Array as an example, we demonstrate that the technique works well on real data despite the presence of classical confusion noise, sidelobe confusion noise, and other systematic errors. We searched for transients lasting between 2 minutes and 3 months. We found no transients and set improved upper limits on the transient surface density at 182 MHz for flux densities between similar to 20 and 200 mJy, providing the best limits to date for hour-and month-long transients.