We present the Strategic Scientific Plan (SSP) for the direction and activities of the Gemini Observatory in the 2020s. The overarching goal is to ensure that Gemini best uses the available resources to serve the needs of its international user community throughout the coming decade. The actionable items fall into three general categories: (1) preserving Gemini's current facilities and strengths; (2) developing instrumentation and software systems, including data pipelines, to enable new scientific capabilities that build on those strengths; (3) strategizing how visiting instruments can deliver additional valuable capabilities. We provide a high-level timeline (schematically illustrated in one figure) for the main developments discussed in this SSP. The schedule is ambitious, but in light of the recent Gemini in the Era of Multi-Messenger Astronomy (GEMMA) award from the NSF, the plan becomes achievable. Lists of milestones are given for gauging progress. As these milestones are reached and new instruments become available, some current instruments will need to be retired; we make recommendations in this regard. The final section concludes by reemphasizing the importance of a strong partnership committed to the needs of all members.

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.

The JAXA Hayabusa-2 mission was approved in 2010 and launched on December 3, 2014. The spacecraft will arrive at the near-Earth asteroid 162173 Ryugu in 2018 where it will perform a survey, land and obtain surface material, then depart in Dec 2019 and return to Earth in Dec 2020. We observed Ryugu with the Herschel Space Observatory in Apr 2012 at far-IR thermal wavelengths, supported by several ground-based observations to obtain optical lightcurves. We reanalysed previously published Subaru-COMICS and AKARI-IRC observations and merged them with a Spitzer-IRS data set. In addition, we used a large set of Spitzer-IRAC observations obtained in the period Jan to May, 2013. The data set includes two complete rotational lightcurves and a series of ten "point-and-shoot" observations. The almost spherical shape of the target together with the insufficient lightcurve quality forced us to combine radiometric and lightcurve inversion techniques in different ways to find the object's key physical and thermal parameters. We find that the solution which best matches our data sets leads to this C class asteroid having a retrograde rotation with a spin-axis orientation of (lambda = 310-340 deg; beta = -40+/-15 deg) in ecliptic coordinates, an effective diameter (of an equal-volume sphere) of 850 to 880 m, a geometric albedo of 0.044 to 0.050 and a thermal inertia in the range 150 to 300 Jm-2s-0.5K-1. Based on estimated thermal conductivities of the top-layer surface in the range 0.1 to 0.6 WK-1m-1, we calculated that the grain sizes are approximately equal to between 1 and 10 mm. The finely constrained values for this asteroid serve as a `design reference model', which is currently used for various planning, operational and modelling purposes by the Hayabusa2 team.

We present a compilation of UBV RIz light curves of 51 type II supernovae discovered during the course of four different surveys during 1986 to 2003: the Cerro Tololo Supernova Survey, the Calan/Tololo Supernova Program (C&T), the Supernova Optical and Infrared Survey (SOIRS), and the Carnegie Type II Supernova Survey (CATS). The photometry is based on template-subtracted images to eliminate any potential host galaxy light contamination, and calibrated from foreground stars. This work presents these photometric data, studies the color evolution using different bands, and explores the relation between the magnitude at maximum brightness and the brightness decline parameter (s) from maximum light through the end of the recombination phase. This parameter is found to be shallower for redder bands and appears to have the best correlation in the B band. In addition, it also correlates with the plateau duration, being thus shorter (longer) for larger (smaller) s values.

We investigated the magnitude-phase relation of (162173) 1999 JU3, a target asteroid for the JAXA Hayabusa 2 sample return mission. We initially employed the international Astronomical Union's H-G formalism but found that it fits less well using a single set of parameters. To improve the inadequate fit, we employed two photometric functions, the Shevchenko and Hapke functions. With the Shevchenko function, we found that the magnitude-phase relation exhibits linear behavior in a wide phase angle range (alpha = 5-75 deg) and shows weak nonlinear opposition brightening at alpha< 5 deg, providing a more reliable absolute magnitude of Hv = 19.25 +- 0.03. The phase slope (0.039 +- 0.001 mag/deg) and opposition effect amplitude (parameterized by the ratio of intensity at alpha=0.3 deg to that at alpha=5 deg, I(0.3)/I(5)=1.31+-0.05) are consistent with those of typical C-type asteroids. We also attempted to determine the parameters for the Hapke model, which are applicable for constructing the surface reflectance map with the Hayabusa 2 onboard cameras. Although we could not constrain the full set of Hapke parameters, we obtained possible values, w=0.041, g=-0.38, B0=1.43, and h=0.050, assuming a surface roughness parameter theta=20 deg. By combining our photometric study with a thermal model of the asteroid (Mueller et al. in preparation), we obtained a geometric albedo of pv = 0.047 +- 0.003, phase integral q = 0.32 +- 0.03, and Bond albedo AB = 0.014 +- 0.002, which are commensurate with the values for common C-type asteroids.

We present an analysis of the diversity of V-band light-curves of hydrogen-rich type II supernovae. Analyzing a sample of 116 supernovae, several magnitude measurements are defined, together with decline rates at different epochs, and time durations of different phases. It is found that magnitudes measured at maximum light correlate more strongly with decline rates than those measured at other epochs: brighter supernovae at maximum generally have faster declining light-curves at all epochs. We find a relation between the decline rate during the 'plateau' phase and peak magnitudes, which has a dispersion of 0.56 magnitudes, offering the prospect of using type II supernovae as purely photometric distance indicators. Our analysis suggests that the type II population spans a continuum from low-luminosity events which have flat light-curves during the 'plateau' stage, through to the brightest events which decline much faster. A large range in optically thick phase durations is observed, implying a range in progenitor envelope masses at the epoch of explosion. During the radioactive tails, we find many supernovae with faster declining light-curves than expected from full trapping of radioactive emission, implying low mass ejecta. It is suggested that the main driver of light-curve diversity is the extent of hydrogen envelopes retained before explosion. Finally, a new classification scheme is introduced where hydrogen-rich events are typed as simply 'SNII' with an s2 value giving the decline rate during the 'plateau' phase, indicating its morphological type.

In classical P-Cygni profiles, theory predicts emission to peak at zero rest velocity. However, supernova spectra exhibit emission that is generally blue shifted. While this characteristic has been reported in many supernovae, it is rarely discussed in any detail. Here we present an analysis of H-alpha emission-peaks using a dataset of 95 type II supernovae, quantifying their strength and time evolution. Using a post-explosion time of 30d, we observe a systematic blueshift of H-alpha emission, with a mean value of -2000 kms-1. This offset is greatest at early times but vanishes as supernovae become nebular. Simulations of Dessart et al. (2013) match the observed behaviour, reproducing both its strength and evolution in time. Such blueshifts are a fundamental feature of supernova spectra as they are intimately tied to the density distribution of ejecta, which falls more rapidly than in stellar winds. This steeper density structure causes line emission/absorption to be much more confined; it also exacerbates the occultation of the receding part of the ejecta, biasing line emission to the blue for a distant observer. We conclude that blue-shifted emission-peak offsets of several thousand kms-1 are a generic property of observations, confirmed by models, of photospheric-phase type II supernovae.

We present a spectroscopic analysis of the H-alpha profiles of hydrogen-rich type II supernovae. A total of 52 type II supernovae having well sampled optical light curves and spectral sequences were analyzed. Concentrating on the H-alpha P-Cygni profile we measure its velocity from the FWHM of emission and the ratio of absorption to emission (a/e) at a common epoch at the start of the recombination phase, and search for correlations between these spectral parameters and photometric properties of the V-band light curves. Testing the strength of various correlations we find that a/e appears to be the dominant spectral parameter in terms of describing the diversity in our measured supernova properties. It is found that supernovae with smaller a/e have higher H-alpha velocities, more rapidly declining light curves from maximum, during the plateau and radioactive tail phase, are brighter at maximum light and have shorter optically thick phase durations. We discuss possible explanations of these results in terms of physical properties of type II supernovae, speculating that the most likely parameters which influence the morphologies of H-alpha profiles are the mass and density profile of the hydrogen envelope, together with additional emission components due to circumstellar interaction.

Cerro Las Campanas located at Las Campanas Observatory (LCO) in Chile has been selected as the site for the Giant Magellan Telescope. We report results obtained since the commencement, in 2005, of a systematic site testing survey of potential GMT sites at LCO. Meteorological (cloud cover, temperature, pressure, wind, and humidity) and DIMM seeing data have been obtained at three potential sites, and are compared with identical data taken at the site of the twin Magellan 6.5m telescopes. In addition, measurements of the turbulence profile of the free-atmosphere above LCO have been collected with a MASS/DIMM. Furthermore, we consider photometric quality, light pollution, and precipitable water vapor (PWV). LCO, and Co. Las Campanas in particular, have dark skies, little or no risk of future light pollution, excellent seeing, moderate winds, PWV adequate for mid-IR astronomy during a reasonable fraction of the nights, and a high fraction of clear nights overall. Finally, Co. Las Campanas meets or exceeds all the defined science requirements.

We present an analysis of the science image quality obtained on the twin 6.5 metre Magellan telescopes over a 1.5 year period, using images of ~10^5 stars. We find that the telescopes generally obtain significantly better image quality than the DIMM-measured seeing. This is qualitatively consistent with expectations for large telescopes, where the wavefront outer scale of the turbulence spectrum plays a significant role. However, the dominant effect is found to be wind speed with Magellan outperforming the DIMMs most markedly when the wind is strongest. Excluding data taken during strong wind conditions (>10 m/s), we find that the Magellan telescopes still significantly outperform the DIMM seeing, and we estimate the site to have L_0 ~ 25 m on average. We also report on the first detection of a negative bias in DIMM data. This is found to occur, as predicted, when the DIMM is affected by certain optical aberrations and the turbulence profile is dominated by the upper layers of the atmosphere.

An extensive dataset for SN 2003hv that covers the flux evolution from maximum light to day +786 is presented. The data are combined with published nebular-phase infrared spectra, and the observations are compared to model light curves and synthetic nebular spectra. SN 2003hv is a normal Type Ia supernova (SN Ia) with photometric and spectroscopic properties consistent with its rarely observed B-band decline-rate parameter, Delta m_15 = 1.61 +- 0.02. The blueshift of the most isolated [Fe II] lines in the nebular-phase optical spectrum appears consistent with those observed in the infrared at similar epochs. At late times there is a prevalent color evolution from the optical toward the near-infrared bands. We present the latest-ever detection of a SN Ia in the near-infrared in Hubble Space Telescope images. The study of the ultraviolet/optical/infrared (UVOIR) light curve reveals that a substantial fraction of the flux is "missing" at late times. Between 300-700 days past maximum brightness, the UVOIR light curve declines linearly following the decay of radioactive Co56, assuming full and instantaneous positron trapping. At 700 days we detect a possible slowdown of the decline in optical bands, mainly in the V band. The data are incompatible with a dramatic infrared catastrophe. However, the idea that an infrared catastrophe occurred in the densest regions before 350 days can explain the missing flux from the UVOIR wavelengths and the flat-topped profiles in the near-infrared. We argue that such a scenario is possible if the ejecta are clumpy. The observations suggest that positrons are most likely trapped in the ejecta.

We obtained optical photometry of SN 2003gs on 49 nights, from 2 to 494 days after T(B_max). We also obtained near-IR photometry on 21 nights. SN 2003gs was the first fast declining Type Ia SN that has been well observed since SN 1999by. While it was subluminous in optical bands compared to more slowly declining Type Ia SNe, it was not subluminous at maximum light in the near-IR bands. There appears to be a bimodal distribution in the near-IR absolute magnitudes of Type Ia SNe at maximum light. Those that peak in the near-IR after T(B_max) are subluminous in the all bands. Those that peak in the near-IR prior to T(B_max), such as SN 2003gs, have effectively the same near-IR absolute magnitudes at maximum light regardless of the decline rate Delta m_15(B). Near-IR spectral evidence suggests that opacities in the outer layers of SN 2003gs are reduced much earlier than for normal Type Ia SNe. That may allow gamma rays that power the luminosity to escape more rapidly and accelerate the decline rate. This conclusion is consistent with the photometric behavior of SN 2003gs in the IR, which indicates a faster than normal decline from approximately normal peak brightness.

Optical and near-infrared photometry and optical spectroscopy are reported for SN 2003bg, starting a few days after explosion and extending for a period of more than 300 days. Our early-time spectra reveal the presence of broad, high-velocity Balmer lines. The nebular-phase spectra, on the other hand, show a remarkable resemblance to those of Type Ib/c supernovae, without clear evidence for hydrogen. Near maximum brightness SN 2003bg displayed a bolometric luminosity comparable to that of other Type I hypernovae unrelated to gamma-ray bursts, implying a rather normal amount of 56Ni production (0.1-0.2 Msun) compared with other such objects. The bolometric light curve of SN 2003bg, on the other hand, is remarkably broad, thus suggesting a relatively large progenitor mass at the moment of explosion. These observations, together with the large value of the kinetic energy of expansion established in the accompanying paper (Mazzali et al. 2009), suggest that SN 2003bg can be regarded as a Type IIb hypernova.

We present infrared photometry of SN 1999em, plus optical photometry, infrared photometry, and optical spectroscopy of SN 2003hn. Both objects were Type II-P supernovae. The V-[RIJHK] color curves of these supernovae evolved in a very similar fashion until the end of plateau phase. This allows us to determine how much more extinction the light of SN 2003hn suffered compared to SN 1999em. Since we have an estimate of the total extinction suffered by SN 1999em from model fits of ground-based and space-based spectra as well as photometry of SN 1999em, we can estimate the total extinction and absolute magnitudes of SN 2003hn with reasonable accuracy. Since the host galaxy of SN 2003hn also produced the Type Ia SN 2001el, we can directly compare the absolute magnitudes of these two SNe of different types.

We observed a stellar occultation by Titan on 2003 November 14 from La Palma Observatory using ULTRACAM with three Sloan filters: u', g', and i' (358, 487, and 758 nm, respectively). The occultation probed latitudes 2 degrees S and 1 degrees N during immersion and emersion, respectively. A prominent central flash was present in only the i' filter, indicating wavelength-dependent atmospheric extinction. We inverted the light curves to obtain six lower-limit temperature profiles between 335 and 485 km (0.04 and 0.003 mb) altitude. The i' profiles agreed with the temperature measured by the Huygens Atmospheric Structure Instrument [Fulchignoni, M. et al., 2005. Nature 438, 785-791] above 415 km (0.01 mb). The profiles obtained from different wavelength filters systematically diverge as altitude decreases, which implies significant extinction in the light curves. Applying an extinction model [Elliot, J.L., Young, L.A., 1992. Astron. J. 103, 991-1015] gave the altitudes of line of sight optical depth equal to unity: 396 +/- 7 km and 401 +/- 20 km (u' immersion and emersion); 354 +/- 7 km and 387 +/- 7 km (g' immersion and emersion); and 336 +/- 5 km and 318 +/- 4 km (i' immersion and emersion). Further analysis showed that the optical depth follows a power law in wavelength with index 1.3 +/- 0.2. We present a new method for determining temperature from scintillation spikes in the occulting body's atmosphere. Temperatures derived with this method are equal to or warmer than those measured by the Huygens Atmospheric Structure Instrument. Using the highly structured, three-peaked central flash, we confirmed the shape of Titan's middle atmosphere using a model originally derived for a previous Titan occultation [Hubbard, W.B. et al., 1993. Astron. Astrophys. 269, 541-563].

We report precipitable water vapor (PWV) measurements made at Las Campanas Observatory using optical spectra of H2O lines obtained with the Magellan echelle spectrograph, and calculated using a robust technique that is accurate to 5-10%. Calibration of the relationship between our PWV measurements and opacity values at 225 GHz was made possible by simultaneous tipping radiometer observations. Based on this calibration, we present Las Campanas Observatory winter-time precipitable water vapor statistics, measured using the tipping radiometer, during a two month campaign. The median value of 2.8 +- 0.3 mm is consistent with that measured at the nearby La Silla Observatory during the VLT site survey. We conclude that, in the Southern hemisphere winter months, we can expect good conditions for infrared observing (<~1.5 mm) approximately 10% of the time at Las Campanas Observatory.

We present results from our two year study of ground-layer turbulence as seen through the 6.5-meter Magellan Telescopes at Las Campanas Observatory. The experiment consists of multiple, moderate resolution, Shack-Hartmann wavefront sensors deployed over a large 16 arcminute field. Over the two years of the experiment, the ground-layer turbulence has been sampled on eleven nights in a variety of seeing and wind conditions. On most nights the ground-layer turbulence contributes 10% to the total visible-band seeing, although a few nights exhibit ground-layer contributions up to 30%. We present the ground-layer turbulence on the sampled nights as well as a demonstration of its strength as a function of field size. This information is combined with data from a MASS-DIMM seeing monitor adjacent to the Magellan Telescopes to infer the annual ground-layer contribution to seeing at Las Campanas.