WMAP 
NineYear Wilkinson Microwave Anisotropy Probe (WMAP) Observations:

WMAP Nineyear Paper Figures, C. Bennett, et al.  
Individual figures are provided for use in talks. Proper display of PNG transparency in PowerPoint requires saving files to your computer before Inserting them. Please acknowledge the WMAP Science Team when using these images. Image Credit: NASA / WMAP Science Team 

Fig.1 The estimated level of artifacts (ξ) that would have occurred in the Kaband map if no processing mask had been used. Banddependent processing masks were used and tailored to minimize these artifacts when converting from timeordered to sky map data. This map is in Galactic coordinates and the high intensity regions arise from observations when one of the beams is near the Galactic center and the processing mask is not used. (See Figure 17 to compare with the analysis sky cuts.) Since bright artifacts originate primarily from beam crossings of bright Galactic plane regions, the nature of the unmasked artifact pattern is similar for all DAs. Although the patterns are similar for all bands, the highest amplitude artifacts occur in K and Ka bands because these have the brightest foregrounds. To prevent significant artifacts, processing masks are constructed for each band by growing the number of pixels in the mask until ξ is sufficiently reduced. The estimated mean residual level of artifacts (ξ) is given in Table 4. We required ξ < 5 μK for all but Kband. Construction of the Kband mask is more complex (see text) yet still achieves ξ < 8 μK. 
PNG (70 kb 1040x664 pixels) PNG (91 kb 2048x1328 pixels) EPS (129 kb) EPSCMYK (708 kb) 

Fig.2 Plots of the maximum and mean magnitude of the estimated map artifacts (ξ) for Kaband versus the number of pixels masked by the processing mask. The vertical line indicates the adopted mask which is the smallest mask for which max(ξ) < 1.15 ξ^{max} _{sat} as described in the text. 
PNG (38 kb 2520x1940 pixels) PNG (115 kb 5040x3880 pixels) EPS (29 kb) EPSCMYK (338 kb) 

Fig.3 Nineyear temperature sky maps in Galactic coordinates shown in a Mollweide projection. Maps have been slightly smoothed with a 0.°2 Gaussian beam. 
PNG (1.2 mb 1040x937 pixels) PNG (4.4 MB 2080x1874 pixels) EPS (921 kb) EPSCMYK (10.5 mb) 

Fig.4 Nineyear Stokes Q polarization sky maps in Galactic coordinates shown in a Mollweide projection. Maps have been smoothed to an effective Gaussian beam of 2.°0. The smooth large angular scale features visible in Wband, and to a lesser extent in Vband, are the result of a pair of modes that are poorly constrained in mapmaking, yet properly deweighted in the analysis. 
PNG (1 mb 1040x937 pixels) PNG (3 MB 2080x1874 pixels) EPS (516 kb) EPSCMYK (11.3 mb) 

Fig.5 Nineyear Stokes U polarization sky maps in Galactic coordinates shown in a Mollweide projection. Maps have been smoothed to an effective Gaussian beam of 2.°0. The smooth large angular scale features visible in Wband, and to a lesser extent in Vband, are the result of a pair of modes that are poorly constrained in mapmaking, yet properly deweighted in the analysis. 
PNG (1 mb 1040x937 pixels) PNG (3.1 MB 2080x1874 pixels) EPS (521 kb) EPSCMYK (11.8 mb) 

Fig.6 Nineyear polarized intensity (P) sky maps in Galactic coordinates shown in a Mollweide projection; P = (Q^{2} + U^{2})^{0.5}, where Q and U are Stokes parameters. Maps have been smoothed to an effective Gaussian beam of 2.°0. Plotted line segments show polarization angles for HEALPix nside = 16 pixels where the signaltonoise exceeds 1. 
PNG (1.1 mb 1040x937 pixels) PNG (3.4 MB 2080x1874 pixels) EPS (1.2 mb) EPSCMYK (12.4 mb) 

Fig.7 Kband images of supernova remnant Tau A (3C 144), at [J2000.0] position (05^{h}34^{m}31^{s}, 22°01') from the first year of WMAP observations. The left panels display the total intensity and right the residuals after removal of a best fit circularly symmetric beam profile. The maps generated with the new partial deconvolution processing (bottom) display significantly reduced beam asymmetry compared with those generated with the standard processing (top). In other words, the apparent asymmetry in Tau A seen in the top left is a result of the asymmetric Kband beam and is not intrinsic to Tau A. The degree of a source's apparent asymmetry is dependent on its sky position and the WMAP frequency at which it is observed: the effect is most pronounced for bright Kband sources at low ecliptic latitudes (Section 4.2). As such, we display the Kband observations of Tau A to demonstrate the effectiveness of the deconvolution in a worst case of beam asymmetry in the normally processed maps. 
PNG (215 kb 1423x1739 pixels) PNG (317 kb 2845x3477 pixels) EPS (184 kb) EPSCMYK (1.2 mb) 

Fig.8 Plots illustrating the formation of the kernel used to generate the symmetrized beam maps for the Q1 DA. The x and y axes are in units of degrees centered on the beam LOS. The zaxis represents weight and panels (a), (e) and (f) use the same scale. (a) The residual (nonaxisymmetric) component of the beam obtained by subtracting the best fit axisymmetric beam from the total beam map. (b) The residual beam after Wiener deconvolution. (c) The cosine apodization function. (d) The convolution kernel used to generate the symmetrized beam maps consisting of the cosine weighted Wiener deconvolved residual map. (e) The convolution kernel reconvolved with the axisymmetric beam. (f) The difference between the residual beam map (a) and the map making kernel convolved with the axisymmetric beam (e). 
PNG (155 kb 2048x1536 pixels) PNG (276 kb 4096x3072 pixels) EPS (301 kb) EPSCMYK (866 kb) 

Fig.9 Verification of effects of asymmetric beams on the power spectrum. Given beam measurements, the formalism in Appendix B of Hinshaw et al. (2007) analytically quantifies the beam asymmetry effect on the power spectrum. This is plotted as a fractional deviation between an ideally deconvolved power spectrum (C^{deconv}_{l} ) and the power spectrum of a normally processed map (C^{np}_{l} ) with no correction for beam asymmetries. These "predictions" of fractional deviations are plotted per DA in the light colored solid lines. The Qband effects become significant at l ~ 400, but Qband is not used in the WMAP cosmological power spectrum. Vband effects become significant at l ~ 1000, however, Vband is deweighted compared to Wband at high l because of its larger beam size. Wband effects from the asymmetric beams can be seen to be . ≤1%. While Hinshaw et al. (2007) provides an analytic prediction, we have explicitly deconvolved the maps in pixel space, allowing for a direct intercomparison of the analytic with the numerical approach. The dark red, green and blue solid lines are the fractional deviations in power spectra for Q, V and Wbands from the directly deconvolved maps. A comparison between the light and dark colored lines per frequency band shows close agreement up to a multipole moment where we expect the spectra derived from the beamsymmetrized maps to break down because the prediction does not account for correlations introduced by the deconvolution. The Qband deviations occur after the window function has dropped below 2.5% and the Vband deviations below 1.5%. The vertical dashed lines indicate where window functions are at 1% of their maximum value. The close agreement between the predictions and explicit deconvolution verifies our understanding of asymmetric beam effects and allows us to conclude that the spectrum from the normally processed (i.e. not deconvolved) maps differs from the ideallydeconvolved spectrum by < 1%. Thus the final WMAP power spectrum is based on the normallyprocessed V and W band maps. 
PNG (47 kb 2048x1503 pixels) PNG (84 kb 4096x3007 pixels) EPS (45 kb) EPSCMYK (369 kb) 

Fig.10 False color image representing the spectral information from multiple WMAP bands. Qband is red, Vband is green, and Wband is blue. In this representation, a CMB thermodynamic spectrum appears as grey. 
PNG (537 kb 1033x529 pixels) PNG (1.6 mb 2066x1058 pixels) EPS (369 kb) EPSCMYK (11.6 mb) 

Fig.11 False color image derived from a combination of WMAP band differences chosen to highlight differing spectral components. Red (WV) highlights regions where thermal emission from dust is highest. Blue (QW) is dominated by freefree emission. Green ((KKa) 1.7(QW)) illustrates contributions from synchrotron and spinning dust. 
PNG (207 kb 514x259 pixels) PNG (698 kb 1029x519 pixels) EPS (214 kb) EPSCMYK (12.9 mb) 

Fig.12 Microwave emission near the Galactic plane is traced by a Kband minus Wband difference map, which eliminates CMB anisotropy. A log scale is used for the color region and blue circles represent the positions of the brightest point sources, as seen by WMAP. 
PNG (169 kb 2052x1046 pixels) PNG (355 kb 4100x2091 pixels) EPS (684 kb) EPSCMYK (1.9 mb) 

Fig.13 Modeling results for Saturn. (Left) Brightness temperatures based on unresolved Saturn observations as a function of ring inclination B are shown in black for each WMAP frequency band. Where there are multiple differencing assemblies per frequency, multiple points are plotted at each inclination. An empirical model including both ring and disk components (see text) is plotted in red. The temperature of the planetary disk predicted by the model occurs at B=0°, when the rings are viewed edgeon. The model is symmetric about B=0°. (Right) Residuals (datamodel) of the model fit to the data are plotted as a function of the ring opening angle. 
PNG (89 kb 1640x2049 pixels) PNG (156 kb 3278x4096 pixels) EPS (142 kb) EPSCMYK (77 kb) 

Fig.14 Saturn model parameters derived from the nineyear analysis. Left: Disk temperatures for 5 WMAP frequencies. Right: Ring system temperatures. Adopted errors for the nineyear analysis have been reduced compared to those in Weiland et al. (2011); errors for T_{disk} are smaller by a factor of 2. 
PNG (29 kb 2048x883 pixels) PNG (62 kb 4096x1766 pixels) EPS (66 kb) EPSCMYK (477 kb) 

Fig.15 Foreground evaluation is generally based on a combination of the data from the five WMAP bands and external observations where the CMB contamination is negligible. The external observations used for foreground fitting and template removal are shown. These provide approximate probes of the synchrotron, freefree, spinning dust, and thermal dust emission. 
PNG (551 kb 1024x812 pixels) PNG (1.7 mb 2048x1623 pixels) EPS (475 kb) EPSCMYK (6.4 mb) 

Fig.16 Polarization template coefficients, scaled to produce model maps in antenna temperature, as a function of frequency. The curves show the predictions of a simple model with synchrotron and thermal dust polarization in which about 2/3 of the dust polarization is traced by the dust template and about 1/3 is traced by the Kband template. 
PNG (28 kb 2048x1463 pixels) PNG (66 kb 4096x2926 pixels) EPS (31 kb) EPSCMYK (338 kb) 

Fig.17 Comparison of nineyear masks to sevenyear masks. At the top KQ75y7 and KQ75y9 are compared, and at the bottom KQ85y7 and KQ85y9. Green regions are masked in both the nineyear and sevenyear masks, yellow regions are newly masked in the nineyear masks, and red regions are masked in the sevenyear masks but no longer in the nineyear masks. 
PNG (148 kb 2108x2344 pixels) PNG (234 kb 4216x4688 pixels) EPS (243 kb) EPSCMYK (1.3 mb) 

Fig.18 The top panel shows spinning dust emissivity spectra predicted by the model of AliHaimoud et al. (2009) and Silsbee et al. (2011) for the nominal physical conditions that they adopted for different ISM environments  cold neutral medium (CNM), warm neutral medium (WNM), warm ionized medium (WIM), molecular cloud (MC), photodissociation region (PDR), reflection nebula (RN), and dark cloud (DC). The spectra were calculated using version 2.01 of the code SpDust provided by the authors, for the case where dust grains are allowed to rotate around nonprincipal axes. The spectra are in units of brightness temperature per H column density. The bottom panel shows the same spectra normalized to a peak of unity and scaled to a common peak frequency (that of the CNM spectrum, 17.8 GHz). The predicted spectral shapes for the different environments are similar. We adopted the CNM case for the shape of the spinning dust spectrum in our foreground fitting. We used this as an externally provided spectral template in our fits, usually with our own arbitrary amplitude and frequency scaling. The fit results in no way imply that the underlying physical mechanisms or the lineofsite conditions have been established. 
PNG (122 kb 2048x2948 pixels) PNG (251 kb 4096x5897 pixels) EPS (82 kb) EPSCMYK (515 kb) 

Fig.19 Parameter maps from the MEM model fit. The top four maps are shown on logarithmic scales and the others are on linear scales. 
PNG (848 kb 1024x1226 pixels) PNG (2.9 mb 2048x2453 pixels) EPS (483 kb) EPSCMYK (8.3 mb) 

Fig.20 Parameter maps from the MCMCg model fit. The top four maps are shown on logarithmic scales and the others are on linear scales. Accurate determination of the CMB close to the Galactic plane is inhibited by CMBforeground covariance. The map for β synchrotron is evaluated at 40 GHz. 
PNG (1.7 mb 1024x1653 pixels) PNG (5.6 mb 2048x3306 pixels) EPS (992 kb) EPSCMYK (13 mb) 

Fig.21 Parameter maps from the Model 9 fit. The top four maps are shown on logarithmic scales and the others are on linear scales. The map for β synchrotron is evaluated at 40 GHz. 
PNG (1.4 mb 1024x1260 pixels) PNG (4.4 mb 2048x2521 pixels) EPS (801 kb) EPSCMYK (10.2 mb) 

Fig.22 Spectra of CMB and foreground anisotropy. The foreground anisotropy results are averages over the three foreground models (MCMCg, MEM, and Model 9). The upper curve for each foreground component shows results for pixels outside of the KQ85 mask, and the lower curve shows results outside of the KQ75 mask. The different foreground models are in good agreement for the total foreground anisotropy. Results for the individual foreground components depend on model assumptions discussed in the text, and typically differ among the three models by 5% to 25%. 
PNG (86 kb 2048x1755 pixels) PNG (240 kb 4096x3511 pixels) EPS (75 kb) EPSCMYK (444 kb) 

Fig.23 Results from foreground degeneracy analysis for sixband Model 9 fitting. The contour plots illustrate the degeneracy between model parameters for a representative single pixel foreground spectrum. Each panel shows the change in Χ^{2} as the selected pair of parameters are varied from their bestfit values while marginalizing over the other parameters. Contours are shown for ΔΧ^{2} values of 0.2, 1, 3, and 10, except values of 0.5, 3, and 10 are used for β_{sync} vs. synchrotron amplitude. There are significant degeneracies between parameter pairs that include either synchrotron amplitude or synchrotron spectral index, except for those that include thermal dust amplitude. 
PNG (85 kb 2048x1888 pixels) PNG (189 kb 4096x3775 pixels) EPS (144 kb) EPSCMYK (741 kb) 

Fig.24 (Left): Thermal dust amplitude at Wband averaged over the MCMCg, MEM and Model 9 fits minus the thermal dust model 8 from Finkbeiner et al. (1999). (Right): Freefree amplitude at Kband averaged over the same three models, minus the freefree template estimated from H_{α} observations. 
PNG (410 kb 1024x395 pixels) PNG (1.4 mb 2048x790 pixels) EPS (278 kb) EPSCMYK (3.3 mb) 

Fig.25 The ratio of Wband predicted thermal dust emission (Finkbeiner et al. (1999) model 8) to the mean over three models (MCMCg, MEM, Model 9) as a function of longitude for b < 5°. Error bars are derived from the rms scatter of the three models about the mean. A line is a plotted at 1.0 to guide the eye. Modeled emission shows systematic variations from the FDS prediction by up to 20%. 
PNG (27 kb 2048x1528 pixels) PNG (63 kb 4096x3056 pixels) EPS (28 kb) EPSCMYK (328 kb) 

Fig.26 Estimates of foreground bias error remaining in the ILC map, on a scale of ±15 μK. Top left: Bias map from the threeyear analysis of Hinshaw et al. (2007). The map is zeroed outside the Kp2 cut. Top right and middle: Bias estimates resulting from the application of the nineyear ILC coefficients to the Galaxy models from MEM, Model 9 and MCMCg analysis. The bias map from the MCMCg analysis is overestimated in the plane (see text). Bottom left: ILC error from foregroundCMB covariance. Within the Kp2 cut, this error and the foreground bias are of comparable magnitude. Bottom right: An estimate of the potential magnitude of ILC foreground bias outside the Kp2 cut, based on the various model results, with heavy weight given to the MCMCg model. Bias errors of 10μK or less are indicated. 
PNG (1.1 mb 1024x992 pixels) PNG (3.8 mb 2048x1984 pixels) EPS (508 kb) EPSCMYK (8.9 mb) 

Fig.27 The top map is the nineyear ILC. The bottom sky map displays the part of the ILC error in each pixel due to foregroundCMB covariance, using the Model 9 foreground estimate from Section 5.3.6. This shows the square root of the diagonal of the covariance matrix, on a linear color scale. Therefore it shows the standard deviation of expected error fluctuations, marginalizing over correlations between pixels. The color scale range was chosen because the r6 ILC map has a CMB standard deviation of 66 μK. Thus, full scale on this map has equal variance with the CMB, and at the halfway point on this color scale the foregroundCMB error variance is down to a quarter of the CMB variance. 
PNG (694 kb 746x1024 pixels) PNG (2.1 mb 1492x2048 pixels) EPS (435 kb) EPSCMYK (4.9 mb) 

Fig.28 The dominant power law in a pixel, combined with information about whether the data in that pixel look like a pure power law, over the WMAP bands. This image was generated by individually specifying the hue, saturation, and value (HSV) for each pixel. The hue, shown in the color scale, describes which power law best fits the data. It is labeled with values of β, where the power law in antenna temperature is T_{A}(ν) ∝ ν^{β}. The saturation describes how well the data fit a power law, so that desaturated (white, gray, black) pixels are not well fit by any power law. Specifically, let n_{A} be a 5vector of the WMAP thermodynamic temperatures, rescaled to be a unit vector, and let n_{p} be a 5vector of the best fit power law in antenna temperature, converted to thermodynamic and then also rescaled to be a unit vector. Then the saturation is n_{A} · n_{p}, which is just the cosine of the angle between these two vectors. The scale is from 0.995 (unsaturated) to 1.0 (completely saturated), so if the two 5vectors are more than 5.73 degrees apart, the pixel is unsaturated. The value in the HSV color space is the magnitude of the data 5vector, so it is the square root of the sum of the squares of the WMAP thermodynamic temperatures, on a scale of 0 to 2 mK. Therefore blacker pixels have less emission in all bands; lighter pixels have more emission. The nineyear ILC was subtracted from the WMAP data, before computing the above image. 
PNG (659 kb 1050x667 pixels) PNG (2.1 mb 2100x1333 pixels) EPS (197 kb) EPSCMYK (10.6 mb) 

Fig.29 Endtoend Monte Carlo pipeline tests. The gray lines are individual l's and the black lines are boxcar smoothed with Δl = 50. In all four cases, the ratio of the Monte Carlo estimated power spectrum and the predicted value is consistent with unity. Top left. Ratio ⟨Ĉ^{×}_{l}⟩_{sig}/C^{fid}_{l} between mean estimated power spectrum of CMBonly simulations and the fiducial input spectrum. Top right. Same as top left panel, but using the autocorrelation estimator Ĉ_{l} instead of the noauto estimator Ĉ^{×}_{l} . Bottom left. Ratio between mean estimated power spectrum of noiseonly simulations and the predicted noise bias, using the autoestimator Ĉ_{l}. Bottom right. Ratio between mean estimated power spectrum of point source simulations and predicted bias. 
PNG (78 kb 2048x1394 pixels) PNG (169 kb 4096x2788 pixels) EPS (162 kb) EPSCMYK (784 kb) 

Fig.30 Binned WMAP7 power spectrum estimates using the optimal pipeline from this paper (left/black error bars), with the estimates from the WMAP7 release (Larson et al. 2011) shown for comparison (right/grey error bars). 
PNG (25 kb 2048x1503 pixels) PNG (61 kb 4096x3007 pixels) EPS (32 kb) EPSCMYK (334 kb) 

Fig.31 Detailed comparison between WMAP7 optimal power spectrum estimator and suboptimal estimator from Larson et al. (2011). Top: Difference (Ĉ^{optimal} _{l}  Ĉ^{subopt} _{l})/Var(Ĉ^{optimal}_{l})^{1/2} between the two estimators in "sigmas", for every l, and boxcarsmoothed with Δl = 10. Bottom: Variance ratio between suboptimal and optimal estimators. 
PNG (85 kb 2048x2779 pixels) PNG (175 kb 4096x5559 pixels) EPS (69 kb) EPSCMYK (472 kb) 

Fig.32 The nineyear WMAP TT angular power spectrum. The WMAP data are in black, with error bars, the best fit model is the red curve, and the smoothed binned cosmic variance curve is the shaded region. The first three acoustic peaks are welldetermined. 
PNG (45 kb 2048x1406 pixels) PNG (99 kb 4096x2812 pixels) EPS (38 kb) EPSCMYK (343 kb) 

Fig.33 The TE spectrum. The WMAP data points and error bars are in black. The red theory curve is fit to the full WMAP data, including the TT angular power spectrum data. Note that the vertical axis on these spectra is (l + 1)C_{l}/(2π) instead of l(l + 1)C_{l}/(2π); this vertical scale differs from that of the TT spectrum plot by a factor of l. The lowest l TE bin where 2 ≤ l ≤ 7 has been adjusted using a pixel likelihood code. 
PNG (28 kb 2048x1115 pixels) PNG (57 kb 4096x2230 pixels) EPS (35 kb) EPSCMYK (344 kb) 

Fig.34 The TB spectrum. The TB spectrum uses the MASTER likelihood code. Note that the vertical axis on these spectra is (l+1)C_{l}/(2π) instead of l(l+1)C_{l}/(2π); this vertical scale differs from that of the TT spectrum plot by a factor of l. 
PNG (24 kb 2048x1146 pixels) PNG (54 kb 4096x2292 pixels) EPS (28 kb) EPSCMYK (326 kb) 

Fig.35 Individual likelihood functions of the low l EE polarized power are shown for l = 2 through 7. When fitting at a particular l, we set C_{l} at all other values of l to the value in the best fit WMAP power spectrum. In addition, at the l in question we set C^{TE} _{l} = 0 to maintain that C^{TE}_{l} ≤ √C^{TT}_{l} C^{EE} _{l} . The black diamonds denote the best fit WMAP EE power spectrum. These likelihood functions include sample variance. 
PNG (52 kb 2048x2064 pixels) PNG (117 kb 4096x4129 pixels) EPS (52 kb) EPSCMYK (410 kb) 

Fig.36 Low l BB spectra. Other C_{l} values are fixed to the best fit WMAP power spectrum. 
PNG (47 kb 2048x2064 pixels) PNG (109 kb 4096x4129 pixels) EPS (52 kb) EPSCMYK (407 kb) 

Fig.37 Top: Even excess ε_{l} in the observed power spectrum, in bins of Δl = 50, compared to the mean and scatter from 512 Monte Carlo realizations. Bottom: ε_{l} as in the top plot, converted to significance units by normalizing to the Monte Carlo scatter in each bin. Only the l = 250299 and l = 300349 bins show a significance greater than 1σ. Black: nineyear results; blue: sevenyear results from Bennett et al. (2011). 
PNG (53 kb 2048x2015 pixels) PNG (104 kb 4096x4031 pixels) EPS (64 kb) EPSCMYK (437 kb) 

Fig.38 The likelihood of the true value of l(l + 1)C^{TT} _{l} /(2π) = 6C^{TT} _{2} /(2π) for l = 2, based on our measured sky. This is computed using the BlackwellRao estimator run on Gibbs samples, and it marginalizes over all other values of C^{TT} _{l} . The maximum likelihood point is shown as the pink line; one and two sigma regions are shown as blue and green lines. The best fit CDM theory spectrum computed on WMAP nineyear data only is shown in red. 
PNG (47 kb 2048x1406 pixels) PNG (94 kb 4096x2812 pixels) EPS (37 kb) EPSCMYK (346 kb) 

Fig.39 The cosmic variance probability distribution for the quadrupole, given the theory power spectrum. This assumes we know l(l + 1)C^{TT} _{l} /(2π) = 6C^{TT} _{2} /(2π) = 1109 μK^{2} (red line) and plots the distribution of quadrupole power values we could measure for random Hubble volumes. Note that 6C^{TT} _{2} /(2π) is the mean of the distribution; due to the skewness of the χ^{2} distribution, the peak of the distribution is substantially lower. One and twosigma regions are shown. The quadrupole cosmic variance distribution has ν = 2l + 1 = 5 degrees of freedom. Assuming f_{sky} ≈ 0.99, we plot a χ^{2} distribution based on ν = (2l + 1)f^{2} _{sky} ≈ 4.9 degrees of freedom. The peak of the distribution is then lower than the mean by a factor of (ν  2)/ν, putting it at 656 μK^{2}. 
PNG (44 kb 2048x1406 pixels) PNG (91 kb 4096x2812 pixels) EPS (32 kb) EPSCMYK (335 kb) 

Fig.40 A test for scaledependent systematics: f^{orth}NL estimates in lbands, with cumulative bestfit value f^{orth} NL = 245 shown by the dotted horizontal line. Each error bar is labeled with the statistical significance of the deviation from the cumulative bestfit value (not the deviation from zero). No evidence for scaledependent systematics is seen. 
PNG (38 kb 2048x1503 pixels) PNG (77 kb 4096x3007 pixels) EPS (38 kb) EPSCMYK (352 kb) 

Fig.41 A visual test for sky location dependent systematics: skymap showing the contribution of different parts of the sky to the f^{orth}_{NL} estimator, in units of "f^{orth}_{NL} per steradian". We do not detect any significant localized features in this map. 
PNG (372 kb 1294x811 pixels) PNG (505 kb 2588x1623 pixels) EPS (174 kb) EPSCMYK (1.2 mb) 

Fig.42 WMAP nineyear constraints on nonGaussianity in singlefield inflation. Upper panel. 68%, 95%, and 99.7% confidence regions in the f^{eq}_{NL}, f^{orth}_{NL} plane, defined by threshold χ^{2} values 2.28, 5.99, 11.62, as appropriate for a χ^{2} random variable with two degrees of freedom. (f^{eq}NL, f^{orth}NL ) = (0, 0) is consistent with the data to within 99% CL. Lower panel. Confidence regions on the dimensionless sound speed c_{s} and interaction coefficient A (defined in Equation (56)), obtained from the top panel via the change of variables in Equation (57). The upper bound on f^{eq}NL gives a lower bound on c_{s}, which is consistent with c_{s} = 1. 
PNG (22 kb 2048x1320 pixels) PNG (56 kb 4096x2640 pixels) EPS (31 kb) EPSCMYK (384 kb) 

Fig.43 Constraints on curvature. Flat universes fall on the Ω_{m} + Ω_{Λ} = 1 line. Allowed regions are shown for WMAP, CMB, and CMB combined with BAO and H_{0} data, all with a hard prior of H_{0} < 100 kms^{1}Mpc^{1}. WMAP data is represented by 290,000 Markov chain points, colored by their value of H_{0}. The WMAP data follow a geometric degeneracy ridge represented by the slightly curved line, a parabola with equation Ω_{Λ}= 0.0620Ω^{2}_{m}  0.825Ω_{m} + 0.947. The most likely point in the WMAPonly chain has Ω_{Λ} = 0.721 and Ω_{m} = 0.279, which is flat to three significant figures, even though this constraint was not enforced. The WMAP data alone require Ω_{Λ} > 0.58 at 68% CL and Ω_{Λ} > 0.22 at 95% CL. The contours show constraints when adding highl CMB data (blue) and BAO and H_{0} data (black). These constraints are consistent with those from WMAP alone, with the tightest constraint being Ω_{tot} = 1.0027^{+0.0038}_{0.0039} (Hinshaw et al. 2012). 
PNG (100 kb 1024x760 pixels) PNG (288 kb 2048x1520 pixels) EPS (194 kb) EPSCMYK (1.9 mb) 

Fig.44 Top  Measurements of the yeartoyear fractional brightness variation of the Galactic plane in WMAP skymaps, obtained by correlating Galactic plane signal in each single year map with Galactic plane signal in the nineyear map. There is a small dependence of these variations on spectral index, which shows that they are caused by variations in effective WMAP band center frequencies over the mission. Bottom  The yeartoyear fractional variation of WMAP band center frequency derived from Galactic plane brightness variations measured for selected spectral index bins. 
PNG (81 kb 2048x2475 pixels) PNG (151 kb 4096x4951 pixels) EPS (17 kb) EPSCMYK (414) 