We study the asymptotic behaviour of needlets-based approximate maximum likelihood estimators for the spectral parameters of Gaussian and isotropic spherical random fields. We prove consistency and asymptotic Gaussianity, in the high-frequency limit, thus generalizing earlier results by Durastanti et al. (2011) based upon standard Fourier analysis on the sphere. The asymptotic results are then illustrated by an extensive Monte Carlo study.

We investigate local variations of the primordial non-Gaussianity parameter f_NL in the WMAP data, looking for possible influence of foreground contamination in the full-sky estimate of f_NL. We first improve the needlet bispectrum estimate in (Rudjord et al. 2009) on the full-sky to f_NL= 73 +/- 31 using the KQ75 mask on the co-added V+W channel. We find no particular values of f_NL estimates close to the galactic plane and conclude that foregrounds are unlikely to affect the estimate of f_NL in the V and W bands even for the smaller KQ85 mask. In the Q band however, we find unexpectedly high values of f_NL in local estimates close to the galactic mask, as well as significant discrepancies between Q band estimates and V/W band estimates. We therefore conclude that the Q band is too contaminated to be used for non-Gaussianity studies even with the larger KQ75 mask. We further noted that the local f_NL estimates on the V+W channel are positive on all equatorial bands from the north to the south pole. The probability for this to happen in a universe with f_NL = 0 is less than one percent.

We consider the statistical analysis of random sections of a spin fibre bundle over the sphere. These may be thought of as random fields that at each point p in $S^2$ take as a value a curve (e.g. an ellipse) living in the tangent plane at that point $T_{p}S^2$, rather than a number as in ordinary situations. The analysis of such fields is strongly motivated by applications, for instance polarization experiments in Cosmology. To investigate such fields, spin needlets were recently introduced by Geller and Marinucci (2008) and Geller et al. (2008). We consider the use of spin needlets for spin angular power spectrum estimation, in the presence of noise and missing observations, and we provide Central Limit Theorem results, in the high frequency sense; we discuss also tests for bias and asymmetries with an asymptotic justification.

We use the full bispectrum of spherical needlets applied to the WMAP data of the cosmic microwave background as an estimator for the primordial non-Gaussianity parameter f_NL. We use needlet scales up to l_max=1000 and the KQ75 galactic cut and find f_NL=84 +/- 40 corrected for point source bias. We also introduce a set of consistency tests to validate our results against the possible influence of foreground residuals or systematic errors. In particular, fluctuations in the value of f_NL obtained from different frequency channels, different masks and different multipoles are tested against simulated maps. All variations in f_NL estimates are found statistically consistent with simulations.

We consider the correlation structure of the random coefficients for a wide class of wavelet systems on the sphere (Mexican needlets) which were recently introduced in the literature by Geller and Mayeli (2007). We provide necessary and sufficient conditions for these coefficients to be asymptotic uncorrelated in the real and in the frequency domain. Here, the asymptotic theory is developed in the high resolution sense. Statistical applications are also discussed, in particular with reference to the analysis of cosmological data.

The purpose of this paper is to join two different threads of the recent literature on random fields on the sphere, namely the statistical analysis of higher order angular power spectra on one hand, and the construction of second-generation wavelets on the sphere on the other. To this aim, we introduce the needlets bispectrum and we derive a number of convergence results. Here, the limit theory is developed in the high resolution sense. The leading motivation of these results is the need for statistical procedures for searching non-Gaussianity in Cosmic Microwave Background radiation.