The NSF-sponsored Undergraduate ALFALFA Team (UAT) promotes long-term collaborative research opportunities for faculty and students from 23 U.S. public and private primarily undergraduate institutions (PUIs) within the context of the extragalactic ALFALFA HI blind legacy survey project. Over twelve project years of partnering with Arecibo and Green Bank Observatories, the UAT has had a demonstrable impact on the health of a legacy astronomy project, science education, and equity/inclusion in astronomy, with successful outcomes for 373 UAT students (39% women; ~30% members of underrepresented groups) and 34 faculty (44% women). The UAT model is adaptable to many large scientific projects and can be supported by relatively modest funding. We recommend that granting agencies identify funding resources to support the model, either as an add-on to legacy grant support or as a stand-alone funding source. This could include encouragement of UAT-like components in large scale projects currently being developed, such as the LSST and TMT. By doing this, we will recognize the high numbers of astronomy research-trained heavy-teaching-load faculty at PUIs as an under-utilized resource of the astronomy community (see also White Paper by Ribaudo et al.). These members of our community have the skills and the strong desire to contribute meaningfully to their field, as well as the ability to encourage and interact closely with many talented and motivated undergraduate students from all backgrounds.

We use the ALFALFA HI survey to examine whether the cold gas reservoirs of galaxies are inhibited or enhanced in large-scale filaments. Our sample includes 9947 late-type galaxies with HI detections, and 4236 late-type galaxies with well-determined HI detection limits that we incorporate using survival analysis statistics. We find that, even at fixed local density and stellar mass, and with group galaxies removed, the HI deficiency of galaxies in the stellar mass range 8.5 <log(M/Mo) < 10.5 decreases with distance from the filament spine, suggesting that galaxies are cut off from their supply of cold gas in this environment. We also find that, at fixed local density and stellar mass, the galaxies that are the most gas-rich are those in small, correlated "tendril" structures within voids: although galaxies in tendrils are in significantly denser environments, on average, than galaxies in voids, they are not redder or more HI deficient. This stands in contrast to the fact that galaxies in tendrils are more massive than those in voids, suggesting a more advanced stage of evolution. Finally, at fixed stellar mass and color, galaxies closer to the filament spine, or in high density environments, are more deficient in HI. This fits a picture where, as galaxies enter denser regions, they first lose HI gas and then redden as star formation is reduced.

We present the HI content of galaxies in nearby groups and clusters as measured by the 70% complete Arecibo Legacy Fast-ALFA (ALFALFA) survey, including constraints from ALFALFA detection limits. Our sample includes 22 systems at distances between 70-160 Mpc over the mass range 12.5<log M/M_sun<15.0, for a total of 1986 late-type galaxies. We find that late-type galaxies in the centers of groups lack HI at fixed stellar mass relative to the regions surrounding them. Larger groups show evidence of a stronger dependence of HI properties on environment, despite a similar dependence of color on environment at fixed stellar mass. We compare several environment variables to determine which is the best predictor of galaxy properties; group-centric distance r and r/R_200 are similarly effective predictors, while local density is slightly more effective and group size and halo mass are slightly less effective. While both central and satellite galaxies in the blue cloud exhibit a significant dependence of HI content on local density, only centrals show a strong dependence on stellar mass, and only satellites show a strong dependence on halo mass. Finally, we see evidence that HI is deficient for blue cloud galaxies in denser environments even when both stellar mass and color are fixed. This is consistent with a picture where HI is removed or destroyed, followed by reddening within the blue cloud. Our results support the existence of pre-processing in isolated groups, along with an additional rapid mechanism for gas removal within larger groups and clusters, perhaps ram-pressure stripping.

We use fractal analysis to systematically study the clustering strength of the distribution of stars, HII regions, molecular gas, and individual giant molecular clouds in M33 over a wide range of spatial scales. We find a clear transition from a scale-free behavior at small spatial scales to a nearly uniform distribution at large scales. The transition region lies in the range 500-1000 pc and it separates the regime of small-scale turbulent motion from that of large-scale galactic dynamics. The three-dimensional fractal dimension of bright young stars and molecular gas at small spatial scales is similar to or less than 1.9 indicating that the interstellar medium in M33 is on average much more fragmented and irregular than the in the Milky Way.

We study the distribution of stars, HII regions, molecular gas, and individual giant molecular clouds in M33 over a wide range of spatial scales. The clustering strength of these components is systematically estimated through the fractal dimension. We find scale-free behavior at small spatial scales and a transition to a larger correlation dimension (consistent with a nearly uniform distribution) at larger scales. The transition region lies in the range 500-1000 pc. This transition defines a characteristic size that separates the regime of small-scale turbulent motion from that of large-scale galactic dynamics. At small spatial scales, bright young stars and molecular gas are distributed with nearly the same three-dimensional fractal dimension (Df <= 1.9), whereas fainter stars and HII regions exhibit higher values (Df = 2.2-2.5). Our results indicate that the interstellar medium in M33 is on average more fragmented and irregular than in the Milky Way.

The autocorrelation function provides an objective test for the existence of special scales in the hierarchical clustering of young stars. We apply this measure to single-star photometry for the brightest main sequence stars in the Small Magellanic Cloud (SMC), the Large Magellanic Cloud (LMC), M33, and M31, using data from the Magellanic Clouds Photometric Survey and the Massey Local Group Survey. Our primary result is the identification of a transition to a higher correlation dimension (weaker clustering) at one kpc in the LMC and M31, and at 300 pc in M33. We suggest that this transition marks the large-scale regime where disk geometry and dynamics set the scale for structure. On smaller scales, the correlation functions for each galaxy are scale-free over at least two orders of magnitude, with a projected correlation dimension varying from 1.0 for M31 to 1.8 for the SMC. This variation is probably caused by a combination of differences in stellar ages and masses, physical environment, and extinction.

We present the correlation dimension of resolved young stars in four actively star-forming dwarf galaxies that are sufficiently resolved and transparent to be modeled as projections of three-dimensional point distributions. We use data in the Hubble Space Telescope archive; photometry for one of them, UGCA 292, is presented here for the first time. We find that there are statistically distinguishable differences in the nature of stellar clustering among the sample galaxies. The young stars of VII Zw 403, the brightest galaxy in the sample, have the highest value for the correlation dimension and also the most dramatic decrease with logarithmic scale, falling from $1.68\pm0.14$ to $0.10\pm0.05$ over less than a factor of ten in $r$. This decrease is consistent with the edge effect produced by a projected Poisson distribution within a 2:2:1 ellipsoid. The young stars in UGC 4483, the faintest galaxy in the sample, exhibit very different behavior, with a constant value of about 0.5 over this same range in $r$, extending nearly to the edge of the distribution. This behavior may indicate either a scale-free distribution with an unusually low correlation dimension, or a two-component (not scale-free) combination of cluster and field stars.