  In this Letter, we present mean optical+NIR color gradient estimates for 5080
early-type galaxies (ETGs) in the grizYJHK wavebands of the Sloan Digital Sky
Survey (SDSS) plus UKIRT Infrared Deep Sky Survey (UKIDSS). The color gradient
is estimated as the logarithmic slope of the radial color profile in ETGs. With
such a large sample size, we study the variation of the mean color gradient as
a function of waveband with unprecedented accuracy. We find that (i) color
gradients are mainly due, on average, to a metallicity variation of about
-0.4dex per decade in galaxy radius; and (ii) a small, but significant,
positive age gradient is present, on average, in ETGs, with the inner stellar
population being slightly younger, by ~0.1dex per radial decade, than the outer
one. Also, we show that the presence of a positive mean age gradient in ETGs,
as found in the present study, implies their effective radius to be smaller at
high z, consistent with observations.
