Turbulence in the interstellar medium (ISM) plays an important role in many physical processes, including forming stars and shaping complex ISM structures. In this work, we investigate the HI turbulent properties of the Small Magellanic Cloud (SMC) to reveal what physical mechanisms drive the turbulence and at what scales. Using the high-resolution HI data of the Galactic ASKAP (GASKAP) survey and multi-point structure functions (SF), we perform a statistical analysis of HI turbulence in 34 subregions of the SMC. Two-point SFs tend to show a linear trend, and their slope values are relatively uniform across the SMC, suggesting that large-scale structures exist and are dominant in the two-point SFs. On the other hand, seven-point SF enables us to probe small-scale turbulence by removing large-scale fluctuations, which is difficult to achieve with the two-point SFs. In the seven-point SFs, we find break features at scales of 34-84 pc, with a median scale of $\sim$50 pc. This result indicates the presence of small-scale turbulent fluctuations in the SMC and quantifies its scale. In addition, we find strong correlations between slope values of the seven-point SFs and the stellar feedback-related quantities (e.g., H$\alpha$ intensities, the number of young stellar objects, and the number of HI shells), suggesting that stellar feedback may affect the small-scale turbulent properties of the HI gas in the SMC. Lastly, estimated sonic Mach numbers across the SMC are subsonic, which is consistent with the fact that the HI gas of the SMC primarily consists of the warm neutral medium.

With the aim of evaluating the roles of the cold neutral medium (CNM) in the cloud-scale baryon cycle, we perform a high-resolution study of the CNM in and around the extreme star-forming region 30 Doradus (30 Dor). For our study, we use Galactic Australian Square Kilometre Array Pathfinder H I Survey data and produce H I emission and absorption cubes on 7 pc scales. To examine the CNM structures toward 30 Dor, we decompose the H I absorption cube into 862 Gaussian components and find that these components are distributed at four velocity ranges (B1, B2, B3, and B4, respectively): 200$-$230 km s$^{-1}$, 230$-$260 km s$^{-1}$, 260$-$277 km s$^{-1}$, and 277$-$300 km s$^{-1}$. We derive line-of-sight average spin temperatures and opacity-corrected total H I column densities and show that the B1$-$B4 structures have systematically different properties, indicating that they are physically distinct. As for the nature of the observed CNM structures, we find that B2 is associated with the main dense structure where ionized, atomic, and molecular gases are concentrated. B3 and B4 trace inflows whose combined mass flux rate of 0.14 $M_{\odot}$ yr$^{-1}$ is comparable to the current star formation rate, while B1 probes outflows with a much lower mass flux rate of 0.007 $M_{\odot}$ yr$^{-1}$. Interestingly, the H I column densities in B1$-$B4 are nearly uniform with a factor of two spatial variations, implying the presence of H I shielding layers for H$_{2}$ formation.