We present details of the Automated Radio Telescope Imaging Pipeline (ARTIP) and results of a sensitive blind search for HI and OH absorbers at $z<0.4$ and $z<0.7$, respectively. ARTIP is written in Python 3.6, extensively uses the Common Astronomy Software Application (CASA) tools and tasks, and is designed to enable the geographically-distributed MeerKAT Absorption Line Survey (MALS) team to collaboratively process large volumes of radio interferometric data. We apply it to the first MALS dataset obtained using the 64-dish MeerKAT radio telescope and 32K channel mode of the correlator. With merely 40 minutes on target, we present the most sensitive spectrum of PKS1830-211 ever obtained and characterize the known HI ($z=0.19$) and OH ($z=0.89$) absorbers. We further demonstrate ARTIP's capabilities to handle realistic observing scenarios by applying it to a sample of 72 bright radio sources observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) to blindly search for HI and OH absorbers. We estimate the numbers of HI and OH absorbers per unit redshift to be $n_{21}(z\sim0.18)<$0.14 and $n_{\rm OH}(z\sim0.40)<$0.12, respectively, and constrain the cold gas covering factor of galaxies at large impact parameters (50 kpc $<\rho<$ 150 kpc) to be less than 0.022. Due to the small redshift path, $\Delta z\sim$13 for HI with column density$>5.4\times10^{19}$ cm$^{-2}$, the survey has probed only the outskirts of star-forming galaxies at $\rho>30$ kpc. MALS with the expected $\Delta z\sim10^{3-4}$ will overcome this limitation and provide stringent constraints on the cold gas fraction of galaxies in diverse environments over $0<z<1.5$.

We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

The Galactic Centre contains a supermassive black hole with a mass of 4 million suns within an environment that differs markedly from that of the Galactic disk. While the black hole is essentially quiescent in the broader context of active galactic nuclei, X-ray observations have provided evidence for energetic outbursts from its surroundings. Also, while the levels of star formation in the Galactic Centre have been approximately constant over the last few hundred Myr, there is evidence of elevated short-duration bursts, strongly influenced by interaction of the black hole with the enhanced gas density present within the ring-like Central Molecular Zone at Galactic longitude |l| < 0.7 degrees and latitude |b| < 0.2 degrees. The inner 200 pc region is characterized by large amounts of warm molecular gas, a high cosmic ray ionization rate, unusual gas chemistry, enhanced synchrotron emission, and a multitude of radio-emitting magnetised filaments, the origin of which has not been established. Here we report radio imaging that reveals bipolar bubbles spanning 1 degree x 3 degrees (140 parsecs x 430 parsecs), extending above and below the Galactic plane and apparently associated with the Galactic Centre. The structure is edge-brightened and bounded, with symmetry implying creation by an energetic event in the Galactic Centre. We estimate the age of the bubbles to be a few million years, with a total energy of 7 x 10^52 ergs. We postulate that the progenitor event was a major contributor to the increased cosmic-ray density in the Galactic Centre, and is in turn the principal source of the relativistic particles required to power the synchrotron emission of the radio filaments within and in the vicinity of the bubble cavities.

We present the results of a monitoring campaign using the KAT-7 and HartRAO 26m telescopes, of hydroxyl, methanol and water vapour masers associated with the high-mass star forming region G9.62+0.20E. Periodic flaring of the main line hydroxyl masers were found, similar to that seen in the 6.7 and 12.2 GHz methanol masers. The 1667 MHz flares are characterized by a rapid decrease in flux density which is coincident with the start of the 12.2 GHz methanol maser flare. The decrease in the OH maser flux density is followed by a slow increase till a maximum is reached after which the maser decays to its pre-flare level. A possible interpretation of the rapid decrease in the maser flux density is presented. Considering the projected separation between the periodic methanol and OH masers, we conclude that the periodic 12.2 methanol masing region is located about 1600 AU deeper into the molecular envelope compared to the location of the periodic OH masers. A single water maser flare was also detected which seems not to be associated with the same event that gives rise to the periodic methanol and OH maser flares.

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

We present early observations of 20 galaxies in the MHONGOOSE survey galaxies using KAT-7, the GBT, and MeerKAT. We present the best calibrators for five of the MHONGOOSE galaxies observed with the KAT-7, and search for signs of gas accretion in the GBT data, down to $3\sigma$ column density levels of $\rm 2.2\times10^{18}\, cm^{-2}$ over a $\rm 20\, km\, s^{-1}$ linewidth, but identify none. Using the KAT-7 and MeerKAT data, we have derived rotation curves and mass models for NGC 3621 and NGC 7424 out to an unprecedented extent. As a precursor to the SKA, the MeerKAT telescope combines both a high spatial resolution and a large field of view, necessary to map the extended neutral hydrogen in local galaxies. The mass models of the two galaxies were constructed for both the Dark Matter (DM) models (the pseudo-isothermal model and the Navarro-Frenk-White model) and MOND. Overall, we find that the DM models provide a better fit than MOND to the galaxies' rotation curves. Furthermore, the pseudo-isothermal model is found to be the most consistent with the observations.

As a product of the maser monitoring program with the 26m telescope of the Hartebeesthoek Radio Astronomy Observatory (HartRAO), we present an unprecedented, contemporaneous flaring event of 10 maser transitions in hydroxyl, methanol, and water that began in 2015 January in the massive star-forming region NGC6334I in the velocity range -10 to -2 km/s. The 6.7 GHz methanol and 22.2 GHz water masers began flaring within 22 days of each other, while the 12.2 GHz methanol and 1665 MHz hydroxyl masers flared 80 and 113 days later respectively. The 1665 MHz, 6.7 GHz, and 22.2 GHz masers have all remained in their flared state for nearly 3 years. The brightest flaring components increased by factors of 66, 21, 26, and 20 in the 12.2 and 6.7 GHz methanol, 1665 MHz hydroxyl and 22.2 GHz water maser transitions respectively; some weaker components increased by up to a factor of 145. We also report new maser emission in the 1720, 6031, and 6035 MHz OH lines and the 23.1 GHz methanol line, along with the detection of only the fifth 4660 MHz OH maser. We note the correlation of this event with the extraordinary (sub)millimeter continuum outburst from the massive protostellar system NGC6334I-MM1 and discuss the implications of the observed time lags between different maser velocity components on the nature of the outburst. Finally, we identify two earlier epoch maser flaring events likely associated with this object, which suggest a recurring accretive phenomenon that generates powerful radiative outbursts.

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation.

The construction of the KAT-7 array in the Karoo region of the Northern Cape in South Africa was intended primarily as an engineering prototype for technologies and techniques applicable to the MeerKAT telescope. This paper looks at the main engineering and scien- tific highlights from this effort, and discusses their applicability to both MeerKAT and other next-generation radio telescopes. In particular we found that the composite dish surface works well, but it becomes complicated to fabricate for a dish lacking circular symmetry; the Stir- ling cycle cryogenic system with ion pump to achieve vacuum works but demands much higher maintenance than an equivalent Gifford-McMahon cycle system; the ROACH (Recon- figurable Open Architecture Computing Hardware)-based correlator with SPEAD (Stream- ing Protocol for Exchanging Astronomical Data) protocol data transfer works very well and KATCP (Karoo Array Telescope Control Protocol) control protocol has proven very flexible and convenient. KAT-7 has also been used for scientific observations where it has a niche in mapping low surface-brightness continuum sources, some extended HI halos and OH masers in star-forming regions. It can also be used to monitor continuum source variability, observe pulsars, and make VLBI observations

A number of mechanisms to understand the periodic class II methanol masers associated with some high-mass star forming regions have been proposed in the past. Two recent proposals, ie. by Parfenov &Sobolev (2014) and Sanna et al. (2015) have been presented in order to explain the periodic masers in sources with light curves similar to the methanol masers in G9.62+0.20E. We evaluate to what extent the proposals and models presented by these authors can explain the light curve of the methanol masers in G9.62+0.20E. It is argued that neither of the proposed mechanisms can reproduce the light curves of the methanol masers in G9.62+0.20E.

Ten new class II methanol masers from the 6.7-GHz Methanol Multibeam survey catalogues III and IV were selected for a monitoring programme at both 6.7 and 12.2 GHz with the 26m Hartebeesthoek Radio Astronomy Observatory (HartRAO) radio telescope for two years and nine months, from August 2012 to May 2015. In the sample, only masers associated with G339.986-0.425 were found to show periodic variability at both 6.7 and 12.2 GHz. The existence of periodic variation was tested with four independent methods. The analytical method gave the best estimation of the period, which was 246 $\pm$ 1 days. The time series of G339.986-0.425 show strong correlations across velocity channels and between the 6.7 and 12.2 GHz masers. The time delay was also measured across channels and shows structure across the spectrum which is continuous between different maser components.

The Antlia Cluster is a nearby, dynamically young structure, and its proximity provides a valuable opportunity for detailed study of galaxy and group accretion onto clusters. We present a deep HI mosaic completed as part of spectral line commissioning of the Karoo Array Telescope (KAT-7), and identify infrared counterparts from the WISE extended source catalog to study neutral atomic gas content and star formation within the cluster. We detect 37 cluster members out to a radius of ~0.9 Mpc with M_HI > 5x10^7 M_Sun. Of these, 35 are new HI detections, 27 do not have previous spectroscopic redshift measurements, and one is the Compton thick Seyfert II, NGC 3281, which we detect in HI absorption. The HI galaxies lie beyond the X-ray emitting region 200 kpc from the cluster center and have experienced ram pressure stripping out to at least 600 kpc. At larger radii, they are distributed asymmetrically suggesting accretion from surrounding filaments. Combining HI with optical redshifts, we perform a detailed dynamical analysis of the internal substructure, identify large infalling groups, and present the first compilation of the large scale distribution of HI, and star forming galaxies within the cluster. We find that elliptical galaxy NGC 3268 is at the center of the oldest substructure and argue that NGC 3258 and its companion population are more recent arrivals. Through the presence of HI and on-going star formation, we rank substructures with respect to their relative time since accretion onto Antlia.

We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with the kinematic distance estimate for this source at 2.5$\pm $0.5 kpc using the latest Solar and Galactic rotation parameters. We find from the Lyman continuum photon flux that the embedded core of the young star is of spectral type B1, demonstrating that luminous 6.7 GHz methanol masers can be associated with high-mass stars towards the lower end of the mass range.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

Eight new class II methanol masers selected from the 6.7 GHz Methanol Multibeam survey catalogues I and II were monitored at 6.7 GHz with the 26m Hartebeesthoek Radio Astronomy Observatory (HartRAO) radio telescope for three years and seven months, from February 2011 to September 2014. The sources were also observed at 12.2 GHz and two were sufficiently bright to permit monitoring. One of the eight sources, namely G358.460-0.391, was found to show periodic variations at 6.7 GHz. The period was determined and tested for significance using the Lomb-Scargle, epoch-folding and Jurkevich methods, and by fitting a simple analytic function. The best estimate for the period of the 6.7 GHz class II methanol maser line associated with G358.460-0.391 is 220.0 $\pm$ 0.2 day.

We present the results of Herschel HOBYS photometric mapping combined with BIMA observations and additional archival data, and perform an in-depth study of the evolutionary phases of the star-forming clumps in W 48A and their surroundings. Age estimates for the compact sources were derived from bolometric luminosities and envelope masses, which were obtained from the dust continuum emission, and agree within an order of magnitude with age estimates from molecular line and radio data. The clumps in W 48A are linearly aligned by age (east-old to west-young): we find a ultra compact (UC) HII region, a young stellar object (YSO) with class II methanol maser emission, a YSO with a massive outflow, and finally the NH_2D prestellar cores from Pillai et al. This remarkable positioning reflects the (star) formation history of the region. We find that it is unlikely that the star formation in the W 48A molecular cloud was triggered by the UCHII region and discuss the Aquila supershell expansion as a mayor influence on the evolution of W 48A. We conclude that the combination of Herschel continuum data with interferometric molecular line and radio continuum data is important to derive trustworthy age estimates and interpret the origin of large scale structures through kinematic information.

We report the results of 10 years of monitoring of six regularly varying 6.7 GHz methanol masers using the Hartebeesthoek 26m telescope. Observations were done at intervals of 1--2 weeks, with faster sampling during flaring episodes. Four of the sources were also monitored at 12.2 GHz and show correlated variations. We find the Lomb-Scargle periodogram to be the most sensitive method to search for periodicity but possibly prone to false detections. Periods range from 132.8 days (with 26 cycles observed) to 509 days (with 7 cycles observed). Five of the sources show arguably periodic variations, while G331.13-0.24 shows strong periodicity in one peak, with large and variable delays in other peaks.

HI observations of the Magellanic-type spiral NGC 3109, obtained with the seven dish Karoo Array Telescope (KAT-7), are used to analyze its mass distribution. Our results are compared to what is obtained using VLA data. KAT-7 is the precursor of the SKA pathfinder MeerKAT, which is under construction. The short baselines and low system temperature of the telescope make it sensitive to large scale low surface brightness emission. The new observations with KAT-7 allow the measurement of the rotation curve of NGC 3109 out to 32', doubling the angular extent of existing measurements. A total HI mass of 4.6 x 10^8 Msol is derived, 40% more than what was detected by the VLA observations. The observationally motivated pseudo-isothermal dark matter (DM) halo model can reproduce very well the observed rotation curve but the cosmologically motivated NFW DM model gives a much poorer fit to the data. While having a more accurate gas distribution has reduced the discrepancy between the observed RC and the MOdified Newtonian Dynamics (MOND) models, this is done at the expense of having to use unrealistic mass-to-light ratios for the stellar disk and/or very large values for the MOND universal constant a0. Different distances or HI contents cannot reconcile MOND with the observed kinematics, in view of the small errors on those two quantities. As for many slowly rotating gas-rich galaxies studied recently, the present result for NGC 3109 continues to pose a serious challenge to the MOND theory.

We present the light curves of the 6.7 and 12.2 GHz methanol masers in the star forming region G9.62+0.20E for a time span of more than 2600 days. The earlier reported period of 244 days is confirmed. The results of monitoring the 107 GHz methanol maser for two flares are also presented. The results show that flaring occurs in all three masing transitions. It is shown that the average flare profiles of the three masing transitions are similar. The 12.2 GHz masers are the most variable of the three masers with the largest relative amplitude having a value of 2.4. The flux densities for the different masing transitions are found to return to the same level during the low phase of the masers, suggesting that the source of the periodic flaring is situated outside the masing region, and that the physical conditions in the masing region are relatively stable. On the basis of the shape of the light curve we excluded stellar pulsations as the underlying mechanism for the periodicity. It is argued that a colliding wind binary can account for the observed periodicity and provide a mechanism to qualitatively explain periodicity in the seed photon flux and/or the pumping radiation field. It is also argued that the dust cooling time is too short to explain the decay time of about 100 days of the maser flare. A further analysis has shown that for the intervals from days 48 to 66 and from days 67 to 135 the decay of the maser light curve can be interpreted as due to the recombination of a thermal hydrogen plasma with densities of approximately $1.6 \times 10^6 \mathrm{cm^{-3}}$ and $6.0 \times 10^5 \mathrm{cm^{-3}}$ respectively.

Time series are presented for the class II methanol maser source G12.89+0.49, which has been monitored for nine years at the Hartebeesthoek Radio Astronomy Observatory. The 12.2 and 6.7 GHz methanol masers were seen to exhibit rapid, correlated variations on timescales of less than a month. Daily monitoring has revealed that the variations have a periodic component with a period of 29.5 days. The period seems to be stable over the 110 cycles spanned by the time series. There are variations from cycle to cycle, with the peak of the flare occurring anywhere within an eleven day window but the minima occur at the same phase of the cycle. Time delays of up to 5.7 days are seen between spectral features at 6.7 GHz and a delay of 1.1 day is seen between the dominant 12.2 GHz spectral feature and its 6.7 GHz counterpart.

(Abridged) Recently, the magnetic field induced Zeeman splitting was measured for the strongest known 6.7 GHz methanol maser, which arises in the massive star forming region G09.62+0.20. This maser is one of a handful of periodically flaring methanol masers. The 100-m Effelsberg telescope was used to monitor the 6.7 GHz methanol masers of G09.62+0.20. With the exception of a two week period during the peak of the maser flare, we measure a constant magnetic field of B_||~11+-2 mG in the two strongest maser components of G09.62+0.20 that are separated by over 200 AU. In the two week period that coincides exactly with the peak of the maser flare of the strongest maser feature, we measure a sharp decrease and possible reversal of the Zeeman splitting. The exact cause of both maser and polarization variability is still unclear, but it could be related to either background amplification of polarized emission or the presence of a massive protostar with a close-by companion. Alternatively, the polarization variability could be caused by non-Zeeman effects related to the radiative transfer of polarized maser emission.

An intensive monitoring program of 54 6.7-GHz methanol maser sources was carried out at the Hartebeesthoek Radio Astronomy Observatory from January 1999 to April 2003. The monitoring program was subsequently continued on 19 sources of interest. Analysis of the resulting time-series stretching over eight years shows that six of the sources are periodic, with periods ranging from 133 days to 504 days. The waveforms in individual sources range from sinusoidal fluctuations to sharp flares and there can be other long term trends in the time-series. The amplitudes of the variations can also change from cycle to cycle. The time-series of the periodic masers will be presented, and possible causes of the variability discussed.

General characteristics of methanol (CH3OH) maser emission are summarized. It is shown that methanol maser sources are concentrated in the spiral arms. Most of the methanol maser sources from the Perseus arm are associated with embedded stellar clusters and a considerable portion is situated close to compact HII regions. Almost 1/3 of the Perseus Arm sources lie at the edges of optically identified HII regions which means that massive star formation in the Perseus Arm is to a great extent triggered by local phenomena. A multiline analysis of the methanol masers allows us to determine the physical parameters in the regions of maser formation. Maser modelling shows that class II methanol masers can be pumped by the radiation of the warm dust as well as by free-free emission of a hypercompact region hcHII with a turnover frequency exceeding 100 GHz. Methanol masers of both classes can reside in the vicinity of hcHIIs. Modelling shows that periodic changes of maser fluxes can be reproduced by variations of the dust temperature by a few percent which may be caused by variations in the brightness of the central young stellar object reflecting the character of the accretion process. Sensitive observations have shown that the masers with low flux densities can still have considerable amplification factors. The analysis of class I maser surveys allows us to identify four distinct regimes that differ by the series of their brightest lines.

The class II methanol maser source G9.62+0.20E undergoes periodic flares at both 6.7 and 12.2 GHz. The flare starting in 2001 October was observed at seven epochs over three months using the VLBA at 12.2 GHz. High angular resolution images (beam size $\sim$ 1.7 x 0.6 mas) were obtained, enabling us to observe changes in 16 individual maser components. It was found that while existing maser spots increased in flux density, no new spots developed and no changes in morphology were observed. This rules out any mechanism which disturbs the masing region itself, implying that the flares are caused by a change in either the seed or pump photon levels. A time delay of 1--2 weeks was observed between groups of maser features. These delays can be explained by light travel time between maser groups. The regularity of the flares can possibly be explained by a binary system.