Resolving the phase structure of neutral hydrogen (HI) is crucial for understanding the life cycle of the interstellar medium (ISM). However, accurate measurements of HI temperature and density are limited by the availability of background continuum sources for measuring HI absorption. Here we test the use of deep learning for extracting HI properties over large areas without optical depth information. We train a 1D convolutional neural network using synthetic observations of 3D numerical simulations of the ISM to predict the fraction of cold neutral medium (f_CNM) and the correction to the optically-thin HI column density for optical depth (R_HI) from $21\rm\,cm$ emission alone. We restrict our analysis to high Galactic latitudes ($|b|>30 deg), where the complexity of spectral line profiles is minimized. We verify that the network accurately predicts f_CNM and R_HI by comparing the results with direct constraints from 21 cm absorption. By applying the network to the GALFA-HI survey, we generate large-area maps of f_CNM and R_HI. Although the overall contribution to the total HI column of cold neutral medium (CNM)-rich structures is small (~5%), we find that these structures are ubiquitous. Our results are consistent with the picture that small-scale structures observed in 21 cm emission aligned with the magnetic field are dominated by CNM. Finally, we demonstrate that the observed correlation between HI column density and dust reddening (E(B-V)) declines with increasing R_HI, indicating that future efforts to quantify foreground Galactic E(B-V) using HI, even at high latitudes, should increase fidelity by accounting for HI phase structure.

We investigate and demonstrate the use of convolutional neural networks (CNNs) for the task of distinguishing between merging and non-merging galaxies in simulated images, and for the first time at high redshifts (i.e. $z=2$). We extract images of merging and non-merging galaxies from the Illustris-1 cosmological simulation and apply observational and experimental noise that mimics that from the Hubble Space Telescope; the data without noise form a "pristine" data set and that with noise form a "noisy" data set. The test set classification accuracy of the CNN is $79\%$ for pristine and $76\%$ for noisy. The CNN outperforms a Random Forest classifier, which was shown to be superior to conventional one- or two-dimensional statistical methods (Concentration, Asymmetry, the Gini, $M_{20}$ statistics etc.), which are commonly used when classifying merging galaxies. We also investigate the selection effects of the classifier with respect to merger state and star formation rate, finding no bias. Finally, we extract Grad-CAMs (Gradient-weighted Class Activation Mapping) from the results to further assess and interrogate the fidelity of the classification model.

A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities match with dominant HI components. We compare the observed radial and tangential velocities of these stars with predictions from the state-of-the-art rotating disk model based on high-resolution 21 cm observations of the SMC from the Australian Square Kilometer Array Pathfinder telescope. We find that the observed kinematics of gas-tracing stars are inconsistent with disk rotation. We conclude that the kinematics of gas in the SMC are more complex than can be inferred from the integrated radial velocity field. As a result of violent tidal interactions with the LMC, non-rotational motions are prevalent throughout the SMC, and it is likely composed of distinct sub-structures overlapping along the line of sight.

This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy.

The spatial distribution of neutral hydrogen (H I) emission is a powerful probe of interstellar medium physics. The small-scale structure in H I channel maps is often assumed to probe the velocity field rather than real density structures. In this work we directly test this assumption, using high-resolution GALFA-H I observations and 50,985 quasar spectra from SDSS. We measure the equivalent widths of interstellar Na I D$_1$ and Na I D$_2$ absorption, and robustly conclude that they depend more strongly on the column density of small-scale structure in H I than on either the large-scale H I structure or the total H I column. This is inconsistent with the hypothesis that small-scale channel map structure is driven by velocity crowding. Instead, the data favor the interpretation that this emission structure predominantly originates in cold, dense interstellar material, consistent with a clumpy cold neutral medium.

The interstellar medium (ISM) of galaxies is composed of a turbulent magnetized plasma. In order to quantitatively measure relevant turbulent parameters of the ISM, a wide variety of statistical techniques and metrics have been developed that are often tested using numerical simulations and analytic formalism. These metrics are typically based on the Fourier power spectrum, which does not capture the Fourier phase information that carries the morphological characteristics of images. In this work we use density slices of magnetohydrodyanmic turbulence simulations to demonstrate that a modern tool, convolutional neural networks, can capture significant information encoded in the Fourier phases. We train the neural network to distinguish between two simulations with different levels of magnetization. We find that, even given a tiny slice of simulation data, a relatively simple network can distinguish sub-Alfv\'enic (strong magnetic field) and super-Alfv\'enic (weak magnetic field) turbulence >98% of the time, even when all spectral amplitude information is stripped from the images. In order to better understand how the neural network is picking out differences betweem the two classes of simulations we apply a neural network analysis method called "saliency maps". The saliency map analysis shows that sharp ridge-like features are a distinguishing morphological characteristic in such simulations. Our analysis provides a way forward for deeper understanding of the relationship between magnetohydrodyanmic turbulence and gas morphology and motivates further applications of neural networks for studies of turbulence. We make publicly available all data and software needed to reproduce our results.

We present a bibliographic analysis of Chandra, Hubble, and Spitzer publications. We find (a) archival data are used in >60% of the publication output and (b) archives for these missions enable a much broader set of institutions and countries to scientifically use data from these missions. Specifically, we find that authors from institutions that have published few papers from a given mission publish 2/3 archival publications, while those with many publications typically have 1/3 archival publications. We also show that countries with lower GDP per capita overwhelmingly produce archival publications, while countries with higher GDP per capital produce guest observer and archival publications in equal amounts. We argue that robust archives are thus not only critical for the scientific productivity of mission data, but also the scientific accessibility of mission data. We argue that the astronomical community should support archives to maximize the overall scientific societal impact of astronomy, and represent an excellent investment in astronomy's future.

Commodity cloud computing, as provided by commercial vendors such as Amazon, Google, and Microsoft, has revolutionized computing in many sectors. With the advent of a new class of big data, public access astronomical facility such as LSST, DKIST, and WFIRST, there exists a real opportunity to combine these missions with cloud computing platforms and fundamentally change the way astronomical data is collected, processed, archived, and curated. Making these changes in a cross-mission, coordinated way can provide unprecedented economies of scale in personnel, data collection and management, archiving, algorithm and software development and, most importantly, science.

Astrophysics lies at the crossroads of big datasets (such as the Large Synoptic Survey Telescope and Gaia), open source software to visualize and interpret high dimensional datasets (such as Glue, WorldWide Telescope, and OpenSpace), and uniquely skilled software engineers who bridge data science and research fields. At the same time, more than 4,000 planetariums across the globe immerse millions of visitors in scientific data. We have identified the potential for critical synergy across data, software, hardware, locations, and content that -- if prioritized over the next decade -- will drive discovery in astronomical research. Planetariums can and should be used for the advancement of scientific research. Current facilities such as the Hayden Planetarium in New York City, Adler Planetarium in Chicago, Morrison Planetarium in San Francisco, the Iziko Planetarium and Digital Dome Research Consortium in Cape Town, and Visualization Center C in Norrkoping are already developing software which ingests catalogs of astronomical and multi-disciplinary data critical for exploration research primarily for the purpose of creating scientific storylines for the general public. We propose a transformative model whereby scientists become the audience and explorers in planetariums, utilizing software for their own investigative purposes. In this manner, research benefits from the authentic and unique experience of data immersion contained in an environment bathed in context and equipped for collaboration. Consequently, in this white paper we argue that over the next decade the research astronomy community should partner with planetariums to create visualization-based research opportunities for the field. Realizing this vision will require new investments in software and human capital.

We investigate the physical properties of structures seen in channel map observations of 21-cm neutral hydrogen (HI) emission. HI intensity maps display prominent linear structures that are well aligned with the ambient magnetic field in the diffuse interstellar medium (ISM). Some literature hold that these structures are "velocity caustics", fluctuations imprinted by the turbulent velocity field, and are not three-dimensional density structures in the ISM. We test this hypothesis by stacking probes of the density field -- broadband far infrared (FIR) emission and the integrated HI column density (NHI) -- at the locations of linear HI intensity structures. We find that the HI intensity features are real density structures and not velocity caustics. We generalize the investigation to all small-scale structure in HI channel maps, and analyze this correlation as a function of velocity channel width, finding no measurable contribution from velocity caustics to the HI channel map emission. Further, we find that small-scale HI channel maps structures have elevated FIR/NHI, implying that this emission originates from a colder, denser phase of the ISM than the surrounding material. The data are consistent with a multi-phase diffuse ISM in which small-scale structures in narrow HI channel maps are preferentially cold neutral medium (CNM) that is anisotropically distributed and aligned with the local magnetic field. The shallow spatial power spectrum (SPS) of narrow HI channels is often attributed to velocity caustics. We conjecture instead that the small-scale structure and narrow linewidths typical of CNM explain the observed relationship between the SPS and channel width.

Galaxies are surrounded by extended atmospheres, which are often called the circumgalactic medium (CGM) and are the least understood part of galactic ecosystems. The CGM serves as a reservoir of both diffuse, metal-poor gas accreted from the intergalactic medium, and metal-rich gas that is either ejected from galaxies by energetic feedback or stripped from infalling satellites. As such, the CGM is empirically multi-phased and complex in dynamics. Significant progress has been made in the past decade or so in observing the cosmic-ray/B-field, as well as various phases of the CGM. But basic questions remain to be answered. First, what are the energy, mass, and metal contents of the CGM? More specifically, how are they spatially distributed and partitioned in the different components? Moreover, how are they linked to properties of host galaxies and their global clustering and intergalactic medium environments? Lastly, what are the origin, state, and life-cycle of the CGM? This question explores the dynamics of the CGM. Here we illustrate how these questions may be addressed with multi-wavelength observations of the CGM.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM) in three dimensions throughout the Milky Way disk and bulge. These developments will enable a fundamentally new understanding of dust properties, including how grains respond to their local environment and how those environments affect dust attenuation of background objects of interest. Distance-resolved maps of dust motion also hold great promise for tracing the flow of interstellar material throughout the Galaxy on a variety of scales, from bar-streaming motions to the collapse and dissolution of individual molecular clouds. These advances require optical and infrared imaging of stars throughout the Galactic midplane, stretching many kiloparsecs from the Sun, matched with very dense spectroscopic coverage to probe the ISM's fine-grained structure.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

Every quasar (QSO) spectrum contains absorption-line signatures from the interstellar medium, disk-halo interface, and circumgalactic medium (CGM) of the Milky Way (MW). We analyze Hubble Space Telescope/Cosmic Origins Spectrograph (COS) spectra of 132 QSOs to study the significance and origin of SiIV absorption at $|v_{\rm LSR}|\leq100$ km/s in the Galactic halo. The gas in the north predominantly falls in at $-50\lesssim v_{\rm LSR}\lesssim 0$ km/s, whereas in the south, no such pattern is observed. The SiIV column density has an average and a standard deviation of $\langle N_{\rm SiIV}\rangle=(3.8\pm1.4)\times10^{13}$ cm$^{-2}$. At $|b|\gtrsim 30$ degree, $N_{\rm SiIV}$ does not significantly correlate with $b$, which cannot be explained by a commonly adopted flat-slab geometry. We propose a two-component model to reconstruct the $N_{\rm SiIV}$-$b$ distribution: a plane-parallel component $N_{\rm DH}^{\perp}$ to account for the MW's disk-halo interface and a global component $N_{\rm G}$ to reproduce the weak dependence on $b$. We find $N_{\rm DH}^{\perp}=1.3^{+4.7}_{-0.7}\times10^{12}$ cm$^{-2}$ and $N_{\rm G}=(3.4\pm0.3)\times10^{13}$ cm$^{-2}$ on the basis of Bayesian analyses and block bootstrapping. The global component is most likely to have a Galactic origin, although its exact location is uncertain. If it were associated with the MW's CGM, we would find $M_{\rm gas, all}\gtrsim4.7\times10^9\ M_{\odot} (\frac{C_f}{1})(\frac{R}{75\ {\rm kpc}})^2 (\frac{f_{\rm SiIV}}{0.3})^{-1}(\frac{Z}{0.3\ Z_{\odot}})^{-1}$ for the cool gas at all velocities in the Galactic halo. Our analyses show that there is likely a considerable amount of gas at $|v_{\rm LSR}|\leq100$ km s$^{-1}$ hidden in the MW's CGM. Along with this work, we make our QSO dataset publicly available as the COS Quasar Database for Galactic Absorption Lines (COS-GAL).

Metal Evolution, Transport, and Abundance in the LMC (METAL) is a large Cycle 24 program on the Hubble Space Telescope aimed at measuring dust extinction properties and interstellar depletions in the Large Magellanic Cloud (LMC) at half-solar metallicity. The 101-orbit program is comprised of COS and STIS spectroscopy toward 33 LMC massive stars between 1150 A and 3180 A, and parallel WFC3 imaging in 7 NUV-NIR filters. The fraction of silicon in the gas-phase (depletion) obtained from the spectroscopy decreases with increasing hydrogen column density. Depletion patterns for silicon differ between the Milky Way, LMC, and SMC, with the silicon depletion level offsetting almost exactly the metallicity differences, leading to constant gas-phase abundances in those galaxies for a given hydrogen column density. The silicon depletion correlates linearly with the absolute-to-selective extinction, R$_V$, indicating a link between gas depletion and dust grain size. Extinction maps are derived from the resolved stellar photometry in the parallel imaging, which can be compared to FIR images from Herschel and Spitzer to estimate the emissivity of dust at LMC metallicity. The full METAL sample of depletions, UV extinction curves, and extinction maps will inform the abundance, size, composition, and optical properties of dust grains in the LMC, comprehensively improve our understanding of dust properties, and the accuracy with which dust-based gas masses, star formation rates and histories in nearby and high-redshift galaxies are estimated. This overview paper describes the goals, design, data reduction, and initial results of the METAL survey.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

The majority of hydrogen in the interstellar medium (ISM) is in atomic form. The transition from atoms to molecules and, in particular, the formation of the H$_2$ molecule, is a key step in cosmic structure formation en route to stars. Quantifying H$_2$ formation in space is difficult, due to the confusion in the emission of atomic hydrogen (HI) and the lack of a H$_2$ signal from the cold ISM. Here we present the discovery of a rare, isolated dark cloud currently undergoing H$_2$ formation, as evidenced by a prominent "ring" of HI self-absorption. Through a combined analysis of HI narrow self-absorption, CO emission, dust emission, and extinction, we directly measured, for the first time, the [HI]/[H$_2$] abundance varying from 2% to 0.2%, within one region. These measured HI abundances are orders of magnitude higher than usually assumed initial conditions for protoplanetary disk models. None of the fast cloud formation model could produce such low atomic hydrogen abundance. We derived a cloud formation timescale of 6$\times$10$^6$ years, consistent with the global Galactic star formation rate, and favoring the classical star formation picture over fast star formation models. Our measurements also help constrain the H$_2$ formation rate, under various ISM conditions.

In this work, we derive a spatially resolved map of the line-of-sight velocity of the interstellar medium and use it, along with a second map of line-of-sight velocity from Paper I of this series, to determine the nature of gaseous spiral structure in the Milky Way. This map is derived from measurements of the 1.527 $\mu$m diffuse interstellar band (DIB) in stellar spectra from the APOGEE survey and covers the nearest 4-5 kpc of the Northern Galactic plane. We cross-check this new DIB-based line-of-sight velocity map with the map derived in Paper I and find that they agree. We then compare these maps with line-of-sight velocity maps derived from simulations of quasi-stationary density wave spiral structure and dynamic, or material, spiral structure in a Milky Way-like galaxy. While none of the maps derived from these simulations is an exact match to the measured velocity field of the Milky Way, the measurements are more consistent with simulations of dynamic spiral structure. In the dynamic spiral structure simulation that best matches the measurements, the Perseus spiral arm is being disrupted.

The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.

The local interstellar medium (ISM) is suffused with "dark" gas, identified by excess infrared and gamma ray emission, yet undetected by standard ISM tracers such as neutral hydrogen (HI) or carbon monoxide emission. Based on observed dust properties from Planck, recent studies have argued that HI mixed with dust is strongly saturated and that dark gas is dominated by optically-thick HI. We test this hypothesis by reproducing this model using data from Planck and new 21 cm emission maps from GALFA-HI -- the first large-area 21cm emission survey with comparable angular resolution to Planck. We compare the results with those from a large sample of HI column densities based on direct observations of HI optical depth, and find that the inferred column density corrections are significantly lower than those inferred by the Planck-based model. Further, we rule out the hypothesis that the pencil-beam HI absorption sight lines preferentially miss opaque "blobs" with small covering fraction, as these structures require densities and pressures which are incompatible with ISM conditions. Our results support the picture that excess dust emission in the local ISM is not dominated by optically-thick HI, but is rather a combination of intrinsic changes in dust grain emissivities and H2 missed by CO observations.

We present an approach for comparing the detections and non-detections of Local Group (LG) dwarf galaxies in large HI surveys to the predictions of a suite of n-body simulations of the LG. This approach depends primarily on a set of empirical scaling relations to connect the simulations to the observations, rather than making strong theoretical assumptions. We then apply this methodology to the Galactic Arecibo L-band Feed Array HI Compact Cloud Catalog (GALFA-HI CCC), and compare it to the ELVIS suite of simulations. This approach reveals a strong tension between the na\"ive results of the model and the observations: while there are no LG dwarfs in the GALFA-HI CCC, the simulations predict $\sim 10$. Applying a simple model of reionization can resolve this tension by preventing low-mass halos from forming gas. However, and if this effect operates as expected, the observations provide a constraint on the mass scale of dwarf galaxy that reionization impacts. Combined with the observed properties of Leo T, the halo virial mass scale at which reionization impacts dwarf galaxy gas content is constrained to be $\sim 10^{8.5} M_\odot$, independent of any assumptions about star formation.

Standardizing and incentivizing the use of digital object identifiers (DOIs) to aggregate and identify both data analyzed and data generated by a research project will advance the field of astronomy to match best practices in other research fields like geosciences and medicine. Increase in the use of DOIs will prepare the discipline for changing expectations among funding agencies and publishers, who increasingly expect accurate and thorough data citation to accompany scientific outputs. The use of DOIs ensures a robust, sustainable, and interoperable approach to data citation in which due credit is given to researchers and institutions who produce and maintain the primary data. We describe in this work the advantages of DOIs for data citation and best practices for integrating a DOI service in an astronomical archive. We report on a pilot project carried out in collaboration with AAS Journals. During the course of the 1.5 year pilot, over 75% of submitting authors opted to use the integrated DOI service to clearly identify data analyzed during their research project when prompted at the time of paper submission.

We present calibrations for star formation rate indicators in the ultraviolet, mid-infrared and radio continuum bands, including one of the first direct calibrations of 150 MHz as a star formation rate indicator. Our calibrations utilize 66 nearby star forming galaxies with Balmer decrement corrected H-alpha luminosities, which span 5 orders of magnitude in star formation rate and have absolute magnitudes of -24<M_r<-12. Most of our photometry and spectrophotometry is measured from the same region of each galaxy, and our spectrophotometry has been validated with SDSS photometry, so our random and systematic errors are small relative to the intrinsic scatter seen in star formation rate indicator calibrations. We find WISE W4 (22.8 micron), Spitzer 24 micron and 1.4 GHz have tight correlations with Balmer decrement corrected H-alpha luminosity, with scatter of only 0.2 dex. Our calibrations are comparable to those from the prior literature for L* galaxies, but for dwarf galaxies our calibrations can give star formation rates that are far greater than those derived from much of the prior literature.

We report the detection of a largely ionized very-high velocity cloud (VHVC; $v_{\rm LSR}\sim-350$ km/s) toward M33 with the Hubble Space Telescope/Cosmic Origin Spectrograph. The VHVC is detected in OI, CII, SiII, and SiIII absorption along five sightlines separated by ~0.06-0.4 degree. On sub-degree scales, the velocities and ionic column densities of the VHVC remain relatively smooth with standard deviations of +/-14 km/s and +/-0.15 dex between the sightlines, respectively. The VHVC has a metallicity of [OI/HI]=-0.56+/-0.17 dex (Z=0.28+/-0.11 Z$_{\odot}$). Despite the position-velocity proximity of the VHVC to the ionized Magellanic Stream, the VHVC's higher metallicity makes it unlikely to be associated with the Stream, highlighting the complex velocity structure of this region of sky. We investigate the VHVC's possible origin by revisiting its surrounding HI environment. We find that the VHVC may be: (1) a MW CGM cloud, (2) related to a nearby HI VHVC -- Wright's Cloud, or (3) connected to M33's northern warp. Furthermore, the VHVC could be a bridge connecting Wright's Cloud and M33's northern warp, which would make it a Magellanic-like structure in the halo of M33.

Near-infrared spectroscopy from APOGEE and wide-field optical photometry from Pan-STARRS1 have recently made possible precise measurements of the shape of the extinction curve for tens of thousands of stars, parameterized by R(V). These measurements revealed structures in R(V) with large angular scales, which are challenging to explain in existing dust paradigms. In this work, we combine three-dimensional maps of dust column density with R(V) measurements to constrain the three-dimensional distribution of R(V) in the Milky Way. We find that variations in R(V) are correlated on kiloparsec scales. In particular, most of the dust within one kiloparsec in the outer Galaxy, including many local molecular clouds (Orion, Taurus, Perseus, California, Cepheus), has a significantly lower R(V) than more distant dust in the Milky Way. These results provide new input to models of dust evolution and processing, and complicate application of locally derived extinction curves to more distant regions of the Milky Way and to other galaxies.

We report the detection of accreting ionized gas at the disk-halo interface of the nearby galaxy M33. We analyze HST/COS absorption-line spectra of seven ultraviolet-bright stars evenly distributed across the disk of M33. We find Si IV absorption components consistently redshifted relative to the bulk M33's ISM absorption along all the sightlines. The Si IV detection indicates an enriched, disk-wide, ionized gas inflow toward the disk. This inflow is most likely multi-phase as the redshifted components can also be observed in ions with lower ionization states (e.g., S II, P II, Fe II, Si II). Kinematic modeling of the inflow is consistent with an accreting layer at the disk-halo interface of M33, which has an accretion velocity of 110$^{+15}_{-20}$ km s$^{-1}$ at a distance of 1.5$^{+1.0}_{-1.0}$ kiloparsec above the disk. The modeling indicates a total mass of $\sim3.9\times10^7$ M$_{\odot}$ for the accreting material at the disk-halo interface on the near side of the M33 disk , with an accretion rate of $\sim2.9$ M$_{\odot}$ yr$^{-1}$. The high accretion rate and the level of metal-enrichment suggest the inflow is likely to be the fall back of M33 gas from a galactic fountain and/or the gas pulled loosed during a close interaction between M31 and M33. Our study of M33 is the first to unambiguously reveal the existence of a disk-wide, ionized gas inflow beyond the Milky Way, providing a better understanding of gas accretion in the vicinity of a galaxy disk.

We have developed a method for deriving the distribution of the Milky Way's interstellar medium as a function of longitude, latitude, distance and line-of-sight velocity. This method takes as input maps of reddening as a function of longitude, latitude, and distance and maps of line emission as a function of longitude, latitude, and line-of-sight velocity. We have applied this method to datasets covering much of the Galactic plane. The output of this method correctly reproduces the line-of-sight velocities of high-mass star forming regions with known distances from Reid et al. (2014) and qualitatively agrees with results from the Milky Way kinematics literature. These maps will be useful for measuring flows of gas around the Milky Way's spiral arms and into and out of giant molecular clouds.

Neutral atomic hydrogen (HI) gas in interstellar space is largely organized into filaments, loops, and shells, the most prominent of which are "supershells". These gigantic structures requiring $\gtrsim 3 \times 10^{52}$ erg to form are generally thought to be produced by either the explosion of multiple supernovae (SNe) in OB associations or alternatively by the impact of high-velocity clouds (HVCs) falling to the Galactic disk. Here we report the detection of a kiloparsec (kpc)-size supershell in the outskirts of the Milky Way with the compact HVC 040+01$-$282 (hereafter CHVC040) at its geometrical center using the "Inner-Galaxy Arecibo L-band Feed Array" HI 21-cm survey data. The morphological and physical properties of both objects suggest that CHVC040, which is either a fragment of a nearby disrupted galaxy or a cloud originated from an intergalactic accreting flow, collided with the disk $\sim 5$ Myrs ago to form the supershell. Our result shows that some compact HVCs can survive their trip through the Galactic halo and inject energy and momentum into the Milky Way disk.

We present medium-resolution, near-ultraviolet VLT/FLAMES observations of the star USNO-A0600-15865535. We adapt a standard method of stellar typing to our measurement of the shape of the Balmer epsilon absorption line to demonstrates that USNO-A0600-15865535 is a blue horizontal branch star, residing in the lower stellar halo at a distance of 4.4 kpc from the Sun. We measure the H & K lines of singly-ionized calcium and find two isolated velocity components, one originating in the disk, and one associated with high-velocity cloud complex WD. This detection demonstrated that complex WD is closer than ~4.4 kpc and is the first distance constraint on the +100 km/s Galactic complex of clouds. We find that Complex WD is not in corotation with the Galactic disk as has been assumed for decades. We examine a number of scenarios, and find that the most likely is that Complex WD was ejected from the solar neighborhood and is only a few kpc from the Sun.

We present a proof-of-concept study of a method to estimate the inclination angle of compact high velocity clouds (CHVCs), i.e. the angle between a CHVC's trajectory and the line-of-sight. The inclination angle is derived from the CHVC's morphology and kinematics. We calibrate the method with numerical simulations, and we apply it to a sample of CHVCs drawn from HIPASS. Implications for CHVC distances are discussed.

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity.

Using high-resolution data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey, we show that linear structure in Galactic neutral hydrogen (HI) correlates with the magnetic field orientation implied by Planck 353 GHz polarized dust emission. The structure of the neutral interstellar medium is more tightly coupled to the magnetic field than previously known. At high Galactic latitudes, where the Planck data are noise-dominated, the HI data provide an independent constraint on the Galactic magnetic field orientation, and hence the local dust polarization angle. We detect strong cross-correlations between template maps constructed from estimates of dust intensity combined with either HI-derived angles, starlight polarization angles, or Planck 353 GHz angles. The HI data thus provide a new tool in the search for inflationary gravitational wave B-mode polarization in the cosmic microwave background, which is currently limited by dust foreground contamination.

We present ultraviolet (UV) follow-up of a sample of potential dwarf galaxy candidates selected for their neutral hydrogen (HI) properties, taking advantage of the low UV background seen by the GALEX satellite and its large and publicly available imaging footprint. The HI clouds, which are drawn from published GALFA-HI and ALFALFA HI survey compact cloud catalogs, are selected to be galaxy candidates based on their spatial compactness and non-association with known high-velocity cloud complexes or Galactic HI emission. Based on a comparison of their UV characteristics to those of known dwarf galaxies, half (48%) of the compact HI clouds have at least one potential stellar counterpart with UV properties similar to those of nearby dwarf galaxies. If galaxies, the star formation rates, HI masses, and star formation efficiencies of these systems follow the trends seen for much larger galaxies. The presence of UV emission is an efficient method to identify the best targets for spectroscopic follow-up, which is necessary to prove that the stars are associated with the compact HI. Further, searches of this nature help to refine the salient HI properties of likely dwarfs (even beyond the Local Group). In particular, HI compact clouds considered to be velocity outliers relative to their neighbor HI clouds have the most significant detection rate of single, appropriate UV counterparts. Correcting for the sky coverage of the two all-Arecibo sky surveys yielding the compact HI clouds, these results may imply the presence of potentially hundreds of new tiny galaxies across the entire sky.

We present a catalog of 1964 isolated, compact neutral hydrogen clouds from the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI DR1). The clouds were identified by a custom machine-vision algorithm utilizing Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds have velocities typically between |VLSR| = 20-400 km/s, linewidths of 2.5-35 km/s, and column densities ranging from 1 - 35 x 10^18 cm^-2. The distances to the clouds in this catalog may cover several orders of magnitude, so the masses may range from less than a Solar mass for clouds within the Galactic disc, to greater than 10^4 Solar Masses for HVCs at the tip of the Magellanic Stream. To search for trends, we separate the catalog into five populations based on position, velocity, and linewidth: high velocity clouds (HVCs); galaxy candidates; cold low velocity clouds (LVCs); warm, low positive-velocity clouds in the third Galactic Quadrant; and the remaining warm LVCs. The observed HVCs are found to be associated with previously-identified HVC complexes. We do not observe a large population of isolated clouds at high velocities as some models predict. We see evidence for distinct histories at low velocities in detecting populations of clouds corotating with the Galactic disc and a set of clouds that is not corotating.

We assess the fraction of the Milky Way's circumgalactic medium (CGM) eluding detection due to its velocity being similar to gas in the disk. This is achieved using synthetic observations of the CGM in a simulated MW-mass galaxy that shows similar CGM kinematics to the MW and external L$\sim$L$_*$ galaxies. As viewed by a mock observer at a location similar to the Sun, only 50$\%$ (by mass) of the gas moves at high velocity ($|v_{\rm LSR}|\geq$100 km s$^{-1}$ or $|v_{\rm DEV}|\geq$50 km s$^{-1}$) in the simulated CGM and would be observable. The low velocity gas is thermodynamically similar to the high velocity gas, indicating the 50$\%$ observable fraction is applicable to spectral lines from the radio to the ultraviolet. We apply the observable mass fraction (50$\%$) to current estimates of the MW's CGM, and find a corrected total mass of 2.8$\times$10$^{8} M_{\odot}$ for gas below 10$^6$K within $\sim15$ kpc (this excludes the Magellanic System). This is less than the total mass of the CGM extending out to $\sim$150 kpc in other L$\sim$L$_*$ galaxies. However, we find similar OVI column densities when the discrepancy in path length between the MW and external galaxies is considered. The coherent spatial and kinematic distribution of low velocity gas in the simulated CGM suggests that current HI observations of the MW's CGM may miss large low velocity HI complexes. In addition, current mass estimates of the MW's CGM based on high-velocity line observations with distance constraints may miss a non-negligible fraction of gas in the outer halo which can be obscured if it moves at a velocity similar to the gas in the lower halo.

As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more detail than anywhere else in the Universe how galaxies acquire fresh gas to fuel their continuing star formation, how they exchange gas with the surrounding intergalactic medium, and turn warm, diffuse gas into molecular clouds and ultimately stars. The $\lambda$21-cm line of atomic hydrogen (HI) is an excellent tracer of these physical processes. With the SKA we will finally have the combination of surface brightness sensitivity, point source sensitivity and angular resolution to transform our understanding of the evolution of gas in the Milky Way, all the way from the halo down to the formation of individual molecular clouds.

Using spectroscopically selected galaxies from the Sloan Digital Sky Survey we present a detection of reddening due to dust in the circumgalactic medium of galaxies. We detect the mean change in the colors of "standard crayons" correlated with the presence of foreground galaxies at z ~0.05 as a function of angular separation. Following Peek & Graves (2010), we create standard crayons using passively evolving galaxies corrected for Milky Way reddening and color-redshift trends, leading to a sample with as little as 2% scatter in color. We devise methods to ameliorate possible systematic effects related to the estimation of colors, and we find an excess reddening induced by foreground galaxies at a level ranging from 10 to 0.5 millimagnitudes on scales ranging from 30 kpc to 1 Mpc. We attribute this effect to a large-scale distribution of dust around galaxies similar to the findings of Menard et al. 2010. We find that circumgalactic reddening is a weak function of stellar mass over the range $6 \times 10^9 M_\odot$ -- $6 \times 10^{10} M_\odot$ and note that this behavior appears to be consistent with recent results on the distribution of metals in the gas phase.

We present observations of a new group of structures in the diffuse Galactic ISM: slender, linear HI features we dub "fibers" that extend for many degrees at high Galactic latitude. To characterize and measure the extent and strength of these fibers, we present the Rolling Hough Transform (RHT), a new machine vision method for parameterizing the coherent linearity of structures in the image plane. With this powerful new tool we show the fibers are oriented along the interstellar magnetic field as probed by starlight polarization. We find that these low column density (N(HI) ~ 5 x 10^18 cm^-2) fiber features are most likely a component of the local cavity wall, about 100 pc away. The HI data we use to demonstrate this alignment at high latitude are from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey and the Parkes Galactic All Sky Survey (GASS). We find better alignment in the higher resolution GALFA-HI data, where the fibers are more visually evident. This trend continues in our investigation of magnetically aligned linear features in the Riegel-Crutcher HI cold cloud, detected in the Southern Galactic Plane Survey (SGPS). We propose an application of the RHT for estimating the field strength in such a cloud, based on the Chandrasekhar-Fermi method. We conclude that data-driven, quantitative studies of ISM morphology can be very powerful predictors of underlying physical quantities.

We have searched for infrared dust emission from subsets of compact, Galactic neutral hydrogen clouds, with the purpose of looking for dust in high-velocity clouds, identifying low-velocity halo clouds, and investigating the cloud populations defined in the GALFA-HI Compact Cloud Catalog. We do not detect dust emission from high-velocity clouds. The lack of dust emission from a group of low-velocity clouds supports the claim that they are low-velocity halo clouds. We detect dust in the remaining low-velocity clouds, indicating a Galactic origin, with a significantly greater dust-to-gas ratio for clouds with linewidths near 15 km/s. We propose that this is due to dust associated with ionized gas.

We present an atlas of 129 spectral energy distributions for nearby galaxies, with wavelength coverage spanning from the UV to the mid-infrared. Our atlas spans a broad range of galaxy types, including ellipticals, spirals, merging galaxies, blue compact dwarfs and luminous infrared galaxies. We have combined ground-based optical drift-scan spectrophotometry with infrared spectroscopy from Spitzer and Akari, with gaps in spectral coverage being filled using MAGPHYS spectral energy distribution models. The spectroscopy and models were normalized, constrained and verified with matched-aperture photometry measured from Swift, GALEX, SDSS, 2MASS, Spitzer and WISE images. The availability of 26 photometric bands allowed us to identify and mitigate systematic errors present in the data. Comparison of our spectral energy distributions with other template libraries and the observed colors of galaxies indicates that we have smaller systematic errors than existing atlases, while spanning a broader range of galaxy types. Relative to the prior literature, our atlas will provide improved K-corrections, photometric redshifts and star-formation rate calibrations.

We present a dust-column--dependent extinction curve parameters for ultraviolet wavelengths at high Galactic latitudes. This extinction function diverges from previous work in that it takes into account the results of Peek & Schiminovich 2013 (Paper I), which demonstrated that there is more reddening in the GALEX bands than would be otherwise expected for E(B-V) < 0.2. We also test the biases in the Planck and SFD extinction maps, and find that the SFD extinction maps are significantly biased at E(B-V) < 0.2. We find that while an extinction function that that takes into account a varying R_FUV with E(B-V) dramatically improves our estimation of FUV-NUV colors, a fit that also includes HI column density dependence is superior. The ultraviolet extinction function we present here follows the model of Fitzpatrick 1999, varying only the amplitude of the FUV rise parameter to be consistent with the data.

We investigate turbulent properties of the non-star-forming, translucent molecular cloud, MBM16 by applying the statistical technique of a two-dimensional spatial power spectrum (SPS) on the neutral hydrogen (HI) observations obtained by the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey. The SPS, calculated over the range of spatial scales from 0.1 to 17 pc, is well represented with a single power-law function, with a slope ranging from -3.3 to -3.7 and being consistent over the velocity range of MBM16 for a fixed velocity channel thickness. However, the slope varies significantly with the velocity slice thickness, suggesting that both velocity and density contribute to HI intensity fluctuations. By using this variation we estimate the slope of 3D density fluctuations in MBM16 to be -3.7\pm0.2. This is significantly steeper than what has been found for HI in the Milky Way plane, the Small Magellanic Cloud, or the Magellanic Bridge, suggesting that interstellar turbulence in MBM16 is driven on scales >17 pc and that the lack of stellar feedback could be responsible for the steep power spectrum.

We search for fast-expanding Hi shells associated with Galactic supernova remnants (SNRs) in the longitude range l \approx 32\arcdeg to 77\arcdeg using 21-cm line data from the Inner-Galaxy Arecibo L-band Feed Array (I-GALFA) Hi survey. Among the 39 known Galactic SNRs in this region, we find such Hi shells in four SNRs: W44, G54.4-0.3, W51C, and CTB 80. All four were previously identified in low-resolution surveys, and three of those (excluding G54.4-0.3) were previously studied with the Arecibo telescope. A remarkable new result, however, is the detection of Hi emission at both very high positive and negative velocities in W44 from the receding and approaching parts of the Hi expanding shell, respectively. This is the first detection of both sides of an expanding shell associated with an SNR in Hi 21-cm emission. The high-resolution I-GALFA survey data also reveal a prominent expanding Hi shell with high circular symmetry associated with G54.4-0.3. We explore the physical characteristics of four SNRs and discuss what differentiates them from other SNRs in the survey area. We conclude that these four SNRs are likely the remnants of core-collapse supernovae interacting with a relatively dense (>~1 cm^-3) ambient medium, and we discuss the visibility of SNRs in the Hi 21-cm line.

In order to study the properties and effects of high Galactic latitude dust we present an analysis of 373,303 galaxies selected from the Galaxy Evolution Explorer (GALEX) All-Sky Survey and Wide-Field Infrared Explorer (WISE) All-Sky Data Release. By examining the variation in aggregate ultraviolet colors and number density of these galaxies we measure the extinction curve at high latitude. We additionally consider a population of spectroscopically selected galaxies from the Sloan Digital Sky Survey (SDSS) to measure extinction in the optical. We find that dust at high latitude is neither quantitatively nor qualitatively consistent with standard reddening laws. Extinction in the FUV and NUV is ~10% and ~35% higher than expected, with significant variation across the sky. We find that no single R_V parameter fits both the optical and ultraviolet extinction at high latitude, and that while both show detectable variation across the sky, these variations are not related. We propose that the overall trends we detect likely stem from an increase in very small silicate grains in the ISM.

Both the Galactic 21-cm line flux from neutral hydrogen (HI) in interstellar medium and the far-infrared (FIR) emission from Galactic dust grains have been used to estimate the strength of Galactic reddening of distant sources. In this work we use a collection of uniform color distant galaxies as color standards to determine whether the HI method or the FIR method is superior. We find that the two methods both produce reasonably accurate maps, but that both show significant errors as compared to the typical color of the background galaxies. We find that a mixture of the FIR and HI maps in roughly a 2-to-1 ratio is clearly superior to either map alone. We recommend that future reddening maps should use both sets of data, and that well-constructed FIR and HI maps should both be vigorously pursued.

Galactic halo gas traces inflowing star formation fuel and feedback from a galaxy's disk and is therefore crucial to our understanding of galaxy evolution. In this review, we summarize the multi-wavelength observational properties and origin models of Galactic and low redshift spiral galaxy halo gas. Galactic halos contain multiphase gas flows that are dominated in mass by the ionized component and extend to large radii. The densest, coldest halo gas observed in neutral hydrogen (HI) is generally closest to the disk (< 20 kpc), and absorption line results indicate warm and warm-hot diffuse halo gas is present throughout a galaxy's halo. The hot halo gas detected is not a significant fraction of a galaxy's baryons. The disk-halo interface is where the multiphase flows are integrated into the star forming disk, and there is evidence for both feedback and fueling at this interface from the temperature and kinematic gradient of the gas and HI structures. The origin and fate of halo gas is considered in the context of cosmological and idealized local simulations. Accretion along cosmic filaments occurs in both a hot (> 10^5.5 K) and cold mode in simulations, with the compressed material close to the disk the coldest and densest, in agreement with observations. There is evidence in halo gas observations for radiative and mechanical feedback mechanisms, including escaping photons from the disk, supernova-driven winds, and a galactic fountain. Satellite accretion also leaves behind abundant halo gas. This satellite gas interacts with the existing halo medium, and much of this gas will become part of the diffuse halo before it can reach the disk. The accretion rate from cold and warm halo gas is generally below a galaxy disk's star formation rate, but gas at the disk-halo interface and stellar feedback may be important additional fuel sources.

A survey of the Milky Way disk and the Magellanic System at the wavelengths of the 21-cm atomic hydrogen (HI) line and three 18-cm lines of the OH molecule will be carried out with the Australian Square Kilometre Array Pathfinder telescope. The survey will study the distribution of HI emission and absorption with unprecedented angular and velocity resolution, as well as molecular line thermal emission, absorption, and maser lines. The area to be covered includes the Galactic plane (|b|< 10deg) at all declinations south of delta = +40deg, spanning longitudes 167deg through 360deg to 79deg at b=0deg, plus the entire area of the Magellanic Stream and Clouds, a total of 13,020 square degrees. The brightness temperature sensitivity will be very good, typically sigma_T ~ 1 K at resolution 30arcsec and 1 km/s. The survey has a wide spectrum of scientific goals, from studies of galaxy evolution to star formation, with particular contributions to understanding stellar wind kinematics, the thermal phases of the interstellar medium, the interaction between gas in the disk and halo, and the dynamical and thermal states of gas at various positions along the Magellanic Stream.

We perform high-resolution hydrodynamic simulations of a Milky Way-mass galaxy in a fully cosmological setting using the adaptive mesh refinement code, Enzo, and study the kinematics of gas in the simulated galactic halo. We find that the gas inflow occurs mostly along filamentary structures in the halo. The warm-hot (10^5 K < T < 10^6 K) and hot (T > 10^6 K) ionized gases are found to dominate the overall mass accretion in the system (with dM/dt = 3-5 M_solar/yr) over a large range of distances, extending from the virial radius to the vicinity of the disk. Most of the inflowing gas (by mass) does not cool, and the small fraction that manages to cool does so primarily close to the galaxy (R <~ 20 kpc), perhaps comprising the neutral gas that may be detectable as, e.g., high-velocity clouds. The neutral clouds are embedded within larger, accreting filamentary flows, and represent only a small fraction of the total mass inflow rate. The inflowing gas has relatively low metallicity (Z/Z_solar < 0.2). The outer layers of the filamentary inflows are heated due to compression as they approach the disk. In addition to the inflow, we find high-velocity, metal-enriched outflows of hot gas driven by supernova feedback. Our results are consistent with observations of halo gas at low z.

To investigate the fundamental principles of H2 formation in a giant molecular cloud (GMC), we derive the HI and H2 surface density (Sigma_HI and Sigma_H2) images of the Perseus molecular cloud on sub-pc scales (~0.4 pc). We use the far-infrared data from the Improved Reprocessing of the IRAS Survey and the V-band extinction image provided by the COMPLETE Survey to estimate the dust column density image of Perseus. In combination with the HI data from the Galactic Arecibo L-band Feed Array HI Survey and an estimate of the local dust-to-gas ratio, we then derive the Sigma_H2 distribution across Perseus. We find a relatively uniform Sigma_HI ~ 6-8 Msun pc^-2 for both dark and star-forming regions, suggesting a minimum HI surface density required to shield H2 against photodissociation. As a result, a remarkably tight and consistent relation is found between Sigma_H2/Sigma_HI and Sigma_HI+Sigma_H2. The transition between the HI- and H2-dominated regions occurs at N(HI)+2N(H2) ~ (8-14) x 10^20 cm^-2. Our findings are consistent with predictions for H2 formation in equilibrium, suggesting that turbulence may not be of primary importance for H2 formation. However, the importance of a warm neutral medium for H2 shielding, an internal radiation field, and the timescale of H2 formation still remain as open questions. We also compare H2 and CO distributions and estimate the fraction of "CO-dark" gas, f_DG ~ 0.3. While significant spatial variations of f_DG are found, we do not find a clear correlation with the mean V-band extinction.

Using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope, we have obtained high-resolution ultraviolet spectra of the C I absorption toward two stars behind the Local Leo Cold Cloud (LLCC). At a distance (~20 pc) that places it well inside the Local Bubble, the LLCC is the nearest example of the coldest known (T~20 K) diffuse interstellar clouds. The STIS measurements of the C I fine-structure excitation toward HD 85259 and HD 83023 indicate that the thermal gas pressure of the LLCC is much greater than that of the warm clouds in the Local Bubble. The mean LLCC pressure measured toward these two stars (60,000 cm$^{-3}$ K) implies an H I density of ~3000 cm$^{-3}$ and a cloud thickness of ~200 AU at the 20 K cloud temperature. Such a thin, cold, dense structure could arise at the collision interface between converging flows of warm gas. However, the measured LLCC pressure is appreciably higher than that expected in the colliding cloud interpretation given the velocity and column density constraints on warm clouds in the HD 85259 and HD 83023 sightlines. Additional STIS measurements of the Zn II, Ni II, and Cr II column densities toward HD 85259 indicate that the LLCC has a modest "warm cloud" dust depletion pattern consistent with its low dust-to-gas ratio determined from H I 21 cm and 100 micron observations. In support of the inferred sheet-like geometry for the LLCC, a multi-epoch comparison of the Na I absorption toward a high-proper-motion background star reveals a 40% column density variation indicative of LLCC Na I structure on a scale of ~50 AU.

We present a multi-wavelength study of the local Leo cold cloud (LLCC), a very nearby, very cold cloud in the interstellar medium. Through stellar absorption studies we find that the LLCC is between 11.3 pc and 24.3 pc away, making it the closest known cold neutral medium cloud and well within the boundaries of the local cavity. Observations of the cloud in the 21-cm HI line reveal that the LLCC is very cold, with temperatures ranging from 15 K to 30 K, and is best fit with a model composed of two colliding components. The cloud has associated 100 micron thermal dust emission, pointing to a somewhat low dust-to-gas ratio of 48 x 10^-22 MJy sr^-1 cm^2. We find that the LLCC is too far away to be generated by the collision among the nearby complex of local interstellar clouds, but that the small relative velocities indicate that the LLCC is somehow related to these clouds. We use the LLCC to conduct a shadowing experiment in 1/4 keV X-rays, allowing us to differentiate between different possible origins for the observed soft X-ray background. We find that a local hot bubble model alone cannot account for the low-latitude soft X-ray background, but that isotropic emission from solar wind charge exchange does reproduce our data. In a combined local hot bubble and solar wind charge exchange scenario, we rule out emission from a local hot bubble with an 1/4 keV emissivity greater than 1.1 Snowdens / pc at 3 sigma, 4 times lower than previous estimates. This result dramatically changes our perspective on our local interstellar medium.

We present the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey, and its first full data release (DR1). GALFA-HI is a high resolution (~ 4'), large area (13000 deg^2), high spectral resolution (0.18 km/s), wide band (-700 < v_LSR < +700 km/s) survey of the Galactic interstellar medium in the 21-cm line hyperfine transition of neutral hydrogen conducted at Arecibo Observatory. Typical noise levels are 80 mK RMS in an integrated 1 km/s channel. GALFA-HI is a dramatic step forward in high-resolution, large-area Galactic HI surveys, and we compare GALFA-HI to past, present, and future Galactic HI surveys. We describe in detail new techniques we have developed to reduce these data in the presence of fixed pattern noise, gain variation, and inconsistent beam shapes, and we show how we have largely mitigated these effects. We present our first full data release, covering 7520 square degrees of sky and representing 3046 hours of integration time, and discuss the details of these data.

Observations from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey of the tail of Complex C are presented and the halo clouds associated with this complex cataloged. The properties of the Complex C clouds are compared to clouds cataloged at the tail of the Magellanic Stream to provide insight into the origin and destruction mechanism of Complex C. Magellanic Stream and Complex C clouds show similarities in their mass distributions (slope = -0.7 and -0.6, respectively) and have a common linewidth of 20 - 30 km/s (indicative of a warm component), which may indicate a common origin and/or physical process breaking down the clouds. The clouds cataloged at the tail of Complex C extend over a mass range of 10^1.1 to 10^4.8 solar masses, sizes of 10^1.2 to 10^2.6 pc, and have a median volume density of 0.065 cm^(-3) and median pressure of (P/k) = 580 K cm^{-3}. We do not see a prominent two-phase structure in Complex C, possibly due to its low metallicity and inefficient cooling compared to other halo clouds. From assuming the Complex C clouds are in pressure equilibrium with a hot halo medium, we find a median halo density of 5.8 x 10^(-4) cm^(-3), which given a constant distance of 10 kpc, is at a z-height of ~3 kpc. Using the same argument for the Stream results in a median halo density of 8.4 x 10^(-5) x (60kpc/d) cm^(-3). These densities are consistent with previous observational constraints and cosmological simulations. We also assess the derived cloud and halo properties with three dimensional grid simulations of halo HI clouds and find the temperature is generally consistent within a factor of 1.5 and the volume densities, pressures and halo densities are consistent within a factor of 3.

The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and an angular resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of many isolated, very compact HI clouds at low Galactic velocities which are distinctly separated from the HI disk emission. In the limited area of ~4600 deg$^2$ searched so far, we have detected 96 such compact clouds. The detected clouds are cold with kinetic temperature less than 300 K. Moreover, they are quite compact and faint, with median values of 5 arcmin in angular size, 0.75 K in peak brightness temperature, and $5\times10^{18}$ cm$^{-2}$ in HI column density. From the modeling of spatial and velocity distributions of the whole compact cloud population, we find that the bulk of clouds are related to the Galactic disk, and are within a few kpc distance. We present properties of the compact clouds sample and discuss various possible scenarios for the origin of this clouds population and its role in the Galactic interstellar medium studies.

The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and a spatial resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of numerous isolated, very compact HI clouds at low Galactic velocities, which are distinctly separated from the HI disk emission. In the limited area of ~4600 deg$^2$ surveyed so far, we have detected 96 of such compact clouds. The detected clouds are cold with a median T$_{k,max}$ (the kinetic temperature in the case in which there is no non-thermal broadening) of 300 K. Moreover, these clouds are quite compact and faint, with median values of 5 arcmin in angular size, 0.75 K in peak brightness temperature, and $5 \times 10^{18}$ cm$^{-2}$ in HI column density. Most of the clouds deviate from Galactic rotation at the 20-30 km/s level, and a significant fraction show evidence for a multiphase medium and velocity gradients. No counterparts for these clouds were found in other wavebands. From the modeling of spatial and velocity distributions of the whole compact cloud population, we find that the bulk of the compact clouds are related to the Galactic disk, and their distances are likely to be in the range of 0.1 to a few kpc. We discuss various possible scenarios for the formation and maintenance of this cloud population and its significance for Galactic ISM studies.

We present corrections to the Schlegel, Finkbeiner, Davis (SFD98) reddening maps over the Sloan Digital Sky Survey northern Galactic cap area. To find these corrections, we employ what we dub the "standard crayon" method, in which we use passively evolving galaxies as color standards by which to measure deviations from the reddening map. We select these passively evolving galaxies spectroscopically, using limits on the H alpha and O II equivalent widths to remove all star-forming galaxies from the SDSS main galaxy catalog. We find that by correcting for known reddening, redshift, color-magnitude relation, and variation of color with environmental density, we can reduce the scatter in color to below 3% in the bulk of the 151,637 galaxies we select. Using these galaxies we construct maps of the deviation from the SFD98 reddening map at 4.5 degree resolution, with 1-sigma error of ~ 1.5 millimagnitudes E(B-V). We find that the SFD98 maps are largely accurate with most of the map having deviations below 3 millimagnitudes E(B-V), though some regions do deviate from SFD98 by as much as 50%. The maximum deviation found is 45 millimagnitudes in E(B-V), and spatial structure of the deviation is strongly correlated with the observed dust temperature, such that SFD98 underpredicts reddening in regions of low dust temperature. Our maps of these deviations, as well as their errors, are made available to the scientific community as supplemental correction to SFD98 at http://www.peekandgraves2010.com.

In this paper we present the results of the statistical analysis of high-latitude HI turbulence in the Milky Way. We have observed HI in the 21 cm line, obtained with the Arecibo L-Band Feed Array (ALFA) receiver at the Arecibo radio telescope. For recovering of velocity statistics we have used the Velocity Coordinate Spectrum (VCS) technique. In our analysis we have used direct fitting of the VCS model, as its asymptotic regimes are questionable for Arecibo's resolution and given the restrictions from thermal smoothing of the turbulent line. We have obtained a velocity spectral index $3.87 \pm 0.11$, an injection scale of $140 \pm 80$ pc, and an HI cold phase temperature of $52 \pm 11$ K. The spectral index is steeper than the Kolmogorov index and can be interpreted as being due to shock-dominated turbulence.

The disruption of the M33 galaxy is evident from its extended gaseous structure. We present new data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey that show the full extent and detailed spatial and kinematic structure of M33's neutral hydrogen. Over 18% of the HI mass of M33 (M_HI(tot) =1.4 x 10^9 Msun) is found beyond the star forming disk as traced in the far-ultraviolet (FUV). The most distinct features are extended warps, an arc from the northern warp to the disk, diffuse gas surrounding the galaxy, and a southern cloud with a filament back to the galaxy. The features extend out to 22 kpc from the galaxy center (18 kpc from the edge of the FUV disk) and the gas is directly connected to M33's disk. Only five discrete clouds (i.e., gas not directly connected to M33 in position-velocity space) are catalogued in the vicinity of M33, and these clouds show similar properties to Galactic and M31 halo clouds. M33's gaseous features most likely originate from the tidal disruption of M33 by M31 1-3 Gyr ago as shown from an orbit analysis which results in a tidal radius < 15 kpc in the majority of M33's possible orbits. M33 is now beyond the disruptive gravitational influence of M31 and the gas appears to be returning to M33's disk and redistributing its star formation fuel. M33's high mean velocity dispersion in the disk (~18.5 km/s) may also be consistent with the previous interaction and high rate of star formation. M33 will either exhaust its star formation fuel in the next few Gyrs or eventually become star formation fuel for M31. The latter represents the accretion of a large gaseous satellite by a spiral galaxy, similar to the Magellanic Clouds' relationship to the Galaxy.

The I-GALFA survey is mapping all the HI in the inner Galactic disk visible to the Arecibo 305m telescope within 10 degrees of the Galactic plane (longitudes of 32 to 77 deg at b=0 deg). The survey, which will obtain 1.3 million independent spectra, became possible with the installation of the 7-beam Arecibo L-Band Feed Array (ALFA) receiver in 2004. ALFA's 3.4 arcmin resolution and tremendous sensitivity offer a great opportunity to observe the fine details of HI in the Galaxy. The I-GALFA survey began in 2008 May and will be completed in 2009 September. The data will be made publicly available when the calibrated and gridded cubes are completed. Further information on the I-GALFA project may be found at www.naic.edu/~galfa.

We discuss the support for the cooling of gas directly at the disk-halo interface as a mechanism to continually fuel a galaxy. This may be an important reservoir as there is not enough cold gas observed in galaxy halos and simulations indicate the existing cold clouds will be rapidly destroyed as they move through the surrounding halo medium. We show possible evidence for a net infall of the WIM layer in the Milky Way, simulation results showing the recooling of warm clouds at the disk-halo interface, and GALFA HI data of small, cold HI clouds that could represent this recooling.

Ongoing star formation in the Milky Way requires continuous gaseous fuel from accretion. Previous work has suggested that the accretion of dwarf galaxies could provide the needed gas for this process. In this work we investigate whether dwarf galaxy accretion is consistent with the radial profile of star formation observed in the Milky Way, which is strongly concentrated toward the center of the gaseous disk of the Galaxy. Using realistic parameters for the Galactic potential, gaseous halo, Galactic disk, velocities of dwarf galaxies, and effects of drag on stripped gas in the halo, we model the distribution of radii at which dwarf galaxies accrete onto the Galactic disk. We also model the radial distribution of the accretion of gas that cools directly out of the halo by examining the results of recent simulations. We show that dwarf galaxies cannot reproduce the concentration of accretion toward the center of the Galaxy required by star formation. We also show that clouds that cool directly from the halo can reproduce this central concentration, and conclude that this difference is largely due to the discrepancy in absolute specific angular momentum between the two mechanisms.

We present 21-cm observations of four Galactic globular clusters, as part of the on-going GALFA-HI Survey at Arecibo. We discovered a peculiar HI cloud in the vicinity of the distant (109 kpc) cluster Pal 4, and discuss its properties and likelihood of association with the cluster. We conclude that an association of the HI cloud and Pal 4 is possible, but that a chance coincidence between Pal 4 and a nearby compact high-velocity cloud cannot be ruled out altogether. New, more stringent upper limits were derived for the other three clusters: M 3, NGC 5466, and Pal 13. We briefly discuss the fate of globular cluster gas and the interaction of compact clouds with the Galactic Halo gas.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

Models that reproduce the observed high-velocity clouds (HVCs) also predict clouds at lower radial velocities that may easily be confused with Galactic disk (|z| < 1 kpc) gas. We describe the first search for these low-velocity halo clouds (LVHCs) using IRAS data and the initial data from the Galactic Arecibo L-band Feed Array survey in HI (GALFA-HI). The technique is based upon the expectation that such clouds should, like HVCs, have very limited infrared thermal dust emission as compared to their HI column density. We describe our 'displacement-map' technique for robustly determining the dust-to-gas ratio of clouds and the associated errors that takes into account the significant scatter in the infrared flux from the Galactic disk gas. We find that there exist lower-velocity clouds that have extremely low dust-to-gas ratios, consistent with being in the Galactic halo - candidate LVHCs. We also confirm the lack of dust in many HVCs with the notable exception of complex M, which we consider to be the first detection of warm dust in HVCs. We do not confirm the previously reported detection of dust in complex C. In addition, we find that most Intermediate- and Low-Velocity clouds that are part of the Galactic disk have a higher 60 micron/100 micron flux ratio than is typically seen in Galactic HI, which is consistent with a previously proposed picture in which fast-moving Galactic clouds have smaller, hotter dust grains.

We explain the entire process by which we conduct the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey. The survey is a high resolution (3.4'), large area (13000 deg^2), large Galactic velocity range (-750 to +750 km/s), high spectral resolution (0.18 km/s) survey of the Galaxy in the 21 cm line hyperfine transition of hydrogen conducted at Arecibo Observatory. We touch on some of the new Galactic science being conducted using the GALFA-HI survey, ranging from High-Velocity Clouds to HI narrow-line self-absorption. We explain the many technical challenges that confront such a survey, including baseline ripple, gain variation and asymmetrical beam shapes. To correct for these systematic effects we use various newly developed methods, which we describe in detail. We also explain the data reduction process step by step, starting with the raw time-ordered data and ending with fully calibrated maps. The effects of each step of the data reduction on the final data product is shown sequentially. We conclude with future directions for the ongoing survey.

Radio observations of gas in the Milky Way and Local Group are vital for understanding how galaxies function as systems. The unique sensitivity of Arecibo's 305m dish, coupled with the 7-beam Arecibo L-Band Feed Array (ALFA), provides an unparalleled tool for investigating the full range of interstellar phenomena traced by the HI 21cm line. The GALFA (Galactic ALFA) HI Survey is mapping the entire Arecibo sky over a velocity range of -700 to +700 km/s with 0.2 km/s velocity channels and an angular resolution of 3.4 arcminutes. We present highlights from the TOGS (Turn on GALFA Survey) portion of GALFA-HI, which is covering thousands of square degrees in commensal drift scan observations with the ALFALFA and AGES extragalactic ALFA surveys. This work is supported in part by the National Astronomy and Ionosphere Center, operated by Cornell University under cooperative agreement with the National Science Foundation.

We report an accurate distance of d = 10+/-2.5kpc to the high-velocity cloud Complex C. Using high signal-to-noise Keck/HIRES spectra of two horizontal-branch stars, we have detected CaII K absorption lines from the cloud. Significant non-detections toward a further 3 stars yield robust lower distance limits. The resulting HI mass of Complex C is 4.9^{+2.8}_{-2.2} x 10^6 Msun; a total mass of 8.2^{+4.6}_{-2.6} x 10^6 Msun is implied, after corrections for helium and ionization. At 10kpc, Complex C has physical dimensions 3x15 kpc, and if it is as thick as it is wide, then the average density is log<n> ~ -2.5. We estimate the contribution of Complex C to the mass influx may be as high as ~0.14 Msun/yr.

Halo clouds have been found about the three largest galaxies of the Local Group and in the halos of nearby spirals. This suggests they are a relatively generic feature of the galaxy evolution process and a source of fuel for galaxy disks. In this review, two main sources of disk star formation fuel, satellite material and clouds condensing from the hot halo medium, are discussed and their contribution to fueling the Galaxy quantified. The origin of the halo gas of M31 and M33 is also discussed.

We present results from neutral hydrogen (HI) observations of the tip of the Magellanic Stream (MS), obtained with the Arecibo telescope as a part of the on-going survey by the Consortium for Galactic studies with the Arecibo L-band Feed Array. We find four large-scale, coherent HI streams, extending continously over a length of 20 degrees, each stream possessing different morphology and velocity gradients. The newly discovered streams provide strong support for the tidal model of the MS formation by Connors et al. (2006), which suggested a spatial and kinematic bifurcation of the MS. The observed morphology and kinematics suggest that three of these streams could be interpreted as a 3-way splitting of the main MS filament, while the fourth stream appears much younger and may have originated from the Magellanic Bridge. We find an extensive population of HI clouds at the tip of the MS. Two thirds of clouds have an angular size in the range 3.5'--10'. We interpret this as being due to thermal instability, which would affect a warm tail of gas trailing through the Galactic halo over a characteristic timescale of a few Myrs to a few hundred Myrs. We show that thermal fragments can survive in the hot halo for a long time, especially if surrounded by a <10^6 K halo gas. If the observed clumpy structure is mainly due to thermal instability, then the tip of the MS is at a distance of ~70 kpc. A significant fraction of HI clouds at the tip of the MS show multi-phase velocity profiles, indicating the co-existence of cooler and warmer gas.

Ongoing accretion onto galactic disks has been recently theorized to progress via the unstable cooling of the baryonic halo into condensed clouds. These clouds have been identified as analogous to the High-Velocity Clouds (HVCs) observed in HI in our Galaxy. Here we compare the distribution of HVCs observed around our own Galaxy and extra-planar gas around the Andromeda galaxy to these possible HVC analogs in a simulation of galaxy formation that naturally generates these condensed clouds. We find a very good correspondence between these observations and the simulation, in terms of number, angular size, velocity distribution, overall flux and flux distribution of the clouds. We show that condensed cloud accretion only accounts for ~ 0.2 M_solar / year of the current overall Galactic accretion in the simulations. We also find that the simulated halo clouds accelerate and become more massive as they fall toward the disk. The parameter space of the simulated clouds is consistent with all of the observed HVC complexes that have distance constraints, except the Magellanic Stream which is known to have a different origin. We also find that nearly half of these simulated halo clouds would be indistinguishable from lower-velocity gas and that this effect is strongest further from the disk of the galaxy, thus indicating a possible missing population of HVCs. These results indicate that the majority of HVCs are consistent with being infalling, condensed clouds that are a remnant of Galaxy formation.

We determine the parallax and proper motion of the flaring, non-thermal radio star GMR A, a member of the Orion Nebula Cluster, using Very Long Baseline Array observations. Based on the parallax, we measure a distance of 389 +24/-21 parsecs to the source. Our measurement places the Orion Nebula Cluster considerably closer than the canonical distance of 480 +/- 80 parsecs determined by Genzel et al. (1981). A change of this magnitude in distance lowers the luminosities of the stars in the cluster by a factor of ~ 1.5. We briefly discuss two effects of this change--an increase in the age spread of the pre-main sequence stars and better agreement between the zero-age main-sequence and the temperatures and luminosities of massive stars.

We present Arecibo L-band Feed Array 21-cm observations of a sub-complex of HVCs at the tip of the Anti-Center Complex. These observations show morphological details that point to interaction with the ambient halo medium and differential drag within the cloud sub-complex. We develop a new technique for measuring cloud distances, which relies upon these observed morphological and kinematic characteristics, and show that it is consistent with H-alpha distances. These results are consistent with distances to HVCs and halo densities derived from models in which HVCs are formed from cooling halo gas.

The high-latitude Galactic H I cloud toward the extragalactic radio source 3C 225 is characterized by very narrow 21 cm emission and absorption indicative of a very low H I spin temperature of about 20 K. Through high-resolution optical spectroscopy, we report the detection of strong, very narrow Na I absorption corresponding to this cloud toward a number of nearby stars. Assuming that the turbulent H I and Na I motions are similar, we derive a cloud temperature of 20 (+6, -8) K (in complete agreement with the 21 cm results) and a line-of-sight turbulent velocity of 0.37+/-0.08 km/s from a comparison of the H I and Na I absorption linewidths. We also place a firm upper limit of 45 pc on the distance of the cloud, which situates it well inside the Local Bubble in this direction and makes it the nearest-known cold diffuse cloud discovered to date.

The consortium for Galactic studies with the Arecibo L-band Feed Array (ALFA) is conducting a neutral hydrogen (HI) survey of the whole Arecibo sky (declination range from -1 to 38 degrees), with high angular (3.5') and velocity resolution (0.2 km/s). The precursor observations with ALFA of a region in the Galactic anti-center reveal numerous isolated, small (a few pc in size), and cold (T_k<400 K) HI clouds at low negative velocities, distinctly separated from the HI disk emission (`low-velocity clouds', LVCs). These clouds are most likely located in the transition region between the Galactic disk and halo (at scale heights of 60--900 pc), yet they have properties of typical cold neutral clouds. LVCs are colder and, most likely, smaller and less massive than Lockman's clouds in the disk/halo interface region of the inner Galaxy. Our observations demonstrate that the cloudy structure of the interface region is most likely a general phenomenon, not restricted to the inner Galaxy. LVCs have sizes and radial velocities in agreement with the expectations for clouds formed in low-temperature fountain flows, although we measure a factor of ten higher HI column densities. Alternatively, LVCs could represent the final stages of the infalling intergalactic material in the on-going construction of the Galaxy. In the same dataset at higher negative velocities, we have discovered a `companion' HI cloud located 50' southwest of CHVC186+19-114. CHVC186+19-114 is a typical compact high velocity cloud (HVC) with a well-defined core/envelope structure. The companion cloud has a diameter of only 7'x9', and is one of the smallest HVCs known, most likely stripped from the main cloud through the interactions with the halo medium.