Time domain and multi-messenger astrophysics are growing and important modes of observational astronomy that will help define astrophysics in the 2020s. Significant effort is being put into developing the components of a follow-up system for dynamically turning survey alerts into data. This system consists of: 1) brokers that will aggregate, classify, and filter alerts; 2) Target Observation Managers (TOMs) for prioritizing targets and managing observations and data; and 3) observatory interfaces, schedulers, and facilities along with data reduction software and science archives. These efforts need continued community support and funding in order to complete and maintain them. Many of the efforts can be community open-source software projects but they will benefit from the leadership of professional software developers. The coordination should be done by institutions that are involved in the follow-up system such as the national observatories (e.g. LSST/Gemini/NOAO Mid-scale/Community Science and Data Center) or a new MMA institute. These tools will help the community to produce the most science from new facilities and will provide new capabilities for all users of the facilities that adopt them.

We describe present and future capabilities of the Mid-Scale Observatories (MSO) of the new national center merging NOAO, Gemini Observatory and LSST Operations. MSO is comprised of Cerro Tololo Interamerican Observatory (CTIO) and the Kitt Peak National Observatory (KPNO). Telescopes at both sites currently operate on a mix of public and private funding. Recent upgrades have equipped the MSO 4-m class telescopes to perform world-class surveys in diverse areas of astrophysics, from dark energy to exoplanets.

We report the discovery of a new Herbig-Haro jet, HH 1165, in SOAR narrow-band imaging of the vicinity of the sigma Orionis cluster. HH 1165 shows a spectacular extended and collimated spatial structure, with a projected length of 0.26 pc, a bent C-shaped morphology, multiple knots, and fragmented bow-shocks at the apparent ends of the flow. The Halpha image shows a bright halo with a clumpy distribution of material seen around the driving source, and curved reflection nebulosity tracing the outflow cavities. The driving source of HH 1165 is a Class I proto-brown dwarf, Mayrit 1701117 (M1701117), with a total (dust+gas) mass of ~36 MJup and a bolometric luminosity of ~0.1 Lsun. High-resolution VLT/UVES spectra of M1701117 show a wealth of emission lines indicative of strong outflow and accretion activity. SOAR/Goodman low-resolution spectra along the jet axis show an asymmetrical morphology for HH 1165. We find a puzzling picture wherein the north-west part exhibits a classical HH jet running into a pre-dominantly neutral medium, while the southern part resembles an externally irradiated jet. The C-shaped bending in HH 1165 may be produced by the combined effects from the massive stars in the ionization front to the east, the sigma Orionis core to the west, and the close proximity to the B2-type star HR 1950. HH 1165 shows all of the signatures to be considered as a scaled-down version of parsec-length HH jets, and can be termed as the first sub-stellar analog of a protostellar HH jet system.

Previous work has shown the Orion Bar to be an interface between ionized and molecular gas, viewed roughly edge on, which is excited by the light from the Trapezium cluster. Much of the emission from any star-forming region will originate from such interfaces, so the Bar serves as a foundation test of any emission model. Here we combine X-ray, optical, IR and radio data sets to derive emission spectra along the transition from H+ to H0 to H2 regions. We then reproduce the spectra of these layers with a simulation that simultaneously accounts for the detailed microphysics of the gas, the grains, and molecules, especially H2 and CO. The magnetic field, observed to be the dominant pressure in another region of the Orion Nebula, is treated as a free parameter, along with the density of cosmic rays. Our model successfully accounts for the optical, IR and radio observations across the Bar by including a significant magnetic pressure and also heating by an excess density of cosmic rays, which we suggest is due to cosmic rays being trapped in the compressed magnetic field. In the Orion Bar, as we had previously found in M17, momentum carried by radiation and winds from the newly formed stars pushes back and compresses the surrounding gas. There is a rough balance between outward momentum in starlight and the total pressure in atomic and molecular gas surrounding the H+ region. If the gas starts out with a weak magnetic field, the starlight from a newly formed cluster will push back the gas and compress the gas, magnetic field, and cosmic rays until magnetic pressure becomes an important factor.

We propose a practical method to calculate Zernike aberrations from analysis of a single long-exposure defocused stellar image. It consists in fitting the aberration coefficients and seeing blur directly to a realistic image binned into detector pixels. This "donut" method is different from curvature sensing in that it does not make the usual approximation of linearity. We calculate the sensitivity of this technique to detector and photon noise and determine optimal parameters for some representative cases. Aliasing of high-order un-modeled aberrations is evaluated and shown to be similar to a low-order Shack-Hartmann sensor. The method has been tested with real data from the SOAR and Blanco 4m telescopes.

We present a proper motion study of the shock waves within the classic stellar jet HH 47 based on Hubble Space Telescope H-alpha and [S II] images of the region taken over two epochs. Individual knots within the jet and in the bow shock/Mach disk working surface of HH 47A move significantly in the five years that separate the images, and the excellent spatial resolution of HST makes it possible to measure the proper motions with enough precision to easily observe differential motions throughout the flow. The bright portion of the jet emerges at 37.5 +/- 2.5 degrees from the plane of the sky with an average velocity of 300 km/s. Dynamical ages of the shock waves in the jet range from a few decades for knots recently ejected by the source to ~ 1300 years for the faint extended bow shock HH 47D. The jet curves, but motions of knots in the jet are directed radially away from the exciting source, and velocity variability in the flow drives the shock waves that heat the jet internally. The jet orientation angle varies with time by about 15 degrees, and currently points to the northwestern portion of a cavity outlined by a reflection nebula, where a quasi-stationary shock deflects the jet. The major working surface HH 47A is more complex than a simple bow shock/Mach disk, and contains numerous clumps that move relative to one another with velocities of ~ +/- 40 km/s. Small clumps or instabilities affect the Mach disk, and dense clumps may move all the way through the working surface to cause the bumpy morphology seen at the bow shock. A localized area between the bow shock and Mach disk varies significantly between the two sets of images.

We present the detailed construction and analysis of the most complete map to date of the circumstellar environment around SN 1987A, using ground and space-based imaging from the past 16 years. PSF-matched difference-imaging analyses of data from 1988 through 1997 reveal material between 1 and 28 ly from the SN. Careful analyses allows the reconstruction of the probable circumstellar environment, revealing a richly-structured bipolar nebula. An outer, double-lobed ``Peanut,'' which is believed to be the contact discontinuity between red supergiant and main sequence winds, is a prolate shell extending 28 ly along the poles and 11 ly near the equator. Napoleon's Hat, previously believed to be an independent structure, is the waist of this Peanut, which is pinched to a radius of 6 ly. Interior to this is a cylindrical hourglass, 1 ly in radius and 4 ly long, which connects to the Peanut by a thick equatorial disk. The nebulae are inclined 41\degr south and 8\degr east of the line of sight, slightly elliptical in cross section, and marginally offset west of the SN. From the hourglass to the large, bipolar lobes, echo fluxes suggest that the gas density drops from 1--3 cm^{-3} to >0.03 cm^{-3}, while the maximum dust-grain size increases from ~0.2 micron to 2 micron, and the Si:C dust ratio decreases. The nebulae have a total mass of ~1.7 Msun. The geometry of the three rings is studied, suggesting the northern and southern rings are located 1.3 and 1.0 ly from the SN, while the equatorial ring is elliptical (b/a < 0.98), and spatially offset in the same direction as the hourglass.

We summarize the analysis of a uniform set of both previously-known and newly-discovered scattered-light echoes, detected within 30" of SN 1987A in ten years of optical imaging, and with which we have constructed the most complete three-dimensional model of the progenitor's circumstellar environment. Surrounding the SN is a richly-structured bipolar nebula. An outer, double-lobed ``peanut,'' which we believe is the contact discontinuity between the red supergiant and main sequence winds, is a prolate shell extending 28 ly along the poles and 11 ly near the equator. Napoleon's Hat, previously believed to be an independent structure, is the waist of this peanut, which is pinched to a radius of 6 ly. Interior, the innermost circumstellar material lies along a cylindrical hourglass, 1 ly in radius and 4 ly long, which connects to the peanut by a thick equatorial disk. The nebulae are inclined 41o south and 8o east of the line of sight, slightly elliptical in cross section, and marginally offset west of the SN. The 3-D geometry of the three circumstellar rings is studied, suggesting the equatorial ring is elliptical (b/a<0.98), and spatially offset in the same direction as the hourglass. Dust-scattering models suggest that between the hourglass and bipolar lobes: the gas density drops from 1--3 cm^{-3} to >0.03 cm^{-3}; the maximum dust-grain size increases from ~0.2 micron to 2 micron; and the Si:C dust ratio decreases. The nebulae have a total mass of ~1.7 Msun, yielding a red-supergiant mass loss around 5*10^{-6} Msun yr^{-1}.

We present new results on the kinematics of the jet HH 110. New proper motion measurements have been calculated from [SII] CCD images obtained with a time baseline of nearly fifteen years. HH 110 proper motions show a strong asymmetry with respect to the outflow axis, with a general trend of pointing towards the west of the axis direction. Spatial velocities have been obtained by combining the proper motions and radial velocities from Fabry-Perot data. Velocities decrease by a factor ~3 over a distance of ~10$^{18}$ cm, much shorter than the distances expected for the braking caused by the jet/environment interaction. Our results show evidence of an anomalously strong interaction between the outflow and the surrounding environment, and are compatible with the scenario in which HH 110 emerges from a deflection in a jet/cloud collision.

We present the results of a spatio-kinematic study of the planetary nebula Menzel 1 using spectro-photometric mapping and a 3-D photoionization code. We create several 2-D emission line images from our long-slit spectra, and use these to derive the line fluxes for 15 lines, the Halpha/Hbeta extinction map, and the [SII] line ratio density map of the nebula. We use our photoionization code constrained by these data to derive the three-dimensional nebular structure and ionizing star parameters of Menzel 1 by simultaneously fitting the integrated line intensities, the density map, and the observed morphologies in several lines, as well as the velocity structure. Using theoretical evolutionary tracks of intermediate and low mass stars, we derive a mass for the central star of 0.63+-0.05 Msolar. We also derive a distance of 1050+_150 pc to Menzel 1.

We present the results of mapping the planetary nebula NGC6369 using multiple long slit spectra taken with the CTIO 1.5m telescope. We create two dimensional emission line images from our spectra, and use these to derive fluxes for 17 lines, the Halpha/Hbeta extinction map, the [SII] line ratio density map, and the [NII] temperature map of the nebula. We use our photoionization code constrained by these data to determine the distance, the ionizing star characteristics, and show that a clumpy hour-glass shape is the most likely three-dimensional structure for NGC6369. Note that our knowledge of the nebular structure eliminates all uncertainties associated with classical distance determinations, and our method can be applied to any spatially resolved emission line nebula. We use the central star, nebular emission line, and optical+IR luminosities to show that NGC6369 is matter bound, as about 70% of the Lyman continuum flux escapes. Using evolutionary tracks from Bloecker(1995) we derive a central star mass of about 0.65 Msolar.

We present ground-based near-infrared imaging and {\em HST} optical imaging and spectroscopy of the interaction between the ejecta of SN 1987A and its equatorial circumstellar ring. This interaction has made a transition, from emission originating in just a few ``hot spots'' at restricted locations in position angle around the ring, to a collision producing optical emission over a nearly continuous distribution, with few breaks larger than 45 degrees. The centroids of the first three spots are measured to move at 2000-3000 km/s, which we interpret as a lower limit of the velocity of the forward blast front. Multi-wavelength light curves of the spots show that they do not evolve uniformly, and change significantly on timescales as short as one month; in particular the first spot shows a significant break in its lightcurve. Implications of observed delays between spots appearances are discussed, which leads to a generalized model of hot spot evolution, and suggests that the early appearance of the first hot spot is explained by its inward radial position and a fairly uniform forward blast wave, rather than extraordinary physical circumstances. Data further suggest that the forward blast is reaching the bulk of the inner ring material to the east, the density of which appears higher than elsewhere in the ring. We study the ring geometry, finding evidence suggestive of an intrinsic ellipticity of 0.95, and find lower and upper distance limits of 47.9 +/- 0.92 kpc and 54.4 +/- 2.1 kpc, respectively.

We present a long-slit spectrum of the Herbig-Haro object HH 111, obtained with SITS on the Hubble Space Telescope. This spectrum has a spectral resolution of approx. 50 km/s, so it gives a good picture of the kinematical properties of the observed object at very high angular resolution. We find that some of the knots along the jet are associated with sudden drops in the radial velocity (modulus), confirming that the emission from the knots is formed in shocks. We interprt the observed position-velocity diagrams in terms of a model of a jet from a variable source, and we attempt to carry out a reconstruction of the ejection variability necessary for reproducing the observed kinematical structure.

We present the discovery of several new regions of interaction between the ejecta and equatorial circumstellar ring of SN1987A, an interaction leading to a much expanded development of the supernova remnant. We also trace the development of the first such "hot spot," discovered in 1997, back to 1995. Later hot spots seem to have emerged by early 1999. We discuss mechanisms for the long delay between the first and later spots.

We present several different measurements of the velocities of structures within the circumstellar envelope of SN 1987A, including the inner, equatorial ring (ER), outer rings (ORs), and the diffuse nebulosity at radii < 5 pc, based on CTIO 4m and HST data. A comparison of STIS and WFPC2 [N II]6583 loci for the rings show that the ER is expanding in radius at 10.5+-0.3 km/s, with the northern OR expanding along the line of sight at about 26 km/s, and for the southern OR, about 23 km/s. Similar results are found with CTIO 4m data. Accounting for inclination, the best fit to all data show both ORs with an expansion from the SN of 26 km/s. The ratio of the ER to OR velocities is nearly equal to the ratio of ER to OR radii, so the rings are roughly homologous, all having kinematic ages corresponding to about 20,000 yr before the SN explosion. This makes previously reported, large compositional differences between the ER and ORs difficult to understand. Additionally, a grid of longslit 4m/echelle spectra centered on the SN shows two velocity components over a region roughly coextensive with the outer circumstellar envelope extending about 5 pc (20 arcsec) from the SN. One component is blueshifted and the other redshifted from the SN centroid by about 10 km/s each. These features may represent a bipolar flow expanding from the SN, in which the ORs are propelled 10-15 km/s faster than that of the surrounding envelope into which they propogate. The kinematic timescale for the entire nebula is at least about 350,000 yr. The kinematics of these various structures constrain possible models for the evolution of the progenitor and its formation of a mass loss nebula.

We reexamine the spectroscopic underpinnings of recent suggestions that [O I] and [Fe II] lines from the Orion H II region are produced in gas where the iron-carrying grains have been destroyed and the electron density is surprisingly high. Our new observations show that previous detections of [O I] 5577 were dominated by telluric emission. Our limits are consistent with a moderate density (10^4 cm^{-3}) photoionized gas. We show that a previously proposed model of the Orion H II region reproduces the observed [O I] and [Fe II] spectrum. These lines are fully consistent with formation in a dusty region of moderate density.