Galaxies in the local Universe are known to follow bimodal distributions in the global stellar populations properties. We analyze the distribution of the local average stellar-population ages of 654,053 sub-galactic regions resolved on ~1-kpc scales in a volume-corrected sample of 394 galaxies, drawn from the CALIFA-DR3 integral-field-spectroscopy survey and complemented by SDSS imaging. We find a bimodal local-age distribution, with an old and a young peak primarily due to regions in early-type galaxies and star-forming regions of spirals, respectively. Within spiral galaxies, the older ages of bulges and inter-arm regions relative to spiral arms support an internal age bimodality. Although regions of higher stellar-mass surface-density, mu*, are typically older, mu* alone does not determine the stellar population age and a bimodal distribution is found at any fixed mu*. We identify an "old ridge" of regions of age ~9 Gyr, independent of mu*, and a "young sequence" of regions with age increasing with mu* from 1-1.5 Gyr to 4-5 Gyr. We interpret the former as regions containing only old stars, and the latter as regions where the relative contamination of old stellar populations by young stars decreases as mu* increases. The reason why this bimodal age distribution is not inconsistent with the unimodal shape of the cosmic-averaged star-formation history is that i) the dominating contribution by young stars biases the age low with respect to the average epoch of star formation, and ii) the use of a single average age per region is unable to represent the full time-extent of the star-formation history of "young-sequence" regions.

While studies of gas-phase metallicity gradients in disc galaxies are common, very little has been done in the acquisition of stellar abundance gradients in the same regions. We present here a comparative study of the stellar metallicity and age distributions in a sample of 62 nearly face-on, spiral galaxies with and without bars, using data from the CALIFA survey. We measure the slopes of the gradients and study their relation with other properties of the galaxies. We find that the mean stellar age and metallicity gradients in the disc are shallow and negative. Furthermore, when normalized to the effective radius of the disc, the slope of the stellar population gradients does not correlate with the mass or with the morphological type of the galaxies. Contrary to this, the values of both age and metallicity at $\sim$2.5 scale-lengths correlate with the central velocity dispersion in a similar manner to the central values of the bulges, although bulges show, on average, older ages and higher metallicities than the discs. One of the goals of the present paper is to test the theoretical prediction that non-linear coupling between the bar and the spiral arms is an efficient mechanism for producing radial migrations across significant distances within discs. The process of radial migration should flatten the stellar metallicity gradient with time and, therefore, we would expect flatter stellar metallicity gradients in barred galaxies. However, we do not find any difference in the metallicity or age gradients in galaxies with without bars. We discuss possible scenarios that can lead to this absence of difference.

In a companion paper we have presented many products derived from the application of the spectral synthesis code STARLIGHT to datacubes from the CALIFA survey, including 2D maps of stellar population properties and 1D averages in the temporal and spatial dimensions. Here we evaluate the uncertainties in these products. Uncertainties due to noise and spectral shape calibration errors and to the synthesis method are investigated by means of a suite of simulations based on 1638 CALIFA spectra for NGC 2916, with perturbations amplitudes gauged in terms of the expected errors. A separate study was conducted to assess uncertainties related to the choice of evolutionary synthesis models. We compare results obtained with the Bruzual & Charlot models, a preliminary update of them, and a combination of spectra derived from the Granada and MILES models. About 100k CALIFA spectra are used in this comparison. Noise and shape-related errors at the level expected for CALIFA propagate to 0.10-0.15 dex uncertainties in stellar masses, mean ages and metallicities. Uncertainties in A_V increase from 0.06 mag in the case of random noise to 0.16 mag for shape errors. Higher order products such as SFHs are more uncertain, but still relatively stable. Due to the large number statistics of datacubes, spatial averaging reduces uncertainties while preserving information on the history and structure of stellar populations. Radial profiles of global properties, as well as SFHs averaged over different regions are much more stable than for individual spaxels. Uncertainties related to the choice of base models are larger than those associated with data and method. Differences in mean age, mass and metallicity are ~ 0.15 to 0.25 dex, and 0.1 mag in A_V. Spectral residuals are ~ 1% on average, but with systematic features of up to 4%. The origin of these features is discussed. (Abridged)

Fossil record methods based on spectral synthesis techniques have matured over the past decade, and their application to integrated galaxy spectra fostered substantial advances on the understanding of galaxies and their evolution. Yet, because of the lack of spatial resolution, these studies are limited to a global view, providing no information about the internal physics of galaxies. Motivated by the CALIFA survey, which is gathering Integral Field Spectroscopy over the full optical extent of 600 galaxies, we have developed an end-to-end pipeline which: (i) partitions the observed data cube into Voronoi zones in order to, when necessary and taking due account of correlated errors, increase the S/N, (ii) extracts spectra, including propagated errors and bad-pixel flags, (iii) feeds the spectra into the STARLIGHT spectral synthesis code, (iv) packs the results for all galaxy zones into a single file, (v) performs a series of post-processing operations, including zone-to-pixel image reconstruction and unpacking the spectral and stellar population properties into multi-dimensional time, metallicity, and spatial coordinates. This paper provides an illustrated description of this whole pipeline and its products. Using data for the nearby spiral NGC 2916 as a show case, we go through each of the steps involved, presenting ways of visualizing and analyzing this manifold. These include 2D maps of properties such as the v-field, stellar extinction, mean ages and metallicities, mass surface densities, star formation rates on different time scales and normalized in different ways, 1D averages in the temporal and spatial dimensions, projections of the stellar light and mass growth (x,y,t) cubes onto radius-age diagrams, etc. The results illustrate the richness of the combination of IFS data with spectral synthesis, providing a glimpse of what is to come from CALIFA and future surveys. (Abridged)

Integral field spectroscopy (IFS) is presented for star forming regions in M33. A central area of 300 x 500 pc^2 and the external HII region IC 132, at a galactocentric distance {\sim} 19arcmin (4.69 kpc) were observed with the Potsdam Multi Aperture Spectrophotometer (PMAS) instrument at the 3.5 m telescope of the Calar Alto Hispano- Alem\'an observatory (CAHA). The spectral coverage goes from 3600 A to 1{\mu}m to include from [OII]{\lambda}3727 A to the near infrared lines required for deriving sulphur electron temperature and abundance diagnostics. Local conditions within individual HII regions are presented in the form of emission line fluxes and physical conditions for each spatial resolution element (spaxel) and for segments with similar H{\alpha} surface brightness. A clear dichotomy is observed when comparing the central to outer disc HII regions. While the external HII region has higher electron temperature plus larger H{\beta} equivalent width, size and excitation, the central region has higher extinction and metal content. The dichotomy extends to the BPT diagnostic diagrams that show two orthogonal broad distributions of points. By comparing with pseudo-3D photoionization models we conclude that the bulk observed differences are probably related to a different ionization parameter and metallicity. Wolf-Rayet features are detected in IC 132, and resolved into two concentrations whose integrated spectra were used to estimate the characteristic number of WR stars. No WR features were detected in the central HII regions despite their higher metallicity.