VERTICO V: The environmentally driven evolution of the inner cold gas discs of Virgo cluster galaxies

Abstract The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (H i) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival H i data, we study the H i and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ( 
$\Sigma_{\star}$
 ), H i ( 
$\Sigma_{\text{H}\,{\small\text{I}}}$
 ), and molecular gas ( 
$\Sigma_{\text{mol}}$
 ) surface densities. Adopting H i deficiency as a measure of environmental impact, we find evidence that, in addition to removing the H i at large radii, the cluster processes also lower the average 
$\Sigma_{\text{H}\,{\small\text{I}}}$
 of the remaining gas even in the central 
$1.2\,$
 kpc. The impact on molecular gas is comparatively weaker than on the H i, and we show that the lower 
$\Sigma_{\text{mol}}$
 gas is removed first. In the most H i-deficient galaxies, however, we find evidence that environmental processes reduce the typical 
$\Sigma_{\text{mol}}$
 of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of 
$\Sigma_{\text{H}\,{\small\text{I}}}$
 or 
$\Sigma_{\text{mol}}$
 in H i-deficient galaxies. Using the ratio of 
$\Sigma_{\text{mol}}$
 -to- 
$\Sigma_{\text{H}\,{\small\text{I}}}$
 in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results.

Indeed, there is significant evidence that hydrodynamical stripping is a dominant mechanism for removing gas from cluster galaxies (see the reviews by Cortese et al., 2021;Boselli et al., 2022).In the canonical outside-in picture of stripping, the circumgalactic medium is first removed and cuts galaxies off from further accretion (Larson et al., 1980), before progressing to strip the ISM directly (e.g., Kenney et al., 2004;Chung et al., 2007;Cortese et al., 2010).Further, cometary tails of stripped gas from the ISM have been observed with multi-wavelength tracers from the X-rays to the radio continuum (e.g.Sivanandam et al., 2014;Jáchym et al., 2014Jáchym et al., , 2017;;Moretti et al., 2020b;Longobardi et al., 2020;Roberts et al., 2021;Ignesti et al., 2022), and reproduced in controlled and cosmological simulations of galaxy evolution (e.g.Tonnesen & Bryan, 2010;Marasco et al., 2016;Yun et al., 2019).
Most of our understanding of how environmental processes impact the cold gas in galaxies comes from observations of the 21-cm emission line from neutral atomic hydrogen (H ), as its large radial extent makes it a sensitive tracer of external environmental mechanisms.In the case of hydrodynamical stripping, galaxies exhibit clear truncation of their H distribution at a radius that depends on the strength of the ram pressure (Gunn & Gott, 1972;Cayatte et al., 1990;Roediger & Hensler, 2005;Tonnesen & Bryan, 2009;Abramson et al., 2011;Yoon et al., 2017;Reynolds et al., 2021Reynolds et al., , 2022)).The impact of these environmental processes on the molecular hydrogen (H 2 ) is less clear; however, truncation of the star formation in stripped galaxies (e.g., Koopmann & Kenney, 2004;Boselli et al., 2006;Cortese et al., 2012;Fossati et al., 2013Fossati et al., , 2018) ) indicates that the H 2 is also affected once stripping proceeds within the stellar disc.Further, the presence of star formation inside this truncation radius means that some gas can survive the stripping process.Whether the environment has also influenced this surviving gas, which would also affect the timescale for star formation to halt (e.g.Boselli et al., 2014;Cortese et al., 2021), remains unknown.As this surviving gas is often dominated by H 2 (Boselli et al., 2014;Cortese et al., 2016;Loni et al., 2021), understanding the impact of environmental processes on the molecular ISM in galaxies is essential to completing our picture of galaxy evolution in dense environments.
Under most ISM conditions, however, H 2 cannot be directly observed as it is a symmetric molecule without a permanent dipole.Instead, most of our understanding of the molecular ISM (from here, 'molecular gas' or subscript 'mol') comes from observing the rotational transitions of CO and its isotopologues (Bolatto et al., 2013;Saintonge et al., 2017).These observations have revealed that molecular gas is not affected by environmental processes as strongly as H .This resistance is primarily because a large fraction of molecular gas exists in giant molecular clouds (GMCs), and it is typically more centrally concentrated where the gravitational restoring force is stronger, making it more resistant to direct stripping than the H (e.g., Engargiola et al., 2003;Bolatto et al., 2008;Yamagami & Fujita, 2011;Stevens et al., 2021).Instead, environmental processes typically lower the molecular gas content in the outskirts of galaxies first (Lee et al., 2020;Brown et al., 2021;Zabel et al., 2022), while significant reduction of the molecular gas mass of galaxies is only observed once galaxies have lost a large fraction of their total H mass (Fumagalli et al., 2009;Boselli et al., 2014).
Although this picture is consistent with the outside-in stripping scenario, there is growing evidence that it is more complicated than just radial truncation.Simulations have shown that the multi-phase nature of the ISM allows the ICM to stream through low-density regions of the surviving gas inside the H truncation radius, allowing environmental processes to affect the cold gas across the disc (Tonnesen & Bryan, 2009;Choi et al., 2022).Disturbed molecular gas morphologies in cluster galaxies can resemble the H (e.g., Lee et al., 2017), and highresolution studies of individual objects have shown GMCs decoupled from large plumes of trailing gas (e.g., Cramer et al., 2020).Further, the mass of molecular gas detected in the tails of stripped gas trailing behind cluster galaxies requires some level of mass-loading from the ISM (e.g., Sivanandam et al., 2010;Moretti et al., 2018;Lee & Chung, 2018;Jáchym et al., 2019), and there is evidence for the direct stripping of cold ISM from dust observations (e.g., Crowl et al., 2005;Longobardi et al., 2020).Alternately, stripping of the molecular gas may not be a dominant process, and the reduction of molecular gas mass could be due to the preferential removal of H across the disc, leaving insufficient fuel to maintain efficient molecule formation (Fumagalli et al., 2009;Tonnesen & Bryan, 2009;Zabel et al., 2022;Villanueva et al., 2022).To further complicate the picture, there is evidence for increased molecular gas formation efficiency due to gas compression during the ICM-ISM interaction (e.g., Tonnesen & Bryan, 2009;Moretti et al., 2020a,b).Clearly, our picture of how the environment impacts the molecular gas in galaxies is incomplete.
In this work, we present a new analysis of the cold gas properties within the stellar components of Virgo cluster galaxies using observations from the Virgo Environment Traced in CO(2-1) (VERTICO, Brown et al., 2021) ALMA ACA a large program.The Virgo cluster is an excellent target for furthera Atacama Large Millimetre/Sub-millimetre Array and Atacama Compact Array ing our understanding of the environment's impact on galaxies' cold gas content.At a distance of 16.5 Mpc (Mei et al., 2007) it is close, allowing for detailed studies of individual galaxies, and it has been observed comprehensively across nearly the entire electromagnetic spectrum (York et al., 2000;Chung et al., 2009;Davies et al., 2010;Urban et al., 2011;Boselli et al., 2011Boselli et al., , 2018a)).The presence of infalling sub-structures also means that there are populations of galaxies at various stages of infall, and thus environmental impact (Gavazzi et al., 1999;Mei et al., 2007;Kim et al., 2016;Yoon et al., 2017).Additionally, the VERTICO sample was selected from the VLA Imaging of Virgo galaxies in Atomic gas survey (VIVA, Chung et al., 2009) of H emission.Together, these data enable us to study how the environmental processes influence the cold gas properties of galaxies in the Virgo cluster on ∼1.2-kpc scales.
We present our study using the framework of spatially resolved scaling relations (e.g., Bigiel et al., 2008;Morselli et al., 2020;Ellison et al., 2021a,b;Pessa et al., 2022;Abdurro'uf et al., 2022a,b;Baker et al., 2022), in particular those between stellar, H , and molecular gas mass surface densities.While azimuthally averaged radial profiles are useful for increasing sensitivity to faint emission, in this work, we focus on the bright emission within the optical components of galaxies.Thus, comparing the physical quantities of kpc-scale regions within galaxies increases our statistics, and reduces the effects of averaging over larger areas that include many different physical regions.For analyses of how the environment impacts the radial profiles of cold gas content in VERTICO galaxies, we refer the reader to the works of Brown et al. (2021); Zabel et al. (2022), andVillanueva et al. (2022).
This paper is structured as follows.In §2 we describe our VERTICO sample, control sample, the homogenisation of the data, and ancillary data products.In §3 we investigate the resolved scaling relations with all galaxies combined and discuss some of the caveats of this approach.We present our main results in §4, where we study how the environment impacts the spatially resolved H and molecular gas content of galaxies across their stellar discs and within the truncation of their gas discs.In §5 we investigate the physical mechanisms responsible for our observations, and in §6 we summarise our results.
As described in Brown et al. (2021), all physical quantities used in this work were computed using a Kroupa (2001) stellar initial mass function, a ΛCDM cosmology with h = 0.7, Ω M = 0.3 and Ω Λ = 0.7, and we adopt a constant distance for Virgo galaxies of 16.5 Mpc (Mei et al., 2007).

VERTICO sample
We selected our main galaxy sample from the VERTICO survey (Brown et al., 2021), which obtained CO(2-1) observations b for 51 Virgo cluster galaxies extracted from the VIVA survey.We refer the reader to Brown et al. (2021) for a detailed description of the parent sample selection and data reduction.In this work, we focus on massive (log(M /M ) ≥ 10) galaxies to minimise metallicity-dependent uncertainties in tracing the molecular gas content of galaxies with CO(2-1) emission, and we enforce an optical inclination cut of i ≤ 70 • to avoid large de-projection uncertainties.We also removed two galaxies that met these criteria, NGC 4772 and NGC 4698, as they are passive, early-type discs that have recently accreted their H gas (Haynes et al., 2000;Cortese & Hughes, 2009;Yoon et al., 2017), and thus their H and molecular gas properties do not reflect the action of the Virgo environment on their cold gas content.Our final VERTICO sample consists of 18 galaxies, and the properties of the galaxies used in this work are summarised in Table 1, along with their sources and derivations that are also described below.
The native spatial resolution of VERTICO data is ∼ 8arcsec, and data products with 9-arcsec and 15-arcsec resolution were also created using the routine of the Common Astronomy Software Applications package (CASA, McMullin et al., 2007).We used the coarser 15-arcsec data products sampled with 8-arcsec spaxels, which are closer to the spatial resolution of the H data.The VIVA H data used in this work were re-imaged by author A. Chung to circularise the synthesised beam and match the VERTICO data products as closely as possible, resulting in datacubes with 8 arcsec spaxels and synthesised beams varying from 15.2-20 arcsec, with an average of 17 arcsec.At the assumed 16.5 Mpc distance to Virgo, the 15/8-arcsec beam/spaxels of the VERTICO data correspond to a 1200/650 pc physical resolution.
Signal masks were created for the CO(2-1) and H datacubes as described in Brown et al. (2021, but see also Sun et al. 2018 andLeroy et al. 2021), and moment 0 maps created for each galaxy.The H mass in each spaxel was calculated assuming optically thin emission (e.g., Meyer et al., 2017), where D L is the luminosity distance to the object (assumed to be 16.5 Mpc for all VERTICO galaxies), z is the redshift, and S H is the velocity-integrated spectral flux density (in Jy km s -1 ).Unless otherwise specified, we do not correct our H masses for the Helium abundance.The velocity-integrated CO(2-1) flux density (S CO ) was used to calculate the total CO(2-1) luminosity using the method of Solomon & Vanden Bout (2005), where ν is the observed frequency.The molecular gas mass was then computed using where R 21 is the integrated flux ratio of the CO(2-1) and CO(1-0) lines and α CO is the molecular gas mass to CO(1-0) luminosity correction factor.We used a constant R 21 = 0.8, as derived by Brown et al. (2021), and a constant α CO = Table 1.Properties of the galaxies used in this work.VERTICO galaxies are contained in the top section, while the HERACLES control sample is in the bottom segment.(Mei et al., 2007;Walter et al., 2008); 3) Position angle (Brown et al., 2021); 4) inclination (Brown et al., 2021); 5) B-band optical radius (Chung et al., 2009;Walter et al., 2008); 6) log 10 integrated stellar mass (Leroy et al., 2019); 7) log 10 integrated star-formation rate (Leroy et al., 2019); 8) H deficiency parameter (Chung et al., 2009, this work); 9) average log Σ at the H -radius Σ H = 1 M pc -2 (this work).

Galaxy
4.35 M pc -2 based on observations of the Milky Way disc by Bolatto et al. (2013), which includes a 36 per cent contribution from Helium.The H and molecular gas mass of spaxels were converted to surface densities Σ H and Σ mol in units of M pc -2 using the 650-pc physical pixel size.
Global stellar masses were drawn from the z = 0 Multiwavelength Galaxy Synthesis (z0MGS, Leroy et al., 2019) sample, a collection of GALEX (Galaxy Evolution Explorer, Martin et al., 2005), and WISE (Wide-field Infrared Survey Explorer, Wright et al., 2010) data for ∼ 11 000 galaxies with D < 50 Mpc and brighter than M B = -18 mag.Stellar-mass estimates were computed from WISE band-1 observations using the WISE band-1 and -3 colour to determine the massto-light ratio, as described in Leroy et al. (2019).The same prescriptions were applied on a pixel-by-pixel basis to create resolved maps of stellar surface density (Σ ) at the same resolution as the H and molecular gas data.Photometric position angles (PA) and inclinations (i) for each galaxy were calculated by fitting Kron apertures to SDSS c r-band images (Brown et al., 2021) and assuming an intrinsic axial ratio of q = 0.2, and B-band 25 th mag arcsec -2 optical radii (R 25 ) were sourced from Chung et al. (2009).All surface densities used in this work were de-projected to face-on values using a factor of cos(i).
To quantify the impact of the environment on VERTICO galaxies, we used the type-independent H deficiency parameter (def-H ) for each galaxy from Chung et al. (2009), which compares the observed H content to the average content in isolated galaxies with the same optical size.It is defined as where Σ H , obs = S H /(2 R 25 ) 2 is the mean surface density within the optical disc.Following Chung et al. (2009), this c Sloan Digital Sky Survey (York et al., 2000;Abazajian et al., 2009) definition uses S H in units of Jy km s -1 and R 25 in arcminutes, and log Σ H = 0.37 is the average type-independent value of log(M H /D 2 25 ) from Haynes & Giovanelli (1984, Table 4, 6.81) after a unit conversion from M kpc -2 .A larger def-H value corresponds to a smaller H mass compared to isolated galaxies, and the typical scatter in the def-H distribution for isolated galaxies is ∼ 0.2 -0.3 dex (Haynes & Giovanelli, 1984;Jones et al., 2018, thus automatically encapsulating the variations in the spatial distribution of H surface density observed in isolated galaxies).In this work, we define 'H -deficient' galaxies as those with def-H ≥ 0.3 dex, namely, a factor of 2 more H -poor than expected.
We are primarily interested in the properties of the gas within the optical disc of galaxies, so we masked each galaxy to only include pixels/spaxels within the B-band R 25 , where the radius of each spaxel was calculated from the photometric PA and inclination.While this allows us to study the impact of the environment across the discs of galaxies, we are also interested in how the environment impacts the surviving cold gas disc.The extent of the H gas is a good tracer of the gas truncation radius, outside of which stripping is usually assumed to be effective.As our analysis uses the relationships between Σ , Σ H , and Σ mol , it is useful to quantify the size of the H reservoir in terms of Σ .
We computed radial profiles of log Σ H and log Σ using 4 arcsec annuli, which we found gave the best consistency with the Chung et al. (2009) VIVA measurements.Defining the H truncation radius (R H ) to be where the radial average Σ H = 1 M pc -2 (e.g., Wang et al., 2014Wang et al., , 2016;;Stevens et al., 2019a), we used the value of log Σ at this radius to determine the average log Σ at R H , log Σ H .In the most H -rich galaxies, R H can fall outside of the B-band R 25 , but this does not affect any of our results as these galaxies are not the focus of our analysis.Further, in all galaxies, less than 5 per cent of spaxels with Σ ≥ log Σ H have radius > R H , meaning trends with higher (lower) Σ trace trends with smaller (larger) galactocentric radius.We also note here that while the H observations cover the optical disc of each galaxy, the size of the VERTICO CO(2-1) observation footprints were based on the R H of each target.Subsequently, there is some variation in the physical coverage of the molecular gas in galaxies (0.44 -1 × R 25 ), and we accounted for this in our analysis of the molecular gas content.
While the VERTICO sample contains some H -normal systems, environmental processes can begin to impact galaxies before they become H -deficient. Thus, we chose to draw our control sample from a different survey to ensure that we had an H -normal reference that was as free from environmental effects as possible.We include the H -normal VERTICO systems in some of our figures (Figs. 1,3,and 6) and discuss them briefly in §4, but the main focus of our analysis is the H -deficient galaxies.

HERACLES control sample
We extracted a control sample of galaxies from the HERA d CO(2-1) Legacy Survey (HERACLES, Leroy et al., 2009).HERACLES targeted the CO(2-1) emission line in 33 galaxies from The H Nearby Galaxy Survey (THINGS, Walter et al., 2008), a sample of nearby (D < 15 Mpc), late-type galaxies spanning from dwarfs to massive spirals.After excluding nondetections and five galaxies that overlap with VERTICO, there are 25 HERACLES galaxies with public data (for more details see Brown et al., 2021).Physical quantities for HERACLES galaxies were calculated using the same methods and sources as described for VERTICO, and B-band R 25 were drawn from Walter et al. (2008).We did not compute log Σ H for HER-ACLES galaxies as the R H was nearly always outside R 25 .
The HERACLES native resolution datacubes have a synthesised beam-width of 13 arcsec, while the THINGS data are more varied depending on the VLA configuration.Before calculating physical quantities, the HERACLES and THINGS data were smoothed using the routine of CASA to match the physical spaxel size and resolution of the VERTICO data, and circularise the synthesised beam.As the distances to HERACLES sample galaxies vary, this left 11 galaxies that could be matched to a 1200-pc physical resolution in both CO(2-1) and H .
To ensure that our HERACLES control sample is a reasonable representation of star-forming galaxies, we restricted this sample to galaxies more star-forming than 1σ below the star-forming main-sequence using the parametrisation from the z0MGS sample (σ = 0.36, with star-formation rates (SFR) from the same catalogue, Leroy et al., 2019).After the stellar mass and inclination selection cuts, our final HERACLES control sample consists of eight galaxies, listed in Table 1.We compared the integrated properties of the final control sample galaxies to a subset of star-forming galaxies from the xCOLD GASS e (Saintonge et al., 2017) sample, selected using the same SFR threshold.While the HERACLES galaxies lie on or slightly above (< 1σ) the locus of the star-forming main-sequence, they span the same range in log M H /M and log M MG /M as the xCOLD GASS sample, indicating that there is no bias toward H -or molecule-rich systems in our control sample.
Last, we briefly mention the detection limits of each dataset after the homogenisation of the spatial resolutions and the masking procedures.The HERACLES:THINGS data are shallower than the VERTICO:VIVA data in both H (Σ H ≥ 1 -3 M pc -2 within B-band R 25 compared to Σ H ≥ 1 M pc -2 ) and molecular gas (Σ mol ≥ 1 -3 M pc -2 compared to Σ mol ≥ 0.4 -1 M pc -2 ).Given that we are interested in environmentdriven reductions in both H and molecular gas surface densities, the fact that the VERTICO:VIVA data are deeper ensures that our findings are not due to differences in sensitivity between the two datasets.The shallower nature of HERACLES means that it can have an increased fraction of non-detected spaxels near the edge of the disc, and we account for this in our analysis.

Combined Σ H -Σ and Σ mol -Σ relations
In the left panel of Focusing on HERACLES galaxies, the rH SM relation of H -normal systems occupies a ∼0.75 dex range in Σ H across the ∼2.5 dex in Σ .H is detected across the whole Σ range, with non-detected spaxels having preferentially lower Σ compared to detections, also shown in the Σ histograms.This distribution of HERACLES points reflects the known result that H surface densities typically saturate around ∼ 10 M pc -2 within the stellar components of galaxies (Martin & Kennicutt, 2001;Krumholz et al., 2009;Bigiel et al., 2008;Bigiel & Blitz, 2012).Focusing now on VERTICO, it is clear that the VERTICO combined rH SM relation inhabits a different region of the parameter space compared to HERACLES.First, there is a decreased number of H -detected spaxels across the range of Σ , but preferentially at lower Σ .This change is also visible in the Σ histograms where both the average VERTICO detected and non-detected spaxels are shifted to higher Σ compared to HERACLES.Second, the typical Σ H of VERTICO galaxies is lower.It is missing the largest Σ H values, there is a cloud of red points at log(Σ /M pc -2 ) > 1.5 with log(Σ H /M pc -2 ) ∼ 0 -0.5 that does not exist in HER-ACLES, and the log Σ H histogram of VERTICO galaxies is broader, extending to lower Σ H .The median detected spaxel has log(Σ H /M pc -2 ) = 0.58 for VERTICO galaxies, while HERACLES has log(Σ H /M pc -2 ) = 0.93, and a two-sample Kolmogorov-Smirnov (K-S) test rejects the null hypothesis that the two Σ H distributions are drawn from the same parent distribution with high significance (p < 10 -4 ).
Given that the VERTICO sample spans a range of def-H values (-0.3 -2.25 dex) and Σ increases toward the centres of galaxies, these changes are broadly consistent with galaxies undergoing ram-pressure stripping from the ICM.Namely, the ISM at large radii is affected first, while lowering the average gas surface density within the central regions where the gas is more gravitationally bound requires stronger stripping that has proceeded within the stellar disc (Roediger & Brüggen, 2007;Fumagalli et al., 2009;Steinhauser et al., 2016;Yoon et al., 2017), and these two effects clearly emerge in this figure.
How are these differences reflected in the molecular gas surface density of galaxies?In the right panel of Fig. 1, we show the resolved molecular gas-stellar mass (rMGSM) relation f for the combined HERACLES and VERTICO samples using the same colouring and non-detection labelling as the rH SM relation.Unlike Σ H , Σ mol shows a positive correlation with Σ , and differences between the combined rMGSM relation of HERACLES and VERTICO are less distinct.Nonf Often called the 'resolved molecular gas main sequence' in other works (Lin et al., 2019;Ellison et al., 2021b;Pessa et al., 2021).
Unweighted median Weighted median 16 th & 84 th percentile 0.0 0.5 PDF 2.1 1.9 Unweighted Weighted are shown with grey points in the background, the median relation in bins of log Σ for each galaxy is shown with thick, coloured lines, and the 16 th -84 th percentiles are shaded in the same colour, which matches the coloured name of each galaxy in the top le corner.Galaxies are coloured from red through to black increasing integrated molecular gas mass and alternate between panels for better visual separation.Third panel: The combined rH SM relation (top) and rMGSM (bottom) are shown with grey points in the background, overlaid are the medians (solid lines) and 16 th -84 th percentiles (dashed lines) treating all points equally (unweighted, blue) and each galaxy equally (weighted, red).The rightmost panels show the density-normalised histograms for the unweighted and weighted treatments.The key message of this Figure is that galaxies lie along well-defined, individual sequences within the rH SM and rMGSM relations.
As the location of an individual sequence is best described by global galaxy properties, we choose to weight each galaxy equally.HERACLES galaxies do not show a large di erence between the weighted and unweighted statistics, as they are restricted to the star-forming regime.
detected spaxels in both samples exist preferentially at lower Σ , although the typical Σ of VERTICO non-detected spaxels is larger than HERACLES.
The clearest difference, however, is that the VERTICO contours extend to lower Σ mol at all Σ , and this is true even at higher Σ where all Σ mol are above the detection limit.This difference is also visible in the VERTICO Σ mol histogram, which extends to lower Σ mol .Quantitatively, the median detected VERTICO spaxel has log(Σ mol /M pc -2 ) = 0.73, lower than the HERACLES value of log(Σ mol /M pc -2 ) = 0.85, and we can exclude that the distributions are drawn from the same parent sample (p < 10 -4 ).Some of these differences will be due to the shallower HERACLES detection limit (causing more non-detections at lower Σ ) and the smaller spatial coverage of VERTICO galaxies with smaller H discs (causing an undersampling of lower Σ ).However, they also suggest that the molecular gas is being lost preferentially at lower Σ (at least to our detection threshold), while the remaining, detectable molecular gas has lower surface density.However, the magnitude of these effects appears to be much smaller, and the affected region(s) are less well defined compared to the rH SM relation.
How, then, do we extract the comparatively weaker signal (if any) of the impact of environmental processes on the rMGSM relation?Using a sample of galaxies from the AL-MaQUEST survey (ALMA-MaNGA QUEnching and STar Formation, Lin et al., 2020) with similar M and SFR ranges to our VERTICO sample, Ellison et al. (2021a,b) showed that, to first order, the scatter within the rMGSM relation is not determined by the properties of individual spaxels.Instead, galaxies form well-defined individual sequences in the parameter space, and their scatter from the median relation is best described with global galaxy properties (e.g., the location of an individual rMGSM is best described by specific SFR, Ellison et al., 2021a).We show these individual sequences for both the rH SM and rMGSM relations of our HERACLES control sample in the left panel of Fig. 2, ordered from red to black by increasing total molecular gas mass.The distribution of HERACLES points is shown in the background, with the median relation (thick coloured lines) and 16 th -84 th percentiles (coloured, shaded regions) for each galaxy in bins of log Σ overlaid.Clearly, the galaxy-galaxy variation is larger than the scatter within individual objects, and this is particularly evident in the higher Σ regions of the rH SM relation, where galaxies might have a central H depression.
As the VERTICO sample contains both H -normal and H -deficient galaxies, combining the entire sample into one relation acts to wash out the signature of environmental effects, particularly in the molecular gas.Thus, to understand how environmental processes impact the resolved properties of galaxies, we need to control for their global properties.Further, we must also ensure that any histograms and statistics we compute are not biased toward larger galaxies and galaxies that are less affected by the environment, as these systems contain more spaxels.To do this, whenever we compute a histogram such as log Σ H or log Σ mol we weight each spaxel by 1/(N sp N G ), where N sp is the number of spaxels contributed by it's host galaxy and N G the number of contributing galaxies.Under this definition, the sum of the weights is always 1, and it ensures that each galaxy contributes to the histogram (and any subsequently computed statistics) equally.All statistics presented in this paper from now on are weighted values.4) and log Σ mol (Fig. 7) for the control sample, and the moderate and large def-H samples, in bins of log Σ .The Σ range used to select spaxels are given in the first column, the second to fourth columns use all spaxels, while the remaining columns use only detected spaxels for each gas phase.Uncertainties are calculated with bootstrap re-sampling using 10 4 iterations where entire galaxies are selected in each subsample to maintain equal weighting.In the middle panels of Fig. 2 we compare the weighted (red) and unweighted (blue) median (solid line), and 16 th and 84 th percentiles (dashed lines) g , of the rH SM relation and rMGSM relation for the HERACLES control sample.The rightmost panels show the weighted and unweighted histograms of log Σ H and log Σ mol across all Σ .The advantage of this weighting is clear in Σ mol at log(Σ /M pc -2 ) 1.5, where the HERACLES data begin to reach their detection threshold.While the unweighted statistics become nondetections, the weighted statistics better trace the relation as the typical galaxy is still detected.However, to avoid biasing our statistics, we will restrict some of our quantitative Σ mol comparisons to log(Σ /M pc -2 ) ≥ 1.25, as this is the typical Σ where the HERACLES galaxies transition to being dominated by non-detections.At these higher Σ , and at all Σ for Σ H , there are negligible differences between the two weighting methods.This is expected as our control sample is defined to be star-forming systems, which exhibit less galaxyto-galaxy variation in their spatially resolved cold gas scaling relations (Ellison et al., 2021a).As we see below, this variation is much larger for VERTICO galaxies.

Impact of environment on the spatially resolved cold gas properties of galaxies
Environmental effects on the H content of galaxies are a wellstudied topic, though many previous works have focused on the gas at large radii due to the extended nature of H .In §4.1, we present an analysis of the H content in the central regions of galaxies using the framework of the spatially resolved scaling relations.Not only does this present a new view of how environmental processes impact the resolved H content of g For N data points x i with weights w i , the weighted percentile p is the data point x j for which the cumulative sum of the weights of the ascending-sorted x i is equal to p. i.e., j i w i ≤ p and N j w i ≤ (1p).To compute the median, p = 0.5.
galaxies, but it is also useful to outline the visualisations and techniques that we apply to the molecular gas in §4.2.

Atomic hydrogen
In Fig. 3, we present the rH SM relations of VERTICO galaxies in individual panels.H -detected spaxels are shown with coloured points and non-detections with downward triangles.The median log Σ H in 0.25 dex bins of Σ (which are combined with adjacent bins if they contain less than 10 points) is shown with white, bordered symbols, where a circle denotes a median detection and a downward triangle denotes that the median is a non-detection.The median (thick, grey line) and 16 th and 84 th percentiles (thin, grey lines) of the HERACLES control sample are shown in the background, and the average log Σ measured at R H ( log Σ H , as calculated in §2.1) is shown with a vertical dashed line.Panels are ordered in rows from left to right by increasing def-H (the numerical value is displayed in the top left of each panel, in-line with the galaxy name), which also determines the colour of the individual points (with a colour bar on the right).Fig. 3 is rich in physical information, and we dissect this in the following analysis.
We first focus on the top row, which contains VERTICO galaxies considered to be H -normal based on their def-H (-0.3 < def-H [dex] < 0.3).Three galaxies are consistent with the HERACLES control sample, while NGC 4567 appears to be a slight outlier, particularly at lower Σ .NGC 4567 is interacting with NGC 4568 (row 2, column 2) and the galaxies cannot be separated in velocity space (Chung et al., 2009), causing an overlap of the rH SM relation of each galaxy and meaning that it has experienced some environmental mechanisms.
The rH SM relations of individual galaxies show two clear trends as a function of increasing def-H that occur distinctly at lower and higher Σ .First, the value of log Σ H (measuring R H , vertical dashed lines), on average, increases with increasing def-H .In other words, the H is stripped (at least to our sensitivity limit) to smaller radii within the disc as def-H increases.Second, inside R H , the individual rH SM relations increasingly fall below the median and then the 16 th (lower) percentile of the HERACLES control sample at fixed Σ .That is, there is an increasing reduction of the typical Σ H within the R H of galaxies as def-H increases.This effect is largest in the highest def-H systems, but several moderate def-H galaxies also begin to show evidence for this reduction (NGC 4321, NGC 4501, NGC 4535).These trends show that the impact of environmental processes on the H is two-fold: not only is the H stripped from the outskirts, but the typical surface density of the remaining gas is also reduced.Additionally, we observe no evidence for an environment-driven enhancement of Σ H , at least at the 1.2-kpc spatial scale that we trace.
To quantify these effects, and because this is a complex parameter space, we use two different strategies for computing and comparing log Σ H histograms of VERTICO galaxies to the control sample.In all comparisons, we split the VERTICO galaxies into two bins of global def-H h , def-H = [0.3,0.7) dex and def-H ≥ 0.7 dex, which we refer to these as the 'moderate' and 'large' def-H samples, respectively.These thresholds were defined as approximately 1 and 2 times the scatter in the def-H of isolated galaxies (σ=0.3 dex), with a slight adjustment to ensure an equal number of galaxies in each sample.Histograms are computed by weighting each galaxy equally (see §3), and we also computed weighted medians for each def-H sample and the control sample, which are presented in Table 2.
In our first comparison, shown in Fig. 4, we divide the rH SM relation into 0.5 dex bins of log Σ and compute the histogram of log Σ H in each bin.The control sample is shown as a grey, solid histogram, the moderate def-H population as a blue, open histogram, and the large def-H as a red, open histogram.All H non-detected spaxels are included in the bin located at log(Σ H /M pc -2 ) = -0.5, and if the value of this bin exceeds the figure limits, its value is written in the top left.
Fig. 4 again highlights the trends seen in Fig. 3. First, both def-H samples show a larger fraction of non-detected H spaxels than the control sample, and the fraction of non-detected spaxels is higher in the larger def-H galaxies.Namely, the radius at which H is below the detection limit gets closer to the galaxy centre as def-H increases.This increased fraction of non-detections is also evident in Table 2, where the weighted median for detections only is often different to the value for all points in both def-H samples.Second, at all Σ , the typical detected Σ H is also lower as def-H increases.Table 2 shows that the moderate def-H sample has a weighted median for detected Σ H spaxels that is 0.34 -0.83 dex lower than the control sample, while it is 0.53 -0.85 dex lower for the large def-H sample, and these differences are greater than the estimated uncertainties.Last, neither def-H sample shows higher h As mentioned in §2.1, our def-H bins are wide enough to fully encapsulate any small, second-order correlation of Σ H and def-H due to the variation of Σ H profiles in galaxies, guaranteeing that our analysis is tracing environmentally-related processes.Non-detected spaxels are included in the bin centered at log(Σ H /M pc -2 ) = -0.5, and if the value of this bin exceeds the y-axis limit then its value is quoted in the top-le of the panel in the same colour as its histogram.At fixed Σ , galaxies with larger def-H have a higher fraction of non-detected spaxels and their detected spaxels have lower Σ H .The di erence is larger at lower Σ .
Σ H than the control sample, highlighting that we observe no environment-driven enhancement of Σ H .We note that there is some overlap between the 1 σ scatter in the Σ H distributions.For example, in the highest Σ bin the 16 th percentile of the control sample (0.42) is smaller than the 84 th percentiles of both the moderate (0.73) and large (0.49) def-H samples.As we do not expect environmental mechanisms to operate across the whole gas disc with the same magnitude (e.g., the centre or the trailing edge will be less affected), some overlap between the distributions is expected.Nevertheless, in all Σ bins we find that the two distributions cannot be drawn from the same sample i (p < 10 -4 ), confirming that we are witnessing the effect of environmental processes already in the moderate def-H regime.
Considering the K-S test results, Fig. 4 also contains interesting insights into how environmental processes change the distribution of Σ H in the remaining gas.While the bins with log(Σ /M pc -2 ) < 1.25 are dominated by non-detections, in the higher log Σ bins, the log Σ H histogram of the moderate def-H galaxies (blue) is similar to the control sample.Namely, there is a peak at higher Σ H and a tail toward lower surface densities.The main difference is a decrease in the average Σ H , while the increase in non-detections comes from the low surface density tail falling below the detection limit.The large def-H galaxies show a larger decrease in Σ H values, but also, the shape of the log Σ H histogram becomes more gaussian rather than skewed.This change in shape suggests that in addition to the decrease in the typical surface density, the densest H regions are no longer present.
In our second comparison, we consider only H within the R H of each galaxy.For both def-H samples, we used all spaxels with log Σ ≥ log Σ H to compute equally-weighted density-normalised log Σ H histograms. Achieving a fair comparison with the control sample required an extra step, however, as the log Σ H of each galaxy is different.For each VERTICO galaxy, we computed a unique, equally-weighted control sample histogram using all spaxels from each control galaxy with log Σ ≥ log Σ H .The control histograms for each galaxy in a given def-H sample were then stacked and re-normalised using the same procedure j .This dependence on log Σ H results in slightly different control samples for each def-H sample k and is visible when comparing the grey histograms in the two panels of Fig. 5. Weighted medians for these histograms are listed in Table 3.
In Fig. 5, we compare histograms of log Σ H for spaxels within R H between the control sample (grey, solid), the moderate def-H sample (top panel, blue) and large def-H sample i As we are using weighted data, we used a weighted K-S test, and our implementation can be found here.
j This was done by introducing an additional factor of 1/D to the weight for each spaxel in each unique control sample histogram, where D is the number of galaxies in the def-H sample.
k We tested a different approach of combining all spaxels from the control galaxies and found the same results for the H using this method.This is because Σ H is essentially flat with Σ , so there is little dependence on the Σ range covered.The same is not true for the molecular gas, so for consistency, we use the same approach for both gas phases.(bottom panel, red).Here, we can see the impact of environmental processes on the remaining H gas that has not been removed.The moderate def-H galaxies have median detected log(Σ H /M pc -2 ) = 0.52 ± 0.07, while their matched control sample has log(Σ H /M pc -2 ) = 0.95 ± 0.02, indicating that the typical Σ H is 0.43 dex lower (Table 3).The difference is larger in the large def-H galaxies, which have a median log(Σ H /M pc -2 ) = 0.28±0.06that is 0.66 dex lower than the median of their matched control sample, log(Σ H /M pc -2 ) = 0.94 ± 0.02.We also see the same changes in the shape of the Σ H histogram described above in Fig. 4, and we can rule out the control and def-H samples being drawn from the same parent distribution (p < 10 -4 ).Thus, when considering the remaining gas inside R H , environmental processes typically lower the average Σ H of moderate def-H galaxies, while large def-H galaxies can no longer support the densest H gas.
Our results so far have provided a unique view of how environmental mechanisms impact the resolved H content of galaxies, and agree with results from previous studies (e.g.Vollmer et al., 2001;Tonnesen & Bryan, 2009;Chung et al., 2009;Yoon et al., 2017;Cortese et al., 2021;Boselli et al., 2022;Lee et al., 2022).However, how these mechanisms impact Table 3. Weighted median log Σ H (Fig. 5) and log Σ mol (Fig. 8) within R H for the moderate and large def-H samples, and their matched control samples.As the control sample histogram is constructed by selecting spaxels with log Σ ≥ log Σ H in each galaxy, the two def-H samples have unique control samples, and thus unique weighted control sample medians.The first column gives the def-H sample, the second and third columns the control sample and def-H sample medians using all spaxels, and the remaining columns use detected spaxels only.Uncertainties are calculated using jackknife re-sampling and rejecting an entire galaxy in each iteration.the molecular gas in galaxies on these spatial scales has so far remained unquantified, and is discussed next.

Molecular hydrogen
In Fig. 6, we present the individual rMGSM relations of VER-TICO galaxies in the same way as the rH SM relations.Focusing first on the top row of H -normal galaxies, aside from NGC 4567 (interacting, as discussed in §4.1), NGC 4651 seems to be molecule-poor for an H -normal galaxy.The kinematic position angles of the H inside and outside the optical disc of NGC 4651 are different (Chung et al., 2009), and it has stellar and H tidal features (Chung et al., 2009;Martínez-Delgado et al., 2010;Morales et al., 2018), suggesting that much of its H disc could have been recently accreted.Thus, the molecular gas content likely does not reflect that of an H -normal galaxy.NGC 4254 appears slightly molecule-rich at fixed Σ .
There is evidence that it recently underwent a fly-by interaction (Vollmer et al., 2005;Haynes et al., 2007;Boselli et al., 2018b), and it has a molecular bar and asymmetric spiral arms containing dense molecular gas (Sofue et al., 2003).These differences between the control sample and NGC 4651 and NGC 4254 highlight how environmental mechanisms can influence galaxies before they become significantly H -poor, and the importance of selecting a control sample from field galaxies.Before dissecting Fig. 6, it is useful to point out some interesting features.First, the median rMGSM relation of four galaxies (NGC 4298, NGC 4689, NGC 4457, NGC 4450) transition from detections to non-detections coincidently with R H , log Σ H , indicating that these galaxies show co-spatial H and molecular gas truncation.Second, several galaxies show changes in the slope of their rMGSM relation.Some of these changes are close to R H and are likely the signature of the environment (e.g.NGC 4298, NGC 4689).Other galaxies show changes in their more central regions, such as the superposition of two slopes in NGC 4535, a flattening and then upturn in NGC 4548, or the zig-zag structure of NGC 4579.These are not solely due to changes in the Σ mol content but also reflect changes in the central stellar structure of galaxies.NGC 4535 and NGC 4548 are barred with central bulges, causing a steeper Σ -Σ mol relation in their centres and bar regions and a flatter relation in the disc.NGC 4579 has a bar and a ring, causing an alternation between elevated Σ mol and Σ in the centre, lowered Σ mol in a higher Σ region between the centre and the ring, and higher Σ mol in a lower Σ region in the ring and disc.Beyond having interesting effects on the shape of the rMGSM relation, this suggests that we must use caution when comparing Σ mol properties at fixed Σ , particularly in the ranges where stellar structures could be important.
Comparing the rMGSM relations of VERTICO galaxies to the control sample, there is less difference as a function of def-H than was visible for the rH SM relation.In the second and third rows, namely at moderate def-H , the rMGSM relations of individual galaxies typically agree well with the median of the control sample, though there are some small signatures of deviation near R H in some galaxies.Instead, it is primarily the largest def-H systems (the bottom two rows) where the individual rMGSM relations begin to deviate below the control sample median and occasionally below the 16 th percentile, at fixed Σ .These differences are largest at lower Σ mol and closer to R H , though there is galaxy-galaxy variation and at high Σ some galaxies are consistent with the control sample (e.g.NGC 4457, NGC 4569), while others show reduced Σ mol (e.g.NGC 4380, NGC 4450).
These results imply that the Σ mol at lower Σ (outer parts of the disc) is more strongly impacted than at higher Σ (the inner regions).Further, as Σ mol is typically lower at lower Σ , the environmental processes preferentially impact the lower surface density molecular gas within galaxies.Interestingly, this can begin even at moderate def-H (e.g.NGC 4298 and its Σ mol slope-change and Σ mol truncation) and effectively changes the slope of the individual rMGSM relations, making them steeper.These results agree with previous studies, which found steeper radial profiles of Σ mol in galaxies with more disturbed H morphology and larger def-H (Fumagalli et al., 2009;Zabel et al., 2022;Villanueva et al., 2022).Additionally, like Σ H we observe no environment-driven elevation of Σ mol at the spatial scale that we are sensitive to.
In Fig. 7 we compare the log Σ mol histograms of VER-TICO galaxies to the control sample in 0.5 dex log Σ bins.We do not expect stellar structures, as mentioned above, to cause any biases in these histograms as they are primarily contained in the log(Σ /M pc -2 ) ≥ 2.25 bin.These histograms reflect the smaller offsets between the individual VERTICO rMGSM relations and the control sample.Moderate def-H galaxies (blue) show little differences from the control sample, except in the lowest Σ bin that has a significant contribution from non-detections.In this bin, the relatively larger number of non-detected spaxels in the control sample is due to the shallower detection limit of HERACLES.Overall, the main difference between the control sample and moderate def-H sample is the presence of a slightly more extended tail toward lower Σ mol .This tail imprints on the median detected Σ mol being up to ∼0.18 dex lower than the control sample, but the distributions peak in similar locations and in the two higher Σ bins the medians are consistent within their estimated uncertainties (Table 2).The presence of this tail, however, is sufficient to rule out the control sample and moderate def-H sample being drawn from the same distribution in all Σ bins (p < 4 × 10 -4 ).
In the large def-H sample (red), the impact of environmental processes is stronger.Non-detections dominate the smallest Σ bin, and despite the VERTICO data being deeper, the relative number of non-detected spaxels is larger than the control sample in this bin.Considering the detected spaxels, the weighted median log Σ mol is 0.34 -0.59 dex lower than the control sample across all log Σ bins, all of which are significant considering their uncertainties (Table 2), and we can exclude that the distributions are drawn from the same parent sample (p < 10 -4 ).
These results support the interpretation that it is primarily the lower surface density molecular gas at lower Σ that is affected by the environment, although as galaxies become significantly H -deficient ( 0.7), the impact is larger and can also affect the gas at higher Σ .Further, in the histograms of Fig. 7 we see no evidence for any environment-driven elevation of Σ mol , although in the largest Σ bin the moderate def-H sample appears to show slightly higher Σ mol .However, the medians are consistent within the uncertainties, and this enhancement disappears if we instead consider spaxels with log(Σ /M pc -2 ) ≥ 2 or log(Σ /M pc -2 ) = [2, 2.5), meaning that it is due to the slightly different Σ coverage of galaxies in this bin.
In Fig. 8, we compare the histograms of log Σ mol within R H for moderate def-H galaxies (top, blue) and large def- .At fixed Σ , moderate def-H galaxies have a similar range of detected Σ mol to HERACLES galaxies, with a slight extension to lower Σ mol , while large def-H galaxies show a lower typical Σ mol .
H galaxies (bottom, red) to their respective control sample log Σ mol histograms l (grey, solid).The Σ mol histogram of moderate def-H galaxies is similar to the control sample, and the median detected Σ mol is 0.15 dex lower (a factor of ∼1.4,Table 3), which is significant considering the uncertainty.Much like Fig. 7, the main difference is the tail toward lower Σ mol values, and there is also a reduction of spaxels with surface densities near the peak of the histogram.The control sample has slightly more non-detected spaxels due to the shallower detection limit, and we verified that the result remains if consider only spaxels with log(Σ /M pc -2 ) ≥ 1.5.We also note that the tail is primarily driven by NGC 4394, which sits below the control sample rMGSM relation at all Σ .We recomputed the log Σ mol histogram with NGC 4394 excluded and found that a smaller tail toward lower surface density Σ mol l As the molecular gas observations of some galaxies do not extend to R H (e.g., NGC 4321, NGC 4535), we picked the minimum Σ of each VERTICO galaxy (used to select control sample spaxels) to be the largest out of log Σ H , log(Σ /M pc -2 ) = 1.25 (see §3), and the minimum log Σ for the galaxy.In this way, only the observed Σ range is used to compute the weighted control sample. .Histograms of log Σ mol within the R H of galaxies, presented in the same way as Fig. 5. Inside the R H of moderate def-H galaxies, environmental processes a ect primarily the molecular gas with low surface density, while in large def-H galaxies even the typical Σ mol is reduced.
still remains.
The log Σ mol histogram of the large def-H galaxies shows a different behaviour.Unlike the moderate def-H galaxies, there is no extended tail of lower Σ mol , but instead, it has a similar shape to the control sample, only shifted to a lower average Σ mol .The weighted median detected spaxel is 0.45 dex lower (a factor of 2.8 or 64 per cent) than the control sample and significant considering the uncertainties (Table 3), reflecting a larger change, which would be expected in larger def-H galaxies.Thus, in the largest def-H galaxies, environmental processes act to lower the typical Σ mol within R H .These results support the interpretation that at moderate def-H , the lower surface density molecular gas within R H is impacted, while at large def-H the typical Σ mol is lowered.
To summarise, we have presented clear evidence that environmental processes impact the molecular gas content of galaxies on 1.2-kpc physical scales, even within the H truncation radius.The effect is primarily on the lower surface density molecular gas and can begin even at moderate def-H , while in the largest def-H systems, the denser molecular gas shows a reduction in the average molecular gas surface density.

Physical drivers of Σ mol reduction
The observed reduction of Σ mol at fixed Σ suggests that cluster environmental processes are impacting the molecular gas in VERTICO galaxies.This reduction could be due to direct stripping of the molecular gas, similar to the H (e.g., Moretti et al., 2018;Cramer et al., 2020).Alternately, it could result from inefficient molecule formation caused by a change in the ISM physical conditions due to the stripping of the H (e.g., Blitz & Rosolowsky, 2006;Krumholz et al., 2009;Fumagalli et al., 2009).In this case, the suppressed molecular gas content is due to consumption by star formation and the lack of cold ISM replenishment (akin to strangulation/starvation). However, for this last scenario to work requires the assumption that the ISM in a galaxy is able to respond to any changes on a timescale equal to, or shorter than, the stripping timescale of the H ( 200 Myr; Vollmer et al., 2001;Roediger & Brüggen, 2007;Tonnesen & Bryan, 2009;Choi et al., 2022).Indeed, this might be expected due to the short lifetimes of GMCs ( 30 Myr, Chevance et al., 2020b).In this case, VERTICO galaxies should follow the scaling relations used to interpret the resolved molecular-to-atomic gas mass ratio (Σ mol /Σ H ) in normal galaxies, which is sensitive to the physical conditions within the ISM.Instead, we show here that these scaling relations break down when applied to VERTICO galaxies, and they cannot describe their Σ mol /Σ H , suggesting that changes to the ISM physical conditions alone cannot explain the observed reduction of Σ mol .
More detailed modelling of the physical state of the cold gas reservoir, and its relationship to Σ mol /Σ H , is outside the scope of this work and will be investigated in future VER-TICO papers.For an analysis of Σ mol /Σ H across the discs of VERTICO galaxies using radial profiles, we refer the reader to Villanueva et al. (2022), who found elevated ratios in galaxies with perturbed and truncated H morphologies.Here, we compare Σ mol /Σ H in individual spaxels to two parameterisations of the ISM physical conditions to try and understand what physics might be responsible for the reduction in molecular gas content of the VERTICO galaxies.To be consistent with previous works, in this section, we correct our H masses for a 36 per cent contribution from Helium.Our analysis is restricted to detected spaxels only, and we also only consider spaxels with log Σ ≥ 1.25 so that we only study the ISM parameterisations where the control sample is dominated by detected spaxels.We also denote the total cold gas surface density as It is also useful to quantify how the removal of H alone manifests in the scaling relations that we will investigate.To do this, we computed the offset between the median log Σ H of each VERTICO galaxy and the median of the control sample at fixed Σ .We then subtracted this offset (as a function of Σ ) from the spaxels of each HERACLES galaxy to create mock 'H -stripped' galaxies.Using these mocks, we re-computed weighted histograms and statistics to estimate how the HER-ACLES control sample would shift if it had experienced H reduction equal to what is observed for each VERTICO galaxy.
We use these mock galaxies in the following analysis to show where deviations from the ISM scaling relations can be explained by H -only removal, and where we also require extra reduction of the molecular gas content.Indeed, it is our ability to perform this kind of spatially resolved analysis that allows us to identify these small differences, which are on the order of 0.1-0.2dex, that would likely be washed out in studies considering only integrated quantities such as global H and molecular gas deficiencies.For conciseness, we only show individual scaling relations for a subset of VERTICO galaxies in this section.Individual relations for all VERTICO galaxies, and the combined VERTICO and HERACLES relations are shown in Appendix 1.
First, we use a proxy for the hydrostatic pressure at the midplane of the disc, which was predicted to determine the molecular gas fraction by Elmegreen (1993, but see also Blitz &Rosolowsky 2006 andLeroy et al. 2008).We adopt the prescription from Blitz & Rosolowsky (2006), which we estimate as Σ 0.5 Σ gas .Σ and Σ gas are the dominant terms in the Blitz & Rosolowsky (2006) prescription, with the stellar scale height and gas velocity dispersion having minor contributions.If this prescription is an accurate predictor of physics determining Σ mol /Σ H , then we might expect that stripping of the ISM (predominantly H ) would reduce the pressure and cause Σ mol /Σ H to decrease correspondingly.
In Fig. 9, we show Σ mol /Σ H as a function of Σ 0.5 Σ gas for two moderate def-H galaxies (left two panels) and two large def-H galaxies (right two panels), presented similarly to Fig. 3 and Fig. 6.We also show, with red lines, how the HERACLES control sample median would shift if it had experienced H reduction equal to what is observed for the displayed VERTICO galaxy.All four VERTICO galaxies have larger Σ mol /Σ H than expected at fixed Σ 0.5 Σ gas showing that, despite a reduction in pressure due to the stripping of the H , the Σ mol /Σ H is larger.Thus, internal hydrostatic pressure cannot explain the Σ mol /Σ H ratio in these VERTICO galaxies, suggesting that stripping of the gas (at least of the H ) is more dominant than the hydrostatic pressure, at least over the timescale traced by VERTICO galaxies.Further, there are regions within Fig. 9 where the mock H -stripped control sample over-estimates the Σ mol /Σ H of VERTICO galaxies.This over-estimation is typically larger at lower pressures (corresponding to larger radii) and in the larger def-H galaxies, which are the regions and galaxies where we observed molecular gas reduction in §4.2.Thus, changes in the hydrostatic pressure plus H removal alone cannot explain the observed molecular gas reduction in these galaxies.
In Fig. 10, we show that the story told by these four galaxies holds for our whole VERTICO sample, using normalised histograms of log(Σ mol /Σ H ) in two bins of log Σ 0.5 Σ gas .The blue and red empty histograms correspond to the moderate and large def-H samples, and the solid grey histogram shows the Σ mol /Σ H histogram for the control sample.The solid blue and red histograms show the mock H -stripped control samples, matched to the moderate and large def-H samples respectively.In both bins, the VERTICO galaxies have elevated Σ mol /Σ H compared to the HERACLES control sample; but not so elevated that they match with the histograms of the mock H -stripped control samples.
Quantitatively, in the lowest Σ 0.5 Σ gas bin, an H -only removal scenario would result in Σ mol /Σ H distributions with medians of 0.24 ± 0.03 dex and 0.65 ± 0.03 dex for the moderate and large def-H populations, which over-estimate their .Ratio of molecular-to-atomic gas mass of individual spaxels as a function of the midplane pressure proxy Σ 0.5 Σgas for four VERTICO galaxies, coloured and ordered by def-H .Galaxies are presented in the same format as in Fig. 3 and Fig. 6, with the le two belonging to the moderate def-H sample and the right two to the large def-H sample.The red lines show how the median and 16 th and 84 th percentiles of the HERACLES control sample would shi if it had experienced the same magnitude of H reduction at fixed Σ as each VERTICO galaxy.At fixed midplane pressure, VERTICO galaxies have a larger Σ mol /Σ H than expected for the same pressure in a control sample galaxy.However, this di erence is not as large as the H -only removal scenario, particularly at lower pressures and larger def-H .observed values of 0.09 ± 0.07 dex and 0.31 ± 0.04 dex by 0.15 dex and 0.34 dex, respectively.In the highest Σ 0.5 Σ gas bin, the large def-H galaxies similarly show 0.2 dex smaller Σ mol /Σ H elevation (median 0.61 ± 0.04 dex) than expected from H removal alone (0.81 ± 0.04 dex).A K-S test between the def-H samples and their respective mock control sam-ples rules out that they are drawn from the same distribution (p < 10 -4 ).Conversely, the moderate def-H sample in this bin matches its mock control sample, which is reflected in the expected and observed Σ mol /Σ H medians of 0.47 ± 0.03 dex and 0.48 ± 0.04 dex, respectively.This match between the mock control sample and the moderate def-H sample agrees with our results in §4.2, where no reduction in Σ mol was observed in these galaxies.Further, we can rule out that the elevated Σ mol /Σ H is due to ram pressure providing additional pressure support and increasing the conversion rate of H to molecular gas, as we found no evidence for elevated molecular gas content in any of our VERTICO galaxies (at least at our 1.2-kpc spatial scale).All these points indicate that some reduction in Σ mol is required to explain the Σ mol /Σ H of VERTICO galaxies, but this parameterisation with hydrostatic pressure alone does not provide a physical explanation.
An alternate scenario based on ISM physics is that, due to the removal of H , there is insufficient fuel to form new molecular gas.We test this hypothesis using the total gas surface density (Σ gas ), suggested by Krumholz et al. (2009, see also Fumagalli et al. 2009) to be the primary regulator of Σ mol /Σ H with metallicity and the interstellar radiation field playing secondary roles.If the stripping of the ISM lowers the total surface density such that the formation of molecular gas is inefficient, we would expect VERTICO galaxies to lie along the same relation between Σ mol /Σ H and Σ gas as the control sample, as claimed in Fumagalli et al. (2009).In Fig. 11, we show Σ mol /Σ H as a function of Σ gas in the same style as shown in Fig. 9 for the pressure relation.The displayed moderate (left two panels) and large (right two panels) def-H galaxies do not lie along the same relation as the control sample.Further, the removal of H alone would over-estimate the observed Σ mol /Σ H at fixed Σ gas , particularly in lower Σ gas regions.This trend remains for the whole sample, and in Fig. 12, we show histograms of Σ mol /Σ H in two bins of log Σ gas in the same style as Fig. 10.In the lower Σ gas bin, the expected median Σ mol /Σ H from H -only removal scenarios, 0.17 ± 0.04 dex and 0.63 ±0.04 dex, over-estimate the observed values .Ratio of molecular-to-atomic gas mass of individual spaxels as a function of Σgas for the same galaxies as Fig. 9.The Σ mol /Σ H of VERTICO galaxies is systematically larger than in control sample galaxies at fixed Σgas, but not as large as the H -only removal scenario shown with the red lines.  of 0.05±0.07dex 0.47±0.03dex for the moderate and large def-H samples by 0.12 dex and 0.16 dex respectively.In the higher Σ gas bin, the elevation in Σ mol /Σ H of the large def-H sample is smaller by 0.11 dex, with an expected median 0.74±0.03dex compared to 0.63 ± 0.04 dex observed, while the moderate def-H galaxies agree more closely with their mock control sample, having 0.38 ± 0.03 dex expected and 0.41 ± 0.05 dex observed.We also can exclude that both def-H samples are drawn from the same parent distribution as their respective mock control samples (p < 10 -4 ).Thus, in VERTICO galaxies it is not Σ gas that sets Σ mol /Σ H , and explaining the observed Σ mol /Σ H requires a reduction in their molecular gas content that is not related to insufficient fuel or an inefficient conversion of H into molecular gas.We must also consider if there are any systematic differences between how molecular gas is traced in VERTICO and HERACLES.As the data have been homogenised and both use the same constant R 21 and α CO to convert CO(2-1) luminosity to molecular gas mass, this scenario is unlikely.Although R 21 and α CO have their own uncertainties, this homogenisation means that our results are actually re-scaled differences between the observed CO(2-1) emission in our samples.Thus, to explain the difference between the typical Σ mol inside the R H of large def-H galaxies would require VERTICO galaxies to preferentially have CO(2-1)-dark gas such that the observed CO(2-1) luminosity under-estimates the actual amount of molecular gas present by a factor of 2.8 (0.45 dex).
Metallicity variations cannot explain a reduction in CO(2-1) emission, as while at smaller metallicities there is more CO(2-1) dark gas, cluster galaxies show similar or elevated metallicity over isolated systems (e.g., Skillman et al., 1996;Gupta et al., 2018;Franchetto et al., 2020).Changes to the UV radiation field due to reduction in star formation also are unlikely to cause any differences in CO(2-1) emission.VERTICO galaxies show reduced SFR surface densities up to ∼1 dex compared to star-forming galaxies at fixed Σ (Brown et al., in prep); however, in the same regions, we observe reduced H and molecular gas densities with similar magnitude, which reduces the gas column's ability to self-shield.Thus, any changes in the UV field are likely balanced by changes in the gas density (Krumholz et al., 2009).VERTICO galaxies also show similar R 21 values to normal galaxies from xCOLD GASS (Brown et al., 2021), suggesting that there is no significant difference in the CO(2-1) excitation that could under-estimate the amount of molecular gas.Additionally, Smith et al. (2014) and Glover & Smith (2016) found that CO(2-1)-dark gas is typically diffuse and low surface density, making it prone to turbulent compression and gravitational collapse and thus less likely to exist in the central regions of galaxies that we have studied.
Further evidence that VERTICO galaxies do not have substantial CO(2-1)-dark gas comes from studies of the Kennicutt-Schmidt (Schmidt, 1959;Kennicutt, 1998) star formation law.If the CO(2-1) emission line underestimated the total molecu-lar gas mass, then galaxies would deviate from the Kennicutt-Schmidt law and appear to have elevated SFR at fixed molecular gas mass.However, VERTICO galaxies show similar galaxyto-galaxy variation in their resolved Kennicutt-Schmidt relations to what is observed in star-forming galaxies (Vollmer et al., 2012;Jiménez-Donaire et al., 2022).This consistency with unperturbed systems implies that the physics of star formation, and by extension the ISM physical conditions giving rise to the CO(2-1) emission, do not show sufficient difference between VERTICO galaxies and the control sample to explain our results.Thus, we need to invoke a process that is unrelated to the ISM physical conditions and can reduce molecular gas content on timescales similar to the H gas. Direct stripping of the molecular gas meets this requirement.
There is evidence that not all molecular gas in galaxies resides in bound GMCs, with 20 -50 per cent (Polk et al., 1988;Wilson & Walker, 1994;Pety et al., 2013;Roman-Duval et al., 2016;Chevance et al., 2020a;Saintonge & Catinella, 2022) existing in a more diffuse phase or a molecular thick disc that has velocity dispersion similar to the H (e.g.Caldú-Primo et al., 2013, 2015;Mogotsi et al., 2016), making it much easier to strip than GMC gas.Our data cannot differentiate between diffuse and dense gas; however, the 64 per cent reduction ( §4.2) in the typical molecular gas content within R H observed in our large def-H galaxies is greater than the 'typical' diffuse gas fractions mentioned above.This is interesting, as there is also evidence that GMCs can be disrupted by the ICM.Cramer et al. (2020) found that only 12 per cent of the molecular gas on the leading edge of NGC 4402 exists in GMCs that have the survived the ongoing stripping, while the remaining gas exists in a trailing plume with lower surface density.The authors suggest that this large fraction of lower surface density molecular gas implies that an indirect stripping process, by which GMCs are converted to lower surface density molecular gas or H before being directly stripped, could be taking place.This interpretation is also supported by simulations, which have shown that the multi-phase structure of the ISM allows stripping to affect gas within the truncation radius (e.g., Quilis et al., 2000;Tonnesen & Bryan, 2009).Further, momentumtransfer from the ICM to the ISM through turbulent mixing can effectively ablate the dense, cold ISM into warmer phases (e.g., Gronke & Oh, 2018;Sparre et al., 2019) and accelerate it out of the galaxy (Tonnesen & Bryan, 2021;Choi et al., 2022).The typical ICM pressures expected for Virgo galaxies are within the ranges used in these simulations (e.g., Vollmer et al., 2001;Roediger & Brüggen, 2007;Chung et al., 2007;Abramson et al., 2011), making stripping of the molecular gas the most feasible explanation for our observed molecular gas reduction.However, future work with higher spatial resolution observations is needed to allow quantification of how the diffuse and dense phases of the ISM are impacted by the Virgo environment, and search for signatures of direct stripping of the molecular gas.

Summary
In this paper, we used spatially resolved scaling relations between Σ , Σ H , and Σ mol to study how cluster environmental processes impact the cold gas content within the optical extent of massive (log(M /M ) ≥ 10) galaxies on 1.2-kpc scales.We presented a new view of how the H gas distribution changes with increasing H deficiency, and we have shown, statistically, how environmental processes impact the spatial distribution of the molecular gas, the distribution of Σ mol , and molecular-toatomic gas mass ratio using observations from the VERTICO survey (Brown et al., 2021).Our main results are as follows: • While the stripping of gas from galaxies impacts the H gas at large radii, it has the additional effect of lowering the average H surface density of the remaining gas within the H -stripping radius, R H .This reduction in surface density increases as the gas removal proceeds inward until the densest H regions can no longer form (Figs. 3, 4, 5).• Some galaxies show co-spatial truncation of their H and molecular gas (Fig. 6), indicating that environmental processes can effectively remove both phases on similar timescales in the same spatial region.• Overall, the molecular gas is less affected than the H , and we find evidence that environmental mechanisms first reduce the lower surface density molecular gas that exists preferentially at lower Σ (and thus typically larger radii) in galaxies (Fig. 7).However, we find that once the removal of H proceeds sufficiently far within the optical disc, the average surface density of the remaining molecular gas within R H is reduced (Fig. 8).• At the spatial scales we are sensitive to, we see no evidence for environment-driven elevation of the H or molecular gas surface density outside the ranges observed in starforming galaxies (Figs. 4, 7).• VERTICO galaxies have systematically larger molecularto-atomic gas mass ratios in individual regions compared to the control sample, but smaller values than expected from the removal of H alone (Figs. 9,11).Using proxies for the physical conditions in the ISM, we showed that changes to the physical state of the gas could not explain the observed molecular-to-atomic gas mass ratios (Figs. 10,12).Thus, ablation and/or stripping of the molecular gas on a similar timescale as the H in a given 1.2-kpc region is the best explanation for the reduction of molecular gas content in VERTICO galaxies.
Our results agree with, and extend upon, previous literature focused on gas stripping and environmental effects in cluster galaxies (Cortese et al., 2021;Boselli et al., 2022).Combined with the recent, complementary studies from the VER-TICO team (Brown et al., 2021;Villanueva et al., 2022;Zabel et al., 2022), a coherent evolutionary scenario for Virgo cluster galaxies begins to emerge.Compared to isolated systems, statistically, Virgo galaxies have normal molecular gas content (Boselli et al., 2014;Cortese et al., 2016) and elevated molecularto-atomic gas mass ratios (Mok et al., 2016;Villanueva et al., 2022) when considering their integrated properties only.Subsequently, there is poor correlation between H and molecular gas deficiency measurements (Zabel et al., 2022), but this does not mean that the molecular gas is not affected by environmental processes.Instead, major changes in the integrated molecular gas mass occur only once a significant fraction of the cold gas within the stellar disc has been affected.
Our analysis reveals how the molecular gas content of galaxies is affected by cluster environmental processes in a spatially-resolved way.We have shown that the reduction of molecular gas begins at the edge of the stellar disc even in moderately H -deficient galaxies, while in the highest Hdeficiency systems the molecular gas becomes significantly reduced across the disc.While this result has been shown using integrated measurements and radial profiles (Villanueva et al., 2022;Zabel et al., 2022), our approach allows us to quantify this effect as a function of Σ with exquisite detail, and has allowed us to disentangle between environmentally-driven and secular/internal processes.
Interestingly, previous works have found evidence for the compression of molecular gas due to the ICM-ISM interaction (e.g., Lee et al., 2017;Moretti et al., 2020a,b), particularly in the early stages of infall.We find no evidence for elevated H or molecular gas surface densities, suggesting that in our sample, which spans a larger (and later) range of infall stages and at the 1.2-kpc scales we are sensitive to, environment-driven elevation of cold gas content is not significant.
Our results also emphasise that the gas removal process is more complex than a simple outside-in stripping scenario, and it is clear that the environmental processes affect both the H and molecular gas within the H truncation radius (see also Lee et al., 2022), which is typically assumed to be the boundary of where stripping is efficient.Even at moderate H deficiency, we find that the average surface density of H remaining within the disc is reduced, indicating that the environment does not just affect the H gas in the outskirts.This result can be understood as the multi-phase ISM enabling the ICM to stream through low-density regions and affect the gas across the disc (e.g., Tonnesen & Bryan, 2009;Lee et al., 2020).However, we have also shown that the molecular gas is affected similarly, although it requires a stronger ICM-ISM interaction (as traced by more H lost). Thus, the ICM can ablate and strip gas across the entire discs of galaxies and does not discriminate which phase is affected so long as the interaction is strong enough (e.g., Tonnesen & Bryan, 2021;Choi et al., 2022).Indeed, this agrees well with observations of larger molecular gas masses in the cometary tails of stripped galaxies when the associated H α emission is brighter or the stripping is more advanced (Moretti et al., 2018).In these cases, it is likely that the stronger stripping has led to a larger mass-loading of the outflowing gas.Future work with larger samples of field and cluster galaxies, but in particular, higher resolution observations that can separate dense and diffuse CO(2-1) emission will be key to understanding the details of gas suppression in cluster galaxies.
Last, these results have interesting implications for the quenching of star formation in cluster galaxies.In particular, the reduction of molecular gas within the inner regions of galaxies will cause star formation to shut down earlier, and with a different distribution, than just radial truncation.Comparison to deep H α observations of VERTICO galaxies from the VESTIGE (A Virgo Environmental Survey Tracing Ionised Gas Emission, Boselli et al., 2018a) survey would enable the study of changes to recent, massive star formation in regions where molecular gas is reduced.Further, future observations with sensitive optical integral-field spectroscopic instruments such as MUSE (Multi Unit Spectroscopic Explorer, Bacon et al., 2017) will allow the study of the physics of galaxy quenching in unprecedented detail, including the additional impact of feedback from active galactic nuclei (e.g., Boselli et al., 2016b;George et al., 2019;Peluso et al., 2022).In particular, employing spectro-photometric fitting techniques (e.g., Fossati et al., 2018) will enable high spatial resolution studies of star formation and star-formation history reconstruction across the discs of Virgo cluster galaxies.These results can then be connected to cold gas reduction on sub-kpc scales where ISM physics operates.

dFigure 1 .
Figure1.Spatially resolved Σ Hand Σ mol -Σ scaling relations for the combined samples.The le panel shows log Σ H vs log Σ , the resolved Hstellar mass scaling relation, and the right panel shows log Σ H vs log Σ , the resolved molecular gas-stellar mass scaling relation.Contours for both samples enclose 5, 16, 50, 84, and 95 per cent of the H or molecular gas detected spaxels, while non-detected spaxels are shown with downward triangles at log(Σ H /M pc -2 ) = log(Σ mol /M pc -2 ) = -0.5.Density-normalised histograms of log Σ , log Σ H , and log Σ mol are shown opposite their respective axes, with filled histograms representing H -or molecular gas-detected spaxels while open histograms represent non-detections.The HERACLES Σ histogram for H non-detections uses 0.25 dex bins compared to the 0.1 dex used for VERTICO as it has fewer points (331 compared to 5717).Compared to HERACLES, the combined VERTICO sample has fewer Σ H detected spaxels, particularly at lower log Σ , and a lower average log Σ H . Similarly, but less clearly, the combined VERTICO sample has fewer Σ mol detections at lower Σ , while the histogram of detected Σ mol has a larger tail toward lower values.
Fig. 1 we plot Σ H against Σ , in other words, the resolved H -stellar mass (rH SM) relation.Detected spaxels for HERACLES (blue) and VERTICO (red) are shown with contours enclosing 5, 16, 50, 84, and 95 per cent of the population, and non-detections are shown at log(Σ H /M pc -2 ) = -0.5 with downward triangles.Densitynormalised histograms of log Σ and log Σ H are shown opposite to their respective axes, computed with 0.25 dex bins.Filled histograms show detected spaxels, while the open histograms on the Σ axis show the histogram of log Σ for H non-detected spaxels.

Figure 2 .
Figure 2.Resolved Σ Hand Σ mol -Σ relations for HERACLES control galaxies.Le two panels: The combined rH SM relation (top) and rMGSM (bottom) are shown with grey points in the background, the median relation in bins of log Σ for each galaxy is shown with thick, coloured lines, and the 16 th -84 th percentiles are shaded in the same colour, which matches the coloured name of each galaxy in the top le corner.Galaxies are coloured from red through to black increasing integrated molecular gas mass and alternate between panels for better visual separation.Third panel: The combined rH SM relation (top) and rMGSM (bottom) are shown with grey points in the background, overlaid are the medians (solid lines) and 16 th -84 th percentiles (dashed lines) treating all points equally (unweighted, blue) and each galaxy equally (weighted, red).The rightmost panels show the density-normalised histograms for the unweighted and weighted treatments.The key message of this Figure is that galaxies lie along well-defined, individual sequences within the rH SM and rMGSM relations.As the location of an individual sequence is best described by global galaxy properties, we choose to weight each galaxy equally.HERACLES galaxies do not show a large di erence between the weighted and unweighted statistics, as they are restricted to the star-forming regime.
Figure3.Individual rH SM relations for VERTICO galaxies, with the galaxy name given in the top-right of each panel.The panels are ordered in rows by def-H (noted in the top-le corner), which also sets the colour of each panel, and a colour bar is displayed on the right-hand side.Solid points indicate detected spaxels, while downward triangles at log(Σ H /M pc -2 ) = -0.5 indicate non-detections.The median rH SM for each galaxy computed in 0.25 dex Σ bins is shown with large white markers, where downward triangles indicate the median is a non-detection, and error bars show the 16 th and 84 th percentiles of the log Σ H distribution in each bin.The HERACLES control sample is shown with thick (median) and thin (16 th and 84 th percentiles), and the vertical dashed, coloured lines show log Σ H , the average log Σ at the H truncation radius (R H ). The two galaxies marked with stars are interacting, and their H observations are partially confused.The impact of the Virgo cluster on the rH SM relation is two-fold.As def-H increases there are more non-detected spaxels and log Σ H moves to larger values; and there is an increasing reduction of the typical Σ H within log Σ H .

Figure 4 .
Figure 4. log Σ H histograms in bins of log Σ .Each panel shows densitynormalised histograms of log Σ H for spaxels in the Σ range quoted in the top right corner.HERACLES control sample galaxies are shown with a grey, filled histogram, VERTICO moderate def-H galaxies with a blue, open histogram, and large def-H galaxies with a red, open histogram.Non-detected spaxels are included in the bin centered at log(Σ H /M pc -2 ) = -0.5, and if the value of this bin exceeds the y-axis limit then its value is quoted in the top-le of the panel in the same colour as its histogram.At fixed Σ , galaxies with larger def-H have a higher fraction of non-detected spaxels and their detected spaxels have lower Σ H .The di erence is larger at lower Σ .
Figure5.Histograms of log Σ H within the R H of galaxies.We compare density-normalised histograms of log Σ H for spaxels with log Σ ≥ log Σ H for galaxies with moderate def-H (top panel; blue, open histogram) and large def-H (bottom panel; red, open histogram) to their matched HERACLES control samples (both panels, grey filled histogram).As def-H increases, the fraction of non-detected spaxels within R H increases, and the typical detected Σ H is lower.In galaxies with the largest def-H the histogram of detected log Σ H is gaussian, rather than skewed, indicating that the densest H is no longer present.Neither def-H sample shows evidence for elevated Σ H compared to the control sample.

Figure 6 .
Figure6.Individual rMGSM relations for VERTICO galaxies, presented in the same way as Fig.3.As def-H increases, the clearest deviation of Σ mol from the HERACLES control sample occurs at lower Σ , typically close to or outside of R H ( log Σ H ).Only some of the galaxies with the largest def-H have reduced Σ mol at higher Σ and within their R H (e.g., NGC 4579 NGC 4380, NGC 4450).

Figure 7 .
Figure7.Histograms of log Σ mol in bins of Σ , presented in the same way as Fig.4.At fixed Σ , moderate def-H galaxies have a similar range of detected Σ mol to HERACLES galaxies, with a slight extension to lower Σ mol , while large def-H galaxies show a lower typical Σ mol .
Figure8.Histograms of log Σ mol within the R H of galaxies, presented in the same way as Fig.5.Inside the R H of moderate def-H galaxies, environmental processes a ect primarily the molecular gas with low surface density, while in large def-H galaxies even the typical Σ mol is reduced.
Figure9.Ratio of molecular-to-atomic gas mass of individual spaxels as a function of the midplane pressure proxy Σ 0.5 Σgas for four VERTICO galaxies, coloured and ordered by def-H .Galaxies are presented in the same format as in Fig.3and Fig.6, with the le two belonging to the moderate def-H sample and the right two to the large def-H sample.The red lines show how the median and 16 th and 84 th percentiles of the HERACLES control sample would shi if it had experienced the same magnitude of H reduction at fixed Σ as each VERTICO galaxy.At fixed midplane pressure, VERTICO galaxies have a larger Σ mol /Σ H than expected for the same pressure in a control sample galaxy.However, this di erence is not as large as the H -only removal scenario, particularly at lower pressures and larger def-H .

Figure 10 .
Figure10.Histograms of log(Σ mol /Σ H ) in bins of midplane pressure, presented in the same way as Fig.4and Fig.7.Additionally, solid blue and red histograms show how the log(Σ mol /Σ H ) histogram of the control sample would move if it had experienced H removal of equal magnitude as the moderate and large def-H samples, respectively.Moderate and large def-H galaxies have systematically larger Σ mol /Σ H than control sample galaxies at fixed Σ 0.5 Σgas.Further, the H -only removal scenario over-estimates the observed Σ mol /Σ H , except in moderate def-H galaxies in the largest Σ 0.5 Σgas bin.
Figure11.Ratio of molecular-to-atomic gas mass of individual spaxels as a function of Σgas for the same galaxies as Fig.9.The Σ mol /Σ H of VERTICO galaxies is systematically larger than in control sample galaxies at fixed Σgas, but not as large as the H -only removal scenario shown with the red lines.

Figure 13 .Figure 14 .
Figure13.Combined versions of the spatially resolved Σ mol /Σ H relations as a function of two proxies for the ISM physical conditions.The le panel shows Σ mol /Σ H vs the midplane-pressure proxy log Σ 0.5 Σgas, and the right panel vs the total gas surface density Σgas.Contours levels for both samples enclose the 5, 16, 50, 84, and 95 per cent of the sample.Density-normalised histograms of log(Σ mol /Σ H ), log Σ 0.5 Σgas, and log Σgas are shown opposite their respective axes.

Table 2 .
Weighted median log Σ H (Fig. 7. Additionally, solid blue and red histograms show how the log(Σ mol /Σ H ) histogram of the control sample would move if it had experienced H removal of equal magnitude as the moderate and large def-H samples, respectively.Moderate and large def-H galaxies have systematically larger Σ mol /Σ H than control sample galaxies at fixed Σ 0.5 Σgas.Further, the H -only removal scenario over-estimates the observed Σ mol /Σ H , except in moderate def-H galaxies in the largest Σ 0.5 Σgas bin.
Figure12.Histograms of log(Σ mol /Σ H ) in bins of total cold gas surface density, presented in the same way as Fig.10.Moderate and large def-H galaxies have systematically larger Σ mol /Σ H than expected if molecular gas formation was set primarily by Σgas.The H -only removal scenario overestimates Σ mol /Σ H , except at large Σgas in moderate def-H galaxies.