Machine learning (automated processes that learn by example in order to classify, predict, discover or generate new data) and artificial intelligence (methods by which a computer makes decisions or discoveries that would usually require human intelligence) are now firmly established in astronomy. Every week, new applications of machine learning and artificial intelligence are added to a growing corpus of work. Random forests, support vector machines, and neural networks (artificial, deep, and convolutional) are now having a genuine impact for applications as diverse as discovering extrasolar planets, transient objects, quasars, and gravitationally-lensed systems, forecasting solar activity, and distinguishing between signals and instrumental effects in gravitational wave astronomy. This review surveys contemporary, published literature on machine learning and artificial intelligence in astronomy and astrophysics. Applications span seven main categories of activity: classification, regression, clustering, forecasting, generation, discovery, and the development of new scientific insight. These categories form the basis of a hierarchy of maturity, as the use of machine learning and artificial intelligence emerges, progresses or becomes established.

We search Dark Energy Survey (DES) Year 3 imaging for galaxy-galaxy strong gravitational lenses using convolutional neural networks, extending previous work with new training sets and covering a wider range of redshifts and colors. We train two neural networks using images of simulated lenses, then use them to score postage stamp images of 7.9 million sources from the Dark Energy Survey chosen to have plausible lens colors based on simulations. We examine 1175 of the highest-scored candidates and identify 152 probable or definite lenses. Examining an additional 20,000 images with lower scores, we identify a further 247 probable or definite candidates. After including 86 candidates discovered in earlier searches using neural networks and 26 candidates discovered through visual inspection of blue-near-red objects in the DES catalog, we present a catalog of 511 lens candidates.

We search Dark Energy Survey (DES) Year 3 imaging data for galaxy-galaxy strong gravitational lenses using convolutional neural networks. We generate 250,000 simulated lenses at redshifts > 0.8 from which we create a data set for training the neural networks with realistic seeing, sky and shot noise. Using the simulations as a guide, we build a catalogue of 1.1 million DES sources with 1.8 < g - i < 5, 0.6 < g -r < 3, r_mag > 19, g_mag > 20 and i_mag > 18.2. We train two ensembles of neural networks on training sets consisting of simulated lenses, simulated non-lenses, and real sources. We use the neural networks to score images of each of the sources in our catalogue with a value from 0 to 1, and select those with scores greater than a chosen threshold for visual inspection, resulting in a candidate set of 7,301 galaxies. During visual inspection we rate 84 as "probably" or "definitely" lenses. Four of these are previously known lenses or lens candidates. We inspect a further 9,428 candidates with a different score threshold, and identify four new candidates. We present 84 new strong lens candidates, selected after a few hours of visual inspection by astronomers. This catalogue contains a comparable number of high-redshift lenses to that predicted by simulations. Based on simulations we estimate our sample to contain most discoverable lenses in this imaging and at this redshift range.

The ability to quickly detect transient sources in optical images and trigger multi-wavelength follow up is key for the discovery of fast transients. These include events rare and difficult to detect such as kilonovae, supernova shock breakout, and "orphan" Gamma-ray Burst afterglows. We present the Mary pipeline, a (mostly) automated tool to discover transients during high-cadenced observations with the Dark Energy Camera (DECam) at CTIO. The observations are part of the "Deeper Wider Faster" program, a multi-facility, multi-wavelength program designed to discover fast transients, including counterparts to Fast Radio Bursts and gravitational waves. Our tests of the Mary pipeline on DECam images return a false positive rate of ~2.2% and a missed fraction of ~3.4% obtained in less than 2 minutes, which proves the pipeline to be suitable for rapid and high-quality transient searches. The pipeline can be adapted to search for transients in data obtained with imagers other than DECam.

We train and apply convolutional neural networks, a machine learning technique developed to learn from and classify image data, to Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) imaging for the identification of potential strong lensing systems. An ensemble of four convolutional neural networks was trained on images of simulated galaxy-galaxy lenses. The training sets consisted of a total of 62,406 simulated lenses and 64,673 non-lens negative examples generated with two different methodologies. The networks were able to learn the features of simulated lenses with accuracy of up to 99.8% and a purity and completeness of 94-100% on a test set of 2000 simulations. An ensemble of trained networks was applied to all of the 171 square degrees of the CFHTLS wide field image data, identifying 18,861 candidates including 63 known and 139 other potential lens candidates. A second search of 1.4 million early type galaxies selected from the survey catalog as potential deflectors, identified 2,465 candidates including 117 previously known lens candidates, 29 confirmed lenses/high-quality lens candidates, 266 novel probable or potential lenses and 2097 candidates we classify as false positives. For the catalog-based search we estimate a completeness of 21-28% with respect to detectable lenses and a purity of 15%, with a false-positive rate of 1 in 671 images tested. We predict a human astronomer reviewing candidates produced by the system would identify ~20 probable lenses and 100 possible lenses per hour in a sample selected by the robot. Convolutional neural networks are therefore a promising tool for use in the search for lenses in current and forthcoming surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope.

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.

Planetary ephemerides have been developed and improved over centuries. They are a fundamental tool for understanding solar system dynamics, and essential for planetary and small body mass determinations, occultation predictions, high-precision tests of general relativity, pulsar timing, and interplanetary spacecraft navigation. This paper presents recent results from a continuing program of high-precision astrometric very-long-baseline interferometry (VLBI) observations of the Cassini spacecraft orbiting Saturn, using the Very Long Baseline Array (VLBA). We have previously shown that VLBA measurements can be combined with spacecraft orbit determinations from Doppler and range tracking and VLBI links to the inertial extragalactic reference frame (ICRF) to provide the most accurate barycentric positions currently available for Saturn. Here we report an additional five years of VLBA observations along with improved phase reference source positions, resulting in an improvement in residuals with respect to the Jet Propulsion Laboratory's dynamical ephemeris.

We study the variation of the shape of the proton intensity-time profiles in simulated gradual Solar Energetic Particle (SEP) events with the relative observer's position in space with respect to the main direction of propagation of an interplanetary (IP) shock. Using a three-dimensional (3D) magnetohydrodynamic (MHD) code to simulate such a shock, we determine the evolution of the downstream-to-upstream ratios of the plasma variables at its front. Under the assumption of an existing relation between the normalized ratio in speed across the shock front and the injection rate of shock-accelerated particles, we model the transport of the particles and we obtain the proton flux profiles to be measured by a grid of 18 virtual observers located at 0.4 and 1.0 AU, with different latitudes and longitudes with respect to the shock nose. The differences among flux profiles are the result of the way each observer establishes a magnetic connection with the shock front, and we find that changes in the observer's latitude may result in intensity changes of up to one order of magnitude at both radial distances considered here. The peak intensity variation with the radial distance for the pair of observers located at the same angular position is also derived. This is the first time that the latitudinal dependence of the peak intensity with the observer's heliocentric radial distance has been quantified within the framework of gradual SEP event simulations.

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.

Most simulations of coronal mass ejections (CMEs) to date either focus on the interplanetary propagation of a giant plasma "blob" without paying too much attention to its origin and to the formation process or they focus on the complex evolution of the coronal magnetic field due to (sub-)photospheric motions which result in an eruption. Here, we present global simulations of CMEs where coronal motions are used to produce a realistic evolution of the coronal magnetic field and cause an eruption. We focus on active region 10069, which produced a number of eruptions in late August 2002, including the August 24, 2002 CME - a fast (~2000 km/s) eruption originating from W81-, as well as a slower eruption on August 22, 2002 (originating from W62). Using a three-dimensional magneto-hydrodynamic (MHD) simulation of these ejections with the Space Weather Modeling Framework (SWMF), we show how a realistic initiation mechanism enables us to study the deflection of the CME in the corona and in the heliosphere. Reconnection of the erupting magnetic field with that of neighboring streamers and active regions modify the solar connectivity of the field lines connecting to Earth and change the expected solar energetic particle fluxes. Comparing the results at 1 AU of our simulations with in situ observations by the ACE spacecraft, we propose an alternate solar origin for the shock wave observed at L1 on August 26.

Eruption of a coronal mass ejection (CME) drags and "opens" the coronal magnetic field, presumably leading to the formation of a large-scale current sheet and the field relaxation by magnetic reconnection. We analyze physical characteristics of ray-like coronal features formed in the aftermath of CMEs, to check if the interpretation of this phenomenon in terms of reconnecting current sheet is consistent with the observations. The study is focused on measurements of the ray width, density excess, and coronal velocity field as a function of the radial distance. The morphology of rays indicates that they occur as a consequence of Petschek-like reconnection in the large scale current sheet formed in the wake of CME. The hypothesis is supported by the flow pattern, often showing outflows along the ray, and sometimes also inflows into the ray. The inferred inflow velocities range from 3 to 30 km s$^{-1}$, consistent with the narrow opening-angle of rays, adding up to a few degrees. The density of rays is an order of magnitude larger than in the ambient corona. The density-excess measurements are compared with the results of the analytical model in which the Petschek-like reconnection geometry is applied to the vertical current sheet, taking into account the decrease of the external coronal density and magnetic field with height. The model results are consistent with the observations, revealing that the main cause of the density excess in rays is a transport of the dense plasma from lower to larger heights by the reconnection outflow.

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 Hubble Space Telescope (HST) imagery and photometry of the active galaxy Markarian 421 and its companion galaxy 14 arcsec to the ENE. The HST images indicate that the companion is a morphological spiral rather than elliptical as previous ground--based imaging has concluded. The companion has a bright, compact nucleus, appearing unresolved in the HST images. This is suggestive of Seyfert activity, or possibly a highly luminous compact star cluster. We also report the results of high dynamic range long-slit spectroscopy with the slit placed to extend across both galaxies and nuclei. We detect no emission lines in the companion nucleus, though there is evidence for recent star formation. Velocities derived from a number of absorption lines visible in both galaxies indicate that the two systems are probably tidally bound and thus in close physical proximity. Using the measured relative velocities, we derive a lower limit on the MKN 421 mass within the companion orbit (R \sim 10 kpc) of 5.9 \times 10^{11} solar masses, and a mass-to-light ratio of >= 17. Our spectroscopy also shows for the first time the presence of H\alpha and [NII] emission lines from the nucleus of MKN 421, providing another example of the appearance of new emission features in the previously featureless spectrum of a classical BL Lac object. We see both broad and narrow line emission, with a velocity dispersion of several thousand km s^{-1} evident in the broad lines.

Summarizes current status of radio interferometry at radio frequencies between Earth-based receivers, for astrometric and geodetic applications. Emphasizes theoretical models of VLBI observables that are required to extract results at the present accuracy levels of 1 cm and 1 nanoradian. Highlights the achievements of VLBI during the past two decades in reference frames, Earth orientation, atmospheric effects on microwave propagation, and relativity.