We introduce a new statistic to measure more accurately the cosmic sound speed of clusters of galaxies at different redshifts. This statistic is evaluated by cross-correlating cosmic microwave background (CMB) fluctuations caused by the Sunyaev-Zel'dovich effect from observed clusters of galaxies with their redshifts. When clusters are distributed in redshift bins of narrow width, one could measure the mean squared cluster peculiar velocity with an error \sigma_{C_S^2}\lsim (300{\rm km/s})^2. This can be done around z>0.3 with clusters of flux above 200 mJy which will be detected by PLANCK, coupled with high resolution microwave images to eliminate the cosmological part of the CMB fluctuations. The latter can be achieved with observations by the planned ALMA array or the NSF South Pole telescope and other surveys. By measuring the cosmic sound speed and the bulk flow in, e.g., 4 spheres of ~ 100h^{-1}Mpc at z=0.3, we could have a direct measurement of the matter density 0.21<\Omega_m<0.47 at 95 % confidence level.

We present a new method for time-efficient and accurate extraction of the power spectrum from future cosmic microwave background (CMB) maps based on properties of peaks and troughs of the Gaussian CMB sky. We construct a statistic describing their angular clustering - analogously to galaxies, the 2-point angular correlation function, $\xi_\nu(\theta)$. We show that for increasing peak threshold, $\nu$, the $\xi_\nu(\theta)$ is strongly amplified and becomes measurable for $\nu\geq$1 on angular scales $\leq 10^\circ$. Its amplitude at every scale depends uniquely on the CMB temperature correlation function, $C(\theta)$, and thus the measured $\xi_\nu$ can be uniquely inverted to obtain $C(\theta)$ and its Legendre transform, the power spectrum of the CMB field. Because in this method the CMB power spectrum is deduced from high peaks/troughs of the CMB field, the procedure takes only $[f(\nu)]^2N^2$ operations where $f(\nu)$ is the fraction of pixels with $|\delta T|\geq\nu$ standard deviations in the map of $N$ pixels and is e.g. 0.045 and 0.01 for $\nu$=2 and 2.5 respectively. We develop theoretical formalism for the method and show with detailed simulations, using MAP mission parameters, that this method allows to determine very accurately the CMB power spectrum from the upcoming CMB maps in only $\sim(10^{-4}-10^{-3})\times N^2$ operations.