A simple cosmological model with only six parameters (matter density, m h 2 , baryon density, b h 2 , Hubble constant, H 0 , amplitude of fluctuations, 8 , optical depth, , and a slope for the scalar perturbation spectrum, n s ) fits not only the 3 year WMAP temperature and polarization data, but also small-scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the bestfit values for cosmological parameters for the power-law flat à cold dark matter (ÃCDM) model are ( m h 2 ; b h 2 ; h; n s ; ; 8 ) ¼ (0:1277 þ0:0080 À0:0079 ;0:02229 AE 0:00073;0:732 þ0:031 À0:032 ;0:958 AE 0:016;0:089 AE 0:030; 0:761 þ0:049 À0:048 ). The 3 year data dramatically shrink the allowed volume in this six-dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power-law spectrum, the WMAP data alone require dark matter and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n s ¼ 1; r ¼ 0). Adding additional data sets improves the constraints on these components and the spectral slope. For power-law models, WMAP data alone puts an improved upper limit on the tensor-to-scalar ratio, r 0:002 < 0:65 (95% CL) and the combination of WMAP and the lensing-normalized SDSS galaxy survey implies r 0:002 < 0:30 (95% CL). Models that suppress largescale power through a running spectral index or a large-scale cutoff in the power spectrum are a better fit to the WMAP and small-scale CMB data than the power-law ÃCDM model; however, the improvement in the fit to the WMAP data is only Á 2 ¼ 3 for 1 extra degree of freedom. Models with a running-spectral index are consistent with a higher amplitude of gravity waves. In a flat universe, the combination of WMAP and the Supernova Legacy Survey (SNLS) data yields a significant constraint on the equation of state of the dark energy, w ¼ À0:967 þ0:073 À0:072 . If we assume w ¼ À1, then the deviations from the critical density, K , are small: the combination of WMAP and the SNLS data implies k ¼ À0:011 AE 0:012. The combination of WMAP 3 year data plus the HST Key Project constraint on H 0 implies k ¼ À0:014 AE 0:017 and à ¼ 0:716 AE 0:055. Even if we do not include the prior that the universe is flat, by combining WMAP, large-scale structure, and supernova data, we can still put a strong constraint on the dark energy equation of state, w ¼ À1:08 AE 0:12. For a flat universe, the combination of WMAP and other astronomical data yield a constraint on the sum of the neutrino masses, P m < 0:66 eV (95%CL). Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps using Minkowski functionals, the bispectrum, trispectrum, and a new statistic designed to detect large-scale anisotropies in the fluctuations. Subject headingg s: cosmic microwave background -cosmology: observations