This paper presents cosmological results based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. Our results are in very good agreement with the 2013 analysis of the Planck nominal-mission temperature data, but with increased precision. The temperature and polarization power spectra are consistent with the standard spatially-flat 6-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations (denoted "base ΛCDM" in this paper). From the Planck temperature data combined with Planck lensing, for this cosmology we find a Hubble constant, H 0 = (67.8 ± 0.9) km s −1 Mpc −1 , a matter density parameter Ω m = 0.308 ± 0.012, and a tilted scalar spectral index with n s = 0.968 ± 0.006, consistent with the 2013 analysis. Note that in this abstract we quote 68% confidence limits on measured parameters and 95% upper limits on other parameters. We present the first results of polarization measurements with the Low Frequency Instrument at large angular scales. Combined with the Planck temperature and lensing data, these measurements give a reionization optical depth of τ = 0.066 ± 0.016, corresponding to a reionization redshift of z re = 8.8+1.7 −1.4 . These results are consistent with those from WMAP polarization measurements cleaned for dust emission using 353-GHz polarization maps from the High Frequency Instrument. We find no evidence for any departure from base ΛCDM in the neutrino sector of the theory; for example, combining Planck observations with other astrophysical data we find N eff = 3.15 ± 0.23 for the effective number of relativistic degrees of freedom, consistent with the value N eff = 3.046 of the Standard Model of particle physics. The sum of neutrino masses is constrained to m ν < 0.23 eV. The spatial curvature of our Universe is found to be very close to zero, with |Ω K | < 0.005. Adding a tensor component as a single-parameter extension to base ΛCDM we find an upper limit on the tensor-to-scalar ratio of r 0.002 < 0.11, consistent with the Planck 2013 results and consistent with the B-mode polarization constraints from a joint analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP B-mode data to our analysis leads to a tighter constraint of r 0.002 < 0.09 and disfavours inflationary models with a V(φ) ∝ φ 2 potential. The addition of Planck polarization data leads to strong constraints on deviations from a purely adiabatic spectrum of fluctuations. We find no evidence for any contribution from isocurvature perturbations or from cosmic defects. Combining Planck data with other astrophysical data, including Type Ia supernovae, the equation of state of dark energy is constrained to w = −1.006 ± 0.045, consistent with the expected Corresponding author: G. Efstathiou, e-mail: gpe@ast.cam.ac.ukArticle published by EDP Sciences A13, page 1 of 63 A&A 594, A13 (2016) value for a cosmological constant. The standard big bang nucleosynthesis predictions for the helium and deuterium abundanc...