In this paper, we propose a new approach to determining cosmological distances to active galactic nuclei (AGN) via light travel-time arguments, which can be extended from nearby sources to very high redshift sources. The key assumption is that the variability seen in AGN is constrained by the speed of light and therefore provides an estimate of the linear size of an emitting region. This can then be compared with the angular size measured with very long baseline interferometry (VLBI) in order to derive a distance. We demonstrate this approach on a specific well studied low redshift (z = 0.0178) source 3C84 (NGC 1275), which is the bright radio core of the Perseus Cluster. We derive an angular diameter distance including statistical errors of $D_{A} = 72^{+5}_{-6}$ Mpc for this source, which is consistent with other distance measurements at this redshift. Possible sources of systematic errors and ways to correct for them are discussed.

We measure the 10 and $18\mu$m silicate features in a sample of 67 local ($z<0.1$) type 1 active galactic nuclei (AGN) with available {\it Spitzer} spectra dominated by non-stellar processes. We find that the $10\mu$m silicate feature peaks at $10.3^{+0.7}_{-0.9}\mu$m with a strength (Si$_{p}$ = ln f$_{p}$(spectrum)/f$_{p}$(continuum)) of $0.11^{+0.15}_{-0.36}$, while the $18\mu$m one peaks at $17.3^{+0.4}_{-0.7}\mu$m with a strength of $0.14^{+0.06}_{-0.06}$. We select from this sample sources with the strongest 10$\mu$m silicate strength ($\sigma_{Si_{10\mu m}}>0.28$, 10 objects). We carry out a detailed modeling of the IRS/{\it Spitzer} spectra by comparing several models that assume different geometries and dust composition: a smooth torus model, two clumpy torus models, a two-phase medium torus model, and a disk+outflow clumpy model. We find that the silicate features are well modeled by the clumpy model of Nenkova et al. 2008, and among all models those including outflows and complex dust composition are the best (Hoenig et al. 2017). We note that even in AGN-dominated galaxies it is usually necessary to add stellar contributions to reproduce the emission at the shortest wavelengths.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

The flat spectrum radio quasar 4C 38.41 showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional monitoring observations at various bands across the electromagnetic spectrum. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components and the evolution of the physical properties seem to be in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1 mG, do not seem to dramatically change during the flares, and are much smaller, by a factor 10,000, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

PKS 1510--089 is a bright and active $\gamma$-ray source that showed strong and complex $\gamma$-ray flares in mid-2015 during which the Major Atmospheric Gamma Imaging Cherenkov telescopes detected variable very high energy (VHE; photon energies $>$100 GeV) emission. We present long-term multi-frequency radio, optical, and $\gamma$-ray light curves of PKS 1510--089 from 2013 to 2018, and results of an analysis of the jet kinematics and linear polarization using 43 GHz Very Long Baseline Array data observed between late 2015 and mid-2017. We find that a strong radio flare trails the $\gamma$-ray flares in 2015, showing an optically thick spectrum at the beginning and becoming optically thin over time. Two laterally separated knots of emission are observed to emerge from the radio core nearly simultaneously during the $\gamma$-ray flares. We detect an edge-brightened linear polarization near the core in the active jet state in 2016, similar to the quiescent jet state in 2008--2013. These observations indicate that the $\gamma$-ray flares may originate from compression of the knots by a standing shock in the core and the jet might consist of multiple complex layers showing time-dependent behavior, rather than of a simple structure of a fast jet spine and a slow jet sheath.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.

A long standing problem in the study of Active Galactic Nuclei (AGNs) is that the observed VLBI core is in fact a blending of the actual AGN core (classically defined by the $\tau=1$ surface) and the upstream regions of the jet or optically thin emitting region flows. This blending may cause some biases towards the observables of the core, such as its flux density, size or brightness temperature, which may lead to misleading interpretation of the derived quantities and physics. We study the effects of such blending under the view of the Korean VLBI Network (KVN) for a sample of AGNs at 43 GHz by comparing their observed properties with observations with the Very Large Baseline Array (VLBA). Our results suggest that the observed core sizes are a factor $\sim11$ larger than these of VLBA, which is similar to the factor expected by considering the different resolutions of the two facilities. We suggest the use of this factor to consider blending effects in KVN measurements. Other parameters, such as flux density or brightness temperature, seem to possess a more complicated dependence.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

Investigating the magnetic field structure in the innermost regions of relativistic jets is fundamental to shed light on the crucial physical processes giving rise to the jet formation, as well as to its extraordinary radiation output up to gamma-ray energies. We study the magnetic field structure of the quasar CTA 102 with 3 and 7 mm-VLBI polarimetric observations, reaching an unprecedented resolution (~50 microarcsec). We also investigate the variability and physical process occurring in the source during the observing period which coincides with a very active state of the source till high-energies. The Faraday rotation analysis between 3 and 7mm shows a gradient in rotation measure with a maximum value of ~6X10^4 rad/m^2 and intrinsic electric vector position angles (EVPAs) oriented around the centroid of the core, suggesting the presence of large-scale helical magnetic fields. Such a magnetic field structure is also visible in 7 mm images when a new superluminal component is crossing the core region. The 7mm EVPAs orientation is different when the component is exiting the core or crossing a stationary feature at ~0.1 mas. The interaction between the superluminal component and a recollimation shock at ~0.1 mas could have triggered the multi-wavelengths flares. The variability Doppler factor associated with such interaction is large enough to explain the high energy emission, as we infer from the analysis of gamma-ray and X-ray data, and it is in agreement with the Doppler factor obtained to explain the extraordinary optical flare by Raiteri et al.(2017).

We report results from a deep polarization imaging of the nearby radio galaxy 3C$\,$84 (NGC$\,$1275). The source was observed with the Global Millimeter VLBI Array (GMVA) at 86$\,$GHz at an ultra-high angular resolution of $50\mu$as (corresponding to 250$R_{s}$). We also add complementary multi-wavelength data from the Very Long Baseline Array (VLBA; 15 & 43$\,$GHz) and from the Atacama Large Millimeter/submillimeter Array (ALMA; 97.5, 233.0, and 343.5$\,$GHz). At 86$\,$GHz, we measure a fractional linear polarization of $\sim2$% in the VLBI core region. The polarization morphology suggests that the emission is associated with an underlying limb-brightened jet. The fractional linear polarization is lower at 43 and 15$\,$GHz ($\sim0.3-0.7$% and $<0.1$%, respectively). This suggests an increasing linear polarization degree towards shorter wavelengths on VLBI scales. We also obtain a large rotation measure (RM) of $\sim10^{5-6}~{\rm rad/m^{2}}$ in the core at $\gtrsim$43$\,$GHz. Moreover, the VLBA 43$\,$GHz observations show a variable RM in the VLBI core region during a small flare in 2015. Faraday depolarization and Faraday conversion in an inhomogeneous and mildly relativistic plasma could explain the observed linear polarization characteristics and the previously measured frequency dependence of the circular polarization. Our Faraday depolarization modeling suggests that the RM most likely originates from an external screen with a highly uniform RM distribution. To explain the large RM value, the uniform RM distribution, and the RM variability, we suggest that the Faraday rotation is caused by a boundary layer in a transversely stratified jet. Based on the RM and the synchrotron spectrum of the core, we provide an estimate for the magnetic field strength and the electron density of the jet plasma.

Located at the center of the Perseus cluster, 3C 84 is an extremely bright and nearby radio galaxy. Because of the strong diffuse thermal emission from the cluster in X-rays, the detailed properties and the origin of a power-law component from the central active galactic nucleus (AGN) remains unclear in the source. We report here the first NuSTAR observations of 3C 84. The source was observed for 24.2 and 32 ks on February 01 and 04, 2018, respectively. NuSTAR observations spectrally decompose the power-law AGN component above 10 keV. The power-law component dominates the spectrum above 20 keV with a photon index $\sim$1.9 and an energy flux F$_{20-30 keV}$ = 1.0 $\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$, corresponding to an isotropic luminosity, $L_{20-30 keV}$ = 7.4$\times$10$^{42}$ erg s$^{-1}$. We discuss possible emitting sites for the power-law component. The expected thermal emission from the accretion disk is not hot enough to account for the hard X-rays detected from the source. Similar X-ray and $\gamma$-ray photon indices and long-term flux variations, the absence of cutoff energy in the hard X-ray spectrum of the source, correlated hard X-ray flux and hardness ratio variations, and the similarity of optical-X-ray slope to blazar rather than Seyfert galaxies supports the hard X-ray power-law component originating from the jet.

Recent Fermi-Large Area Telescope (LAT) light curves indicate an active $\gamma$-ray state spanning about five months from 2016 June to 2016 October in the BL Lac object 1749+096 (OT 081). During this period, we find two notable $\gamma$-ray events: an exceptionally strong outburst followed by a significant enhancement (local peak). In this study, we analyze multi-waveband light curves (radio, optical, X-ray, and $\gamma$-ray) plus very-long baseline interferometry (VLBI) data to investigate the nature of the $\gamma$-ray events. The $\gamma$-ray outburst coincides with flux maxima at longer wavelengths. We find a spectral hardening of the $\gamma$-ray photon index during the $\gamma$-ray outburst. The photon index shows a transition from a softer-when-brighter to a harder-when-brighter trend at around 1.8 $\times$ $10^{-7}$ ph cm$^{-2}$ s$^{-1}$. We see indication that both the $\gamma$-ray outburst and the subsequent enhancement precede the propagation of a polarized knot in a region near the VLBI core. The highest polarized intensity, 230\,mJy, and an electric vector position angle rotation, by $\sim$32$^{\circ}$, are detected about 12 days after the $\gamma$-ray outburst. We conclude that both $\gamma$-ray events are caused by the propagation of a disturbance in the mm-wave core.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We study the linear polarization of the radio cores of eight blazars simultaneously at 22, 43, and 86 GHz with observations obtained by the Korean VLBI Network (KVN) in three epochs between late 2016 and early 2017 in the frame of the Plasma-physics of Active Galactic Nuclei (PAGaN) project. We investigate the Faraday rotation measure (RM) of the cores; the RM is expected to increase with observing frequency if core positions depend on frequency due to synchrotron self-absorption. We find a systematic increase of RMs at higher observing frequencies in our targets. The RM--$\nu$ relations follow power-laws with indices distributed around 2, indicating conically expanding outflows serving as Faraday rotating media. Comparing our KVN data with contemporaneous optical polarization data from the Steward Observatory for a few sources, we find indication that the increase of RM with frequency saturates at frequencies of a few hundreds GHz. This suggests that blazar cores are physical structures rather than simple $\tau=1$ surfaces. A single region, e.g. a recollimation shock, might dominate the jet emission downstream of the jet launching region. We detect a sign change in the observed RMs of CTA 102 on a time scale of $\approx$1 month, which might be related to new superluminal components emerging from its core undergoing acceleration/deceleration and/or bending. We see indication for quasars having higher core RMs than BL Lac objects, which could be due to denser inflows/outflows in quasars.

The flat spectrum radio quasar 1633+382 (4C~38.41) showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional radio monitoring observations with the OVRO 40m telescope, the Boston University VLBI program, and the Submillimeter Array. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components. Analysis of the spectral energy distribution indicates a higher turnover frequency after the flaring events. The evolution of the flare in the turnover frequency - turnover flux density plane probes the adiabatic losses in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1mG, do not seem to dramatically change during the flares, and are much smaller, by a factor $10^4$, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

Understanding the launching, acceleration, and collimation of jets powered by active galactic nuclei remains an outstanding problem in relativistic astrophysics. This is partly because observational tests of jet formation models suffer from the limited angular resolution of ground-based very long baseline interferometry that has thus far been able to probe the transverse jet structure in the acceleration and collimation zone of only two sources. Here we report radio interferometric observations of 3C 84 (NGC 1275), the central galaxy of the Perseus cluster, made with an array including the orbiting radio telescope of the RadioAstron mission. The obtained image transversely resolves the edge-brightened jet in 3C 84 only 30 microarcseconds from the core, which is ten times closer to the central engine than what has been possible in previous ground-based observations, and it allows us to measure the jet collimation profile from ~ 100 to ~10000 gravitational radii from the black hole. The previously found, almost cylindrical jet profile on scales larger than a few thousand r_g is now seen to continue at least down to a few hundred r_g from the black hole and we find a broad jet with a transverse radius larger than about 250 r_g at only 350 r_g from the core. If the bright outer jet layer is launched by the black hole ergosphere, it has to rapidly expand laterally on scales smaller than 100 r_g. If this is not the case, then this jet sheath is likely launched from the accretion disk.

3C 84 (NGC 1275) is a well-studied mis-aligned Active Galactic Nucleus (AGN), which has been active in Gamma rays since at least 2008. We have monitored the source at four wavelengths (14 mm, 7 mm, 3 mm and 2 mm) using the Korean VLBI network (KVN) since 2013 as part of the interferometric monitoring of $\gamma$-ray bright AGN (iMOGABA) program. 3C 84 exhibits bright radio emission both near the central supermassive black hole (SMBH) feature known as C1 and from a moving feature located to the south known as C3. Other facilities have also detected these short-term variations above a slowly rising trend at shorter wavelengths, such as in Gamma ray and 1 mm total intensity light-curves. We find that the variations in the $\gamma$ rays and 1 mm total intensity light-curves are correlated, with the $\gamma$ rays leading and lagging the radio emission. Analysis of the 2 mm KVN data shows that both the Gamma rays and 1 mm total intensity short-term variations are better correlated with the SMBH region than C3, likely placing the short-term variations in C1. We interpret the emission as being due to the random alignment of spatially separated emission regions. We place the slowly rising trend in C3, consistent with previous results. Additionally, we report that since mid-2015, a large mm-wave radio flare has been occurring in C3, with a large Gamma ray flare coincident with the onset of this flare at all radio wavelengths.

3C 84 is a nearby Active Galactic Nucleus (AGN) that is unique in that is believed that we are observing near the true jet launching region - unlike blazars. The source is active in Gamma rays and has been detected with Fermi since its launch in 2008, including being detected at TeV energies with other instruments. Due to the relative proximity of the source (z=0.018), it provides a unique opportunity to pinpoint the location of the $\gamma$-ray emission by combining the Gamma ray data with very long baseline inteferometry (VLBI) data. A study using the Korean VLBI network (KVN) showed that the Gamma rays occur in both downstream jet emission and the region near where the jet is launched. Further analysis of the kinematics using Wavelet Image Segmentation and Evaluation (WISE) algorithm, which uses 2-dimensional cross-correlations to statistically derive the kinematics of high-resolution 7 mm VLBA data show that the Gamma ray emission is caused by a fast-travelling shock catching a slower moving shock and then interacting with the external medium, in behaviour reminiscent of a long duration gamma-ray burst (GRB). This could explain why such high energy flaring is seen in such low Doppler boosted sources. Finally, we show some early results from a study of the jet launching region using the Global mm-VLBI Array (GMVA). The nucleus appears to have a consistent double nuclear structure that is likely too broad to be the true jet base.

We study the centimeter- to millimeter-wavelength synchrotron spectrum of the core of the radio galaxy M87 at $\lesssim0.8\,{\rm mas}~\sim110R_{s}$ spatial scales using four years of fully simultaneous, multi-frequency VLBI data obtained by the Korean VLBI Network (KVN). We find a core spectral index $\alpha$ of $\gtrsim-0.37$ ($S\propto \nu^{+\alpha}$) between 22GHz and 129GHz. By combining resolution-matched flux measurements from the Very Long Baseline Array (VLBA) at 15GHz and taking the Event Horizon Telescope (EHT) 230GHz core flux measurements in epochs 2009 and 2012 as lower limits, we find evidence of a nearly flat core spectrum across 15GHz and 129GHz, which could naturally connect the 230GHz VLBI core flux. The extremely flat spectrum is a strong indication that the jet base does not consist of a simple homogeneous plasma, but of inhomogeneous multi-energy components, with at least one component with the turn-over frequency $\gtrsim100$GHz. The spectral shape can be qualitatively explained if both the strongly (compact, optically thick at $>$100GHz) and the relatively weakly magnetized (more extended, optically thin at $<$100GHz) plasma components are colocated in the footprint of the relativistic jet.

Atmospheric propagation effects at millimeter wavelengths can significantly alter the phases of radio signals and reduce the coherence time, putting tight constraints on high frequency Very Long Baseline Interferometry (VLBI) observations. In previous works, it has been shown that non-dispersive (e.g. tropospheric) effects can be calibrated with the frequency phase transfer (FPT) technique. The coherence time can thus be significantly extended. Ionospheric effects, which can still be significant, remain however uncalibrated after FPT as well as the instrumental effects. In this work, we implement a further phase transfer between two FPT residuals (i.e. so-called FPT-square) to calibrate the ionospheric effects based on their frequency dependence. We show that after FPT-square, the coherence time at 3 mm can be further extended beyond 8~hours, and the residual phase errors can be sufficiently canceled by applying the calibration of another source, which can have a large angular separation from the target (>20 deg) and significant temporal gaps. Calibrations for all-sky distributed sources with a few calibrators are also possible after FPT-square. One of the strengths and uniqueness of this calibration strategy is the suitability for high-frequency all-sky survey observations including very weak sources. We discuss the introduction of a pulse calibration system in the future to calibrate the remaining instrumental effects and allowing the possibility of imaging the source structure at high frequencies with FPT-square, where all phases are fully calibrated without involving any additional sources.

We present multi-frequency simultaneous VLBI radio observations of the flat spectrum radio quasar 1633+382 (4C~38.41) as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program combined with additional observations in radio, optical, X-rays and $\gamma-$rays carried out between the period 2012 March - 2015 August. The monitoring of this source reveals a significant long-lived increase in its activity since approximately two years in the radio bands, which correlates with a similar increase on all other bands from sub-millimeter to $\gamma-$rays. A significant correlation is also found between radio fluxes and simultaneous spectral indices during this period. The study of the discrete correlation function (DCF) indicates time lags smaller than the $\sim40$ days uncertainties among both radio bands and also high-energy bands, and a time lag of $\sim$70 days, with $\gamma-$rays leading radio. We interpret that the high-energy and radio fluxes are arising from different emitting regions, located at $1\pm12$ and $40\pm13$ pc from the central engine respectively.

We present total and linearly polarized 3 mm Global mm-VLBI Array images of a sample of blazars and radio galaxies from the VLBA-BU-BLAZAR 7 mm monitoring program designed to probe the innermost regions of active galactic nuclei (AGN) jets and locate the sites of gamma-ray emission observed by the Fermi-LAT. The lower opacity at 3 mm and improved angular resolution, on the order of 50 microarcseconds, allow us to distinguish features in the jet not visible in the 7 mm VLBA data. We also compare two different methods used for the calibration of instrumental polarisation and we analyze the resulting images for some of the sources in the sample.

We present a study of the inexplicit connection between radio jet activity and gamma-ray emission of BL Lacertae (BL Lac; 2200+420). We analyze the long-term millimeter activity of BL Lac via interferometric observations with the Korean VLBI Network (KVN) obtained at 22, 43, 86, and 129 GHz simultaneously over three years (from January 2013 to March 2016); during this time, two gamma-ray outbursts (in November 2013 and March 2015) can be seen in gamma-ray light curves obtained from Fermi observations. The KVN radio core is optically thick at least up to 86 GHz; there is indication that it might be optically thin at higher frequencies. To first order, the radio light curves decay exponentially over the time span covered by our observations, with decay timescales of 411+/-85 days, 352+/-79 days, 310+/-57 days, and 283+/-55 days at 22, 43, 86, and 129 GHz, respectively. Assuming synchrotron cooling, a cooling time of around one year is consistent with magnetic field strengths B~2microT and electron Lorentz factors gamma~10,000. Taking into account that our formal measurement errors include intrinsic variability and thus over-estimate the statistical uncertainties, we find that the decay timescale tau scales with frequency nu like tau~nu^{-0.2}. This relation is much shallower than the one expected from opacity effects (core shift), but in agreement with the (sub-)mm radio core being a standing recollimation shock. We do not find convincing radio flux counterparts to the gamma$ray outbursts. The spectral evolution is consistent with the `generalized shock model' of Valtaoja et al. (1992). A temporary increase in the core opacity and the emergence of a knot around the time of the second gamma-ray event indicate that this gamma-ray outburst might be an `orphan' flare powered by the `ring of fire' mechanism.

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

We present the results of very long baseline interferometry (VLBI) observations of gamma-ray bright blazar S5 0716+714 using the Korean VLBI Network (KVN) at the 22, 43, 86, and 129 GHz bands, as part of the Interferometric Monitoring of Gamma-ray Bright AGNs (iMOGABA) KVN key science program. Observations were conducted in 29 sessions from January 16, 2013 to March 1, 2016, with the source being detected and imaged at all available frequencies. In all epochs, the source was compact on the milliarcsecond (mas) scale, yielding a compact VLBI core dominating the synchrotron emission on these scales. Based on the multi-wavelength data between 15 GHz (Owens Valley Radio Observatory) and 230 GHz (Submillimeter Array), we found that the source shows multiple prominent enhancements of the flux density at the centimeter (cm) and millimeter (mm) wavelengths, with mm enhancements leading cm enhancements by -16$\pm$8 days. The turnover frequency was found to vary between 21 to 69GHz during our observations. By assuming a synchrotron self-absorption model for the relativistic jet emission in S5 0716+714, we found the magnetic field strength in the mas emission region to be $\le$5 mG during the observing period, yielding a weighted mean of 1.0$\pm$0.6 mG for higher turnover frequencies (e.g., >45 GHz).

The mis-aligned active galactic nuclei (AGN), 3C\,84, has been observed approximately monthly at four frequencies simultaneously (22\,GHz -- 129\,GHz) on the Korean VLBI Network (KVN) since 2013 as part of the interferometric monitoring of $\gamma$-ray bright AGN (iMOGABA) program. 3C\,84 is known to have to sites of bright radio emission, with the first thought to be near the central super-massive black hole and the second in a slowly moving emission feature, several milliarcseconds south of the core. Analysis of this data suggests that the highly variable Gamma-ray emission originates near the SMBH and the slowly rising $\gamma$-ray emission originates in the slow moving feature.

We present results of single-epoch very long baseline interferometry (VLBI) observations of gamma-ray bright active galactic nuclei (AGNs) using the Korean VLBI Network (KVN) at 22, 43, 86, and 129~GHz bands, which are part of a KVN key science program, Interferometric Monitoring of Gamma-ray Bright AGNs (iMOGABA). We selected a total of 34 radio-loud AGNs of which 30 sources are gamma-ray bright AGNs with flux densities of $>6\times10^{-10}$~ph~cm$^{-2}$~s$^{-1}$. Single-epoch multi-frequency VLBI observations of the target sources were conducted during a 24-hr session on 2013 November 19 and 20. All observed sources were detected and imaged at all frequency bands with or without a frequency phase transfer technique which enabled the imaging of 12 faint sources at 129~GHz, except for one source. Many of the target sources are resolved on milliarcsecond scales, yielding a core-jet structure with the VLBI core dominating the synchrotron emission on the milliarcsecond scale. CLEAN flux densities of the target sources are 0.43-28~Jy, 0.32-21~Jy, 0.18-11~Jy, and 0.35-8.0~Jy in the 22, 43, 86, and 129~GHz bands, respectively.

Powerful radio galaxies are often associated with gas-rich galaxy mergers. These mergers may provide the fuel to trigger starburst and active galactic nuclear (AGN) activity. In this Research Note, we study the host galaxies of three seemingly young or re-started radio sources that drive fast outflows of cool neutral hydrogen (HI) gas, namely 3C 293, 3C 305 and 4C 12.50 (PKS 1345+12). Our aim is to link the feedback processes in the central kpc-scale region with new information on the distribution of stars and gas at scales of the galaxy. For this, we use deep optical V-band imaging of the host galaxies, complemented with HI emission-line observations to study their gaseous environments. We find prominent optical tidal features in all three radio galaxies, which confirm previous claims that 3C 293, 3C 305 and 4C 12.50 have been involved in a recent galaxy merger or interaction. Our data show the complex morphology of the host galaxies, and identify the companion galaxies that are likely involved in the merger or interaction. The radio sources appear to be (re-)triggered at a different stage of the merger; 4C 12.50 is a pre-coalescent and possibly multiple merger, 3C 293 is a post-coalescent merger that is undergoing a minor interaction with a close satellite galaxy, while 3C 305 appears to be shaped by an interaction with a gas-rich companion. For 3C 293 and 3C 305, we do not detect HI beyond the inner ~30-45 kpc region, which shows that the bulk of the cold gas is concentrated within the host-galaxy, rather than along the widespread tidal features.

The high radio frequency polarization imaging of non-thermal emission from active galactic nuclei (AGN) is a direct way to probe the magnetic field strength and structure in the immediate vicinity of supermassive black holes (SMBHs) and is crucial in testing the jet-launching scenario. To explore the the magnetic field configuration at the base of jets in blazars, we took advantage of the full polarization capabilities of the Global Millimeter VLBI Array (GMVA). With an angular resolution of $\sim$50 micro-arcseconds ($\mu$as) at 86 GHz, one could resolve scales up to $\sim$450 gravitational radii (for a 10$^9$ solar mass black hole at a redshift of 0.1). We present here the preliminary results of our study on the blazar BL~Lac. Our results suggest that on sub-mas scales the core and the central jet of BL Lac are significantly polarized with two distinct regions of polarized intensity. We also noted a great morphological similarity between the 7mm/3mm VLBI images at very similar angular resolution.

The calibration of Very Long Baseline Interferometry (VLBI) data has long been a time consuming process. The Korean VLBI Network (KVN) is a simple array consisting of three identical antennas. Because four frequencies are observed simultaneously, phase solutions can be transferred from lower frequencies to higher frequencies in order to improve phase coherence and hence sensitivity at higher frequencies. Due to the homogeneous nature of the array, the KVN is also well suited for automatic calibration. In this paper we describe the automatic calibration of single-polarisation KVN data using the KVN Pipeline and comparing the results against VLBI data that has been manually reduced. We find that the pipelined data using phase transfer produces better results than a manually reduced dataset not using the phase transfer. Additionally we compared the pipeline results with a manually reduced phase-transferred dataset and found the results to be identical.

The Gamma-ray BL Lac object OJ 287 is known to exhibit inner-parsec "jet-wobbling", high degrees of variability at all wavelengths and quasi-stationary features including an apparent (~100 deg) position angle change in projection on the sky plane. Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps together with cm/mm total intensity and Gamma-ray observations from Fermi/LAT from 2008-2014, we aimed to determine the location of Gamma-ray emission and to explain the inner-mas structural changes. Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and allow us to observe above the synchrotron self-absorption peak frequency. The jet was spectrally decomposed at multiple locations along the jet. From this we derived estimates of the magnetic field. How the field decreases down the jet allowed an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics we attempt to locate the site of Gamma-ray activity, radio flares and spectral changes. Strong Gamma-ray flares appeared to originate from either the "core" region, a downstream stationary feature, or both, with Gamma-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be >1.6 G in the "core" and >0.4 G in the downstream quasi-stationary feature. We therefore found upper limits on the location of the "core" as >6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss.

One of the most intriguing and challenging quests of current astrophysics is to understand the physical conditions and processes responsible for production of high-energy particles, and emission of \gamma-rays. A combination of high-resolution Very Long Baseline Interferometry (VLBI) images with broadband flux variability measurements is a unique way to probe the emission mechanisms at the bases of jets. Our analysis of \gamma-ray flux variability observed by the Fermi-LAT (Large Area Telescope) along with the parsec-scale jet kinematics suggests that the $\gamma$-ray emission in blazar S5 0716+714 has a significant correlation with the mm-VLBI core flux and the orientation of jet outflow on parsec scales. These results indicate that the inner jet morphology has a tight connection with the observed $\gamma$-ray flares. An overview of our current understanding on high-energy radiation processes, their origin, and location is presented here.

Context. Blazars are among the most energetic objects in the Universe. In 2008 August, Fermi/LAT detected the blazar PKS 1502+106 showing a rapid and strong gamma-ray outburst followed by high and variable flux over the next months. This activity at high energies triggered an intensive multi-wavelength campaign covering also the radio, optical, UV, and X-ray bands indicating that the flare was accompanied by a simultaneous outburst at optical/UV/X-rays and a delayed outburst at radio bands. Aims: In the current work we explore the phenomenology and physical conditions within the ultra-relativistic jet of the gamma-ray blazar PKS 1502+106. Additionally, we address the question of the spatial localization of the MeV/GeV-emitting region of the source. Methods: We utilize ultra-high angular resolution mm-VLBI observations at 43 and 86 GHz complemented by VLBI observations at 15 GHz. We also employ single-dish radio data from the F-GAMMA program at frequencies matching the VLBI monitoring. Results: PKS 1502+106 shows a compact core-jet morphology and fast superluminal motion with apparent speeds in the range 5--22 c. Estimation of Doppler factors along the jet yield values between ~7 up to ~50. This Doppler factor gradient implies an accelerating jet. The viewing angle towards the source differs between the inner and outer jet, with the former at ~3 degrees and the latter at ~1 degree, after the jet bends towards the observer beyond 1 mas. The de-projected opening angle of the ultra-fast, magnetically-dominated jet is found to be (3.8 +/- 0.5) degrees. A single jet component can be associated with the pronounced flare both at high-energies and in radio bands. Finally, the gamma-ray emission region is localized at less than 5.9 pc away from the jet base.

We report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

The Global mm-VLBI Array (GMVA) is a network of 14 3\,mm and 7\,mm capable telescopes spanning Europe and the United States, with planned extensions to Asia. The array is capable of sensitive maps with angular resolution often exceeding 50\,$\mu$as. Using the GMVA, a large sample of prominent $\gamma$-ray blazars have been observed approximately 6 monthly from later 2008 until now. Combining 3\,mm maps from the GMVA with near-in-time 7\,mm maps from the VLBA-BU-BLAZAR program and 2\,cm maps from the MOJAVE program, we determine the sub-pc morphology and high frequency spectral structure of $\gamma$-ray blazars. The magnetic field strength can be estimated at different locations along the jet under the assumption of equipartition between magnetic field and relativistic particle energies. Making assumptions on the jet magnetic field configuration (e.g. poloidal or toroidal), we can estimate the separation of the mm-wave "core" and the jet base, and estimate the strength of the magnetic field there. The results of this analysis show that on average, the magnetic field strength decreases with a power-law $B \propto r^{-n}$, $n=0.3 \pm 0.2$. This suggests that on average, the mm-wave "core" is $\sim 1-3$\,pc downstream of the de-projected jet apex and that the magnetic field strength is of the order $B_{\rm{apex}} \sim 5-20$\,kG, broadly consistent with the predictions of magnetic jet launching (e.g. via magnetically arrested disks (MAD)).

Blazars are among the most variable objects in the universe. They feature energetic jets of plasma that launch from the cores of these active galactic nuclei (AGN), triggering activity from radio up to gamma-ray energies. Spatial localization of the region of their MeV/GeV emission is a key question in understanding the blazar phenomenon. The flat spectrum radio quasar (FSRQ) PKS 1502+106 has exhibited extreme and correlated, radio and high-energy activity that triggered intense monitoring by the Fermi-GST AGN Multi-frequency Monitoring Alliance (F-GAMMA) program and the Global Millimeter VLBI Array (GMVA) down to $\lambda$3 mm (or 86 GHz), enabling the sharpest view to date towards this extreme object. Here, we report on preliminary results of our study of the gamma-ray loud blazar PKS 1502+106, combining VLBI and single dish data. We deduce the critical aspect angle towards the source to be $\theta_{\rm c} = 2.6^{\circ}$, calculate the apparent and intrinsic opening angles and constrain the distance of the 86 GHz core from the base of the conical jet, directly from mm-VLBI but also through a single dish relative timing analysis. Finally, we conclude that gamma rays from PKS 1502+106 originate from a region between ~1-16 pc away from the base of the hypothesized conical jet, well beyond the bulk of broad-line region (BLR) material of the source.

The radio galaxy 3C84 is a representative of gamma-ray-bright misaligned active galactic nuclei (AGN) and one of the best laboratories to study the radio properties of subparsec scale jets. We discuss here the past and present activity of the nuclear region within the central 1pc and the properties of subparsec-sized components C1, C2 and C3. We compare these results with the high resolution space-VLBI image at 5GHz obtained with the RadioAstron satellite and we shortly discuss the possible correlation of radio emission with the gamma-ray emission.

We present a high-frequency very long baseline interferometry (VLBI) kinematical study of the BL Lac object S5 0716+714 over the time period of September 2008 to October 2010. The aim of the study is to investigate the relation of the jet kinematics to the observed broadband flux variability. We find significant non-radial motions in the jet outflow of the source. In the radial direction, the highest measured apparent speed is \sim37 c, which is exceptionally high, especially for a BL Lac object. Patterns in the jet flow reveal a roughly stationary feature \sim0.15 mas downstream of the core. The long-term fits to the component trajectories reveal acceleration in the sub-mas region of the jet. The measured brightness temperature, T_{B}, follows a continuous trend of decline with distance, T_B \propto r_{jet}^{-(2.36\pm0.41)}, which suggests a gradient in Doppler factor along the jet axis. Our analysis suggest that a moving disturbance (or a shock wave) from the base of the jet produces the high-energy (optical to \gamma-ray) variations upstream of the 7 mm core, and then later causes an outburst in the core. Repetitive optical/\gamma-ray flares and the curved trajectories of the associated components suggest that the shock front propagates along a bent trajectory or helical path. Sharper \gamma-ray flares could be related to the passage of moving disturbances through the stationary feature. Our analysis suggests that the \gamma-ray and radio emission regions have different Doppler factors.

The analysis of $\gamma$-ray flux variability along with the parsec-scale jet kinematics suggests that the high-energy radiation in the BL Lac object S5 0716+714 has a significant correlation with the mm-VLBI core flux density and with the local orientation of the inner jet flow. For the first time in any blazar, we report a significant correlation between the $\gamma$-ray flux variations and the variations in the local orientation of the jet outflow (position angle). We find that the $\gamma$-ray flux variations lead the 7~mm VLBI core flux variations by 82$\pm$32~days, which suggests that the high-energy emission in S5 0716+714 is coming from a region located 3.8$\pm$1.9~parsecs closer to the central black hole than the "core" seen on the mm-VLBI images. The results imply a strong physical and casual connection between $\gamma$-ray emission and the inner jet morphology in the source.

We have used the broadband backend available at the ATCA to study the fast interstellar scintillation of quasar PKS 1257-326, resolving the core shift as a function of frequency on scales less than 10 microarcseconds. In this short paper we discuss the jet direction implied from the microarcsecond-scale core shift in PKS 1257-326.

Using millimeter-very long baseline interferometry (VLBI) observations of the BL Lac object S5 0716+714 from August 2008 to September 2013, we investigate variations in the core flux density and orientation of the sub-parsec scale jet i.e. position angle. The gamma-ray data obtained by the Fermi-LAT (Large Area Telescope) are used to investigate the high-energy flux variations over the same time period. For the first time in any blazar, we report a significant correlation between the gamma-ray flux variations and the position angle (PA) variations in the VLBI jet. The cross-correlation analysis also indicates a positive correlation such that the mm-VLBI core flux density variations are delayed with respect to the gamma-ray flux by 82$\pm$32 days. This suggests that the high-energy emission is coming from a region located $\geq$(3.8$\pm$1.9) parsecs upstream of the mm-VLBI core (closer to the central black hole). These results imply that the observed inner jet morphology has a strong connection with the observed gamma-ray flares.

The Global mm-VLBI Array (GMVA) is the highest angular resolution imaging interferometer currently available as a common user facility. It is capable of angular resolutions on the order of 40 microarcseconds. Currently 14 stations in the United States and Europe participate in global 3 mm VLBI observations. The GMVA is used for continuum and spectroscopic imaging, probing the central regions of active galaxies and the origin of jets as these regions are typically not observable at longer wavelengths due to synchrotron self-absorption. In early 2012, fringes were detected to the three stations of the Korean VLBI Network (KVN), opening the possibility of extending the baseline coverage of the VLBI array to the East. In these proceedings, we will present recent images from a monitoring program of gamma-ray blazars using the GMVA, including the sources 3C454.3 and 0235+164, and an update of its current status and abilities.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

Global VLBI imaging at millimeter and sub-millimeter wavelength overcomes the opacity barrier of synchrotron self-absorption in AGN and opens the direct view into sub-pc scale regions not accessible before. Since AGN variability is more pronounced at short millimeter wavelength, mm-VLBI can reveal structural changes in very early stages after outbursts. When combined with observations at longer wavelength, global 3mm and 1mm VLBI adds very detailed information. This helps to determine fundamental physical properties at the jet base, and in the vicinity of super-massive black holes at the center of AGN. Here we present new results from multi-frequency mm-VLBI imaging of OJ287 during a major outburst. We also report on a successful 1.3mm VLBI experiment with the APEX telescope in Chile. This observation sets a new record in angular resolution. It also opens the path towards future mm-VLBI with ALMA, which aims at the mapping of the black hole event horizon in nearby galaxies, and the study of the roots of jets in AGN.

We describe a new tool for studying the structure and physical characteristics of ultracompact AGN jets and their surroundings with microarcsecond precision. This tool is based on the frequency dependence of the light curves observed for intra-day variable radio sources, where the variability is caused by interstellar scintillation. We apply this method to PKS1257-326 to resolve the core-shift as a function of frequency on scales well below ~12 microarcseconds. We find that the frequency dependence of the position of the scintillating component is r \propto \nu^{-0.1 \pm 0.24} (99% confidence interval) and the frequency dependence of the size of the scintillating component is d \propto \nu^{-0.64 \pm 0.006}. Together, these results imply that the jet opening angle increases with distance along the jet: d \propto r^{n_d}$ with n_d > 1.8. We show that the flaring of the jet, and flat frequency dependence of the core position is broadly consistent with a model in which the jet is hydrostatically confined and traversing a steep pressure gradient in the confining medium with p \propto r^{-n_p} and n_p > 7. Such steep pressure gradients have previously been suggested based on VLBI studies of the frequency dependent core shifts in AGN.