Astrometry, the measurement of positions and motions of the stars, is one of the oldest disciplines in Astronomy, extending back at least as far as Hipparchus' discovery of the precession of Earth's axes in 190 BCE by comparing his catalog with those of his predecessors. Astrometry is fundamental to Astronomy, and critical to many aspects of Astrophysics and Geodesy. In order to understand our planet's and solar system's context within their surroundings, we must be able to to define, quantify, study, refine, and maintain an inertial frame of reference relative to which all positions and motions can be unambiguously and self-consistently described. It is only by using this inertial reference frame that we are able to disentangle our observations of the motions of celestial objects from our own complex path around our star, and its path through the galaxy, and the local group. Every aspect of each area outlined in the call for scientific frontiers in astronomy in the era of the 2020-2030 timeframe will depend on the quality of the inertial reference frame. In this white paper, we propose support for development of radio Very Long Baseline Interferometry (VLBI) capabilities, including the Next Generation Very Large Array (ngVLA), a radio astronomy observatory that will not only support development of a next generation reference frame of unprecedented accuracy, but that will also serve as a highly capable astronomical instrument in its own right. Much like its predecessors, the Very Long Baseline Array (VLBA) and other VLBI telescopes, the proposed ngVLA will provide the foundation for the next three decades for the fundamental reference frame, benefitting astronomy, astrophysics, and geodesy alike.

We investigate a sample of 3412 {\it International Celestial Reference Frame} (ICRF3) extragalactic radio-loud sources with accurate positions determined by VLBI in the S/X band, mostly active galactic nuclei (AGN) and quasars, which are cross-matched with optical sources in the second Gaia data release (Gaia DR2). The main goal of this study is to determine a core sample of astrometric objects that define the mutual orientation of the two fundamental reference frames, the Gaia (optical) and the ICRF3 (radio) frames. The distribution of normalized offsets between the VLBI sources and their optical counterparts is non-Rayleigh, with a deficit around the modal value and a tail extending beyond the 3-sigma confidence level. A few filters are applied to the sample in order to discard double cross-matches, confusion sources, and Gaia astrometric solutions of doubtful quality. {\it Panoramic Survey Telescope and Rapid Response System} (Pan-STARRS) and {\it Dark Energy Survey} (DES) stacked multi-color images are used to further deselect objects that are less suitable for precision astrometry, such as extended galaxies, double and multiple sources, and obvious misidentifications. After this cleaning, 2642 quasars remain, of which 20\% still have normalized offset magnitudes exceeding 3, or 99\% confidence level. We publish a list of 2118 radio-loud quasars of prime astrometric quality. The observed dependence of binned median offset on redshift shows the expected decline at small redshifts, but also an unexpected rise at $z\sim 1.6$, which may be attributed to the emergence of the C IV emission line in the Gaia's $G$ band. The Gaia DR2 parallax zero-point is found to be color-dependent, suggesting an uncorrected instrumental calibration effect.

The VLBI USNO 2016A (U16A) solution is part of a work-in-progress effort by USNO towards the preparation of the ICRF3. Most of the astrometric improvement with respect to the ICRF2 is due to the re-observation of the VCS sources. Our objective in this paper is to assess U16A's astrometry. A comparison with ICRF2 shows statistically significant offsets of size 0.1 mas between the two solutions. While Gaia DR1 positions are not precise enough to resolve these offsets, they are found to be significantly closer to U16A than ICRF2. In particular, the trend for typically larger errors for Southern sources in VLBI solutions are decreased in U16A. Overall, the VLBI-Gaia offsets are reduced by 21%. The U16A list includes 718 sources not previously included in ICRF2. Twenty of those new sources have statistically significant radio-optical offsets. In two-thirds of the cases, these offsets can be explained from PanSTARRS images.

The goal of this presentation is to report the latest progress in creation of the next generation of VLBI-based International Celestial Reference Frame, ICRF3. Two main directions of ICRF3 development are improvement of the S/X-band frame and extension of the ICRF to higher frequencies. Another important task of this work is the preparation for comparison of ICRF3 with the new generation optical frame GCRF expected by the end of the decade as a result of the Gaia mission.

We present an all-sky sample of ~ 1.4 million AGNs meeting a two color infrared photometric selection criteria for AGNs as applied to sources from the Wide-Field Infrared Survey Explorer final catalog release (AllWISE). We assess the spatial distribution and optical properties of our sample and find that the results are consistent with expectations for AGNs. These sources have a mean density of ~ 38 AGNs per square degree on the sky, and their apparent magnitude distribution peaks at g ~ 20, extending to objects as faint as g ~ 26. We test the AGN selection criteria against a large sample of optically-identified stars and determine the "leakage" (that is, the probability that a star detected in an optical survey will be misidentified as a QSO in our sample) rate to be < 4.0 x 10^-5. We conclude that our sample contains almost no optically-identified stars (< 0.041%), making this sample highly promising for future celestial reference frame work by significantly increasing the number of all-sky, compact extragalactic objects. We further compare our sample to catalogs of known AGNs/QSOs and find a completeness value of > 84% (that is, the probability of correctly identifying a known AGN/QSO is at least 84%) for AGNs brighter than a limiting magnitude of R < 19. Our sample includes approximately 1.1 million previously uncatalogued AGNs.

(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.

Quasometry is precision measurement of celestial positions and apparent motion of very distant extragalactic objects, such as quasars, galactic nuclei, and QSOs. We use this term to identify a specific area of research, the methodology of which differs from that of general astrometry. The main purpose of quasometry is to link the sub-milliarcsecond radio frame (ICRF) with the existing and emerging optical reference frames of similar accuracy, constructed by astrometric satellites. Some of the main difficulties in achieving this goal are discussed, e.g., the extended structures of quasar hosts, apparent motion on the sky, optical variability, galactic companions, faintness. Besides the strategic purpose, quasometry is undoubtedly useful for global astrometric surveys, as it helps to verify or even correct the resulting reference frames. There are two options of using measurements of distant quasars in a global astrometric solution: 1) hard constraints embedded in the fabric of observational equations; 2) {\it a posteriori} fitting of zonal errors. The relative benefits and shortcoming of the two options are reviewed. A relatively small set of about 200 carefully selected reference quasars can go a long way in improving the astrometric value of a space mission, if they are sufficiently bright, stable, fairly uniformly distributed on the sky, and are defining sources in the ICRF. We present an ongoing program at the USNO to construct a quality set of optical quasars with the required properties and to enhance the ICRF with new sources in the areas where known, well-observed quasars are scarce.

We present a detailed multi-epoch analysis of 31 potential southern hemisphere radio calibrators that were originally observed as part of a program to maintain the International Celestial Reference Frame (ICRF). At radio wavelengths, the primary calibrators are Active Galactic Nuclei (AGN), powerful radio emitters which exist at the centre of most galaxies. These are known to vary at all wavelengths at which they have been observed. By determining the amount of radio source structure and variability of these AGN, we determine their suitability as phase calibrators for long baseline radio interferometry at 2.3 GHz. For this purpose, we have used a set of complementary metrics to classify these 31 southern sources into five categories pertaining to their suitability as VLBI calibrators. We find that all of the sources in our sample would be good interferometric calibrators and almost ninety per cent would be very good calibrators.

Four close radio sources in the International Celestial Reference Frame (ICRF) catalog were observed using phase referencing with the VLBA at 43, 23 and 8.6 GHz, and with VERA at 23 GHz over a one year period. The goal was to determine the stability of the radio cores, and to assess structure effects associated with positions in the ICRF. Although the four sources were compact at 8.6 GHz, the VLBA images at 43 GHz with 0.3-mas resolution showed that all were composed of several components. A component in each source was identified as the radio core using some or all of the following emission properties: compactness, spectral index, location at the end of the extended emission region, and stationary in the sky. Over the observing period, the relative positions between the four radio cores were constant to 0.02 mas---the phase referencing positional accuracy obtained at 23 and 43 GHz among the sources---suggesting that once a radio core is identified, it remains stationary in the sky to this accuracy. Other radio components in two of the four sources had detectable motion in the radio jet direction. Comparison of the 23 and 43 GHz VLBA images with the VLBA 8.6 GHz images and the ICRF positions suggests that some ICRF positions are dominated by a moving jet component; hence, they can be displaced up to 0.5 mas from the radio core, and may also reflect the motion of the jet component. Future astrometric efforts to determine a more accurate quasar reference frame at 23 and 43 GHz and from the VLBI2010 project are discussed, and supporting VLBA or EVN observations of ICRF sources at 43 GHz are recommended in order to determine the internal structure of the sources. A future collaboration between the radio (ICRF) and the optical frame of Gaia is discussed.

We introduce the TANAMI program (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) which is monitoring an initial sample of 43 extragalactic jets located south of -30 degrees declination at 8.4 GHz and 22 GHz since 2007. All aspects of the program are discussed. First epoch results at 8.4 GHz are presented along with physical parameters derived therefrom. We present first epoch images for 43 sources, some observed for the first time at milliarcsecond resolution. Parameters of these images as well as physical parameters derived from them are also presented and discussed. These and subsequent images from the TANAMI survey are available at http://pulsar.sternwarte.uni-erlangen.de/tanami/ We obtain reliable, high dynamic range images of the southern hemisphere AGN. All the quasars and BL Lac objects in the sample have a single-sided radio morphology. Galaxies are either double-sided, single-sided or irregular. About 28% of the TANAMI sample has been detected by LAT during its first three months of operations. Initial analysis suggests that when galaxies are excluded, sources detected by LAT have larger opening angles than those not detected by LAT. Brightness temperatures of LAT detections and non-detections seem to have similar distributions. The redshift distributions of the TANAMI sample and sub-samples are similar to those seen for the bright gamma-ray AGN seen by LAT and EGRET but none of the sources with a redshift above 1.8 have been detected by LAT.

We have measured the sub-milli-arcsecond structure of 274 extragalactic sources at 24 and 43 GHz in order to assess their astrometric suitability for use in a high frequency celestial reference frame (CRF). Ten sessions of observations with the Very Long Baseline Array have been conducted over the course of $\sim$5 years, with a total of 1339 images produced for the 274 sources. There are several quantities that can be used to characterize the impact of intrinsic source structure on astrometric observations including the source flux density, the flux density variability, the source structure index, the source compactness, and the compactness variability. A detailed analysis of these imaging quantities shows that (1) our selection of compact sources from 8.4 GHz catalogs yielded sources with flux densities, averaged over the sessions in which each source was observed, of about 1 Jy at both 24 and 43 GHz, (2) on average the source flux densities at 24 GHz varied by 20%-25% relative to their mean values, with variations in the session-to-session flux density scale being less than 10%, (3) sources were found to be more compact with less intrinsic structure at higher frequencies, and (4) variations of the core radio emission relative to the total flux density of the source are less than 8% on average at 24 GHz. We conclude that the reduction in the effects due to source structure gained by observing at higher frequencies will result in an improved CRF and a pool of high-quality fiducial reference points for use in spacecraft navigation over the next decade.

We present a summary of the observation strategy of TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry), a monitoring program to study the parsec-scale structure and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Australian Long Baseline Array (LBA) and the trans-oceanic antennas Hartebeesthoek, TIGO, and O'Higgins. TANAMI is focusing on extragalactic sources south of -30 degrees declination with observations at 8.4 GHz and 22 GHz every ~2 months at milliarcsecond resolution. The initial TANAMI sample of 43 sources has been defined before the launch of the Fermi Gamma Ray Space Telescope to include the most promising candidates for bright gamma-ray emission to be detected with its Large Area Telescope (LAT). Since November 2008, we have been adding new sources to the sample, which now includes all known radio- and gamma-ray bright AGN of the Southern Hemisphere. The combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN and the nature of relativistic jets.

The TANAMI (Tracking AGN with Austral Milliarcsecond Interferometry) and associated programs provide comprehensive radio monitoring of extragalactic gamma-ray sources south of declination -30 degrees. Joint quasi-simultaneous observations between the Fermi Gamma-ray Space Telescope and ground based observatories allow us to discriminate between competing theoretical blazar emission models. High resolution VLBI observations are the only way to spatially resolve the sub-parsec level emission regions where the high-energy radiation originates. The gap from radio to gamma-ray energies is spanned with near simultaneous data from the Swift satellite and ground based optical observatories. We present early results from the TANAMI program in the context of this panchromatic suite of observations.

TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) is a monitoring program to study the parsec-scale structures and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Long Baseline Array and associated telescopes. Extragalactic jets south of -30 degrees declination are observed at 8.4 GHz and 22 GHz every two months at milliarcsecond resolution. The initial TANAMI sample is a hybrid radio and gamma-ray selected sample since the combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN.

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

We analyze a sample of 58 multi-wavelength, Very Long Baseline Array observations of active galactic nuclei (AGN) to determine their scattering properties. Approximately 75% of the sample consists of AGN that exhibit centimeter-wavelength intraday variability (interstellar scintillation) while the other 25% do not show intraday variability. We find that interstellar scattering is measurable for most of these AGN, and the typical broadening diameter is 2 mas at 1 GHz. We find that the scintillating AGN are typically at lower Galactic latitudes than the non-scintillating AGN, consistent with the scenario that intraday variability is a propagation effect from the Galactic interstellar medium. The magnitude of the inferred interstellar broadening measured toward the scintillating AGN, when scaled to higher frequencies, is comparable to the diameters inferred from analyses of the light curves for the more well-known intraday variable sources. However, we find no difference in the amount of scattering measured toward the scintillating versus non-scintillating AGN. A consistent picture is one in which the scintillation results from localized regions ("clumps") distributed throughout the Galactic disk, but which individually make little contribution to the angular broadening. Of the 58 AGN observed, 37 (64%) have measured redshifts. At best, a marginal trend is found for scintillating (non-scintillating) AGN to have smaller (larger) angular diameters at higher redshifts. We also use our observations to try to constrain the possibility of intergalactic scattering. While broadly consistent with the scenario of a highly turbulent intergalactic medium, our observations do not place significant constraints on its properties.

(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 used the Very Large Array (VLA), linked with the Pie Town Very Long Baseline Array antenna, to determine astrometric positions of 46 radio stars in the International Celestial Reference Frame (ICRF). Positions were obtained in the ICRF directly through phase referencing of the stars to nearby ICRF quasars whose positions are accurate at the 0.25 mas level. Radio star positions are estimated to be accurate at the 10 mas level, with position errors approaching a few milli-arcseconds for some of the stars observed. Our measured positions were combined with previous measurements taken from as early as 1978 to obtain proper motion estimates for all 46 stars with average uncertainties of ~1.7 mas/yr. We compared our radio star positions and proper motions with the Hipparcos Catalogue data, and find consistency in the reference frames produced by each data set on the 1-sigma level, with errors of ~2.7 mas per axis for the reference frame orientation angles at our mean epoch of 2003.78. No significant spin is found between our radio data frame and the Hipparcos Celestial Reference Frame (HCRF) with largest rotation rates of +0.55 and -0.41 mas/yr around the x and z axes, respectively, with 1-sigma errors of 0.36 mas/yr. Thus, our results are consistent with a non-rotating Hipparcos frame with respect to the ICRF.

The current realization of the International Celestial Reference Frame (ICRF) comprises a total of 717 extragalactic radio sources distributed over the entire sky. An observing program has been developed to densify the ICRF in the northern sky using the European VLBI network (EVN) and other radio telescopes in Spitsbergen, Canada and USA. Altogether, 150 new sources selected from the Jodrell Bank-VLA Astrometric Survey were observed during three such EVN+ experiments conducted in 2000, 2002 and 2003. The sources were selected on the basis of their sky location in order to fill the "empty" regions of the frame. A secondary criterion was based on source compactness to limit structural effects in the astrometric measurements. All 150 new sources have been successfully detected and the precision of the estimated coordinates in right ascension and declination is better than 1 milliarcsecond (mas) for most of them. A comparison with the astrometric positions from the Very Long baseline Array Calibrator Survey for 129 common sources indicates agreement within 2 mas for 80% of the sources.

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.

We have used the Very Large Array, linked with the Pie Town Very Long Baseline Array antenna, to determine astrometric positions of 19 radio stars in the International Celestial Reference Frame (ICRF). The positions of these stars were directly linked to the positions of distant quasars through phase referencing observations. The positions of the ICRF quasars are known to 0.25 mas, thus providing an absolute reference at the angular resolution of our radio observations. Average values for the errors in our derived positions for all sources were 13 mas and 16 mas in R.A. and declination respectively, with accuracies approaching 1-2 mas for some of the stars observed. Differences between the ICRF positions of the 38 quasars, and those measured from our observations showed no systematic offsets, with mean values of -0.3 mas in R.A. and -1.0 mas in declination. Standard deviations of the quasar position differences of 17 mas and 11 mas in R.A. and declination respectively, are consistent with the mean position errors determined for the stars. Our measured positions were combined with previous Very Large Array measurements taken from 1978-1995 to determine the proper motions of 15 of the stars in our list. With mean errors of approximately 1.6 mas/yr, the accuracies of our proper motions approach those derived from Hipparcos, and for a few of the stars in our program, are better than the Hipparcos values. Comparing the positions of our radio stars with the Hipparcos catalog, we find that at the epoch of our observations, the two frames are aligned to within formal errors of approximately 3 mas. This result confirms that the Hipparcos frame is inertial at the expected level.

We present multi-frequency observations with the Very Long Baseline Array of the compact double radio source B2 2050+36. Our observations are at 0.33, 0.61, 1.67, 2.3, and 8.4 GHz, with the 0.61 GHz observations forming the third epoch of observation of this source at that frequency. At 0.61 GHz, the structure of B2 2050+36 is dominated by two components 56 mas apart. Within the uncertainties of the various measurements, this separation has remained unchanged for the past 16 years. Any differential image wander caused by refractive interstellar scattering is less than 4 mas. Both the lack of differential image wander and the frequency dependence of the angular diameter of B2 2050+36 below 1 GHz indicate that the electron density power spectrum along this line of sight has a spectral index near the Kolmogorov value, with a value of 4 being highly unlikely. We conclude that diffractive scattering dominates along this line of sight; results in the literature indicate that this conclusion also holds true for the line of sight to the pulsar PSR B2020+28 (8.7 deg. from B2 2050+36). Comparison of our 1.67 GHz observations with those obtained 21 years previously place a limit on the projected linear separation velocity of the two components of c.

We review observational constraints on the structures responsible for extreme scattering events, focussing on a series of observations of the quasar PKS 1741-038. VLA observations were conducted to search for changes in the rotation measure and H I absorption during the ESE, while VLBI observations sought ESE-induced changes in the source's image. No RM changes were found implying B_|| < 12 mG, and no H I opacity changes were found implying N(H I) < 6.4 x 10^{17} cm^{-2}. No multiple imaging was observed, but the diameter of the source increased by 0.7 mas, contrary to what is predicted by simple refractive lens modeling of ESEs. We summarize what these limits imply about the structure responsible for this ESE.

(Abridged) We report multi-epoch VLBI observations of the source PKS 1741-038 as it underwent an extreme scattering event. Observations at four epochs were obtained, and images were produced at three of these. During the event the source consisted of a dominant, compact component, essentially identical to the structure seen outside the event. However, the source's diameter increased slightly at 13 cm during the ESE. An increase in the source's diameter is inconsistent with a simple refractive model. We also see no evidence for ESE-induced substructure within the source or the formation of multiple images, as would occur in a strongly refractive lens. However, a model in which the decrease in flux density during the ESE occurs solely because of stochastic broadening within the lens requires a larger broadening diameter during the ESE than is observed. Thus, the ESE toward 1741-038 involved both stochastic broadening and refractive defocussing within the lens. If the structure responsible for the ESE has a size of order 1 AU, the level of scattering within an ESE lens may be a factor of 10^7 larger than that in the ambient medium. A filamentary structure could reduce the difference between the strength of scattering in the lens and ambient medium, but we conclude that, if ESEs arise from filamentary structures, they occur when the filamentary structures are seen lengthwise. We predict the amount of pulse broadening that would result from a comparable lens passing in front of a pulsar. The pulse broadening would be no more than 1.1 microseconds, consistent with the lack of pulse broadening detected during ESEs toward the pulsars PSR B1937+21 and PSR J1643-1224.

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 consider the geometrical optics for refraction of a distant radio source by an interstellar plasma lens, with application to a lens with a Gaussian electron column density profile. The refractive properties of the lens are specified completely by a dimensionless parameter, alpha, which is a function of the wavelength of observation, the lens' electron column density, the lens-observer distance, and the transverse diameter of the lens. Relative motion of the observer and lens produces modulations in the source's light curve. Plasma lenses are diverging so the light curve displays a minimum, when the lens is on-axis, surrounded by enhancements above the unlensed flux density. Lensing can also produce caustics, multiple imaging, and angular position wander of the background source. If caustics are formed, the separation of the outer caustics can constrain alpha, while the separation of the inner caustics can constrain the size of the lens. We apply our analysis to 0954+654, a source for which we can identify caustics in its light curve, and 1741-038, for which polarization observations were obtained during and after the scattering event. We find general agreement between modelled and observed light curves at 2.25 GHz, but poor agreement at 8.1 GHz. The discrepancies may result from a combination of lens substructure or anisotropic shape, a lens that only grazes the source, or unresolved source substructure. Our analysis places the following constraints on the lenses: Toward 0954+654 (1741-038) the lens was 0.38 AU (0.065 AU) in diameter, with a peak column density of 0.24 pc cm^{-3} (1E-4 pc cm^{-3}) and an electron density of 1E5 cm^{-3} (300 cm^{-3}). The angular wander caused by the lens was 250 mas (0.4 mas) at 2.25 GHz. For 1741-038, we place an upper limit of 100 mG on the lens' magnetic field.

The radio observations are based on dual frequency Mark--III VLBI data from 1979 until the end of 1993. Currently 436 sources define the celestial inertial frame. A comparison to the JPL94R01 catalog is made. Position vs. time plots of selected sources show significant apparent non linear proper motions. The status of the project and recent results from optical observations are presented.