Braving the Storm: Quantifying Disk-wide Ionized Outflows in the Large Magellanic Cloud with ULLYSES

The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas (T ~ 10^4-5 K) in HI, SII, SiIV, and CIV using 110 stellar sight lines from the Hubble Space Telescope's Ultraviolet Legacy Library of Young Stars as Essential Standards program. We develop a continuum fitting algorithm based on the concept of Gaussian process regression and identify reliable LMC interstellar absorption over v_helio = 175-375 km s-1. Our analyses show disk-wide ionized outflows in Si iv and C iv across the LMC with bulk velocities of divided by v out, bulk divided by similar to 20-60 km s-1, which indicates that most of the outflowing mass is gravitationally bound. The outflows' column densities correlate with the LMC's star formation rate surface densities (Sigma_SFR), and the outflows with higher Sigma SFR tend to be more ionized. Considering outflows from both sides of the LMC as traced by C iv, we conservatively estimate a total outflow rate of Mout greater than or similar to 0.03 Msun yr-1 and a mass-loading factor of eta greater than or similar to 0.15. We compare the LMC's outflows with those detected in starburst galaxies and simulation predictions, and find a universal scaling relation of divided by v_out,bulk divided by proportional to Sigma_SFR~0.23 over a wide range of star-forming conditions (Sigma_SFR~10^-4.5-10^2 Msun yr-1 kpc-2). Lastly, we find that the outflows are corotating with the LMC's young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way.