Search

We present an analytical model for the evolution of extended active galactic nuclei (AGNs) throughout their full lifecycle, including the initial jet expansion, lobe formation, and eventual remnant phases. A particular focus of our contribution is on the early jet expansion phase, which is traditionally not well captured in analytical models. We implement this model within the Radio AGN in Semi-Analytic Environments (RAiSE) framework, and find that the predicted radio source dynamics are in good agreement with hydrodynamic simulations of both low-powered Fanaroff-Riley Type-I and high-powered Type-II radio lobes. We construct synthetic synchrotron surface brightness images by complementing the original RAiSE model with the magnetic field and shock-acceleration histories of a set of Lagrangian tracer particles taken from an existing hydrodynamic simulation. We show that a single set of particles is sufficient for an accurate description of the dynamics and observable features of Fanaroff-Riley Type-II radio lobes with very different jet parameters and ambient density profile normalisations. Our new model predicts that the lobes of young (< 10 Myr) sources will be both longer and brighter than expected at the same age from existing analytical models which lack a jet-dominated expansion phase; this finding has important implications for interpretation of radio galaxy observations. The RAiSE code, written in Python, is publicly available on GitHub and PyPI., 22 pages, 11 figures, 1 table; accepted in MNRAS

Quantifying the energetics and lifetimes of remnant radio-loud active galactic nuclei (AGNs) is much more challenging than for active sources due to the added complexity of accurately determining the time since the central black hole switched off. Independent spectral modelling of remnant lobes enables the derivation of the remnant ratio, $R_\mathrm{rem}$, (i.e. `off-time/source age'), thus reducing this complexity back to that of active sources however, the requirement of high-frequency ($\gtrsim5\,$GHz) coverage makes the application of this technique over large-area radio surveys difficult. In this work we propose a new method, which relies on the observed brightness of backflow of Fanaroff-Riley type~II lobes, combined with the \emph{Radio AGN in Semi-Analytic Environments} (RAiSE) code, to measure the duration of the remnant phase. Sensitive radio observations of the remnant radio galaxy J2253-34 are obtained to provide a robust comparison of this technique with the canonical spectral analysis and modelling methods. We find that the remnant lifetimes modelled by each method are consistent; spectral modelling yields $R_\mathrm{rem} = 0.23\pm0.02$, compared to $R_\mathrm{rem} = 0.26\pm0.02$ from our new method. We examine the viability of applying our proposed technique to low-frequency radio surveys using mock radio source populations, and examine whether the technique is sensitive to any intrinsic properties of radio AGNs. Our results show that the technique can be used to robustly classify active and remnant populations, with the most confident predictions for the remnant ratio, and thus off-time, in the longest-lived radio sources ($>50~$Myr) and those at higher redshifts ($z > 0.1$)., 20 pages, 9 figures, 8 tables. Accepted for publication to MNRAS

We present a method for applying spatially resolved adiabatic and radiative loss processes to synthetic radio emission from hydrodynamic simulations of radio sources from active galactic nuclei (AGN). Lagrangian tracer particles, each representing an ensemble of electrons, are injected into simulations and the position, grid pressure, and time since the last strong shock are recorded. These quantities are used to track the losses of the electron packet through the radio source in a manner similar to the Radio AGN in Semi-analytic Environments (RAiSE) formalism, which uses global source properties to calculate the emissivity of each particle ex-situ. Freedom in the choice of observing parameters, including redshift, is provided through the post-processing nature of this approach. We apply this framework to simulations of jets in different environments, including asymmetric ones. We find a strong dependence of radio source properties on frequency and redshift, in good agreement with observations and previous modelling work. There is a strong evolution of radio spectra with redshift due to the more prominent inverse-Compton losses at high redshift. Radio sources in denser environments have flatter spectral indices, suggesting that spectral index asymmetry may be a useful environment tracer. We simulate intermediate Mach number jets that disrupt before reaching the tip of the lobe, and find that these retain an edge-brightened Fanaroff-Riley Type II morphology, with the most prominent emission remaining near the tip of the lobes for all environments and redshifts we study., 17 pages, 15 figures. Accepted for publication in MNRAS

We present the discovery of another Odd Radio Circle (ORC) with the Australian Square Kilometre Array Pathfinder (ASKAP) at 944 MHz. The observed radio ring, ORC J0102-2450, has a diameter of ~70 arcsec or 300 kpc, if associated with the central elliptical galaxy DES J010224.33-245039.5 (z ~ 0.27). Considering the overall radio morphology (circular ring and core) and lack of ring emission at non-radio wavelengths, we investigate if ORC J0102-2450 could be the relic lobe of a giant radio galaxy seen end-on or the result of a giant blast wave. We also explore possible interaction scenarios, for example, with the companion galaxy, DES J010226.15-245104.9, located in or projected onto the south-eastern part of the ring. We encourage the search for further ORCs in radio surveys to study their properties and origin., MNRAS Letters, accepted (5 pages, 3 figures)

In this paper, we present the identification of five previously unknown giant radio galaxies (GRGs) using Data Release 1 of the Radio Galaxy Zoo citizen science project and a selection method appropriate to the training and validation of deep learning algorithms for new radio surveys. We associate one of these new GRGs with the brightest cluster galaxy (BCG) in the galaxy cluster GMBCG J251.67741+36.45295 and use literature data to identify a further 13 previously known GRGs as BCG candidates, increasing the number of known BCG GRGs by $\gt 60$ per cent. By examining local galaxy number densities for the number of all known BCG GRGs, we suggest that the existence of this growing number implies that GRGs are able to reside in the centres of rich (∼1014 M⊙) galaxy clusters and challenges the hypothesis that GRGs grow to such sizes only in locally underdense environments.

Radio galaxies are known to go through cycles of activity, where phases of apparent quiescence can be followed by repeated activity of the central supermassive black hole. A better understanding of this cycle is crucial for ascertaining the energetic impact that the jets have on the host galaxy, but little is known about it. We used deep LOFAR images at 150 MHz of the Lockman Hole extragalactic field to select a sample of 158 radio sources with sizes $> 60^{\prime\prime}$ in different phases of their jet life cycle. Using a variety of criteria (e.g. core prominence combined with low-surface brightness of the extended emission and steep spectrum of the central region) we selected a subsample of candidate restarted radio galaxies representing between 13% and 15% of the 158 sources of the main sample. We compare their properties to the rest of the sample, which consists of remnant candidates and active radio galaxies. Optical identifications and characterisations of the host galaxies indicate similar properties for candidate restarted, remnant, and active radio galaxies, suggesting that they all come from the same parent population. The fraction of restarted radio galaxies is slightly higher with respect to remnants, suggesting that the restarted phase can often follow after a relatively short remnant phase (the duration of the remnant phase being a few times 10$^{7}$ years). This confirms that the remnant and restarted phases are integral parts of the life cycle of massive elliptical galaxies. A preliminary investigation does not suggest a strong dependence of this cycle on the environment surrounding any given galaxy., 16 pages, 9 figures, 4 tables, A&A accepted, companion paper to Shabala et al., also in today arXiv posting

This paper demonstrates a novel and efficient unsupervised clustering method with the combination of a Self-Organising Map (SOM) and a convolutional autoencoder. The rapidly increasing volume of radio-astronomical data has increased demand for machine learning methods as solutions to classification and outlier detection. Major astronomical discoveries are unplanned and found in the unexpected, making unsupervised machine learning highly desirable by operating without assumptions and labelled training data. Our approach shows SOM training time is drastically reduced and high-level features can be clustered by training on auto-encoded feature vectors instead of raw images. Our results demonstrate this method is capable of accurately separating outliers on a SOM with neighbourhood similarity and K-means clustering of radio-astronomical features complexity. We present this method as a powerful new approach to data exploration by providing a detailed understanding of the morphology and relationships of Radio Galaxy Zoo (RGZ) dataset image features which can be applied to new radio survey data., Accepted in Publications of the Astronomical Society of the Pacific, special issue on Machine Intelligence in Astronomy and Astrophysics. 23 pages, 8 full-page colour figures

Gaseous halos play a key role for understanding inflow, feedback and the overall baryon budget in galaxies. Literature models predict transitions of the state of the gaseous halo between cold and hot accretion, winds, fountains and hydrostatic halos at certain galaxy masses. Since luminosities of radio AGN are sensitive to halo densities, any significant transition would be expected to show up in the radio luminosities of large samples of galaxies. The Low Frequency Array (LOFAR) Two Metre Sky Survey (LoTSS) has indeed identified a galaxy stellar mass scale, $10^{11} M_\odot$ , above which the radio luminosities increase disproportionately. Here, we investigate, if radio luminosities of galaxies, especially the marked rise at galaxy masses around $10^{11} M_\odot$, can be explained with standard assumptions on jet powers, scaling between black hole-mass and galaxy mass and gaseous halos. We developed models for the radio luminosity of radio AGN in halos under infall, galactic wind and hydrostatic conditions based on observational data and theoretical constraints, and compared it to LoTSS data for a large sample of galaxies in the mass rangebetween $10^{8.5} M_\odot$ and $10^{12} M_\odot$. Assuming the same characteristic upper limit to jet powers as is known from high galaxy masses to hold at all masses, we find that the maximum radio luminosities for the hydrostatic gas halos fit well with the upper envelope of the distribution of the LOFAR data. The marked rise in radio luminosity at $10^{11} M_\odot$ is matched in our model, and is related to significant change in halo gas density around this galaxy mass, which is a consequence of the lower cooling rates at higher virial temperature. Wind and infall models overpredict the radio luminosities at small galaxy masses and have no particular steepening of the run of the radio luminosities predicted at any galaxy mass. [...], 6 pages, 3 figures, accepted by A&A, comments welcome

We describe a physical model of the outflows produced as a result of gas accretion onto a black hole, and the resultant changes to star formation rates and efficiencies in galaxies, using the Radio-SAGE semi-analytic galaxy formation model. We show that the ratio of outflow rate to SFR of galaxies is mainly driven by black hole mass and virial halo mass, and show that the SFR is higher than the outflow rate at low black hole masses. The model consistently reproduces the observed evolution of star formation rate density from z = 6 to z = 0, as well as the trend of the stellar mass - halo mass relations. We show the characteristic growth of massive galaxies influenced by AGN feedback at different redshifts. We find feedback to be prevalent in the most massive galaxy halos, inhibiting the cooling catastrophe., 15 pages, 8 figures. Accepted for publication in MNRAS

We investigate the role of environment on radio galaxy properties by constructing a sample of large ($\gtrsim100$~kpc), nearby ($z, 18 pages, 31 figures, 1 table; accepted in MNRAS

The upcoming next-generation large area radio continuum surveys can expect tens of millions of radio sources, rendering the traditional method for radio morphology classification through visual inspection unfeasible. We present ClaRAN - Classifying Radio sources Automatically with Neural networks - a proof-of-concept radio source morphology classifier based upon the Faster Region-based Convolutional Neutral Networks (Faster R-CNN) method. Specifically, we train and test ClaRAN on the FIRST and WISE images from the Radio Galaxy Zoo Data Release 1 catalogue. ClaRAN provides end users with automated identification of radio source morphology classifications from a simple input of a radio image and a counterpart infrared image of the same region. ClaRAN is the first open-source, end-to-end radio source morphology classifier that is capable of locating and associating discrete and extended components of radio sources in a fast (< 200 milliseconds per image) and accurate (>= 90 %) fashion. Future work will improve ClaRAN's relatively lower success rates in dealing with multi-source fields and will enable ClaRAN to identify sources on much larger fields without loss in classification accuracy., 22 pages, 16 figures, Accepted in Monthly Notices of the Royal Astronomical Society

We consider the problem of determining the host galaxies of radio sources by cross-identification. This has traditionally been done manually, which will be intractable for wide-area radio surveys like the Evolutionary Map of the Universe (EMU). Automated cross-identification will be critical for these future surveys, and machine learning may provide the tools to develop such methods. We apply a standard approach from computer vision to cross-identification, introducing one possible way of automating this problem, and explore the pros and cons of this approach. We apply our method to the 1.4 GHz Australian Telescope Large Area Survey (ATLAS) observations of the Chandra Deep Field South (CDFS) and the ESO Large Area ISO Survey South 1 (ELAIS-S1) fields by cross-identifying them with the Spitzer Wide-area Infrared Extragalactic (SWIRE) survey. We train our method with two sets of data: expert cross-identifications of CDFS from the initial ATLAS data release and crowdsourced cross-identifications of CDFS from Radio Galaxy Zoo. We found that a simple strategy of cross-identifying a radio component with the nearest galaxy performs comparably to our more complex methods, though our estimated best-case performance is near 100 per cent. ATLAS contains 87 complex radio sources that have been cross-identified by experts, so there are not enough complex examples to learn how to cross-identify them accurately. Much larger datasets are therefore required for training methods like ours. We also show that training our method on Radio Galaxy Zoo cross-identifications gives comparable results to training on expert cross-identifications, demonstrating the value of crowdsourced training data.

Kinetic jet power estimates based exclusively on observed monochromatic radio luminosities are highly uncertain due to confounding variables and a lack of knowledge about some aspects of the physics of active galactic nuclei (AGNs). We propose a new methodology to calculate the jet powers of the largest, most powerful radio sources based on combinations of their size, lobe luminosity and shape of their radio spectrum; this approach avoids the uncertainties encountered by previous relationships. The outputs of our model are calibrated using hydrodynamical simulations and tested against independent X-ray inverse-Compton measurements. The jet powers and lobe magnetic field strengths of radio sources are found to be recovered using solely the lobe luminosity and spectral curvature, enabling the intrinsic properties of unresolved high-redshift sources to be inferred. By contrast, the radio source ages cannot be estimated without knowledge of the lobe volumes. The monochromatic lobe luminosity alone is incapable of accurately estimating the jet power or source age without knowledge of the lobe magnetic field strength and size respectively. We find that, on average, the lobes of the 3C radio sources have magnetic field strengths approximately a factor three lower than the equipartition value, inconsistent with equal energy in the particles and the fields at the 5$\sigma$ level. The particle content of 3C radio lobes is discussed in the context of complementary observations; we do not find evidence favouring an energetically-dominant proton population., Comment: 21 pages, 16 figures, 3 tables; accepted in MNRAS

We study the mutual alignment of radio sources within two surveys, Faint Images of the Radio Sky at Twenty-centimetres (FIRST) and TIFR GMRT Sky Survey (TGSS). This is done by producing two position angle catalogues containing the preferential directions of respectively 30 059 and 11 674 extended sources distributed over more than 7000 and 17 000 deg2. The identification of the sources in the FIRST sample was performed in advance by volunteers of the Radio Galaxy Zoo (RGZ) project, while for the TGSS sample it is the result of an automated process presented here. After taking into account systematic effects, marginal evidence of a local alignment on scales smaller than 2.5 deg is found in the FIRST sample. The probability of this happening by chance is found to be less than 2 per cent. Further study suggests that on scales up to 1.5 deg the alignment is maximal. For one third of the sources, the RGZ volunteers identified an optical counterpart. Assuming a flat Λ cold dark matter cosmology with Ωm=0.31,ΩΛ=0.69 Ωm=0.31,ΩΛ=0.69 , we convert the maximum angular scale on which alignment is seen into a physical scale in the range [19, 38] Mpc h−170 h70−1 . This result supports recent evidence reported by Taylor and Jagannathan of radio jet alignment in the 1.4 deg2 ELAIS N1 field observed with the Giant Metrewave Radio Telescope. The TGSS sample is found to be too sparsely populated to manifest a similar signal.

In this work we present IRAM-30m telescope observations of a sample of bulge-dominated galaxies with large dust lanes, which have had a recent minor merger. We find these galaxies are very gas rich, with H2 masses between 4x10^8 and 2x10^10 Msun. We use these molecular gas masses, combined with atomic gas masses from an accompanying paper, to calculate gas-to-dust and gas-to-stellar mass ratios. The gas-to-dust ratios of our sample objects vary widely (between ~50 and 750), suggesting many objects have low gas-phase metallicities, and thus that the gas has been accreted through a recent merger with a lower mass companion. We calculate the implied minor companion masses and gas fractions, finding a median predicted stellar mass ratio of ~40:1. The minor companion likely had masses between ~10^7 - 10^10 Msun. The implied merger mass ratios are consistent with the expectation for low redshift gas-rich mergers from simulations. We then go on to present evidence that (no matter which star-formation rate indicator is used) our sample objects have very low star-formation efficiencies (star-formation rate per unit gas mass), lower even than the early-type galaxies from ATLAS3D which already show a suppression. This suggests that minor mergers can actually suppress star-formation activity. We discuss mechanisms that could cause such a suppression, include dynamical effects induced by the minor merger., Comment: 12 pages, 5 figures, accepted to MNRAS

Radio AGN in the nearby Universe are more likely to be found in galaxies with early-type morphology, the detection rate in spiral or late-type galaxies (LTGs) being around an order of magnitude lower. We combine the mJy Imaging VLBA Exploration at 20cm (mJIVE-20) survey with the Sloan Digital Sky Survey (SDSS), to study the relatively rare population of AGN in LTGs that have nuclear radio luminosities similar to that in their early-type counterparts. The LTG AGN population is preferentially hosted by galaxies that have high stellar masses (M* > 10^10.8 MSun), red colours and low star-formation rates, with little dependence on the detailed morphology or local environment of the host LTG. The merger fraction in the LTG AGN is around 4 times higher than that in the general LTG population, indicating that merging is an important trigger for radio AGN in these systems. The red colours of our systems extend recent work which indicates that merger-triggered AGN in the nearby Universe appear after the peak of the associated starburst, implying that they do not strongly regulate star formation. Finally, we find that in systems where parsec-scale jets are clearly observed in our VLBI images, the jets are perpendicular to the major axis of the galaxy, indicating alignment between the accretion disc and the host galaxy stellar disc., Submitted to MNRAS

We explore the processes that trigger local AGN and the role of these AGN in regulating star formation, using ~350 nearby galaxies observed by the mJy Imaging VLBA Exploration at 20cm (mJIVE) survey. The >10^7 K brightness temperature required for an mJIVE detection cannot be achieved via star formation alone, allowing us to unambiguously detect nearby radio AGN and study their role in galaxy evolution. Radio AGN are an order of magnitude more common in early-type galaxies (ETGs) than in their late-type counterparts. The VLBI-detected ETGs in this study have a similar stellar mass distribution to their undetected counterparts, are typically not the central galaxies of clusters and exhibit merger fractions that are significantly higher than in the average ETG. This suggests that these radio AGN (which have VLBI luminosities >10^22 W Hz^-1) are primarily fuelled by mergers, and not by internal stellar mass loss or cooling flows. Our radio AGN are a factor of ~3 times more likely to reside in the UV-optical red sequence than the average ETG. Furthermore, typical AGN lifetimes (a few 10^7 yr) are much shorter than the transit times from blue cloud to red sequence (~1.5 Gyr). This indicates that the AGN are not triggered promptly and appear several dynamical timescales into the associated star formation episode, implying that they typically couple only to residual gas, at a point where star formation has already declined significantly. While evidence for AGN feedback is strong in systems where the black hole is fed by the cooling of hot gas, AGN triggered by mergers appear not to strongly regulate the associated star formation. The inability of the AGN to rapidly quench merger-driven star formation is likely to make merging the dominant mode of star formation in nearby ETGs, in line with evidence for minor mergers being the primary driver of stellar mass growth in these systems., Submitted to MNRAS

We use multi-wavelength (0.12 - 500 micron) photometry from Herschel-ATLAS, WISE, UKIDSS, SDSS and GALEX, to study 23 nearby spheroidal galaxies with prominent dust lanes (DLSGs). DLSGs are considered to be remnants of recent minor mergers, making them ideal laboratories for studying both the interstellar medium (ISM) of spheroids and minor-merger-driven star formation in the nearby Universe. The DLSGs exhibit star formation rates (SFRs) between 0.01 and 10 MSun yr^-1, with a median of 0.26 MSun yr^-1 (a factor of 3.5 greater than the average SG). The median dust mass, dust-to-stellar mass ratio and dust temperature in these galaxies are around 10^7.6 MSun yr^-1, ~0.05% and ~19.5 K respectively. The dust masses are at least a factor of 50 greater than that expected from stellar mass loss and, like the SFRs, show no correlation with galaxy luminosity, suggesting that both the ISM and the star formation have external drivers. Adopting literature gas-to-dust ratios and star formation histories derived from fits to the panchromatic photometry, we estimate that the median current and initial gas-to-stellar mass ratios in these systems are ~4% and ~7% respectively. If, as indicated by recent work, minor mergers that drive star formation in spheroids with (NUV-r)>3.8 (the colour range of our DLSGs) have stellar mass ratios between 1:6 and 1:10, then the satellite gas fractions are likely >50%., MNRAS in press

We examine the relationship between source radio luminosity and kinetic power in Active Galactic Nucleus (AGN) jets. We show that neglecting various loss processes can introduce a systematic bias in the jet powers inferred from radio luminosities for a sample of radio galaxies. This bias can be corrected for by considering source size as well as radio luminosity; effectively the source size acts as a proxy for source age. Based on a sample of FR-II radio sources with jet powers derived from the measured hotspot parameters, we empirically determine a new expression for jet power that accounts for the source size, Q_jet / 10^{36} W = 1.5 (L_151 / 10^{27} W/Hz)^{0.8} (1+z)^{1.0} (D / kpc)^{0.58 \pm 0.17}, where D is source size and L_151 the 151 MHz radio luminosity. By comparing a flux-limited and volume-limited sample, we show that any derived radio luminosity - jet power relation depends sensitively on sample properties, in particular the source size distribution and the size-luminosity correlation inherent in the sample. Such bias will affect the accuracy of the kinetic luminosity function derived from lobe radio luminosities and should be treated with caution., 8 pages, 5 figures, accepted for publication in ApJ

The AuScope geodetic Very Long Baseline Interferometry array consists of three new 12 m radio telescopes and a correlation facility in Australia. The telescopes at Hobart (Tasmania), Katherine (Northern Territory) and Yarragadee (Western Australia) are co-located with other space geodetic techniques including Global Navigation Satellite Systems (GNSS) and gravity infrastructure, and in the case of Yarragadee, Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) facilities. The correlation facility is based in Perth (Western Australia). This new facility will make significant contributions to improving the densification of the International Celestial Reference Frame in the Southern Hemisphere, and subsequently enhance the International Terrestrial Reference Frame through the ability to detect and mitigate systematic error. This, combined with the simultaneous densification of the GNSS network across Australia will enable the improved measurement of intraplate deformation across the Australian tectonic plate. In this paper we present a description of this new infrastructure and present some initial results, including telescope performance measurements and positions of the telescopes in the International Terrestrial Reference Frame. We show that this array is already capable of achieving centimetre precision over typical long-baselines and that network and reference source systematic effects must be further improved to reach the ambitious goals of VLBI2010., Journal of Geodesy, March 2013

We report on the variability of 443 flat spectrum, compact radio sources monitored using the VLA for 3 days in 4 epochs at ~ 4 month intervals at 5 GHz as part of the Micro-Arcsecond Scintillation-Induced Variability (MASIV) survey. Over half of these sources exhibited 2-10% rms variations on timescales over 2 days. We analyzed the variations by two independent methods, and find that the rms variability amplitudes of the sources correlate with the emission measure in the ionized Interstellar Medium along their respective lines of sight. We thus link the variations with interstellar scintillation of components of these sources, with some (unknown) fraction of the total flux density contained within a compact region of angular diameter in the range 10-50 micro-arcseconds. We also find that the variations decrease for high mean flux density sources and, most importantly, for high redshift sources. The decrease in variability is probably due either to an increase in the apparent diameter of the source, or a decrease in the flux density of the compact fraction beyond z ~ 2. Here we present a statistical analysis of these results, and a future paper will the discuss the cosmological implications in detail., Comment: 62 pages, 13 figures. Accepted for publication in the Astrophysical Journal

We present observations of the H91$\alpha$ recombination line emission towards a sample of nine HII regions associated with 6.7-GHz methanol masers, and report arcsecond-scale emission around compact cores. We derive physical parameters for our sources, and find that although simple hydrostatic models of region evolution reproduce the observed region sizes, they significantly underestimate emission measures. We argue that these findings are consistent with young source ages in our sample, and can be explained by existence of density gradients in the ionised gas., Comment: 11 pages, 6 figures; accepted for publication in MNRAS

The AUSTRAL observing program was started in 2011, performing geodetic and astrometric very long baseline interferometry (VLBI) sessions using the new Australian AuScope VLBI antennas at Hobart, Katherine, and Yarragadee, with contribution from the Warkworth (New Zealand) 12 m and Hartebeesthoek (South Africa) 15 m antennas to make a southern hemisphere array of telescopes with similar design and capability. Designed in the style of the next-generation VLBI system, these small and fast antennas allow for a new way of observing, comprising higher data rates and more observations than the standard observing sessions coordinated by the International VLBI Service for Geodesy and Astrometry (IVS). In this contribution, the continuous development of the AUSTRAL sessions is described, leading to an improvement of the results in terms of baseline length repeatabilities by a factor of two since the start of this program. The focus is on the scheduling strategy and increased number of observations, aspects of automated operation, and data logistics, as well as results of the 151 AUSTRAL sessions performed so far. The high number of the AUSTRAL sessions makes them an important contributor to VLBI end-products, such as the terrestrial and celestial reference frames and Earth orientation parameters. We compare AUSTRAL results with other IVS sessions and discuss their suitability for the determination of baselines, station coordinates, source coordinates, and Earth orientation parameters.