REVIEW 2 major objections 4 minor 54 references
Resolved and Integrated Stellar Masses in the SDSS-IV/MaNGA Survey, Paper II: Applications of PCA-based stellar mass estimates
T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper shows that stellar masses from galaxy-coadded (unresolved) spectra are systematically low, with the largest deficits in dusty, edge-on galaxies, and identifies differential dust attenuation as the dominant cause.
desk verdict A workmanlike MaNGA application paper with a useful catalog; the main new claim about dust-driven luminosity-weighting bias is built on a clearly flagged proxy that deserves a direct check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central working object is the PCA-based i-band stellar mass-to-light ratio, $\log \Upsilon^*_i$: a basis set of six eigenspectra, trained on a library of about 40,000 composite stellar population spectra, onto which each observed spectrum is projected to recover a mass-to-light ratio with covariate uncertainties. Because it is evaluated spaxel-by-spaxel across MaNGA datacubes, it yields resolved stellar mass maps that power every comparison in the paper: luminosity-weighting those maps reproduces the effect of an unresolved galaxy spectrum; summing inside and outside the IFU sets up the aperture corrections; and deprojected radial profiles can be placed against kinematic mass surface densities. The luminosity-weighted combination rule is therefore the load-bearing operation for the paper's main negative finding about spatial binning.
What would settle it
Coadd the MaNGA spaxel spectra for a sample of galaxies into one spectrum per galaxy, run the same PCA fit on the coadded spectrum, and compare the resulting stellar mass with the sum of the resolved PCA masses; if the coadded-fit masses do not show a deficit that correlates with $\sigma_{\tau_V}$ and inclination, the luminosity-weighting proxy—not unresolved spectroscopy itself—is creating the reported bias.
Extended reading notes
Core claim
The central claim is that stellar mass is not invariant under spatial binning: when per-spaxel PCA mass-to-light estimates are combined by luminosity weighting to mimic a single galaxy spectrum, the resulting total stellar mass is lower than the sum over the resolved spaxels. The deficit is small on average but strongly structured—edge-on disks and galaxies with prominent dust lanes show the largest shortfalls, while face-on, low-dust galaxies are nearly unbiased. An ordinary least-squares fit of the deficit to average V-band optical depth, its spaxel-to-spaxel dispersion, and central specific star-formation rate finds that the dispersion term, $\sigma_{\tau_V}$, is the dominant driver, with confidence intervals excluding zero; central star formation matters only at roughly a tenth of the dust term's strength. The paper interprets this as evidence that differential dust attenuation, rather than simple outshining by young stars, is the main mechanism behind unresolved-spectrum mass deficits. The same mass estimates also appear heavier than photometric catalog masses by about 0.1–0.15 dex at red colors, and the comparison with DiskMass dynamical surface densities can be brought into qualitative agreement by reducing the adopted disk scale-heights by a factor of about 1.5.
Load-bearing premise
The analysis assumes that the stellar mass-to-light ratio a real coadded galaxy spectrum would produce equals the luminosity-weighted average of the per-spaxel PCA mass-to-light ratios, but the paper never coadds the MaNGA spectra and reruns the fit; if spectral fitting responds nonlinearly to signal-to-noise or mixed stellar populations, the claimed unresolved-spectrum bias could be misstated.
Editorial extensions
If this is right
- Single-fiber and low-spatial-resolution surveys will tend to underestimate total stellar masses, with the deficit growing for edge-on and dusty galaxies.
- The correction for unresolved-spectrum masses must include differential dust attenuation; average attenuation alone and star-formation rate are not enough.
- The color-based mass-to-light ('CMLR') aperture-correction method is preferred over the 'ring' method, and the two agree to better than 0.1 dex in total mass for most galaxies.
- If the PCA/SPS masses are close to correct, the DiskMass scale-heights would need to shrink by about a factor of 1.5 to reconcile the two mass scales.
- The soon-to-be-released SDSS value-added catalog will provide both resolved mass maps and aperture-corrected total stellar masses for MaNGA galaxies.
Reading between the lines
- A testable extension the paper leaves implicit: the fitted regression could serve as an empirical correction for unresolved surveys, if the spatial scatter in optical depth $\sigma_{\tau_V}$ can be estimated from multiband imaging or attenuation maps.
- Because the test uses luminosity-weighted resolved mass-to-light ratios rather than a real coadded refit, the true bias of actual coadded spectra could be larger or smaller; a direct coadd-and-refit experiment on the same datacubes would settle the size of the effect.
- If the disk scale-height tension survives in a larger sample, disks would be closer to maximal than the DiskMass fiducial assumptions imply, shifting dark-matter halo fits for disk galaxies.
- The ~0.12 dex scatter between measured attenuation and axis ratio implies that inclination-based dust corrections in unresolved surveys introduce galaxy-to-galaxy scatter at about that level.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper (Paper II of a series) presents applications of PCA-based stellar mass-to-light ratio estimates to resolved MaNGA spectroscopy. It compares PCA-derived stellar mass surface densities with DiskMass dynamical surface densities for three galaxies, finding that moderate overestimates of disk scale-heights (by a factor of roughly 1.5) would reconcile the two for two of the three galaxies. It constructs a total stellar mass catalog using two aperture-correction methods (the 'ring' method and a color-based CMLR method), recommends the CMLR method, and compares the resulting masses with NSA and JHU-MPA photometric masses, finding typical discrepancies of about 0.1 dex. The paper's most novel result is an analysis of luminosity-weighting biases: comparing IFU-summed masses with masses obtained by luminosity-weighting the resolved mass-to-light ratios, the authors find mass deficits that are largest for edge-on systems and systems with dust lanes, and a multivariate regression attributes the deficit primarily to the spatial dispersion of V-band optical depth (differential attenuation) rather than to average attenuation or central star formation. The paper also describes a planned SDSS Value-Added Catalog.
Significance. If validated, this is a useful contribution for the MaNGA community and for resolved-spectroscopy mass estimation more broadly. The external comparisons to kinematical measurements and to independently derived photometric masses provide helpful anchoring for the PCA method, and the paper is transparent about the main caveats, including the DiskMass scale-height systematics and the absence of uncertainties in the regression predictors. The luminosity-weighting analysis is the most novel element, and the paper correctly identifies a potential source of systematic error in unresolved or spatially binned mass estimates. However, the central claim about differential dust attenuation rests on a proxy whose equivalence to a true coadded-spectrum measurement is not directly demonstrated; this is a load-bearing assumption that must be validated or substantially de-emphasized before the main result can be fully accepted.
major comments (2)
- [Section 4.3, Eq. (4), Table 2] The response variable in the central regression is log(M*/M_LW), but M_LW is not obtained by coadding the MaNGA spaxel spectra and rerunning the PCA fitting, as the 'galaxy-coadded spectrum' framing implies. It is instead the total IFU luminosity multiplied by the luminosity-weighted average of the resolved PCA mass-to-light ratios. This replacement is exact only if the PCA-based M/L estimator is a linear functional of the input spectrum; PCA projection with covariate uncertainties followed by a nonlinear mapping from PCA scores to M/L is not such a functional. If the nonlinearity has a signal-to-noise dependence or correlates with dust geometry, the measured mass deficit, and hence the conclusion that sigma_tauV dominates, could be an artifact of the proxy. The external support cited (Ibarra-Medel et al. 2019) uses simulations with a different spatial-binning scheme and does not validate this specific proxy. I request a direct validation: coadd the spectra (or realistic synthetic spectra) and rerun the full PCA fit, or alternatively demonstrate with synthetic MaNGA-like observations that the luminosity-weighted proxy reproduces the coadded-spectrum result to within the claimed 0.05–0.1 dex.
- [Section 4.3, Table 2] The ordinary least-squares fit in Eq. (4) uses predictors (tauV, sigma_tauV, eta) that are derived from the same PCA fits that produce the response variable, and the uncertainties on these predictors are not propagated into the reported confidence intervals or p-values. The authors note this caveat, but because the claim that differential attenuation is the dominant factor is one of the paper's main new results, the significance should be checked for robustness to correlated errors, for example by Monte Carlo sampling the PCA posteriors or by performing a split-sample test by inclination and star-formation activity. As presented, the p-values in Table 2 likely overstate the confidence in the beta1 coefficient.
minor comments (4)
- [Abstract] The abstract contains a typo ('galaxys' should be 'galaxy's'); similar minor typographical errors appear elsewhere in the text.
- [References] The two Gallazzi et al. entries for 2005a and 2005b are identical in the reference list; one of these citations is likely intended to refer to a different article.
- [Abstract and Table 1] The abstract's statement that the DiskMass tension 'could be resolved if the disk scale-heights were overestimated by a factor of ~1.5' is only supported by two of the three galaxies (f_hz = 1.83 and 1.68), while UGC3997 gives f_hz = 0.90. The abstract should make this sample size and heterogeneity explicit, or the claim should be softened.
- [Section 4.3, Table 3] Several of the twenty galaxies with the largest mass deficits are flagged as containing foreground stars, overlapping companions, or poor fits, so the visual impression from Figure 11 may overstate the dust-lane/edge-on connection. The claim is better supported by the regression than by the table, and the text should clarify that the table is illustrative rather than definitive evidence.
Circularity Check
No significant circularity: the Paper II mass estimates are benchmarked against external dynamical and photometric data, and the luminosity-weighting mass deficit is an explicitly constructed diagnostic rather than a hidden prediction.
full rationale
The paper's derivation chain is not circular. The PCA-based mass-to-light estimates are inherited from Paper I, but Paper II does not re-derive them from its own conclusions; it applies them to MaNGA data and validates them against external benchmarks: DiskMass dynamical surface densities (Section 3), NSA and JHU-MPA photometric masses (Section 4.2), and an independent simulation study by Ibarra-Medel et al. (2019) for spatial-coadding effects. The scale-height factor fhz is fit to the DMS comparison, but the authors explicitly state that 'claims regarding specific values of fhz are outside the scope of this paper,' so it is presented as an exploratory diagnostic, not as a predicted quantity. The Section 4.3 mass deficit is also not circular: the paper explicitly defines M_LW as the mass obtained by 'flux-weighting the corresponding, spatially-resolved mass-to-light ratios, and multiplying by the total luminosity in the IFU.' This is a transparent construction of an implicitly luminosity-weighted diagnostic, not a claim that a true coadded-spectrum fit was run. The later OLS regression of that diagnostic on tauV dispersion is an internal correlation with acknowledged limitations ('we do not include uncertainties on tauV, sigma_tauV, or log eta'); concerns about proxy fidelity are validity caveats rather than definitional circularity. Self-citations to Paper I and Chen et al. (2012) supply the fitting method, but they are not load-bearing in the sense of forcing the new results: the core comparisons are external, and the mass-deficit analysis is an internal comparison of two explicitly defined combinations of the same resolved maps. No equation reduces by construction to its own inputs, and no fitted parameter is renamed as a prediction.
Assumptions & free parameters
free parameters (2)
- fhz =
0.90, 1.83, 1.68 for UGC3997, UGC4107, UGC4368
- OLS regression coefficients alpha, beta0, beta1, beta2 =
alpha ~ 0.02, beta0 ~ -0.003, beta1 ~ 0.05, beta2 ~ 0.003 (Table 2)
assumptions (5)
- domain assumption The Bruzual & Charlot (2003) stellar population synthesis models and the 40,000-SFH training library from Paper I accurately represent real galaxy optical spectra.
- domain assumption A Chabrier or Kroupa stellar initial mass function applies to all galaxies.
- domain assumption For the three DiskMass comparison galaxies, the stellar disk dominates the vertical dynamical mass, so the stellar mass surface density should be slightly below the dynamical mass surface density.
- ad hoc to paper The luminosity-weighted combination of resolved mass-to-light ratios is equivalent to fitting a single coadded galaxy spectrum.
- domain assumption The outermost 0.5 Re ring's median mass-to-light ratio is representative of the galaxy outskirts for the ring aperture correction, and the color-mass-to-light relation from Paper I applies to the missing flux for the CMLR correction.
Cite this review
Pith. "Pith review of Resolved and Integrated Stellar Masses in the SDSS-IV/MaNGA Survey, Paper II: Applications of PCA-based stellar mass estimates." pith.science (2026). https://pith.science/paper/2VZCMFOD
@misc{pith2026190802331,
author = {Pith},
title = {Pith review of: Resolved and Integrated Stellar Masses in the SDSS-IV/MaNGA Survey, Paper II: Applications of PCA-based stellar mass estimates},
year = {2026},
howpublished = {\url{https://pith.science/paper/2VZCMFOD}},
note = {Machine review of arXiv:1908.02331}
}
abstract
A galaxy's stellar mass is one of its most fundamental properties, but it remains challenging to measure reliably. With the advent of very large optical spectroscopic surveys, efficient methods that can make use of low signal-to-noise spectra are needed. With this in mind, we created a new software package for estimating effective stellar mass-to-light ratios $\log \Upsilon^*$ that uses principal component analysis(PCA) basis set to optimize the comparison between observed spectra and a large library of stellar population synthesis models. In Paper I, we showed that a with a set of six PCA basis vectors we could faithfully represent most optical spectra from the Mapping Nearby Galaxies at APO (MaNGA) survey;and we tested the accuracy of our M/L estimates using synthetic spectra. Here, we explore sources of systematic error in our mass measurements by comparing our new measurements to data from the literature. We compare our stellar mass surface density estimates to kinematics-derived dynamical mass surface density measurements from the DiskMass Survey and find some tension between the two which could be resolved if the disk scale-heights used in the kinematic analysis were overestimated by a factor of $\sim 1.5$. We formulate an aperture-corrected stellar mass catalog for the MaNGA survey, and compare to previous stellar mass estimates based on multi-band optical photometry, finding typical discrepancies of 0.1 dex. Using the spatially resolved MaNGA data, we evaluate the impact of estimating total stellar masses from spatially unresolved spectra, and we explore how the biases that result from unresolved spectra depend upon the galaxy's dust extinction and star formation rate. Finally, we describe a SDSS Value-Added Catalog which will include both spatially resolved and total (aperture-corrected) stellar masses for MaNGA galaxies.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[1]
2011, ApJS, 193, 29
Aihara, H., Allende Prieto, C., An, D., et al. 2011, ApJS, 193, 29
2011
-
[2]
Aniyan, S., Freeman, K. C., Gerhard, O. E., Arnaboldi, M., & Flynn, C. 2016, MNRAS, 456, 1484
work page 2016
-
[3]
Aniyan, S., Freeman, K. C., Arnaboldi, M., et al. 2018, ArXiv e-prints, arXiv:1802.00465 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33
arXiv 2018
-
[4]
K., Balogh, M
Baldry, I. K., Balogh, M. L., Bower, R. G., et al. 2006, MNRAS, 373, 469
2006
-
[5]
L., Baldry, I
Balogh, M. L., Baldry, I. K., Nichol, R., et al. 2004, ApJL, 615, L101
2004
-
[6]
F., & de Jong, R
Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212
2001
-
[7]
F., McIntosh, D
Bell, E. F., McIntosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289
2003
-
[8]
Bershady, M. A., Martinsson, T. P. K., Verheijen, M. A. W., et al. 2011, ApJL, 739, L47
work page 2011
Show all 54 references
-
[9]
2011, AJ, 142, 31
Price-Whelan, A. 2011, AJ, 142, 31
2011
-
[10]
R., & Roweis, S
Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 —. 2017, kcorrect: Calculate K-corrections between observed and desired bandpasses, Astrophysics Source Code Library, , , ascl:1701.010
2007
-
[11]
R., Bershady, M
Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, AJ, 154, 28
2017
-
[12]
F., Weijmans, A.-M., van den Bosch, R., et al
Boardman, N. F., Weijmans, A.-M., van den Bosch, R., et al. 2017, MNRAS, 471, 4005
2017
-
[13]
2003, MNRAS, 344, 1000
Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000
2003
-
[14]
A., & Law, D
Bundy, K., Bershady, M. A., & Law, D. R. e. a. 2015, ApJ, 798, 7
2015
-
[15]
2013, Star Formation Rate Indicators, ed
Calzetti, D. 2013, Star Formation Rate Indicators, ed. J. Falc´ on-Barroso & J. H. Knapen, 419
2013
-
[16]
2003, MNRAS, 342, 345
Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345
2003
-
[17]
2003, PASP, 115, 763
Chabrier, G. 2003, PASP, 115, 763
2003
-
[18]
Charlot, S., & Fall, S. M. 2000, ApJ, 539, 718
2000
-
[19]
A., et al
Chen, Y.-M., Kauffmann, G., Tremonti, C. A., et al. 2012, MNRAS, 421, 314
2012
-
[20]
S., Schlegel, D
Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, AJ, 145, 10
2013
-
[21]
2010, AJ, 139, 1628
Doi, M., Tanaka, M., Fukugita, M., et al. 2010, AJ, 139, 1628
2010
-
[22]
A., et al
Drory, N., MacDonald, N., Bershady, M. A., et al. 2015, AJ, 149, 77
2015
-
[23]
Eigenbrot, A., & Bershady, M. A. 2018, ApJ, 853, 114
2018
-
[24]
E., et al
Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748
1996
-
[25]
Gallazzi, A., & Bell, E. F. 2009, ApJS, 185, 253
2009
-
[26]
Gallazzi, A., Charlot, S., Brinchmann, J., & White, S. D. M. 2006, MNRAS, 370, 1106
2006
-
[27]
Gallazzi, A., Charlot, S., Brinchmann, J., White, S. D. M., & Tremonti, C. A. 2005b, MNRAS, 362, 41 Gonz´ alez Delgado, R. M., P´ erez, E., Cid Fernandes, R., et al. 2014, A&A, 562, A47
2014
-
[28]
E., Siegmund, W
Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332
2006
-
[29]
Hessman, F. V. 2017, MNRAS, 469, 1147
2017
-
[30]
Hunter, J. D. 2007, Computing In Science & Engineering, 9, 90
2007
-
[31]
J., Avila-Reese, V., S´ anchez, S
Ibarra-Medel, H. J., Avila-Reese, V., S´ anchez, S. F., Gonz´ alez-Samaniego, A., & Rodr´ ıguez-Puebla, A. 2019, MNRAS, 483, 4525 Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003a, MNRAS, 341, 33 Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003b, MNRAS, 346, 1055
2019
-
[32]
J., Dopita, M
Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121 La Barbera, F., Vazdekis, A., Ferreras, I., et al. 2016, MNRAS, 457, 1468
2001
-
[33]
R., Yan, R., Bershady, M
Law, D. R., Yan, R., Bershady, M. A., et al. 2015, AJ, 150, 19
2015
-
[34]
R., Cherinka, B., Yan, R., et al
Law, D. R., Cherinka, B., Yan, R., et al. 2016, AJ, 152, 83
2016
-
[35]
2017, ApJ, 838, 77
Li, H., Ge, J., Mao, S., et al. 2017, ApJ, 838, 77
2017
-
[36]
2019, arXiv e-prints, arXiv:1903.09282 Mart´ ın-Navarro, I., La Barbera, F., Vazdekis, A., Falc´ on-Barroso, J., & Ferreras, I
Li, R., Li, H., Shao, S., et al. 2019, arXiv e-prints, arXiv:1903.09282 Mart´ ın-Navarro, I., La Barbera, F., Vazdekis, A., Falc´ on-Barroso, J., & Ferreras, I. 2015, MNRAS, 447, 1033 Mart´ ınez-Garc´ ıa, E. E., Gonz´ alez-L´ opezlira, R. A., Gladis
2019 arXiv
-
[37]
Magris, C., & Bruzual, A. G. 2017, ApJ, 835, 93
2017
-
[38]
S., & de Blok, W
McGaugh, S. S., & de Blok, W. J. G. 1997, ApJ, 481, 689 18 Pace et al
1997
-
[39]
J., Tremonti, C., Chen, Y., et al
Pace, Z. J., Tremonti, C., Chen, Y., et al. 2019, ApJ
2019
-
[40]
M., Charlot, S., et al
Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267
2007
-
[41]
Schechtman-Rook, A., & Bershady, M. A. 2013, ApJ, 773, 45 —. 2014, ApJ, 795, 136
2013
-
[42]
2010, in 9th Python in Science Conference
Seabold, S., & Perktold, J. 2010, in 9th Python in Science Conference
2010
-
[43]
A., Gunn, J
Smee, S. A., Gunn, J. E., Uomoto, A., et al. 2013, AJ, 146, 32
2013
-
[44]
2015, MNRAS, 452, 235
Sorba, R., & Sawicki, M. 2015, MNRAS, 452, 235
2015
-
[45]
A., Bershady, M
Swaters, R. A., Bershady, M. A., Martinsson, T. P. K., et al. 2014, ApJL, 797, L28
2014
-
[46]
Tinsley, B. M. 1972, A&A, 20, 383 —. 1973, ApJ, 186, 35
1972
-
[47]
R., Cardone, V
Tortora, C., Napolitano, N. R., Cardone, V. F., et al. 2010, MNRAS, 407, 144
2010
-
[48]
R., Romanowsky, A
Tortora, C., Napolitano, N. R., Romanowsky, A. J., et al. 2011, MNRAS, 418, 1557
2011
-
[49]
A., Heckman, T
Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898
2004
-
[50]
A., Bundy, K., Diamond-Stanic, A
Wake, D. A., Bundy, K., Diamond-Stanic, A. M., et al. 2017, AJ, 154, 86
2017
-
[51]
B., Bershady, M
Westfall, K. B., Bershady, M. A., & Verheijen, M. A. W. 2011, ApJS, 193, 21
2011
-
[52]
B., Cappellari, M., Bershady, M
Westfall, K. B., Cappellari, M., Bershady, M. A., et al. 2019, arXiv e-prints, arXiv:1901.00856
2019 arXiv
-
[53]
G., Adelman, J., Anderson, Jr., J
York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579
2000
-
[54]
2009, MNRAS, 400, 1181
Zibetti, S., Charlot, S., & Rix, H.-W. 2009, MNRAS, 400, 1181
2009
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.