A Review of Recent Observations of Galactic Winds Driven by Star Formation

TL;DR

This review summarizes recent observational advances in understanding galactic winds driven by star formation, highlighting their complex structures, multiphase nature, and persistent scaling relationships across cosmic time.

Contribution

It compiles and analyzes recent observational techniques and findings, emphasizing the multiphase structure and widespread occurrence of galactic winds since 2010.

Findings

Galactic winds are common in star-forming galaxies across cosmic history.

Multiphase structures of winds are revealed by advanced imaging spectroscopy.

Scaling relationships between winds and galaxy properties are consistent over time.

Abstract

Galaxy-scale outflows of gas, or galactic winds (GWs), driven by energy from star formation are a pivotal mechanism for regulation of star formation in the current model of galaxy evolution. Observations of this phenomenon have proliferated through the wide application of old techniques on large samples of galaxies, the development of new methods, and advances in telescopes and instrumentation. I review the diverse portfolio of direct observations of stellar GWs since 2010. Maturing measurements of the ionized and neutral gas properties of nearby winds have been joined by exciting new probes of molecular gas and dust. Low-$z$ techniques have been newly applied in large numbers at high $z$. The explosion of optical and near-infrared 3D imaging spectroscopy has revealed the complex, multiphase structure of nearby GWs. These observations point to stellar GWs being a common feature of rapidly star-forming galaxies throughout at least the second half of cosmic history, and suggest that scaling relationships between outflow and galaxy properties persist over this period. The simple model of a modest-velocity, biconical flow of multiphase gas and dust perpendicular to galaxy disks continues to be a robust descriptor of these flows.

Tables

Table 1. Table 1: Log-Log Fits of Outflow vs. Galaxy Properties.

Axes	Tracer	IP (eV)	$N$	Range	Slope	$p$	Reference
$v_{50\%}$ vs. SFR	Si 2	16.3	48	${10}^{-1}-{10}^2$ M⊙ yr-1	0.22 $\pm$ 0.04	0.001	(Chisholm et al., 2015)
$v_{50\%}$ vs. SFR	Na 1	5.1	41	${10}^{-1}-{10}^3$ M⊙ yr-1	0.35 $\pm$ 0.06	$\cdots$	(Martin, 2005)
$v_{50\%}$ vs. SFR	Na 1	5.1	13	${10}^{0.8}-{10}^{2.2}$ M⊙ yr-1	0.15 $\pm$ 0.06 c	$\cdots$	(Cazzoli et al., 2016)
$v_{50\%}$ vs. SFR	Na 1	5.1	13	${10}^{0.8}-{10}^{2.2}$ M⊙ yr-1	0.30 $\pm$ 0.05 c	$\cdots$	(Cazzoli et al., 2016)
$v_{90\%}$ vs. SFR	Si 2	16.3	48	${10}^{-1}-{10}^2$ M⊙ yr-1	0.08 $\pm$ 0.02	0.002	(Chisholm et al., 2015)
$v_{90\%}$ vs. SFR	Na 1	5.1	59	${10}^{-1}-{10}^3$ M⊙ yr-1	0.21 $\pm$ 0.04	0.001	(Rupke et al., 2005)
$v_{90\%}$ vs. SFR	H 1, N 2	13.6–29.6	48	${10}^{0.7}-{10}^{2.6}$ M⊙ yr-1	0.24 $\pm$ 0.05	0.001	(Arribas et al., 2014)
$v_{\max }$ vs. SFR	Si 2, C 2, Mg 2	15.0–24.4	48	${10}^{-2}-{10}^3$ M⊙ yr-1	0.32 $\pm$ 0.02	0.0001	(Heckman and Borthakur, 2016)
$v_{50\%}$ vs. $M_{\ast }$	Si 2	16.3	48	${10}^9-{10}^{11.5}$ $M_{\odot }$	0.20 $\pm$ 0.05	0.002	(Chisholm et al., 2015)
$v_{90\%}$ vs. $M_{\ast }$	Si 2	16.3	48	${10}^9-{10}^{11.5}$ $M_{\odot }$	0.12 $\pm$ 0.03	0.003	(Chisholm et al., 2015)
$v_{90\%}$ vs. $M_{\ast }$b	Na 1	5.1	52	$\cdots$	0.28 $\pm$ 0.08	0.001	(Rupke et al., 2005)
$v_{50\%}$ vs. $v_{\mathrm{circ}}$a	Si 2	16.3	48	${10}^{1.8}-{10}^{2.5}$ km s-1	0.87 $\pm$ 0.17	0.002	(Chisholm et al., 2015)
$v_{90\%}$ vs. $v_{\mathrm{circ}}$	Na 1	5.1	20	${10}^{1.4}-{10}^{2.7}$ km s-1	0.85 $\pm$ 0.15	0.001	(Rupke et al., 2005)
$v_{90\%}$ vs. $v_{\mathrm{circ}}$a	Si 2	16.3	48	${10}^{1.8}-{10}^{2.5}$ km s-1	0.44 $\pm$ 0.09	0.003	(Chisholm et al., 2015)
$v_{\max }$ vs. $v_{\mathrm{circ}}$a	Si 2, C 2, Mg 2	15.0–24.4	48	${10}^{1.3}-{10}^{2.5}$ km s-1	1.16 $\pm$ 0.37	0.0001	(Heckman and Borthakur, 2016)
$\eta$ vs. $M_{\ast }$	O 1, Si 2–Si 4	13.6–45.1	7	${10}^7-{10}^{11}$ $M_{\odot }$	$-$0.43 $\pm$ 0.07	0.001	(Chisholm et al., 2017)
$\eta$ vs. $M_{\ast }$	H 1, N 2	13.6–29.6	33	${10}^{9.6}-{10}^{11.2}$ $M_{\odot }$	$-$0.43	$\cdots$	(Arribas et al., 2014)
$\eta$ vs. $M_{\ast }$b	Na 1	5.1	42	${10}^{10}-{10}^{11}$ $M_{\odot }$	$-$0.95 $\pm$ 0.20	0.006	(Rupke et al., 2005)
$\eta$ vs. $v_{\mathrm{circ}}$a	O 1, Si 2–Si 4	13.6–45.1	7	$\cdots$	$-$1.56 $\pm$ 0.25	0.001	(Chisholm et al., 2017)

Full text

\pubvolume

6 \issuenum4 \articlenumber138

\historyReceived: 14 September 2018; Accepted: 5 December 2018; Published: 9 December 2018 \updatesyes

\TitleA Review of Recent Observations of Galactic Winds Driven by Star Formation \AuthorDavid S. N. Rupke\orcidA \AuthorNamesDavid S. N. Rupke

1 Introduction

Significant amounts of gas in galaxies move in outward radial trajectories due to energy imparted by star formation. This energy originates from a combination of phenomena rooted in stellar processes: radiation, winds, explosive events, and cosmic rays. These gas outflows powered by star formation, or stellar galactic winds (GWs), continue to be a dynamic topic of observational research in the era of large galaxy surveys and multi-messenger astronomy. Outflows are challenging to characterize largely because of two factors: the large contrast between the outflow and an underlying galaxy disk; and the complex, multiphase structure of outflows. High quality data in many gas tracers are essential to adequately quantify the mass, momentum, and energy budget of the wind.

Advances in observations of GWs driven by stellar processes in this decade have come from the use of new observational techniques and the application of old techniques to much larger samples. Both have been aided by new-and-improved telescopes or instrumentation. In the first category fall molecular gas transitions newly applied to study outflows (notably the hydroxyl molecule), mid- and far-infrared (MIR/FIR) imaging of dust, and resonant-line emission in the rest-frame ultraviolet (UV). In the second category are surveys with integral field spectrographs (IFSs), with multi-object long-slit spectrographs, and with the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST).

Excellent and thorough reviews from previous decades of theory and observations of GWs provide in-depth discussions of the quantities of observational interest and the astrophysics of GWs (e.g., (Heckman et al., 1993; Veilleux et al., 2005)). They also show the trajectory and progress of the field. The scope of this review is narrower. I synthesize observational results from the current decade (approximately 2010 through the present). I focus on direct measures of outflows and do not discuss studies that infer the presence or properties of GWs by studying other phenomena.

As an example of an indirect measurement of the presence or properties of GWs, studies of the mass-metallicity relation constrain how efficiently GWs eject metals from galaxy disks compared to metal production and reaccretion (e.g., (Zahid et al., 2014)). A second example is the significant reservoirs of highly ionized carbon and oxygen that preferentially arise in the circumgalactic media of actively star-forming galaxies (Tumlinson et al., 2011; Borthakur et al., 2013; Bordoloi et al., 2014; Borthakur et al., 2016; Heckman et al., 2017). A logical source of these metals is stellar GWs. Third, the hot halos that appear to be a common feature of star-forming galaxies (Li and Wang, 2013) may be produced by stellar GWs. Finally, cosmological simulations typically use numerical prescriptions for the unresolved physics of stellar GWs and compare simulated galaxy properties with observed galaxy properties (like the galaxy mass function) to constrain the nature of GWs (e.g., (Oppenheimer et al., 2010; Hopkins et al., 2014; Anglés-Alcázar et al., 2014)). Such indirect measurements are essential for a complete picture of the relationship of GWs to their surrounding environments. However, at present, interpretations of these measurements based on stellar GWs typically compete with other physical models and are rarely definitive.

Stellar processes are not the only possible drivers of GWs. Evidence continues to accumulate that actively accreting supermassive black holes (active galactic nuclei, or AGN) are also important in powering GWs in galaxies with above-average masses (see reviews in (Fabian, 2012; Harrison, 2017)). However, it can be difficult to distinguish whether the AGN is energetically important for the GW if significant star formation is also present, except in the most powerful AGN. This difficulty is due to the possible power sources being cospatial at low resolution and to the uncertain duty cycle of AGN. Consequently, in this review I focus on studies of purely star-forming systems, or at least those where star formation clearly dominates the galaxy’s luminosity. AGN with low-to-moderate Eddington ratios are almost certainly present in many galaxies classified as purely star-forming (Aird et al., 2012; Jones et al., 2016), but their energetic importance for GWs remains unquantified.

I organized this review at the highest level by redshift. This is useful for two reasons. First, star formation and galaxy properties at redshifts above unity differ significantly from those in the local universe. The global star formation rate (SFR) peaked at $z\sim 1.9$ (e.g., (Madau and Dickinson, 2014)) and gas mass fraction increases with increasing redshift (e.g., (Tacconi et al., 2013)). Galaxies may grow from the inside out, meaning that disks are more extended compared to stars at high $z$ (van Dokkum et al., 2015; Nelson et al., 2016). Second, low-$z$ winds are much better characterized because many more techniques can be brought to bear. Some techniques are in practice easier to apply to high-redshift galaxies, but on balance this is not the case.

2 Winds Driven by Star Formation at Low redshift

Most data on nearby stellar GWs in the first 20 years of earnest work came from the optical and X-ray bands. These studies fell largely into two classes: (1) long-slit spectroscopic and/or narrowband imaging surveys (e.g., (Lehnert and Heckman, 1996; Martin, 1998; Rupke et al., 2005)) or (2) case studies using optical 3D spectroscopy, X-ray imaging spectroscopy, or narrowband imaging (e.g., (Veilleux et al., 1994; Lehnert et al., 1999; Veilleux and Rupke, 2002; Strickland and Heckman, 2009)). This work was primarily focused on starburst galaxies, which lie above the star-forming main sequence in SFR. High-mass starbursts are typically luminous in the infrared (IR) and classified as either luminous or ultraluminous infrared galaxies. These LIRGs and ULIRGs are defined to have $L_{\mathrm{IR}}>10^{11}$ $L_{\odot}$ and $L_{\mathrm{IR}}>10^{12}$ $L_{\odot}$, which corresponds to $\mathrm{SFR}>10$ M⊙ yr*-1* and $\mathrm{SFR}>100$ M⊙ yr*-1* if all the IR luminosity is powered by star formation. Low-mass (dwarf) starbursts are luminous in the UV and optical bands.

New observing capabilities, surveys, and archival databases have allowed probes of more physical conditions, broadening the picture of winds sketched first by small studies of the warm and hot ionized and cool neutral phases. They have also extended the study of winds to main-sequence galaxies.

This section is divided into four parts. First, I discuss new UV surveys of nearby starbursts. Second, I visit the use of the Sloan Digital Sky Survey (SDSS) for outflow studies. Third, I summarize results from the widening use of integral field spectroscopy. I end by surveying significant advances in revealing and quantifying the molecular gas and dust in GWs using ground-based submillimeter and space-based FIR telescopes.

2.1 Ultraviolet Surveys

The rest-frame UV contains many interstellar absorption lines and Ly$\alpha$, the brightest UV emission line in unobscured starbursts. Previous UV instruments have been used to study individual systems or small samples, but the sensitivity of COS has enabled population studies of low-$z$, UV-bright starbursts at high signal-to-noise (S/N). These studies are large enough to study correlations between the properties of GWs and galaxy properties.

Two studies find a correlation between the velocity of the ionized outflow and basic properties of the galaxy and starburst itself: stellar mass $M_{*}$ and star formation rate SFR (Chisholm et al., 2015; Heckman and Borthakur, 2016). Table 1 compares these correlations and the ionization states probed. The first (Chisholm et al., 2015) uses COS data on 48 starbursts at $z<0.26$ taken from the literature. The second (Heckman and Borthakur, 2016) combines COS and Far-Ultraviolet Spectroscopic Explorer (FUSE) data on 39 starbursts at $z<0.25$ (Heckman et al., 2015) plus MMT or Keck spectra of 9 compact starbursts at $z=0.4-0.7$ (Sell et al., 2014; Geach et al., 2014). Notably, the studies differ in their definition of outflow velocity. The first (Chisholm et al., 2015) uses the 50% and 90% points in the cumulative velocity distribution (CVD; starting from the red side of the line). The second (Heckman and Borthakur, 2016) uses either the 98% point of the CVD (in the case of data taken from secondary sources (Sell et al., 2014; Geach et al., 2014)) or the authors’ own velocity measurements ($v_{\mathrm{max}}$). I label the correlations from this second study with $v_{\mathrm{max}}$ for simplicity.

While both studies find correlations, the ranges of galaxy properties in one are wider (Heckman and Borthakur, 2016). This study also finds steeper slopes in lines fitted to the data; the slopes from the two studies differ at $>1\sigma$. The second study does not specify a fitting method (Heckman and Borthakur, 2016) , while the first uses a method that accounts for $x$ and $y$ errors, outliers, and scatter (Chisholm et al., 2015). Neither study finds a significant correlation between outflow velocity and specific star formation rate sSFR $\equiv\mathrm{SFR}/M_{*}$. One also finds no correlation with star formation rate surface density $\Sigma_{\mathrm{SFR}}$, but their sample has a range $\Sigma_{\mathrm{SFR}}=10^{-2}-10^{1}$ M⊙ yr*-1* kpc*-2* (Chisholm et al., 2015). The second study adds three orders of magnitude to the upper end of this range and do find a correlation, which they parameterize as $v_{\mathrm{max}}=3296/[(\Sigma_{\mathrm{SFR}}/1307.9)^{-0.34}+(\Sigma_{\mathrm{SFR}}/1307.9)^{0.15}]$ at a significance of $p<0.0001$ (Heckman and Borthakur, 2016).

The mass outflow rate $\dot{M}$ and mass outflow rate normalized to the star formation rate ($\eta\equiv\dot{M}/\mathrm{SFR}$)111This quantity has an uncertain physical interpretation, though at face value it might quantify the capability of a wind to act as negative feedback on star formation. The logic is: if $\eta>1$, then more gas is leaving the region than is forming stars. Thus, the outflow is going to deplete gas more quickly than star formation (maybe leading to fewer stars in the end). It is sometimes referred to as the mass-loading efficiency or factor, perhaps suggesting that it measures the amount of gas that is “loaded” into the wind as it moves through the galaxy. That is, it is the ratio of the mass of gas that eventually emerges from star-forming regions to the mass of outflowing gas that is initially produced by star formation through, e.g., stellar winds and supernovae. However, it is not a direct measure of this. The term mass-loading originated with studies of how much cool, ambient gas is mixed into a hot wind phase in the model of a wind driven by a hot fluid. were computed for subsamples of these two larger samples (Heckman et al., 2015; Chisholm et al., 2017). One study estimates the wind column density using multiple ions and a stacked spectrum; they apply this single column density to the entire sample of 39 galaxies (Heckman et al., 2015). Their mass outflow rate is then a product of the outflow velocity and starburst radius times a constant. They find a correlation between $\dot{M}$ and SFR that is near-linear (though the slope is unquantified) and which appears to be driven largely by the correlation between $\dot{M}$ and $v_{\mathrm{max}}$. They also find that $\eta$ is inversely correlated with SFR and $M_{*}$ (or $v_{\mathrm{circ}}$).

The other study finds a similar result for the inverse correlation of $\eta$ and galaxy mass ((Chisholm et al., 2017); Figure 1; Table 1). While based on a much smaller sample size, this work uses high signal-to-noise spectra, velocity-resolved optical depths and covering factors, and photoionization modeling to estimate the ionization state, density, metallicity, and inner radius of the wind (assuming a model relating velocity and radius). The measured slope ($-$0.43 $\pm$ 0.07) matches predictions from numerical simulations, which are in the range $-$0.35 to $-$0.50 at low masses (Hopkins et al., 2012; Lagos et al., 2013; Muratov et al., 2015; Mitra et al., 2015).

The velocity measurements in these UV surveys employ lines from ions of relatively low ionization potential (Chisholm et al., 2015; Heckman and Borthakur, 2016), while the estimated wind masses are inferred from a larger range of ions (Heckman et al., 2015; Chisholm et al., 2017). The resulting fits to GW vs. galaxy properties are in overall agreement with similar fits made a decade earlier using the cool, neutral gas phase (Table 1; (Rupke et al., 2005; Martin, 2005)). The spread in fitted slopes illustrates the systematic uncertainties in this enterprise, but a conservative synthesis points to velocity depending on SFR and $M_{*}$ as $v\sim\mathrm{SFR}^{0.1-0.3}$ and $v\sim M_{*}^{0.1-0.3}$. The dependence on circular velocity is thus steeper (close to linear). The upper end of the range of slopes of $v$ vs. SFR matches a prediction from simulated absorption lines in 3D hydrodynamic simulations (Tanner et al., 2017). There are fewer measurements of how $\eta$ depends on various properties, and the systematic uncertainties in $\eta$ are higher than in velocity. However, as discussed above, published fits are surprisingly consistent with predictions from numerical simulations (Mitra et al., 2015; Muratov et al., 2015), including at high masses where the slope steepens (Rupke et al., 2005; Mitra et al., 2015).

This begs the question of what drives the scatter in these relationships. It is worth noting that the pursuit of scaling relations relies on reducing the velocity field of a galactic outflow to a single parameter in a single tracer. Outflows are unlike galaxies in that they are inherently a violent, non-equilibrium process whose structural properties are governed by time-varying power sources, gas hydrodynamics, and radiation transfer rather than gravitational processes. Their properties thus almost certainly depend on multiple host parameters simultaneously. The fact that outflow properties scale at all with bulk galaxy properties such as SFR is not surprising (more energy means higher-velocity gas). There are clearly important details beyond this, however, that are harder to quantify and may reflect the unique structural properties or history of particular galaxies (Is it a merger? Did the outflow just turn on?) or other hidden scalings that are orthogonal to those with SFR and mass. One solution might be to combine multiple galaxy or wind properties into single parameters in search of an “outflow fundamental plane,” but combining, e.g., two parameters into one has yet to noticeably improve matters (Heckman and Borthakur, 2016; Chisholm et al., 2017).

The imprint of GWs in the UV is found not only absorption lines but also in Ly emission lines and the escape of ionizing radiation beyond the Lyman limit. Extended Ly$\alpha$ emission and Lyman continuum (LyC) escape may be caused in part by outflows that create low-density holes in a galaxy’s ISM.

Recent studies of Ly$\alpha$ and absorption lines in samples of order tens of galaxies have quantified the impact of outflows on the properties of Ly$\alpha$ and the escape of LyC. Population studies of Ly$\alpha$ in nearby galaxies are, however, more complicated to interpret than studies of absorption lines because of the resonant emission and absorption behavior of Ly$\alpha$ and its sensitivity to dust. As a result, these studies have not reached firm conclusions on how these tracers reflect GW properties. There is some suggestion that outflow velocity is correlated with Ly$\alpha$ escape fraction (Wofford et al., 2013; Martin et al., 2015; Rivera-Thorsen et al., 2015; Henry et al., 2015; Alexandroff et al., 2015; Chisholm et al., 2017; Jaskot et al., 2017). However, this correlation has significant scatter, and cases exist of low outflow velocity and high LyC escape (Jaskot et al., 2017). Ly$\alpha$ and outflow velocity are also both correlated with SFR and/or $\Sigma_{\mathrm{SFR}}$ (Chisholm et al., 2017; Alexandroff et al., 2015). In one sample, galaxies with and without escaping LyC have similar outflow velocities (Chisholm et al., 2017). The simple presence of outflows may be a necessary but not a sufficient condition for Ly$\alpha$ and/or LyC escape (Rivera-Thorsen et al., 2015); other factors such as the outflow acceleration, low H 1 column density, or low metallicity may also be required for Ly$\alpha$ or LyC to exit the dense regions of a galaxy (Martin et al., 2015; Henry et al., 2015; Chisholm et al., 2017).

2.2 Single-Aperture Surveys

Searches for outflows in single-aperture optical spectroscopic surveys like SDSS are hampered by the spatially unresolved bright host galaxy. Outflow features in emission lines are typically overwhelmed by ionized gas emission from the star-forming host disk, which in galaxies with modest GWs lies in the same velocity range as the outflow. Interstellar absorption lines can be similarly dominated by stellar absorption features.

Both limitations can be overcome with data of high enough S/N and high-quality modeling of the underlying stellar continuum. Stacking of many spectra is commonly used to achieve the required S/N, and stellar models that match the spectral resolution of SDSS are mature. (Even if these models may not uniquely constrain the star formation history or stellar population properties, they fit the data very well.)

Three studies of GWs in star-forming galaxies have stacked SDSS DR7 data. The authors of the first of these stack 105 high-mass ($M_{*}\sim 10^{10}-10^{11}$ $M_{\odot}$) galaxies in bins of various physical parameters to study the properties of cool, neutral outflows using the interstellar Na 1 D doublet (Chen et al., 2010). They detect this resonant line in both absorption and emission. They find that the velocity and equivalent width of absorption scale with inclination, such that face-on galaxies are observed to have faster, higher equivalent width outflows, consistent with the model of minor-axis outflows. Outflow equivalent width and linewidth also correlate with $\Sigma_{\mathrm{SFR}}$ (over the range $10^{-2.5}-10^{-0.5}$ M⊙ yr*-1* kpc*-2*) and $M_{*}$ (over the small range $10^{10.3}-10^{11.2}$ $M_{\odot}$). However, systemic interstellar absorption also scales with $\Sigma_{\mathrm{SFR}}$ and $M_{*}$, so the nature of these correlations is uncertain.

The second study uses 200,000 galaxies to stack by SFR and $M_{*}$ over a range of galaxy properties comparable to the ranges in samples of individual galaxies ((Cicone et al., 2016); Table 1; Sections 2.1 and 2.3): $\mathrm{SFR}=10^{-2.7}-10^{2.3}$ M⊙ yr*-1* and $M_{*}=10^{7.3}-10^{11.8}$ $M_{\odot}$. They fit [O 3] 5007 Å, H$\alpha$, and [N 2] 6548, 6583 Å using a high S/N instrumental profile and extract the line-of-sight velocity distribution (LOSVD) of both the ionized gas and stars. They use excess blueshifted emission-line gas at the extreme end of the gas LOSVD to measure outflow velocity, with the stellar LOSVD serving as a reference. For star-forming galaxies selected by line ratio (Kauffmann et al., 2003), these authors find that $v_{\mathrm{out}}$ correlates significantly with SFR and sSFR. $v_{\mathrm{out}}$ does not correlate with $M_{*}$, though it does with stellar velocity dispersion, suggesting a discrepancy in how these quantities are measured. This method is most sensitive to higher-velocity outflows, and outflows are detected primarily at $\mathrm{SFR}>1$ M⊙ yr*-1* and $\mathrm{sSFR}>10^{-9}$ yr*-1*. They are also preferentially observed in galaxies with SFR values that put them above the main sequence.

Finally, a third study employs a similar technique over a larger sample of 600,000 galaxies in bins of SFR and $M_{*}$ (Concas et al., 2017). These authors apply two methods for parameterizing outflows: (1) they fit the [O 3] emission line in each stack using two Gaussians; and (2) they use the observed [O 3] profile as the gas LOSVD and calculate measures of line width and asymmetry. They then remove the instrumental resolution in quadrature. They find no significant detection of starburst-driven winds in star-forming galaxies. However, they select star-forming galaxies using a more restrictive line-ratio criterion (Stasińska et al., 2006) than the ionized gas study that does detect stellar GWs (Cicone et al., 2016). They argue that a more permissive selection (Kauffmann et al., 2003) allows low-luminosity AGN to contaminate a sample of star-forming galaxies; these AGN could in turn produce outflows in galaxies that are supposedly purely star-forming. Low-luminosity AGN exist throughout the star-forming sequence with a range of contributions to [O 3] (Aird et al., 2012; Jones et al., 2016). It is certainly plausible that star formation, rather than low-luminosity AGN, powers the outflows detected in galaxies whose line emission is dominated by star formation (Cicone et al., 2016), but separating potential contributions from low-luminosity AGN requires further work.

These stacking studies add weight to the conclusions drawn from studies of smaller samples and individual galaxies by probing more statistically complete samples of galaxies across a wider range of basic galaxy properties. There is agreement between two of the stacking studies discussed above that ionized and/or neutral winds are a common feature of galaxies of modest to high SFR across the mass spectrum (Chen et al., 2010; Cicone et al., 2016), consistent with results discussed elsewhere in this review. The disagreement with the third study (Concas et al., 2017) on this point is puzzling, though the difficulty of separating AGN from star formation as a possible energy source in galaxies with supermassive black holes accreting at lower Eddington ratios is a valid concern. Furthermore, the correlations between wind and galaxy properties are akin to those seen using UV and other optical studies (Section 2.1; Table 1), and the connection to inclination angle is consistent with a wide range of observational results.

2.3 IFS Results

Though they cannot yet match the numbers of galaxies observed in the SDSS, spatially resolved spectroscopic surveys with optical IFSs promise to significantly improve the detectability and characterization of GWs. This is true both at ground-based resolution and with enhanced spatial resolution from adaptive optics; the latter is important for probing the nuclear regions of galaxies where GWs emerge from their power sources. Studies of one or a few galaxies with stellar GWs published in the current decade are numerous (Rich et al., 2010; Fogarty et al., 2012; Rupke and Veilleux, 2013; Vogt et al., 2013; Ho et al., 2014; Medling et al., 2015; Martín-Fernández et al., 2016; López-Cobá et al., 2017; Cresci et al., 2017; Cortijo-Ferrero et al., 2017). These detailed case studies highlight the ubiquity and complex, multiphase nature (ionized, neutral, dusty) of these winds (Figure 2) and the power of IFS for leveraging the combination of spectral and spatial information to separate outflows from their hosts.

Small surveys using instruments that target one galaxy at a time have characterized GWs in galaxies with the highest star formation rates. The authors of a study of the ionized gas in 27 LIRGs (Rich et al., 2014, 2015) found frequent shock ionization accompanied by high gas velocity dispersions, a conclusion supported by an earlier long-slit study of ULIRGs (Soto et al., 2012). The line ratios are consistent with models of slow shocks, correlate with $\sigma_{\mathrm{gas}}$, and may be caused in part by GWs (Rich et al., 2014, 2015).

The authors of a larger study of the ionized and cool, neutral gas in $\sim$50 LIRGs and ULIRGs (Arribas et al., 2014; Cazzoli et al., 2016) concur that high ionized gas linewidths are shock-powered. To detect GWs, they spatially integrated their spectra (Arribas et al., 2014; Cazzoli et al., 2016) and/or fit the spatially resolved data (Cazzoli et al., 2016). In the former case, they corrected the velocity of each spatial position for gravitational motions and then summed over the field of view. Fits to the wind velocity vs. SFR give best-fit slopes of 0.24 $\pm$ 0.05 (ionized gas), 0.15 $\pm$ 0.06 (integrated neutral gas), and 0.30 $\pm$ 0.05 (spatially resolved neutral gas). These are consistent with the single-aperture fits discussed above (Section 2.1; Table 1). $\eta$ is between 0.1 and 1 in the neutral and ionized gas, with average values similar to those previously estimated in LIRGs and ULIRGs (Rupke et al., 2005; Rupke and Veilleux, 2013). The study of integrated spectra (Arribas et al., 2014) finds evidence that $\eta$ decreases with increasing dynamical mass over the range $10^{9.6}-10^{11.2}$ $M_{\odot}$, with a log-log slope of $-$0.43 that is identical to that found in the UV over a much wider mass range ((Chisholm et al., 2017); Section 2.1; Table 1). (A steeper slope in this relationship at high masses, as previously measured (Rupke et al., 2005), is predicted by some simulations (Lagos et al., 2013; Mitra et al., 2015).) Finally, these authors fit a relationship between $\eta$ and $\Sigma_{\mathrm{SFR}}$ with a log-log slope of 0.17 (Arribas et al., 2014).

Two major (of order thousands of galaxies) IFS surveys have been ongoing for several years: the SAMI (Sydney-Australian-Astronomical-Observatory Multi-object Integral-Field Spectrograph) Galaxy Survey (Bryant et al., 2015) and MaNGa (Mapping Nearby Galaxies at Apache Point Observatory) (Bundy et al., 2015). These surveys employ instruments that target multiple galaxies at once, each with a single, small IFS. Though limited in spatial resolution (of order 2*′′*), the ability of these surveys to integrate to high S/N, collect large samples of galaxies across a wide parameter range, and separate the outflow from the host galaxy using both spatial and spectral information with IFS promises more accurate statements about the ubiquity and property of stellar GWs over a wider range of galaxy types.

The first result from these surveys relies on a sample of 40 edge-on SAMI disk galaxies (Ho et al., 2016). These authors use minor-axis kinematic asymmetries to detect winds and find them down to very low $\Sigma_{\mathrm{SFR}}$ (10*-3* M⊙ yr*-1* kpc*-2*). The asymmetry increases with $\Sigma_{\mathrm{SFR}}$ and is associated with recent bursts of star formation. The connection of the asymmetry to GWs and its association with high $\Sigma_{\mathrm{SFR}}$ and high-temperature extraplanar gas is bolstered by comparison to simulations (Tescari et al., 2018).

2.4 Molecular Gas and Dust

The most rapid progress in the study of stellar GWs has come from molecular gas and dust continuum measurements. The direct dust measurements have been made in the MIR/FIR using space-based telescopes: AKARI, the Spitzer Space Telescope, and the Herschel Space Observatory. The molecular gas measurements have been led by ground-based interferometers such as the Atacama Large Millimeter/submillimeter Array (ALMA), SubMillimeter Array (SMA), Institut de Radioastronomie Millimétrique (IRAM) Plateau de Bure interferometer (PdBI), and Nobeyama Millimeter Array (NMA), though Herschel has also played an important role. While sample sizes are still small and most studies have focused on individual galaxies in the starburst regime, results point to the ubiquity of dusty molecular gas entrained in stellar GWs.

Recent interferometric measurements of CO transitions have uncovered outflowing cold molecular gas in approximately 11 nearby LIRGs and ULIRGs whose energy output is dominated by star formation (Tsai et al., 2009; Sakamoto et al., 2014; Cicone et al., 2014; García-Burillo et al., 2015; Pereira-Santaella et al., 2016; Sakamoto et al., 2017; Barcos-Muñoz et al., 2018; Falstad et al., 2018). Molecular outflows in four very nearby star-forming galaxies with lower SFR have also been studied in detail: NGC 3628, M82, NGC 253 (Figure 3), and NGC 1808 (Tsai et al., 2012; Bolatto et al., 2013; Leroy et al., 2015; Salak et al., 2016). M82 had been the only previous starburst galaxy known to host a molecular outflow (García-Burillo et al., 2001; Walter et al., 2002).

These molecular outflows are typically confined to the inner kpc of the galaxy (in radius). They have modest outflow velocities in nearby disks (tens to a few hundred km s*-1*) but have a broader range of peak velocities in the more luminous starbursting disks and mergers (300–800 km s*-1*). Similarly, the outflow rates appear to scale with SFR and/or $\Sigma_{\mathrm{SFR}}$. The nearby, low-luminosity starbursts have $\dot{M}$ values of a few to a few tens of M⊙ yr*-1*, while the LIRGs and ULIRGs have a broader range of estimated $\dot{M}$: several to several hundred M⊙ yr*-1* (perhaps even $\sim$1000 M⊙ yr*-1* in F17208$-$0014; (García-Burillo et al., 2015)). Mass outflow rates of these magnitudes are similar to the star formation rates in these systems, and equal to or larger than the mass outflow rates estimated in other gas phases. The larger mass outflow rates in the molecular phase are due primarily to larger total outflowing gas masses.

Dense gas, as probed by characteristic molecules such as HCN and HCO+, is also a common feature of these outflows. It is present in the extended outflows of M82, NGC 253, NGC 1808, and Arp 220, as observed in emission using interferometry (Salas et al., 2014; Kepley et al., 2014; Walter et al., 2017; Salak et al., 2018; Barcos-Muñoz et al., 2018). Outflowing dense molecular gas has also been observed in absorption in the inner regions of the highly obscured binary Arp 220 (Sakamoto et al., 2009; Tunnard et al., 2015; Zschaechner et al., 2016; Martín et al., 2016).

The ubiquity of starburst-driven molecular outflows was first discovered using FIR OH and H2O absorption lines with Herschel (Rangwala et al., 2011; Sturm et al., 2011; Spoon et al., 2013; Veilleux et al., 2013; González-Alfonso et al., 2017). The outflow properties inferred from these spatially unresolved absorption lines and detailed radiation transfer models are comparable to measurements from molecular emission-line interferometry (González-Alfonso et al., 2017): massive, high-velocity but compact flows.

The cold molecular and cool atomic gas share similar dynamics. Global OH properties correlate with those measured from the 158 $\upmu$m [C 2] emission line (Janssen et al., 2016), which traces photo-dissociation regions (PDRs). The 158 $\upmu$m [C 2] line may be a key technique for tracing outflows at high $z$ as it moves into the submillimeter band. OH properties may also correlate with Na 1 D atomic gas measurements (Rupke et al., 2005, 2005; Martin, 2005; Rupke and Veilleux, 2013; Veilleux et al., 2013; Janssen et al., 2016; González-Alfonso et al., 2017), though there is significant scatter in this correlation. On the spatially resolved level, FIR atomic fine structure lines are seen in the GWs of M82 (Contursi et al., 2013), NGC 2146 (Kreckel et al., 2014), and NGC 253 (Krips et al., 2016). In M82 and NGC 253 these atomic PDR tracers correlate well with the molecular gas in terms of morphology and dynamics (Contursi et al., 2013; Krips et al., 2016).

Direct observations of atomic H 1 in outflows may become possible with the next generation of sensitive radio arrays. New Jansky Very Large Array H 1 data on M82 are consistent with a deceleration of the outflow as it moves upward into the halo (Martini et al., 2018), suggestive of a galactic fountain (Leroy et al., 2015).

Warm molecular gas does not contribute significantly to the mass budget of outflows, but the strength of its emission lines in the accessible NIR band makes it a useful tracer of the extent and physical state of molecular gas. Deep near-infrared (NIR) imaging and spectroscopy shows that shocked, warm H2 in the M82 outflow extends several kpc from the disk along the minor axis (Veilleux et al., 2009; Beirão et al., 2015). Outflowing warm molecular gas has been detected in several nearby, starburst-dominated ULIRGs using the 2.12 $\upmu$m H2 1-0 S(1) transition (Emonts et al., 2014; Medling et al., 2015; Emonts et al., 2017). These outflows are compact ($<$2 kpc), with velocities of a few hundred km s*-1*, and the relationship between the cold and warm phases requires further study. Indirect evidence from the global properties of a large sample of ULIRGs implicates shocks in GWs as the origin of excess warm H2 in this population (Hill and Zakamska, 2014).

Dust and molecular gas are abundant under similar physical conditions. Observations of filamentary dust structures along disk minor axes in some wind systems and the correlation of these dust features with outflowing atomic gas columns are indirect evidence for outflowing dust in stellar GWs (Figure 2; e.g., (Phillips, 1993; Chen et al., 2010; Rupke and Veilleux, 2013)). Further evidence in the form of UV reflection nebula and polarized line emission strongly suggests the presence of dusty winds in M82 and NGC 253 (Hoopes et al., 2005; Yoshida et al., 2011; Hutton et al., 2014). Modeling of UV reflection in the M82 wind indicates smaller average grain sizes in the M82 wind compared to the disk (Hutton et al., 2014, 2015).

Observations of thermal emission from dust in galaxies where the disk and outflow can be spatially separated strengthen the case that stellar GWs are dusty in at least some cases. Studies detect 1–4$\times$106 $M_{\odot}$ of dust in M82 (traced at 7–500 $\upmu$m) extending far from the disk plane (Kaneda et al., 2010; Roussel et al., 2010), as well as significant amounts of cold dust emplaced by tidal interaction with M81 (Leeuw and Robson, 2009).

Roughly 106 $M_{\odot}$ of warmer dust (traced at 70–160 $\upmu$m) is also found in the outflows of NGC 253 and NGC 4136 (Kaneda et al., 2009; Meléndez et al., 2015) in structures that correlate with ionized gas emission at other wavelengths. The extended cold dust in NGC 4631, however, may be of tidal origin (Meléndez et al., 2015), as in the case of M82. In the dwarf galaxy NGC 1569, 70–500 $\upmu$m imaging indicates a large reservoir (a few times 105 $M_{\odot}$) of circumgalactic dust, perhaps deposited by its starburst-driven GW (McCormick et al., 2018). At shorter wavelengths, GW features in NGC 1569 are seen in warm dust emission (Suzuki et al., 2018). Contrary to these results, no dust is detected at 37$\upmu$m in the NGC 2146 wind (Kreckel et al., 2014). FIR observations of dwarf galaxies with outflows observed at other wavelengths paint a nuanced picture of dusty outflows in systems with lower SFR (McCormick et al., 2018); few systems show abundant circumgalactic dust or strong correlations with multiwavelength signatures of GWs.

Finally, polycyclic aromatic hydrocarbons (PAHs) are prominent features of GWs in those galaxies in which they have been detected. As alternately large molecules or small dust grains, PAHs bridge the gap between molecular gas and dust and are luminous features in rest-frame 3–20 $\upmu$m spectra. An enormous PAH nebula extends from the M82 disk (Engelbracht et al., 2006) and shows evidence for grain shattering from PAH line ratios (Yamagishi et al., 2012; Beirão et al., 2015). PAHs have been found in GWs in other galaxies (Onaka et al., 2010; McCormick et al., 2013), and the amount of extraplanar PAH emission in galaxies may correlate with $\Sigma_{\mathrm{SFR}}$ (McCormick et al., 2013).

The warm and hot ionized gas traced by UV absorption lines, optical emission lines, and X-rays are prevalent in GWs (Sections 2.1– 2.3). Quantifying the physical relationship between the molecular and ionized gas phases thus appears to be within reach. In a few nearby systems this connection has been made, but there has been no attempt for most galaxies (probably because the ability to detect the molecular phase of the outflow is still relatively new). There is a need for sensitive, multiphase studies at the highest spatial and spectral resolution to connect these outflow components.

M82 and NGC 253 both have deep, high-spatial-resolution maps of warm ionized gas (H$\alpha$), hot ionized gas (soft and hard X-rays), and molecular gas. These two exemplars paint a picture of an innermost hot, ionized wind fluid (Strickland and Heckman, 2009) that entrains warm and hot ionized material which are spatially correlated with each other and surround the hot fluid, perhaps as a shell (Shopbell and Bland-Hawthorn, 1998; Strickland et al., 2000). The cold molecular gas in turn envelops these ionized phases as it is entrained from the disk (Figure 3; (Bolatto et al., 2013; Leroy et al., 2015)). The relationship between the warm molecular and ionized phases, however, is likely more complex. In examples of minor-axis outflows where both are spatially resolved, there is coarse-grained spatial correlation (Arp 220 (Lockhart et al., 2015; Barcos-Muñoz et al., 2018), F08572$+$3915 (Rupke and Veilleux, 2013, 2013), and M82 (Veilleux et al., 2009)), but in other cases there is not (NGC 4945; (Marconi et al., 2000)). At very high spatial resolution, apparent correlation may break down (Veilleux et al., 2009). Note that the outflows in some of these galaxies may be AGN-driven, but they are included here because the sample size of galaxies with resolved molecular and ionized outflows is unfortunately small.

3 The Nearest Galactic Wind

The unexpected recent discovery of the so-called Fermi bubbles–large, diffuse $\gamma$-ray structures that form a bipolar shape above and below the Galactic plane–has revived interest in characterizing the GW in the Milky Way (Su et al., 2010; Ackermann et al., 2014). Previous data provided strong evidence of its existence (see references in (Veilleux et al., 2005)), but the picture-perfect morphology of the Fermi bubbles make it an almost inarguable fact.

This review is concerned with observations of starburst-driven GWs, and the Milky Way’s GW may be driven by star formation (e.g., (Crocker et al., 2015)). It also may be powered by the Galactic nuclear black hole during a previous accretion episode (see, e.g., the recent review of relevant data and models in (Yang et al., 2018)). However, for completeness we note some recent observations, since the Milky Way is an excellent laboratory for studying a GW at high sensitivity and spatial resolution in what is a “typical” galaxy in the Local Universe. Since its discovery, the Fermi bubbles have notably been found to also contain a magnetized radio plasma (Carretti et al., 2013) and have been connected to the previously known microwave “haze” that was re-observed with Planck (Planck Collaboration et al., 2013). The bubbles also appear to host neutral gas clouds moving up to several hundred km s*-1* (McClure-Griffiths et al., 2013; Lockman and McClure-Griffiths, 2016; Di Teodoro et al., 2018), as well as higher-ionization species observed in absorption (Fang and Jiang, 2014; Bordoloi et al., 2017; Savage et al., 2017; Karim et al., 2018). These neutral and ionized clouds may lie along filaments swept up by the bubble along its edges, though their contribution to the structure and mass/energy budget of the outflow are not yet clear.

A recent claim has also been made for a GW in a satellite of the Milky Way. This potential outflow in the Large Magellanic Cloud was extrapolated from absorption-line measurements of a single line of sight through its disk (Barger et al., 2016). The LMC contains 30 Doradus, which is undergoing a significant star formation episode.

4 Winds Driven by Star Formation at High redshift

The presence and ubiquity of stellar GWs outside of the local universe was first evident in blueshifted rest-frame UV absorption lines and complex Ly$\alpha$ profiles in Lyman-break galaxies (LBGs) (e.g., (Pettini et al., 2000, 2001; Steidel et al., 2010)). Information on stellar GWs at high $z$ comes predominantly from single-aperture, down-the-barrel spectroscopy of rest-frame UV lines. (“Down-the-barrel” refers to sightlines toward the galaxy itself.) However, newer instruments and techniques have in the past decade opened other avenues to study high $z$ winds. These include rest-frame optical measurements of emission lines with wide-field, NIR, multi-object spectrographs; multiplexed or wide-field IFS instruments such as the K-band Multi-Object Spectrograph (KMOS) and the Multi-Unit Spectroscopic Explorer (MUSE) that enable multi-object and/or spatially resolved measurements; NIR, adaptive optics IFS to achieve high spatial resolution; FIR and submillimeter observations that probe molecular and atomic gas transitions; and transverse-sightline spectroscopic surveys. These new techniques have established the ubiquity and properties of stellar GWs in new galaxy populations and constrained the redshift evolution of their bulk properties.

Deep spectroscopy of strong rest-frame optical emission lines (mainly H$\alpha$ and [N 2]) at $z\sim 2$ reveals broad wings that arise primarily from bright star-forming regions (Genzel et al., 2011; Le Tiran et al., 2011; Newman et al., 2012a, b). These broad wings, which appear to extend over several kpc and strengthen with increasing $\Sigma_{\mathrm{SFR}}$, have been interpreted as evidence of stellar feedback. Similar wings are found in massive, compact star-forming galaxies over a wider wavelength range (Wisnioski et al., 2018).

Larger spectroscopic surveys (up to $\sim$500 galaxies) of star-forming galaxies at $z=0.3-2$ use the low-ionization metal lines Mg 1, Mg 2, and Fe 2 to probe GWs (Weiner et al., 2009; Rubin et al., 2010; Banerji et al., 2011; Kornei et al., 2012; Erb et al., 2012; Martin et al., 2012; Bradshaw et al., 2013; Bordoloi et al., 2014; Rubin et al., 2014; Finley et al., 2017). Resonant transitions from outflows in absorption (blueshifted) and emission (redshifted), as well as corresponding non-resonant transitions, constrain basic outflow properties such as velocity and ionization state and, potentially, more complex structural parameters (Figure 4; (Prochaska et al., 2011)). These high-$z$, low-ionization outflows are broadly consistent with those found at low $z$: they have modest velocities (up to a few hundred km s*-1* on average); their properties (velocity and equivalent width) correlate with SFR, $\Sigma_{\mathrm{SFR}}$, and $M_{*}$; their detection rates in absorption average a few tens of percent, indicative primarily of the wind geometry (a high frequency of occurrence of non-spherical winds); they have estimated mass-loss rates of order the star formation rate, though with considerable uncertainty; and they are preferentially found in face-on galaxies (or their properties are more extreme in face-on galaxies), consistent with minor-axis flows. The velocities of these winds may increase with increasing $z$ for galaxies of given SFR, possibly due to increasing star formation rate surface density (Erb et al., 2012; Williams et al., 2015; Sugahara et al., 2017; Du et al., 2018), and with increasing ionization potential (Du et al., 2016, 2018; Chisholm et al., 2018). Extended, scattered emission (both resonant and non-resonant) from low-ionization species has been probed in detail in a handful of systems (Rubin et al., 2011; Kornei et al., 2013; Bordoloi et al., 2016; Finley et al., 2017) and is now recognized as a common feature of star-forming galaxies at these epochs (Coil et al., 2011; Martin et al., 2013; Finley et al., 2017). Pure resonant emission is seen at low SFR, transitioning to P-Cygni or pure-absorption resonant profiles plus non-resonant emission at high SFR, indicative of an increasing signature of outflowing gas (Finley et al., 2017). In a handful of gravitationally lensed systems, absorption-line analyses of more ions indicate $\alpha$ enhancement and more robustly constrain mass-loss rates (Pettini et al., 2002; Jones et al., 2018).

At moderate redshifts ($z\sim 0.7-0.8$), the presence of low-velocity GWs in poststarburst galaxies points to the possibility that these winds help to quench star formation (Coil et al., 2011). Very high-velocity winds (1000 km s*-1*) found in bluer, rarer post-starbursts (Tremonti et al., 2007) appear to be driven by very compact starbursts rather than AGN activity (Diamond-Stanic et al., 2012; Sell et al., 2014).

The down-the-barrel technique, which probes low-ionization outflows that emerge along the line-of-sight toward the galaxy disk, is complemented by transverse sightlines through galaxy halos toward background quasars or galaxies. LBGs at $z=2-3$ host both low- and high-ionization gas out to radii $\sim$100 kpc; because this gas is outflowing over a large solid angle (as inferred from down-the-barrel observations), it is plausibly doing so at large radii (Steidel et al., 2010). In other samples, Mg 2 and O 6 absorbers show a preference for alignment with galaxy major or minor axes, and the minor-axis gas is preferentially found in blue galaxies. This absorber alignment is suggestive of major-axis inflow and minor-axis outflow, and the connection to blue galaxies points to star formation as the power source (Kacprzak et al., 2012, 2015). Simple geometric outflow models of halo absorbers yield outflow properties that are consistent with those seen at low $z$ (Bouché et al., 2012; Schroetter et al., 2015; Muzahid et al., 2015).

The prominence of Ly$\alpha$ in the spectra of high $z$ galaxies makes it a tempting target for parameterizing outflows (e.g., (Erb, 2015)). However, as mentioned above (Section 2.1), radiation transfer effects make it an ambiguous indicator. Redshifted Ly$\alpha$ does typically accompany blueshifted low-ionization lines (e.g., (Steidel et al., 2010)). Star forming galaxies also show an increase in the velocity of Ly$\alpha$ as SFR and Ly$\alpha$ equivalent width increase (Hashimoto et al., 2013; Shibuya et al., 2014; Erb et al., 2014). However, whether this is due solely to changing outflow properties or instead to an increase in gas near the systemic redshift is unclear (Erb et al., 2014).

Finally, a handful of molecular gas detections of outflows at moderate-to-high $z$ are emerging. CO has been imaged in two high-velocity, apparently stellar GWs in post-starbursts at $z\sim 0.7$ (Geach et al., 2014, 2018). A long Herschel integration allowed detection of an OH outflow in absorption in a $z=2.3$ ULIRG (George et al., 2014). A serendipitous ALMA discovery of extremely broad CH+ in several $z\sim 2.5$ ULIRGs points to turbulent outflows (Falgarone et al., 2017). Tentative detections of broad, faint [C 2] line wings at $z=5.5$ hint at the possibility of stellar GWs in modestly star-forming galaxies at this epoch (Gallerani et al., 2018). Finally, an OH outflow exists in a gravitationally lensed, dusty galaxy at $z=5.3$ (Spilker et al., 2018).

5 Summary

We can say with reasonable certainty that GWs driven by energy from stellar processes are a common feature of galaxies with moderate-to-high star formation rates and/or surface densities out to $z\sim 2-3$. Stacking analyses of large rest-frame UV and optical spectroscopic surveys have established that the average star-forming galaxy has an ionized and/or neutral wind whose velocity scales with star formation rate, stellar mass, and possibly $\Sigma_{\mathrm{SFR}}$. At low redshift, these winds are most prominent in starburst galaxies that lie above the galaxy main sequence; at higher redshift, where galaxies on the main sequence have higher SFR, the situation may be different. Hints exist that GWs are common but simply hard to detect even in galaxies with low $\Sigma_{\mathrm{SFR}}$ (for instance, the low-surface-brightness Milky Way GW).

The correlations between outflow and galaxy properties found in some of the first large surveys of stellar GWs have been verified and refined by larger and more diverse samples and different gas probes. Besides serving merely as input to parameterizations of outflows in numerical simulations, measurements of GWs can now be compared to the predicted properties of GWs from simulations that better implement the physics of stellar feedback.

Detailed, multiwavelength studies of star-forming galaxies continue to reveal new layers of GWs. Most notable is that stellar GWs entrain large quantities of molecular gas, including dense clumps, and loft dust and soot (PAHs) far above the galactic disk. The promising technique of combining resonant-line absorption and emission with non-resonant re-emission channels has been successfully used to detect winds at high $z$ and may prove a powerful probe of GW structure and extent when widely deployed. Observations of a wider range of galaxies besides the usual suspects (e.g., M82) with 3D imaging spectroscopy shows that a complex, multiphase structure of filaments of dusty ionized and neutral gas collimated along the minor axis is a common feature of GWs. Transverse-sightline spectroscopy and correlations with galaxy inclination at a variety of redshifts bolster this picture. Finally, increasingly in-depth studies of local galaxies with extended Ly$\alpha$ and LyC may eventually put meaningful constraints on how outflows contribute to reionization and help interpret high-$z$ observations of Ly$\alpha$.

Future progress will occur on a variety of fronts. At low $z$, a new generation of ongoing multi-object IFS surveys (SAMI Galaxy Survey, MaNGA) will soon produce results on thousands of nearby galaxies. Future, much larger IFS surveys are being planned (using, e.g., Hector; (Bryant et al., 2016)). Sensitive, wide-field IFS instruments on large telescopes (such as MUSE and KCWI, the Keck Cosmic Web Imager) will probe the full extent of GWs in nearby, well-resolved targets and enable efficient, spatially resolved characterization of many galaxies at once in high-$z$ deep fields. Next generation multi-object spectroscopy surveys (e.g., the Dark Energy Spectroscopic Instrument, or DESI, Survey and 4MOST, the 4-metre Multi-Object Spectroscopic Telescope) will increase the fidelity of stacking analyses over a wider range of redshift and galaxy properties.

Continued measurements with ALMA, particularly at high spatial resolution, will provide more detailed understanding of the structure and chemistry of molecular gas in outflows. ALMA will also certainly expand on its currently short list of detections of high-$z$ stellar GWs. The high resolution and sensitivity of the James Webb Space Telescope (JWST) in the MIR will undoubtedly produce useful measurements of molecular gas and dust in stellar GWs, as well. However, the spatial resolution and sensitivity of JWST is likely to provide the most dramatic advances in measuring the properties of outflows in high-$z$ galaxies by characterizing them in individual main-sequence galaxies at $z\sim 2-3$ and detecting them at very high $z$, where their impact could be especially significant but where measurements currently do not exist.

Finally, we note two areas of study that have seen little recent progress, but whose prospects should eventually rise. Measurements of the radio-emitting plasma in GWs are very rare except for a few recent detections (Kepley et al., 2010; Leslie et al., 2017; Ramírez-Olivencia et al., 2018). The next generation of wide-field radio arrays may make this a growth area. The field of X-rays studies of GWs has also lain fallow, with a few exceptions (e.g., (O’Sullivan et al., 2014; McQuinn et al., 2018)). The hottest gas phase of GWs, which may drive the outflows in starburst galaxies, has proven extremely difficult to detect except in the nearest cases (Strickland and Heckman, 2009). More sensitive X-ray telescopes in the coming two decades will eventually lead to a better characterization of this pivotal component.

Acknowledgements.

The author thanks the referees for their feedback. \fundingThis research received no external funding. The author is supported by the J. Lester Crain chair at Rhodes College. \conflictsofinterestThe author declares no conflict of interest. \abbreviationsThe following abbreviations are used in this manuscript:

[TABLE]

\reftitleReferences