Distribution Of Galaxies Research Articles

Insights into the 21 cm field from the vanishing cross-power spectrum at the epoch of reionization

ABSTRACT The early stages of the epoch of reionization, probed by the 21 cm line, are sensitive to the detailed properties and formation histories of the first galaxies. We use 21cmfast and a simple, self-consistent galaxy model to examine the redshift evolution of the large-scale cross-power spectrum between the 21 cm field and line-emitting galaxies. A key transition in redshift occurs when the 21 cm field shifts from being positively correlated with the galaxy distribution to being negatively correlated. Importantly, this transition redshift is insensitive to the properties of the galaxy tracers but depends sensitively on the thermal and ionization histories traced through the 21 cm field. Specifically, we show that the transition occurs when both ionization fluctuations dominate over 21 cm spin temperature fluctuations and when the average spin temperature exceeds the temperature of the cosmic microwave background. We illustrate this with three different 21 cm models which have largely the same neutral fraction evolution but different heating histories. We find that the transition redshift has a scale dependence, and that this can help disentangle the relative importance of heating and ionization fluctuations. The best prospects for constraining the transition redshift occur in scenarios with late X-ray heating, where the transition occurs at redshifts as low as z ∼ 6–8. In our models, this requires high-redshift galaxy surveys with sensitivities of $\sim 10^{-18}~\rm erg\,s^{-1}\,cm^{-2}$ for optical lines and $\sim 10^{-19}~\rm erg\,s^{-1}\,cm^{-2}$ for far-infrared lines. Future measurements of the transition redshift can help discriminate between 21 cm models and will benefit from reduced systematics.

White Dwarfs in Dwarf Spheroidal Galaxies: A New Class of Compact-Dark-Matter Detectors.

Recent surveys have discovered a population of faint supernovae, known as Ca-rich gap transients, inferred to originate from explosive ignitions of white dwarfs. In addition to their unique spectra and luminosities, these supernovae have an unusual spatial distribution and are predominantly found at large distances from their presumed host galaxies. We show that the locations of Ca-rich gap transients are well matched to the distribution of dwarf spheroidal galaxies surrounding large galaxies, in a scenario where dark matter interactions induce thermonuclear explosions among low-mass white dwarfs that may be otherwise difficult to ignite with standard stellar or binary evolution mechanisms. A plausible candidate to explain the observed event rate are primordial black holes with masses above 10^{21} grams.

Towards quantifying the impact of triaxiality on optical signatures of galaxy clusters: weak lensing and galaxy distributions

ABSTRACT We present observational evidence of the impact of triaxiality on radial profiles that extend to 40 Mpc from galaxy cluster centres in optical measurements. We perform a stacked profile analysis from a sample of thousands of nearly relaxed galaxy clusters from public data releases of the Dark Energy Survey and the Dark Energy Camera Legacy Survey. Using the central galaxy elliptical orientation angle as a proxy for galaxy cluster orientation, we measure cluster weak lensing and excess galaxy density axis-aligned profiles, extracted along the central galaxy’s major or minor axes on the plane of the sky. Our measurements show a ≳ 2σ–3σ difference per radial bin between the normalized axis-aligned profiles. The profile difference between each axis-aligned profile and the azimuthally averaged profile ($\sim \pm 10\,\rm per\ cent-20~{{\ \rm per\ cent}}$ along major/minor axis) appears inside the clusters (∼0.4 Mpc) and extends to the large-scale structure regime (∼10–20 Mpc). The magnitude of the difference appears to be relatively insensitive to cluster richness and redshift, and extends further out in the weak lensing surface mass density than in the galaxy overdensity. Looking forward, this measurement can easily be applied to other observational or simulation data sets and can inform the systematics in cluster mass modelling related to triaxiality. We expect imminent upcoming wide-area deep surveys, such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, to improve our quantification of optical signatures of cluster triaxiality.

The DESI one-per cent survey: exploring the halo occupation distribution of luminous red galaxies and quasi-stellar objects with AbacusSummit

ABSTRACT We present the first comprehensive halo occupation distribution (HOD) analysis of the Dark Energy Spectroscopic Instrument (DESI) One-Percent Survey luminous red galaxy (LRG) and Quasi Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 h−1 Mpc in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of $f_\mathrm{sat} = 11\pm 1~{y{\ \mathrm{per\,cent}}}$, a mean halo mass of $\log _{10}\overline{M}_h/M_\odot =13.40^{+0.02}_{-0.02}$, and a linear bias of $b_\mathrm{lin} = 1.93_{-0.04}^{+0.06}$. For LRGs in 0.6 < z < 0.8, we find $f_\mathrm{sat}=14\pm 1~{{\ \mathrm{per\,cent}}}$, $\log _{10}\overline{M}_h/M_\odot =13.24^{+0.02}_{-0.02}$, and $b_\mathrm{lin}=2.08_{-0.03}^{+0.03}$. For QSOs, we infer $f_\mathrm{sat}=3^{+8}_{-2}\mathrm{per\,cent}$, $\log _{10}\overline{M}_h/M_\odot = 12.65^{+0.09}_{-0.04}$, and $b_\mathrm{lin} = 2.63_{-0.26}^{+0.37}$ in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high fidelity galaxy mocks, forming the basis of systematic tests for DESI Y1 cosmological analyses. We also study the redshift-evolution of the DESI LRG sample from z = 0.4 up to z = 1.1, revealling significant and interesting trends in mean halo mass, linear bias, and satellite fraction.

TESTING AN ENTROPY ESTIMATOR RELATED TO THE DYNAMICAL STATE OF GALAXY CLUSTERS

We propose the entropy estimator HZ, calculated from global dynamical parameters, in an attempt to capture the degree of evolution of galaxy systems. We assume that the observed (spatial and velocity) distributions of member galaxies in these systems evolve over time towards states of higher dynamical relaxation (higher entropy), becoming more random and homogeneous in virial equilibrium. Thus, the HZ-entropy should correspond to the gravitacional assembly state of the systems. This was probed in a sample of 70 well sampled clusters in the Local Universe whose gravitational assembly state, classified from optical and X-ray analysis of substructures, shows clear statistical correlation with HZ. This estimator was also tested on a sample of clusters (halos) from the IllustrisTNG simulations, obtaining results in agreement with the observational ones.

Evolution of cosmic filaments in the MTNG simulation

We present a study of the evolution of cosmic filaments across redshift with an emphasis on some important properties: filament lengths, growth rates, and radial profiles of galaxy densities. Following an observation-driven approach, we built cosmic filament catalogues at z = 0, 1, 2, 3, and 4 from the galaxy distributions of the large hydro-dynamical run of the MilleniumTNG project. We employed the extensively used DisPerSE cosmic web finder code, for which we provide a user-friendly guide, including the details of a physics-driven calibration procedure, with the hope of helping future users. We performed the first statistical measurements of the evolution of connectivity in a large-scale simulation, finding that the connectivity of cosmic nodes (defined as the number of filaments attached) globally decreases from early to late times. The study of cosmic filaments in proper coordinates reveals that filaments grow in length and radial extent, as expected from large-scale structures in an expanding Universe. But the most interesting results arise once the Hubble flow is factored out. We find remarkably stable comoving filament length functions and over-density profiles, showing only little evolution of the total population of filaments in the past ∼12.25 Gyr. However, by tracking the spatial evolution of individual structures, we demonstrate that filaments of different lengths actually follow different evolutionary paths. While short filaments preferentially contract, long filaments expand along their longitudinal direction with growth rates that are the highest in the early, matter-dominated Universe. Filament diversity at a fixed redshift is also shown by the different (∼5σ) density values between the shortest and longest filaments. Our results hint that cosmic filaments can be used as additional probes for dark energy, but further theoretical work is still needed.

Overabundance of orphan galaxies in the UniverseMachine empirical galaxy formation model

ABSTRACT Orphan galaxies that have lost a large fraction of the dark matter subhaloes have often been invoked in semi-analytical as well as empirical models of galaxy formation. However, observations have not yet constrained the distribution of orphan galaxies in such models so far. We run a mock cluster finder that mimics the optical cluster finding technique of the redMaPPer algorithm on a catalogue of galaxies with quenched star formation from one such empirical model, the UniverseMachine, and obtain the prevalence of orphan galaxies in these clusters as a function of their cluster-centric distance. We compare the fraction of orphan galaxies with the upper limits derived based on our prior observations of the weak lensing signals around satellite galaxies from SDSS redMaPPer clusters. Although the orphan fraction from the UniverseMachine is marginally consistent with the upper limits in the innermost regions of galaxy clusters spanning $[0.1, 0.3]\, {h^{-1}\,\rm Mpc}$, we observe that the orphan fractions substantially violate the upper limits in the outer regions of galaxy clusters beyond $0.3\,{h^{-1}\,\rm Mpc}$. We discuss the reasons, plausible improvements to the model, and how observations can be used to constrain such models further.

Project Dinos I: A joint lensing–dynamics constraint on the deviation from the power law in the mass profile of massive ellipticals

ABSTRACT The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing–dynamics analysis of a large statistical sample with 77 galaxy–galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy–galaxy lenses. We find that the Sloan Lens ACS Survey lens galaxies with a mean redshift of 0.2 are consistent with the power-law profile within 1.1σ (2.8σ) and the Strong Lensing Legacy Survey lens galaxies with a mean redshift of 0.6 are consistent within 0.8σ (2.1σ), for a spatially constant (Osipkov–Merritt) stellar anisotropy profile. We adopted the spatially constant anisotropy profile as our baseline choice based on previous dynamical observables of local ellipticals. However, spatially resolved stellar kinematics of lens galaxies are necessary to differentiate between the two anisotropy models. Future studies will use our lens models to constrain the mass distribution individually in the dark matter and baryonic components.

Simulating emission line galaxies for the next generation of large-scale structure surveys

ABSTRACT We investigate emission line galaxies across cosmic time by combining the modified L-Galaxies semi-analytical galaxy formation model with the JiuTian cosmological simulation. We improve the tidal disruption model of satellite galaxies in L-Galaxies to address the time dependence problem. We utilize the public code cloudy to compute emission line ratios for a grid of H ii region models. The emission line models assume the same initial mass function as that used to generate the spectral energy distribution of semi-analytical galaxies, ensuring a coherent treatment for modelling the full galaxy spectrum. By incorporating these emission line ratios with galaxy properties, we reproduce observed luminosity functions for H α, H β, [O ii], and [O iii] in the local Universe and at high redshifts. We also find good agreement between model predictions and observations for autocorrelation and cross-correlation functions of [O ii]-selected galaxies, as well as their luminosity dependence. The bias of emission line galaxies depends on both luminosity and redshift. At lower redshifts, it remains constant with increasing luminosity up to around $\sim 10^{42.5}\, {\rm erg\, s^{-1}}$ and then rises steeply for higher luminosities. The transition luminosity increases with redshift and becomes insignificant above z = 1.5. Generally, galaxy bias shows an increasing trend with redshift. However, for luminous galaxies, the bias is higher at low redshifts, as the strong luminosity dependence observed at low redshifts diminishes at higher redshifts. We provide a fitting formula for the bias of emission line galaxies as a function of luminosity and redshift, which can be utilized for large-scale structure studies with future galaxy surveys.

An improved halo occupation distribution prescription from UNITsim Hα emission-line galaxies: conformity and modified radial profile

ABSTRACT Emission-line galaxies (ELGs) are targeted by the new generation of spectroscopic surveys to make unprecedented measurements in cosmology from their distribution. Accurately interpreting these data require understanding the imprints imposed by the physics of galaxy formation and evolution on galaxy clustering. In this work, we utilize a semi-analytical model of galaxy formation (sage) to explore the necessary components for accurately reproducing the clustering of ELGs. We focus on developing a halo occupation distribution (HOD) prescription able to reproduce the clustering of sage galaxies. Typically, HOD models assume that satellite and central galaxies of a given type are independent events. We investigate the need for conformity, i.e. whether the average satellite occupation depends on the existence of a central galaxy of a given type. Incorporating conformity into HOD models is crucial for reproducing the clustering in the reference galaxy sample. Another aspect we investigate is the radial distribution of satellite galaxies within haloes. The traditional density profile models, Navarro–Frenk–White (NFW) and Einasto profiles, fail to accurately replicate the small-scale clustering measured for sage satellite galaxies. To overcome this limitation, we propose a generalization of the NFW profile, thereby enhancing our understanding of galaxy clustering.

Deprojection and stellar dynamical modelling of boxy/peanut bars in edge-on discs

ABSTRACT We present a new method to infer the 3D luminosity distributions of edge-on barred galaxies with boxy-peanut/X (BP/X) shaped structures from their 2D surface brightness distributions. Our method relies on forward modelling of newly introduced parametric 3D density distributions for the BP/X bar, disc and other components using an existing image fitting software package (imfit). We validate our method using an N-body simulation of a barred disc galaxy with a moderately strong BP/X shape. For fixed orientation angles, the derived 3D BP/X-shaped density distribution is shown to yield a gravitational potential that is accurate to at least 5 per cent and forces that are accurate to at least 15 per cent, with average errors being $\sim 1.5~{{\ \rm per\ cent}}$ for both. When additional quantities of interest, such as the orientation of the bar to the line of sight, its pattern speed, and the stellar mass-to-light ratio are unknown they can be recovered to high accuracy by providing the parametric density distribution to the Schwarzschild modelling code FORSTAND. We also explore the ability of our models to recover the mass of the central supermassive black hole. This method is the first to be able to accurately recover both the orientation of the bar to the line of sight and its pattern speed when the disc is perfectly edge-on.

Lopsided Galaxy Distributions in Massive TNG300 Clusters

Studies of isolated central galaxies and their satellites have shown that the spatial distributions of the satellites are distinctly “lopsided” with respect to the locations of the central galaxies. Here we extend this type of analysis to larger systems by analyzing the lopsidedness of 280 massive galaxy clusters in the IllustrisTNG300 Λ Cold Dark Matter simulation. Using a pairwise clustering statistic, we compute the probability that pairs of cluster galaxies are separated by a given polar angle difference, Δϕ, in the plane of the sky. Relative to the location of the central cluster galaxy, we find a statistically significant excess of galaxy pairs that are located on the same side of the cluster. The lopsidedness of the galaxy distribution is most pronounced for large clusters and for pairs of intrinsically red galaxies. The results summarized here were presented at the 243rd meeting of the American Astronomical Society.

Dark matter effects explanation with the torsion in the Minkowski space

Investigating rotation curves and the Tully–Fisher ratio within galaxies represents a central theme of extensive research and scientific interest. Despite several theoretical models, a comprehensive explanation of the observed correlation between galaxy types and their rotation curves remains elusive. This study endeavors to bridge this knowledge gap by delving into the discernible connection between the presence of dark matter and galaxy classification. By meticulously examining the gravitational field’s dependency on its source’s point symmetry, we introduce a novel theoretical framework that offers a coherent rationale for these empirical findings. Our proposed model explains the appearance of dark matter as a direct consequence of the reduction of point symmetry in gravitational systems. Neither arbitrary systems with a high mass density nor a perfectly spherically symmetric mass distribution give the observable effects of dark matter. Special attention was paid to the axial symmetry scenario as a reasonable approach for modeling the mass distribution in most galaxies. We thoroughly analyzed, showing strong agreement with experimental observations for dwarf, Sb, and Scd galaxies. Thus, our study provides a compelling theoretical foundation for elucidating the intricate interplay between galaxy types, rotation curves, and the presence of dark matter, shedding new light on the dynamics of the cosmos.

The phase-space distribution of the M 81 satellite system

The spatial distribution of dwarf galaxies around their host galaxies is a critical test for the standard model of cosmology because it probes the dynamics of dark matter halos and is independent of the internal baryonic processes of galaxies. Comoving planes of satellites have been found around the Milky Way, the Andromeda galaxy, and the nearby Cen A galaxy, which seems to be at odds with the standard model of galaxy formation. Another nearby galaxy group, with a putative flattened distribution of dwarf galaxies, is the M 81 group. In this paper, we present a quantitative analysis of the distribution of the M 81 satellites using a Hough transform to detect linear structures. Using this method, we confirm a flattened distribution of the dwarf galaxies. Depending on the morphological type, we find a minor-to-major axis ratio of the satellite distribution of 0.5 (all types) or 0.3 (dSph), which is in line with previous results for the M 81 group. Comparing the orientation of this flattened structure in 3D with the surrounding large-scale matter distribution, we find a strong alignment with the local sheet and the planes of satellites around the Andromeda galaxy and Cen A. Furthermore, the satellite system seems to be lopsided. Employing line-of-sight velocities for a subsample of the dwarfs, we find no signal of corotation. Comparing the flattening and motion of the M 81 dwarf galaxy system with TNG50 of the IllustrisTNG suite we find good agreement between observations and simulations, but caution that i) velocity information of half of the satellite population is still missing, ii) current velocities mainly come from dwarf irregulars clustered around NGC 3077, which may indicate an infall of a dwarf galaxy group, and iii) some of the dwarfs in our sample may be tidal dwarf galaxies. From the missing velocities, we predict that the observed frequency within IllustrisTNG may still range between 2 to 29%. Any final conclusions about the agreement or disagreement with cosmological models needs to wait for a more complete picture of the dwarf galaxy system.

Testing CCC+TL Cosmology with Observed Baryon Acoustic Oscillation Features

The primary purpose of this paper is to see how well a recently proposed new model fits (a) the position of the baryon acoustic oscillation (BAO) features observed in the large-scale distribution of galaxies and (b) the angular size measured for the sound horizon due to BAO imprinted in the cosmic microwave background (CMB) anisotropy. The new model is a hybrid model that combines the tired light (TL) theory with a variant of the ΛCDM model in which the cosmological constant is replaced with a covarying coupling constants’ (CCC) parameter α. This model, dubbed the CCC+TL model, can fit the Type Ia supernovae Pantheon+ data as accurately as the ΛCDM model, and also fit the angular size of cosmic dawn galaxies observed by the James Webb Space Telescope, which is in tension with the ΛCDM model. The results we obtained are 151.0 (±5.1) Mpc for the absolute BAO scale at the current epoch, and the angular size of the sound horizon θ sh = 0.°60, matching Planck’s observations at the surface of the last scattering when the baryon density is set to 100% of the matter density and ∣α∣ is increased by 5.6%. It remains to be seen if the new model is consistent with the CMB power spectrum, the Big Bang nucleosynthesis of light elements, and other critical observations.

The origin of lopsided satellite galaxy distribution around isolated systems in MillenniumTNG

ABSTRACT Dwarf satellites in galaxy groups are distributed in an anisotropic and asymmetric manner, which is called the ‘lopsided satellite distribution’. This lopsided signal has been observed not only in galaxy pairs but also in isolated systems. However, the physical origin of the lopsided signal in isolated systems is still unknown. In this work, we investigate this in the state-of-the-art hydrodynamical simulation of the MillenniumTNG Project by tracing each system back to high redshift. We find that the lopsided signal is dominated by satellites located in the outer regions of the halo and is also dominated by recently accreted satellites. The lopsided signal originates from the anisotropic accretion of galaxies from the surrounding large-scale structure and that, after accretion, the non-linear evolution of satellites inside the dark matter halo weakens the lopsidedness. The signal decreases as cosmic time passes because of a competition between anisotropic accretion and internal evolution within dark matter haloes. Our findings provide a useful perspective for the study of galaxy evolution, especially for the origin of the spatial satellite galaxy distributions.

Most of the photons that reionized the Universe came from dwarf galaxies.

The identification of sources driving cosmic reionization, a major phase transition from neutral hydrogen to ionized plasma around 600-800 Myr after the Big Bang1-3, has been a matter of debate4. Some models suggest that high ionizing emissivity and escape fractions (fesc) from quasars support their role in driving cosmic reionization5,6. Others propose that the high fesc values from bright galaxies generate sufficient ionizing radiation to drive this process7. Finally, a few studies suggest that the number density of faint galaxies, when combined with a stellar-mass-dependent model of ionizing efficiency and fesc, can effectively dominate cosmic reionization8,9. However, so far, comprehensive spectroscopic studies of low-mass galaxies have not been done because of their extreme faintness. Here we report an analysis of eight ultra-faint galaxies (in a very small field) during the epoch of reionization with absolute magnitudes between MUV ≈ -17 mag and -15 mag (down to 0.005L⋆ (refs. 10,11)). We find that faint galaxies during the first thousand million years of the Universe produce ionizing photons with log[ξion (Hz erg-1)] = 25.80 ± 0.14, a factor of 4 higher than commonly assumed values12. If this field is representative of the large-scale distribution of faint galaxies, the rate of ionizing photons exceeds that needed for reionization, even for escape fractions of the order of 5%.

Restoration of the Tully–Fisher Relation by Statistical Rectification

I employ the Lucy rectification algorithm to recover the inclination-corrected distribution of local disk galaxies in the plane of absolute magnitude (M i ) and H i velocity width (W 20). By considering the inclination angle as a random variable with a known probability distribution, the novel approach eliminates one major source of uncertainty in studies of the Tully–Fisher relation: inclination angle estimation from axial ratio. Leveraging the statistical strength derived from the entire sample of 28,264 H i-selected disk galaxies at z < 0.06 from the Arecibo Legacy Fast ALFA survey, I show that the restored distribution follows a sharp correlation that is approximately a power law between −16 > M i > −22: Mi=M0−2.5β[log(W20/250km/s)] , with M 0 = −19.77± 0.04 and β = 4.39 ± 0.06. At the brighter end (M i < −22), the slope of the correlation decreases to β ≈ 3.3, confirming previous results. Because the method accounts for measurement errors, the intrinsic dispersion of the correlation is directly measured: σ(logW20)≈0.06 dex between −17 > M i > −23, while σ(M i ) decreases from ∼0.8 in slow rotators to ∼0.4 in fast rotators. The statistical rectification method holds significant potential, especially in the studies of intermediate-to-high-redshift samples, where limited spatial resolution hinders precise measurements of inclination angles.

Topological bias: how haloes trace structural patterns in the cosmic web

ABSTRACT We trace the connectivity of the cosmic web as defined by haloes in the Planck-Millennium simulation using a persistence and Betti curve analysis. We normalize clustering up to the second-order correlation function and use our systematic topological analysis to correlate local information and properties of haloes with their multiscale geometrical environment of the cosmic web (elongated filamentary bridges and sheetlike walls). We capture the multiscale topology traced by the halo distribution through filtrations of the corresponding Delaunay tessellation. The resulting nested alpha shapes are sensitive to the local density, perfectly outline the local geometry, and contain the complete information on the multiscale topology. We find a remarkable linear relationship between halo masses and topology: haloes of different mass trace environments with different topological signature. This is a topological bias, an environmental structure bias independent of the halo clustering bias associated with the two-point correlation function. This mass-dependent linear scaling relation allows us to take clustering into account and determine the overall connectivity from a limited sample of galaxies. The presence of topological bias has major implications for the study of voids and filaments in the observed distribution of galaxies. The (infra)structure and shape of these key cosmic web components will strongly depend on the underlying galaxy sample. Their use as cosmological probes, with their properties influenced by cosmological parameters, will have to account for the subtleties of topological bias. This is of particular relevance with the large upcoming galaxy surveys such as DESI, Euclid, and the Vera Rubin telescope surveys.

A geostatistical analysis of multiscale metallicity variations in galaxies – III. Spatial resolution and data quality limits

ABSTRACT Geostatistical methods are powerful tools for understanding the spatial structure of the metallicity distribution of galaxies, and enable construction of accurate predictive models of the 2D metallicity distribution. However, so far these methods have only been applied to very high spatial resolution metallicity maps, leaving it uncertain if they will work on lower quality data. In this study, we apply geostatistical techniques to high-resolution spectroscopic maps of three local galaxies convolved to eight different spatial resolutions ranging from ∼40 pc to ∼1 kpc per pixel. We fit a geostatistical model to the data at all resolutions, and find that for metallicity maps where small-scale structure is visible by eye (with ≳10 resolution elements per Re), all parameters, including the metallicity correlation scale, can be recovered accurately. At all resolutions tested, we find that point metallicity predictions from such a geostatistical model outperform a circularly symmetric metallicity gradient model. We also explore dependence on the number of data points, and find that N ≳ 100 spatially resolved metallicity values are sufficient to train a geostatistical model that yields more accurate metallicity predictions than a radial gradient model. Finally, we investigate the potential detrimental effects of having spaxels smaller than an individual H ii region by repeating our analysis with metallicities integrated over H ii regions. We see that spaxel-based measurements have more noise, as expected, but the underlying spatial metallicity distribution can be recovered regardless of whether spaxels or integrated regions are used.