We present robust statistical estimates of the accuracy of early-type galaxy stellar masses derived from spectral energy distribution (SED) fitting as functions of various empirical and theoretical assumptions. Using large samples consisting of 40,000 galaxies from the Sloan Digital Sky Survey, of which 5,000 are also in the UKIRT Infrared Deep Sky Survey, with spectroscopic redshifts in the range 0.05 \leq z \leq 0.095, we test the reliability of some commonly used stellar population models and extinction laws for computing stellar masses. Spectroscopic ages (t), metallicities (Z), and extinctions (A) are also computed from fits to SDSS spectra using various population models. These constraints are used in additional tests to estimate the systematic errors in the stellar masses derived from SED fitting, where t, Z, and A are typically left as free parameters. We find reasonable agreement in mass estimates among stellar population models, with variation of the IMF and extinction law yielding systematic biases on the mass of nearly a factor of 2, in agreement with other studies. Removing the near-infrared bands changes the statistical bias in mass by only 0.06 dex, adding uncertainties of 0.1 dex at the 95% CL. In contrast, we find that removing an ultraviolet band is more critical, introducing 2{\sigma} uncertainties of 0.15 dex. Finally, we find that stellar masses are less affected by absence of metallicity and/or dust extinction knowledge. However, there is a definite systematic offset in the mass estimate when the stellar population age is unknown, up to a factor of 2.5 for very old (12 Gyr) stellar populations. We present the stellar masses for our sample, corrected for the measured systematic biases due to photometrically determined ages, finding that age errors produce lower stellar masses by 0.15 dex, with errors of 0.02 dex at the 95% CL for the median stellar age subsample.

We present an analysis of stellar population gradients in 4,546 Early-Type Galaxies with photometry in $grizYHJK$ along with optical spectroscopy. A new approach is described which utilizes color information to constrain age and metallicity gradients. Defining an effective color gradient, $\nabla_{\star}$, which incorporates all of the available color indices, we investigate how $\nabla_{\star}$ varies with galaxy mass proxies, i.e. velocity dispersion, stellar (M_star) and dynamical (M_dyn) masses, as well as age, metallicity, and alpha/Fe. ETGs with M_dyn larger than 8.5 x 10^10, M_odot have increasing age gradients and decreasing metallicity gradients wrt mass, metallicity, and enhancement. We find that velocity dispersion and alpha/Fe are the main drivers of these correlations. ETGs with 2.5 x 10^10 M_odot =< M_dyn =< 8.5 x 10^10 M_odot, show no correlation of age, metallicity, and color gradients wrt mass, although color gradients still correlate with stellar population parameters, and these correlations are independent of each other. In both mass regimes, the striking anti-correlation between color gradient and alpha-enhancement is significant at \sim 4sigma, and results from the fact that metallicity gradient decreases with alpha/Fe. This anti-correlation may reflect the fact that star formation and metallicity enrichment are regulated by the interplay between the energy input from supernovae, and the temperature and pressure of the hot X-ray gas in ETGs. For all mass ranges, positive age gradients are associated with old galaxies (>5-7 Gyr). For galaxies younger than \sim 5 Gyr, mostly at low-mass, the age gradient tends to be anti-correlated with the Age parameter, with more positive gradients at younger ages.