As some of the only Lyman continuum (LyC) emitters at z~0, Green Pea (GP) galaxies are possible analogs of the sources that reionized the universe. We present HST COS spectra of 13 of the most highly ionized GPs, with [O III]/[O II]=6-35, and investigate correlations between Ly-alpha, galaxy properties, and low-ionization UV lines. Galaxies with high [O III]/[O II] have higher H-alpha equivalent widths (EWs), and high intrinsic Ly-alpha production may explain the prevalence of high Ly-alpha EWs among GPs. While Ly-alpha escape fraction is closely linked to low gas covering fractions, implying a clumpy gas geometry, narrow Ly-alpha velocity peak separation (delta_v,LyA) correlates with the ionization state, suggesting a density-bounded geometry. We therefore suggest that delta_v,LyA may trace the residual transparency of low-column-density pathways. Metallicity is associated with both [O III]/[O II] and delta_v,LyA. This trend may result from catastrophic cooling around low-metallicity star clusters, which generates a compact geometry of dense clouds within a low-density inter-clump medium. We find that the relative strength of low-ionization UV emission to absorption correlates with Ly-alpha emission strength and is related to Ly-alpha profile shape. However, as expected for optically thin objects, the GPs with the lowest delta_v,LyA show both weak low-ionization emission and weak absorption. The strengths of the low-ionization absorption and emission lines in a stacked spectrum do not correspond to any individual spectrum. Galaxies with high [O III]/[O II] contain a high fraction of LyC emitter candidates, but [O III]/[O II] alone is an insufficient diagnostic of LyC escape.

Mechanisms regulating the escape of Ly$\alpha$ photons and ionizing radiation remain poorly understood. To study these processes we analyze VLA 21cm observations of one Green Pea (GP), J160810+352809 (hereafter J1608), and HST COS spectra of 17 GP galaxies at $z<0.2$. All are highly ionized: J1608 has the highest [O III] $\lambda5007$/[O II] $\lambda3727$ for star-forming galaxies in SDSS, and the 17 GPs have [O III]/[O II] $\geq6.6$. We set an upper limit on J1608's HI mass of $\log M_{HI}/M_\odot=8.4$, near or below average compared to similar mass dwarf galaxies. In the COS sample, eight GPs show Ly$\alpha$ absorption components, six of which also have Ly$\alpha$ emission. The HI column densities derived from Ly$\alpha$ absorption are high, $\log N_{HI}/$cm$^{-2}=19-21$, well above the LyC optically thick limit. Using low-ionization absorption lines, we measure covering fractions ($f_{\mbox{cov}}$) of $0.1-1$, and find that $f_{\mbox{cov}}$ strongly anti-correlates with Ly$\alpha$ escape fraction. Low covering fractions may facilitate Ly$\alpha$ and LyC escape through dense neutral regions. GPs with $f_{\mbox{cov}}\sim1$ all have low neutral gas velocities, while GPs with lower $f_{\mbox{cov}}=0.2-0.6$ have a larger range of velocities. Conventional mechanical feedback may help establish low $f_{\mbox{cov}}$ in some cases, whereas other processes may be important for GPs with low velocities. Finally, we compare $f_{\mbox{cov}}$ with proposed indicators of LyC escape. Ionizing photon escape likely depends on a combination of neutral gas geometry and kinematics, complicating the use of emission-line diagnostics for identifying LyC emitters.

By clearing neutral gas away from a young starburst, superwinds may regulate the escape of Lyman continuum (LyC) photons from star-forming galaxies. However, models predict that superwinds may not launch in the most extreme, compact starbursts. We explore the role of outflows in generating low optical depths in the Green Peas (GPs), the only known star-forming population with several confirmed and candidate LyC-leaking galaxies. With Hubble Space Telescope UV spectra of 25 low-redshift GPs, including new observations of 13 of the most highly ionized GPs, we compare the kinematics of UV absorption lines with indirect HI optical depth diagnostics: Ly-alpha escape fraction, Ly-alpha peak separation, or low-ionization absorption line equivalent width. The data suggest that high ionization kinematics tracing superwind activity may correlate with low optical depth in some objects. However, the most extreme GPs, including many of the best candidate LyC emitters with weak low-ionization absorption and strong, narrow Ly-alpha profiles, show the lowest velocities. These results are consistent with models for suppressed superwinds, which suggests that outflows may not be the only cause of LyC escape from galaxies.