Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed.

PG 1553+113 is the first blazar showing an approximately two-year quasi-periodic pattern in its gamma-ray light curve. Such quasi-periodicity might have a geometrical origin, possibly related to the precessing nature of the jet, or could be intrinsic to the source and related to pulsational accretion flow instabilities. By means of a ~2yr very long baseline array (VLBA) monitoring at 15, 24, and 43 GHz we investigate the source pc-scale properties during an entire cycle of gamma-ray activity in the period 2015-2017. In contrast to the well-defined periodicity in the gamma-ray emission, at radio frequencies no clear periodic pattern can be recognized. The jet position angle, constrained by means of the total intensity ridge line, varies across the different observing epochs in the range 40-60 deg. We also investigate the time evolution of the source polarization properties, including the rotation measure. The brightness temperature is found to decrease as the frequency increases with an intrinsic value of ~1.5 x 10^10 K and the estimated Doppler factor is ~1.4.

We present 88 multi-epoch Very Long Baseline Array (VLBA) images (most at an observing frequency of 8 GHz), of 20 TeV blazars, all of the HBL class, that have not been previously studied at multiple epochs on the parsec scale. From these 20 sources, we analyze the apparent speeds of 43 jet components that are all detected at four or more epochs. As has been found for other TeV HBLs, the apparent speeds of these components are relatively slow. About two-thirds of the components have an apparent speed that is consistent (within two-sigma) with no motion; and some of these components may be stationary patterns whose apparent speed does not relate to the underlying bulk flow speed. In addition, a superluminal tail to the apparent speed distribution of the TeV HBLs is detected for the first time, with eight components in seven sources having a two-sigma lower limit on the apparent speed exceeding 1c. We combine the data from these 20 sources with an additional 18 sources from the literature to analyze the complete apparent speed distribution of all 38 TeV HBLs that have been studied with VLBI at multiple epochs. The highest two-sigma apparent speed lower limit considering all sources is 3.6c. This suggests that bulk Lorentz factors of up to about 4, but probably not much higher, exist in the parsec-scale radio emitting regions of these sources, consistent with estimates obtained in the radio by other means such as brightness temperatures. This can be reconciled with the high Lorentz factors estimated from the high-energy data if the jet has velocity structures consisting of different emission regions with different Lorentz factors. In particular, we analyze the current apparent speed data for the TeV HBLs in the context of a model with a fast central spine and a slower outer layer.

We present new multi-epoch Very Long Baseline Array (VLBA) observations of a set of TeV blazars drawn from our VLBA program to monitor all TeV-detected high-frequency peaked BL Lac objects (HBLs) at parsec scales. Most of these sources are faint in the radio, so they have not been well observed with VLBI by other surveys. Our previous measurements of apparent jet speeds in TeV HBLs showed apparent jet speeds that were subluminal or barely superluminal, suggesting jets with velocity structures at the parsec-scale. Here we present apparent jet speed measurements for eight new TeV HBLs, which for the first time show a superluminal tail to the apparent speed distribution for the TeV HBLs.

In this work we present a polarization analysis at radio frequencies of Markarian 421, one of the closest (z=0.03) TeV blazars. The observations were obtained, both in total and in polarized intensity, with the Very Long Baseline Array (VLBA) at 15, 24, and 43 GHz throughout 2011, with one observation per month (for a total of twelve epochs). We investigate the magnetic field topology and the polarization structure on parsec scale and their evolution with time. We detect polarized emission both in the core and in the jet region, and it varies with frequency, location and time. In the core region we measure a mean fractional polarization of about 1-2%, with a peak of about 4% in March at 43 GHz; the polarization angle is almost stable at 43 GHz, but it shows significant variability in the range 114-173 deg at 15 GHz. In the jet region the polarization properties show a more stable behavior; the fractional polarization is about 16% and the polarization angle is nearly perpendicular to the jet axis. The higher EVPA variability observed at 15 GHz is due both to a variable Faraday rotation effect and to opacity. The residual variability observed in the intrinsic polarization angle, together with the low degree of polarization in the core region, could be explained with the presence of a blend of variable cross-polarized subcomponents within the beam.

We present Very Long Baseline Array (VLBA) images of 20 TeV blazars (HBLs) not previously well-studied on the parsec scale. Observations were made between August and December 2013, at a frequency of 8.4 GHz. These observations represent the first epoch of a VLBA monitoring campaign on these blazars, and they significantly increase the fraction of TeV HBLs studied with high-resolution imaging. The peak VLBI flux densities of these sources range from about 10 to 100 mJy/beam, and parsec-scale jet structure is detected in all sources. About half of the VLBI cores are resolved, with brightness temperature upper limits of a few times 10^10 K, and we find that a brightness temperature of about 2x10^10 K is consistent with the VLBI data for all but one of the sources. Such brightness temperatures do not require any relativistic beaming to reduce the observed value below commonly invoked intrinsic limits; however, the lack of detection of counter-jets does place a modest limit on the bulk Lorentz factor of greater than about 2. These data are thus consistent with a picture where weak-jet sources like the TeV HBLs develop significant velocity structures on parsec-scales. We also extend consideration to the full sample of TeV HBLs by combining the new VLBI data with VLBI and gamma-ray data from the literature. By comparing measured VLBI and TeV fluxes to samples with intrinsically uncorrelated luminosities generated by Monte Carlo simulations, we find a marginally significant correlation between the VLBI and TeV fluxes for the full TeV HBL sample.

We report the results obtained for the AGN Markarian 421 by model-fitting the data in the visibility plane, studing the proper motion of jet components, the light curve, and the spectral index of the jet features. We compare the radio data with optical light curves obtained at the Steward Observatory, considering also the optical polarization information. Mrk 421 has a bright nucleus and a one-sided jet extending towards the north-west for a few parsecs. The model-fits show that brightness distribution is well described using 6-7 circular Gaussian components, four of which are reliably identified at all epochs; all components are effectively stationary except for component D, at ~0.4 mas from the core, whose motion is however subluminal. Analysis of the light curve shows two different states, with the source being brighter and more variable in the first half of 2011 than in the second half. The highest flux density is reached in February. A comparison with the optical data reveals an increase of the V magnitude and of the fractional polarization simultaneous with the enhancement of the radio activity.

We expand our previous studies of the parsec-scale structure of TeV blazars by presenting first-epoch images from VLBA monitoring of ten newer TeV HBLs. All ten sources were successfully detected and imaged, and all showed a one-sided core-jet structure. Many display a morphology common to TeV HBLs: a short, collimated jet followed by a transition to low surface brightness extended emission with a much broader opening angle. The newly detected TeV HBLs tend to be fainter in the radio; the median core flux density was 22 mJy, and the median brightness temperature was 8x10^9 K. The brightness temperatures are well below the equipartition limit, and thus the VLBI cores do not require strong beaming, consistent with the modest values of Doppler and Lorentz factors found in the VLBI jets of TeV HBLs by other studies, and contrasting with the strong beaming generally required by the TeV emission. We study the full sample of TeV HBLs that have been observed with VLBI, and find a correlation between TeV flux and VLBI core brightness temperature, suggesting different but correlated beaming factors for the TeV and radio emission. We present a discussion of these observations in the context of velocity structures in the jets of the TeV HBLs.

(Abridged) We analyze blazar jet apparent speeds and accelerations from the RDV series of astrometric and geodetic VLBI experiments. From these experiments, we have produced and analyzed 2753 global VLBI images of 68 sources at 8 GHz with a median beam size of 0.9 milliarcseconds (mas), and a median of 43 epochs per source. From this sample, we analyze the motions of 225 jet components in 66 sources. The distribution of the fastest measured apparent speed in each source has a median of 8.3c and a maximum of 44c. Sources in the 2FGL Fermi LAT catalog display higher apparent speeds than those that have not been detected. On average, components farther from the core in a given source have significantly higher apparent speeds than components closer to the core. We measure accelerations of components in orthogonal directions parallel and perpendicular to their average velocity vector. Parallel accelerations have significantly larger magnitudes than perpendicular accelerations, implying observed accelerations are predominantly due to changes in the Lorentz factor (bulk or pattern) rather than projection effects from jet bending. Positive parallel accelerations are significantly more common than negative ones, so the Lorentz factor (bulk or pattern)tends to increase on the scales observed here. Observed parallel accelerations correspond to modest source frame increases in the bulk or pattern Lorentz factor.

We report on our observations of the parsec-scale radio jet structures of five blazars that have been detected by ground-based TeV gamma-ray telescopes. These five blazars all belong to the class of high-frequency peaked BL Lac objects (HBLs), which are the most common blazar type detected at the TeV energy range. Because of their relative faintness in the radio, these HBLs are not well represented in other radio blazar surveys. Our observations consist of five epochs of Very Long Baseline Array (VLBA) imaging from 2006 to 2009, of each of the five blazars 1ES 1101-232, Markarian 180, 1ES 1218+304, PG 1553+113, and H 2356-309, at frequencies from 5 to 22 GHz. Fundamental jet properties, including the apparent jet speeds, that can be measured from these multi-epoch series of VLBA images are presented and compared with other gamma-ray blazars. Confirming prior work, we find that the TeV HBLs have significantly slower apparent jet speeds than radio-selected blazars. Together with other radio properties of the HBL class, this suggests modest Lorentz factors in their parsec-scale radio jets. This in turn suggests some form of Lorentz factor gradient in these jets, since they are likely to have high Lorentz factors in their TeV-emitting regions. The study presented here approximately doubles the number of TeV HBLs with multi-epoch parsec-scale kinematic measurements.

We present the results from an ultra-high-resolution 7mm Very Long Baseline Array (VLBA) study of the relativistic jet in the BL Lacertae object OJ287 from 1995 to 2011 containing 136 total intensity images. Analysis of the image sequence reveals a sharp jet-position-angle swing by >100 deg. during [2004,2006], as viewed in the plane of the sky, that we interpret as the crossing of the jet from one side of the line of sight to the other during a softer and longer term swing of the inner jet. Modulating such long term swing, our images also show for the first time a prominent erratic wobbling behavior of the innermost ~0.4mas of the jet with fluctuations in position angle of up to ~40 deg. over time scales ~2yr. This is accompanied by highly superluminal motions along non-radial trajectories, which reflect the remarkable non-ballistic nature of the jet plasma on these scales. The erratic nature and short time scales of the observed behavior rules out scenarios such as binary black hole systems, accretion disk precession, and interaction with the ambient medium as possible origins of the phenomenon on the scales probed by our observations, although such processes may cause longer-term modulation of the jet direction. We propose that variable asymmetric injection of the jet flow; perhaps related to turbulence in the accretion disk; coupled with hydrodynamic instabilities, leads to the non-ballistic dynamics that cause the observed non-periodic changes in the direction of the inner jet.

We present 23 new VLBA images of the six established TeV blazars Markarian 421, Markarian 501, H 1426+428, 1ES 1959+650, PKS 2155-304, and 1ES 2344+514, obtained from 2005 to 2009. Most images were obtained at 43 GHz, and they reveal the parsec-scale structures of three of these sources (1ES 1959+650, PKS 2155-304, and 1ES 2344+514) at factors of two to three higher resolution than has previously been attained. Most of the remaining images map the linear polarization structures at a lower frequency of 22 GHz. We discuss the transverse structures of the jets as revealed by the high-frequency and polarimetric imaging. The transverse structures include significant limb-brightening in Mrk 421, and spine-sheath structures in the electric vector position angle (EVPA) and fractional polarization distributions in Mrk 421, Mrk 501, and 1ES 1959+650. We use new measured component positions to update measured apparent jet speeds, in many cases significantly reducing the statistical error over previously published results. With the increased resolution at 43 GHz, we detect new components within 0.1-0.2 mas of the core in most of these sources. No motion is apparent in these new components over the time span of our observations, and we place upper limits on the apparent speeds of the components near the core of less than 2c. From those limits, we conclude that Gamma2 < Gamma1^{1/2} at about 10^5 Schwarzschild radii, where Gamma1 and Gamma2 are the bulk Lorentz factors in the TeV-emitting and 43 GHz-emitting regions, respectively, assuming that their velocity vectors are aligned.

We present three 43 GHz images and a single 86 GHz image of Markarian 501 from VLBA observations in 2005. The 86 GHz image shows a partially resolved core with a flux density of about 200 mJy and a size of about 300 Schwarzschild radii, similar to recent results by Giroletti et al. Extreme limb-brightening is found in the inner parsec of the jet in the 43 GHz images, providing strong observational support for a `spine-layer' structure at this distance from the core. The jet is well resolved transverse to its axis, allowing Gaussian model components to be fit to each limb of the jet. The spine-layer brightness ratio and relative sizes, the jet opening angle, and a tentative detection of superluminal motion in the layer are all discussed.

We present Very Long Baseline Array (VLBA) observations of the TeV blazars H 1426+428, 1ES 1959+650, and PKS 2155-304 obtained during the years 2001 through 2004. We observed H 1426+428 at four epochs at 8 GHz, and found that its parsec-scale structure consisted of a ~17 mJy core and a single ~3 mJy jet component with an apparent speed of 2.09 +/- 0.53c. The blazar 1ES 1959+650 was observed at three epochs at frequencies of 15 and 22 GHz. Spectral index information from these dual-frequency observations was used to definitively identify the core of the parsec-scale structure. PKS 2155-304 was observed at a single epoch at 15 GHz with dual-circular polarization, and we present the first VLBI polarimetry image of this source. For 1ES 1959+650 and PKS 2155-304, the current observations are combined with the VLBA observations from our earlier paper to yield improved apparent speed measurements for these sources with greatly reduced measurement errors. The new apparent speed measured for component C2 in 1ES 1959+650 is 0.00 +/- 0.04c (stationary), and the new apparent speed measured for component C1 in PKS 2155-304 is 0.93 +/- 0.31c. We combine the new apparent speed measurements from this paper with the apparent speeds measured in TeV blazar jets from our earlier papers to form a current set of apparent speed measurements in TeV HBLs. The mean peak apparent pattern speed in the jets of the TeV HBLs is about 1c. We conclude the paper with a detailed discussion of the interpretation of the collected VLBA data on TeV blazars in the context of current theoretical models for the parsec-scale structure of TeV blazar jets.

In February 1997, the Japanese radio astronomy satellite HALCA was launched to provide the space-bourne element for the VLBI Space Observatory Programme (VSOP) mission. Approximately twenty-five percent of the mission time was dedicated to the VSOP Survey of bright compact Active Galactic Nuclei (AGN) at 5 GHz. This paper, the fifth in the series, presents images and models for the remaining 140 sources not included in Paper III, which contained 102 sources. For most sources, the plots of the uv-coverage, the visibility amplitude versus uv-distance, and the high resolution image are presented. Model fit parameters to the major radio components are determined, and the brightness temperature of the core component for each source is calculated. The brightness temperature distributions for all of the sources in the VSOP AGN survey are discussed.

(Abridged) We present the results of an analysis of relativistic jet apparent speeds from VLBI images in the Radio Reference Frame Image Database (RRFID). The images are snapshot VLBI images at 8 and 2 GHz using the VLBA plus up to ten additional antennas that provide global VLBI coverage. We have analyzed the 8 GHz images from the first five years of the database (1994-1998), for all sources observed at three or more epochs during this time range. This subset comprises 966 images of 87 sources. The sources in this subset have an average of 11 epochs of observation over the years 1994-1998, with the best-observed sources having 19 epochs. About half of the sources in this RRFID kinematic survey have not been previously studied with multi-epoch VLBI observations. We have measured apparent speeds for a total of 184 jet components in 77 sources, of which the best-measured 94 component speeds in 54 sources are used in the final analysis. The apparent speed distribution shows a peak at low apparent speeds (consistent with stationary components), a tail extending out to apparent speeds of about 30c, and a mean apparent speed of 3.6c. A total of 36 of the sources in this paper are also included in the 2 cm VLBA survey by Kellermann et al., with similar angular resolution, sensitivity, and time range. For those sources, we present a detailed component-by-component comparison of the apparent speeds measured by the 2 cm survey and those measured in this paper. Many of the independent apparent speed measurements agree very well, but for approximately 25% of the components we find significant differences in the apparent speeds measured by the two surveys. The leading cause of these discrepancies are differences in how the two surveys have identified jet components from epoch-to-epoch.

We have measured peak apparent speeds of 25.6 +/- 7.0c, 25.6 +/- 4.4c, and 28.2 +/- 6.6c in the jets of 0235+164, 0827+243, and 1406-076, respectively, based on six epochs of high-sensitivity VLBA observations at 22 and 43 GHz during 2002 and 2003 (H_{0}=71 km/(s Mpc), Omega_{m}=0.27, and Omega_{Lambda}=0.73). These blazars had been identified as potentially having apparent speeds exceeding 40c in an earlier VLBA survey of EGRET blazars by Jorstad et al. We therefore confirm (with high confidence in 0827+243, and lower confidence in 0235+164 and 1406-076) the presence of highly relativistic pattern speeds in these three jets, although not at the >40c levels reported by Jorstad et al. The lower limit to the bulk Lorentz factor implied by the observed apparent speeds is Gamma>25-30 in these three sources, if the pattern speeds are equal to or slower than the bulk flow speed.

We present four high dynamic range, dual-circular polarization, Very Long Baseline Array (VLBA) observations at 22 GHz of Markarian 421, taken throughout the year following the source's unprecedented gamma-ray high state in early 2001. These four new VLBA observations are combined with data from our earlier 1999 paper and archival VLBA data-sets that have become available since 1999 to produce a combined 28 epoch VLBA data-set on Mrk 421 spanning the years 1994 to 2002. No new component associated with the 2001 flares was seen on the total intensity images, but the combined data-set allowed precise measurements of the apparent speeds of the existing components. The peak measured apparent speed was for component C5, which has an apparent speed of 0.1 +/- 0.02 c (H_0=71 km s^{-1} Mpc^{-1}, Omega_m=0.27, and Omega_Lambda=0.73). No counterjet is seen with a limit on the jet to counterjet brightness ratio J >~ 100. These observed VLBI properties of Markarian 421 are consistent with a jet with a bulk Lorentz factor of about 2 and an angle to the line-of-sight of about 1 degree, suggesting a jet that decelerates between the gamma-ray producing region and the parsec scale. The VLBI core and inner jet (component C7) have fractional polarizations of about 5%, and an electric vector position angle (EVPA) aligned with the jet axis. Component C5 (at 1.5 mas from the core) has a higher fractional polarization of about 15%, and an EVPA nearly orthogonal to the jet axis. Significant variability is detected in the EVPA of component C6, which at two of the four epochs shows an EVPA aligned with the jet axis, possibly a sign of propagating disturbances that are only visible on the polarization images. If these propagating disturbances are linked to the 2001 gamma-ray high state, then their inferred apparent speed is between 1 and 3 c.

We present current results from our ongoing project to study the parsec-scale relativistic jet kinematics of sources in the U.S. Naval Observatory's Radio Reference Frame Image Database (RRFID). The RRFID consists of snapshot observations using the VLBA plus up to 9 additional antennas at 8 and 2 GHz. The Image Database currently contains about 3000 images of 450 sources from 1994 to 2004, with some sources having images at 20 epochs or more. We have now completed analysis of the 8 GHz images for all sources observed at 3 or more epochs from 1994 to 1998. The completed analysis comprises 966 images of 87 sources, or an average of 11 epochs per source. Apparent jet speeds have been measured for these sources, and the resulting speed distribution has been compared with results obtained by other large VLBI surveys. The measured apparent speed distribution agrees with those found by the 2 cm survey and Caltech-Jodrell Bank (CJ) survey; however, when a source-by-source comparison is done with the 2 cm survey results, significant disagreement is found in the apparent speed measurements for a number of sources. This disagreement can be traced in most cases to either an insufficient time baseline for the current RRFID results, or to apparent component mis-identification in the 2 cm survey results caused by insufficient time sampling. These results emphasize the need for long time baselines and dense time sampling for multi-epoch monitoring of relativistic jets.

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsecond angular resolution using an orbiting 8-m telescope, HALCA and global arrays of Earth-based telescopes. Approximately 25% of the observing time has been devoted to a survey of compact AGN at 5 GHz which are stronger than 1 Jy -- the VSOP AGN Survey. This paper, the second in a series, describes the data calibration, source detection, self-calibration, imaging and modeling, and gives examples illustrating the problems specific to space VLBI. The VSOP Survey web-site which contains all results and calibrated data is described.

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsec resolution using an orbiting 8 m telescope, HALCA, along with global arrays of Earth-based telescopes. Approximately 25% of the observing time is devoted to a survey of compact AGN which are stronger than 1 Jy at 5 GHz-the VSOP AGN Survey. This paper, the third in the series, presents the results from the analysis of the first 102 Survey sources. We present high resolution images and plots of visibility amplitude versus projected baseline length. In addition, model-fit parameters to the primary radio components are listed, and from these the angular size and brightness temperature for the radio cores are calculated. For those sources for which we were able to determine the source frame core brightness temperature, a significant fraction (53 out of 98) have a source frame core brightness temperature in excess of 10^12 K. The maximum source frame core brightness temperature we observed was 1.2 X 10^13 K. Explaining a brightness temperature this high requires an extreme amount of relativistic Doppler beaming. Since the maximum brightness temperature one is able to determine using only ground-based arrays is of the order of 10^12 K, our results confirm the necessity of using space VLBI to explore the extremely high brightness temperature regime.

This review paper discusses the past and present observations of gamma-ray blazars with both space and ground-based gamma-ray telescopes (such as EGRET and current TeV telescopes), and with the VLBA. I also discuss the more sensitive gamma-ray telescopes such as AGILE and GLAST that will become available over the next several years, and how observations with these telescopes, together with improved VLBA observations, can increase our understanding of blazar physics.

As part of our study of the VLBI properties of TeV blazar jets, we present here a series of high-resolution 15 GHz Very Long Baseline Array (VLBA) images of the parsec-scale jets of the TeV blazars 1ES 1959+650, PKS 2155-304, and 1ES 2344+514, with linear resolutions of about 0.5 pc. Each of these sources was observed with the VLBA at three or four epochs during 1999 and 2000. There is a notable lack of any strong moving components on the VLBI images (in contrast to the rapid superluminal motions seen in EGRET blazars), and the structure of the VLBI jet can be modeled either as a series of stationary Gaussian components, or as a smooth power law for two of the sources (PKS 2155-304 and 1ES 2344+514). The low apparent speeds, together with beaming indicators such as the brightness temperature of the VLBI core, imply only modest Doppler boosting of the VLBI radio emission, and only modest bulk Lorentz factors; in contrast to the more extreme values of these parameters invoked to explain the high-energy emission. The fact that no moving shocks or plasmoids are seen on the parsec-scale suggests that the shocks or plasmoids that are assumed to be responsible for the high-energy flares must dissipate before they separate from the core on the VLBI images. This requires the loss of a substantial amount of bulk kinetic energy on parsec scales, and implies a higher efficiency than is typically assumed for internal shock scenarios.

We present the results of a high angular resolution study of the BL Lac object Markarian 501 in the radio band. We consider data taken at 14 different epochs, ranging between 1.6 GHz and 22 GHz in frequency, and including new Space VLBI observations obtained on 2001 March 5 and 6 at 1.6 and 5 GHz. We study the kinematics of the parsec-scale jet and estimate its bulk velocity and orientation with respect to the line of sight. Limb brightened structure in the jet is clearly visible in our data and we discuss its possible origin in terms of velocity gradients in the jet. Quasi-simultaneous multi-wavelength observations allow us to map the spectral index distribution and to compare it to the jet morphology. Finally, we estimate the physical parameters of the parsec-scale jet.

We have calculated inverse-Compton Doppler factors for 3C 279 using the collection of VLBI data recently published by us, and the collection of multiwavelength spectra recently published by Hartman et al. From the Doppler factor and superluminal apparent speed, we then calculate the Lorentz factor and angle to the line-of-sight of the parsec-scale relativistic jet. We model the jet components as homogeneous spheres and the VLBI core as an unresolved inhomogeneous conical jet. The conical-jet model can be made to match both the observed X-ray emission and the VLBI properties of the core with a suitable choice of Doppler factor, implying the core makes a significant contribution to the X-ray emission. The parameters of the conical models indicate the jet is particle dominated at the radii that produce significant emission, and is not in equipartition. At the inner radius of the conical jet the magnetic field is of order 0.1 G and the relativistic-particle number density is of order 10 cm^{-3}. When all components are included in the calculation, then on average the core produces about half of the X-rays, with the other half being split between the long-lived component C4 and the brightest inner-jet component. We calculate an average speed and angle to the line-of-sight for the region of the jet interior to 1 mas of v=0.992c (gamma=8) and 4 degrees, and an average speed and angle to the line-of-sight for C4 (at a distance from the core of 3 mas) of v=0.997c (gamma=13) and 2 degrees. These values imply average Doppler factors of delta=12 for the inner jet, and delta=21 for C4.

We have combined Very Long Baseline Interferometry (VLBI) data from several programs in order to resolve differences in reported parsec-scale jet speeds for the TeV gamma-ray source Mkn501. Data from the VLBA 2cm survey, and 8 and 15 GHz data from the Radio Reference Frame Image Database, have been combined with data from a 5 GHz VSOP Space VLBI observation to determine the apparent motions of jet components in this source. The combined data set consists of 12 observations between April 1995 and July 1999. Four jet components are detected at most epochs, all of which are clearly sub-luminal (i.e. with apparent speeds < c), and two of which appear stationary. The established TeV gamma-ray sources Mkn 501 and Mkn 421 thus both have sub-luminal parsec-scale jets, in contrast to the apparently super-luminal jets of the majority of >100 MeV sources detected by EGRET. No new VLBI component has emerged from the core following the extended TeV high-state in 1997, again in contrast to the general behavior of >100 MeV gamma-ray sources.

NGC4261 (3C270) is an elliptical galaxy containing a 300-parsec scale nuclear disk of gas and dust imaged by the Hubble Space Telescope (HST), around a central supermassive black hole. Previous VLBI observations of NGC4261 revealed a gap in emission in the radio counterjet, presumably due to free-free absorption in the inner parsec of the accretion disk. Here we present three 8 GHz VLBA observations of NGC4261 that allow us to monitor the location and depth of the gap, and check for motions in the jet and counterjet. The separation between the brightest peak and the gap is stable, with an upper limit to its motion of 0.01 c, supporting the interpretation of the gap as absorption by an accretion disk rather than an intrinsic jet feature. These observations span a time of order that required for orbiting material in the disk to transit the counterjet, so we are able to search for density changes (clumps) in the disk by monitoring the optical depth of the gap. The optical depth of the gap is stable to within 20% over 5 years at tau = 1.1 +/- 0.1, corresponding to an electron density in the disk that is constant to within 10%. We measure an apparent speed in the jet of 0.52 +/- 0.07 c. An apparent speed could not be measured for the counterjet due to a lack of identifiable features. From the apparent jet speed and the jet to counterjet brightness ratio, we calculate the viewing angle of the jet to be 63 +/- 3 degrees and its intrinsic speed to be 0.46 +/- 0.02 c. From the inclination and position angles of the parsec-scale radio jet and outer HST disk rotation axis we calculate a difference between the parsec-scale radio jet and outer HST disk rotation axis of 12 +/- 2 degrees.

We present images from a space-VLBI survey using the facilities of the VLBI Space Observatory Programme (VSOP), drawing our sample from the well-studied Pearson-Readhead survey of extragalactic radio sources. Our survey has taken advantage of long space-VLBI baselines and large arrays of ground antennas, such as the Very Long Baseline Array and European VLBI Network, to obtain high resolution images of 27 active galactic nuclei, and to measure the core brightness temperatures of these sources more accurately than is possible from the ground. A detailed analysis of the source properties is given in accompanying papers. We have also performed an extensive series of simulations to investigate the errors in VSOP images caused by the relatively large holes in the (u,v) plane when sources are observed near the orbit normal direction. We find that while the nominal dynamic range (defined as the ratio of map peak to off-source error) often exceeds 1000:1, the true dynamic range (map peak to on-source error) is only about 30:1 for relatively complex core-jet sources. For sources dominated by a strong point source, this value rises to approximately 100:1. We find the true dynamic range to be a relatively weak function of the difference in position angle (PA) between the jet PA and (u,v) coverage major axis PA. For low signal-to-noise regions typically located down the jet away from the core, large errors can occur, causing spurious features in VSOP images that should be interpreted with caution.

The physical conditions in the inner parsec of accretion disks believed to orbit the central black holes in active galactic nuclei can be probed by imaging the absorption (by ionized gas in the disk) of background emission from a radio counterjet. We report high angular resolution VLBI observations of the nearby (about 40 Mpc) radio galaxy NGC 4261 that confirm free-free absorption of radio emission from a counterjet by a geometrically thin, nearly edge-on disk at 1.6, 4.8, and 8.4 GHz. The angular width and depth of the absorption appears to increase with decreasing frequency, as expected. We derive an average electron density of ~10E4 per cc at a disk radius of about 0.2 pc, assuming that the inner disk inclination and opening angles are the same as at larger radii. Pressure balance between the thermal gas and the magnetic field in the disk implies an average field strength of 0.1 milligauss at a radius of 0.2 pc. These are the closest-in free-free absorption measurements to date of the conditions in an extragalactic accretion disk orbiting a black hole with a well-determined mass. If a standard advection-dominated accretion flow exists in the disk center, then the transition between thin and thick disk regions must occur at a radius less than 0.2 pc (4000 Schwarzschild radii).

We present results of long-term high-frequency VLBI monitoring of the relativistic jet in 3C279, consisting of 18 epochs at 22 GHz from 1991 to 1997 and 10 epochs at 43 GHz from 1995 to 1997. Three major results of this study are: apparent speeds measured for six superluminal components range from 4.8 to 7.5 c (H_{0}=70 km s^{-1} Mpc^{-1}, q_{0}=0.1), variations in the total radio flux are due primarily to changes in the VLBI core flux, and the uniform-sphere brightness temperature of the VLBI core is about 1 x 10^{13} K at 22 GHz after 1995, one of the highest direct estimates of a brightness temperature. If the variability brightness temperature measured for 3C279 by Lahteenmaki & Valtaoja is an actual value and not a lower limit, then the rest-frame brightness temperature of 3C279 is quite high and limited by inverse Compton effects rather than equipartition. The parsec-scale morphology of 3C279 consists of a bright, compact VLBI core, a jet component (C4) that moved from about 2 mas to about 3.5 mas from the core during the course of our monitoring, and an inner jet that extends from the core to a stationary component, C5, at about 1 mas from the core. Components in the inner jet are relatively short-lived, and fade by the time they reach about 1 mas from the core. The components have different speeds and position angles from each other, but these differences do not match the differences predicted by the precession model of Abraham & Carrara. Although VLBI components were born about six months prior to each of the two observed gamma-ray high states, the sparseness of the gamma-ray data prevents a statistical analysis of possible correlations.

The quasar PKS 0637-753, the first celestial X-ray target of the Chandra X-ray Observatory, has revealed asymmetric X-ray structure extending from 3 to 12 arcsec west of the quasar, coincident with the inner portion of the jet previously detected in a 4.8 GHz radio image (Tingay et al. 1998). At a redshift of z=0.651, the jet is the largest (~100 kpc) and most luminous (~10^{44.6} ergs/s) of the few so far detected in X-rays. This letter presents a high resolution X-ray image of the jet, from 42 ks of data when PKS 0637-753 was on-axis and ACIS-S was near the optimum focus. For the inner portion of the radio jet, the X-ray morphology closely matches that of new ATCA radio images at 4.8 and 8.6 GHz. Observations of the parsec scale core using the VSOP space VLBI mission show structure aligned with the X-ray jet, placing important constraints on the X-ray source models. HST images show that there are three small knots coincident with the peak radio and X-ray emission. Two of these are resolved, which we use to estimate the sizes of the X-ray and radio knots. The outer portion of the radio jet, and a radio component to the east, show no X-ray emission to a limit of about 100 times lower flux. The X-ray emission is difficult to explain with models that successfully account for extra-nuclear X-ray/radio structures in other active galaxies. We think the most plausible is a synchrotron self-Compton (SSC) model, but this would imply extreme departures from the conventional minimum-energy and/or homogeneity assumptions. We also rule out synchrotron or thermal bremsstrahlung models for the jet X-rays, unless multicomponent or ad hoc geometries are invoked.

The physical conditions in the inner parsec of accretion disks believed to orbit the central black holes in active galactic nuclei can be probed by imaging the absorption of background radio emission by ionized gas in the disk. We report high angular resolution observations of the nearby galaxy NGC 4261 which show evidence for free-free absorption by a thin, nearly edge-on disk at several frequencies. The angular width, and probably the depth, of the absorption appears to increase with decreasing frequency, as expected. Because free-free absorption is much larger at lower frequencies, the longest possible baselines are needed to provide adequate angular resolution; observing at higher frequencies to improve resolution will not help.

We present the first VLBI Space Observatory Programme (VSOP) observations of the gamma-ray blazar 3C 279 at 1.6 and 5 GHz. The combination of the VSOP and VLBA-only images at these two frequencies maps the jet structure on scales from 1 to 100 mas. On small angular scales the structure is dominated by the quasar core and the bright secondary component `C4' located 3 milliarcseconds from the core (at this epoch). On larger angular scales the structure is dominated by a jet extending to the southwest, which at the largest scale seen in these images connects with the smallest scale structure seen in VLA images. We have exploited two of the main strengths of VSOP: the ability to obtain matched-resolution images to ground-based images at higher frequencies and the ability to measure high brightness temperatures. A spectral index map was made by combining the VSOP 1.6 GHz image with a matched-resolution VLBA-only image at 5 GHz from our VSOP observation on the following day. The spectral index map shows the core to have a highly inverted spectrum, with some areas having a spectral index approaching the limiting value for synchrotron self-absorbed radiation of 2.5. Gaussian model fits to the VSOP visibilities revealed high brightness temperatures (>10^{12} K) that are difficult to measure with ground-only arrays. An extensive error analysis was performed on the brightness temperature measurements. Most components did not have measurable brightness temperature upper limits, but lower limits were measured as high as 5x10^{12} K. This lower limit is significantly above both the nominal inverse Compton and equipartition brightness temperature limits. The derived Doppler factor, Lorentz factor, and angle to the line-of-sight in the case of the equipartition limit are at the upper end of the range of expected values for EGRET blazars.

The structure of AGN accretion disks on sub-parsec scales can be probed through free-free absorption of synchrotron emission from the base of symmetric radio jets. We present new VLBA observations of the nearby FR-I radio galaxy NGC 4261 at 22 and 43 GHz, and combine these with previous VLBA observations at 1.6 and 8.4 GHz to map absorption caused by an inner accretion disk. Assuming the disk is geometrically and optically thin and composed of a uniform 10^4 K plasma, the average electron density in the inner 0.1 pc is 10^3 - 10^8 cm^-3. Equating thermal gas pressure and magnetic field strength gives a disk magnetic field of 10^-4 - 10^-2 Gauss at 0.1 pc. The jet opening angle is between 0.3 and 20 degrees during the first 0.2 pc of the jet, and must be less than 5 degrees during the first 0.8 pc. We include an appendix containing expressions for a simple, optically thin, gas pressure dominated accretion disk model which may be applicable to other galaxies in addition to NGC 4261.

We present thirty VLBI images of the TeV blazar Markarian 421 (1101+384) at fifteen epochs spanning the time range from 1994 to 1997, and at six different frequencies from 2.3 to 43 GHz. The imaged observations include a high-resolution 5 GHz VLBI Space Observatory Programme (VSOP) observation with the HALCA satellite on 1997 November 14; full-track VLBA observations from 1994 April, 1996 November, and 1997 May at frequencies between 5 and 43 GHz; six epochs of VLBA snapshot observations at frequencies between 2 and 15 GHz from Radio Reference Frame studies; and five geodetic VLBI observations at 2 and 8 GHz from the archive of the Washington VLBI Correlator Facility located at the U.S. Naval Observatory. The dense time coverage of the images allows us to unambiguously track components in the parsec-scale jet over the observed time range. We measure the speeds of three inner jet components located between 0.5 and 5 mas from the core (0.3 to 3 pc projected linear distance) to be 0.19 +/- 0.27, 0.30 +/- 0.07, and -0.07 +/- 0.07 c (H_{0}=65 km s^{-1} Mpc^{-1}). If the sole 43 GHz image is excluded, all measured speeds are consistent with no motion. These speeds differ from tentative superluminal speeds measured by Zhang & B\aa\aa th from three epochs of data from the early 1980's. Possible interpretations of these subluminal speeds in terms of the high Doppler factor demanded by the TeV variability of this source are discussed.

We present 27 geodetic VLBI maps of OJ 287 obtained from the archive of the Washington correlator. The observations presented here were made between 1990 October and 1996 December. During this period a sequence of six superluminal components has been identified. We measured the proper motion of these components to be approximately 0.5 mas/yr, which is about twice as high as that seen in previous VLBI observations. These results imply a higher component ejection rate than previously observed, in good agreement with the observed occurrences of radio outbursts. We have examined a possible connection between VLBI components and optical flares in the framework of a binary black hole system.

We present VLBI observations of the EGRET quasars 0202+149, CTA 26, and 1606+106, as well as additional analysis of VLBI observations of 1156+295 presented in Piner & Kingham (1997b). We have produced 8 and 2 GHz VLBI images at 11 epochs, 8 epochs, and 12 epochs, spanning the years 1989 to 1996, of 0202+149, CTA 26, and 1606+106 respectively. The VLBI data have been taken from the Washington VLBI correlator's geodetic database. We have measured the apparent velocities of the jet components and find that CTA 26 and 1606+106 are superluminal sources, with average apparent speeds of 8.9 and 2.9 h^{-1}c respectively (H_{0}=100h km s^{-1} Mpc^{-1}, q_{0}=0.5). The components in 0202+149 are stationary, and we identify this source as a compact F double. These sources all have apparently bent jets, and we detected non-radial motion of components in CTA 26 and 1156+295. We have not yet detected any components emerging subsequent to the gamma-ray flares in CTA 26, 1156+295, and 1606+106, and we derive lower limits on the ejection times of any such components. The misalignment angle distribution of the EGRET sources is compared to the distribution for blazars as a whole, and we find that EGRET sources belong preferentially to neither the aligned nor the misaligned population. We also compare the average values for the apparent velocities and Doppler beaming factors for the EGRET and non-EGRET blazars, and find no significant differences. We thus find no indication, within the measurement errors, that EGRET blazars are any more strongly beamed than their counterparts which have not been detected in gamma-rays.

As part of an ongoing effort to observe high energy gamma-ray blazars with VLBI, we have produced 8 and 2 GHz VLBI images, at ten epochs spanning the years 1988 to 1996, of the quasar 1156+295. The VLBI data have been taken from the Washington VLBI correlator's geodetic database. We have detected detected four components and have measured their apparent speeds to be 8.8 +/- 2.3, 5.3 +/- 1.1, 5.5 +/- 0.9, and 3.5 +/- 1.2 h^{-1}c from the outermost component inwards. (H_{0}=100h km/(s Mpc), q_{0}=0.5 throughout paper). These velocities contradict a previously published very high superluminal velocity of 26 h^{-1}c for this source.