Helical distributed chaos in magnetic field has been studied using results of direct numerical simulations (dominated by magnetic helicity), of a laboratory experiment with plasma wind tunnel and of solar wind measurements (dominated by combined magnetic and cross helicity effects). The solar wind measurements, used for the spectral computations, were produced by Helios-1 and Ulysses missions for low and high heliolatitudes respectively and for high solar wind speed at $0.4 < R <4.5$ AU (where $R$ is distance from the Sun).

It is shown that the long-term solar activity (represented by the dynamically covered monthly time series of the sunspot number for period 1750-2005yy) exhibits spectral properties of the Hamiltonian distributed chaos with spontaneously broken time translational symmetry.

It is shown, by the means of distributed chaos approach and using the experimental data, that at very large Rayleigh number $Ra > 10^{14}$ and Prandtl number $Pr \sim 1$ the Rayleigh-B\'{e}nard turbulence can undergo a transition related to spontaneous breaking of the fundamental Lagrangian relabeling symmetry. Due to the Noether's theorem helicity plays central role in this process. After the transition the temperature spectrum has a stretched exponential form $E (k) \propto \exp(-k/k_{\beta})^{\beta}$ with $\beta =2/5$ both at the cell midplain and at the near-wall (low boundary) regions. There is a similarity between this phenomenon and the effects of polymer additives.

It is shown, using results of recent direct numerical simulations, laboratory experiments and atmospheric measurements, that buoyancy driven turbulence exhibits a broad diversity of the types of distributed chaos with its stretched exponential spectrum $\exp(-k/k_{\beta})^{\beta}$. The distributed chaos with $\beta = 1/3$ (determined by the helicity correlation integral) is the most common feature of the stably stratified turbulence (due to the strong helical waves presence). These waves mostly dominate spectral properties of the vertical component of velocity field, while the horizontal component is dominated by the diffusive processes both for the weak and strong stable stratification ($\beta =2/3$). For the last case influence of the low boundary can overcome the wave effects and result in $\beta =1/2$ for the vertical component of the velocity field (the spontaneous breaking of the space translational symmetry - homogeneity). For the unstably stratified turbulence in the Rayleigh-Taylor mixing zone the diffusive processes ($\beta =2/3$) are the most common dominating processes in the anisotropic chaotic mixing of the two fluids under buoyancy forces. The distributed chaos in Rayleigh-B\'{e}nard turbulent convection in an upright cell is determined by the strong confinement conditions. That is: the spontaneous breaking of the space translational symmetry (homogeneity) by the finite boundaries ($\beta = 1/2$) or by the non-perfect orientation of the cell along the buoyancy direction ($\beta =4/7$). In all types of turbulence appearance of an inertial range of scales results in deformation of the distributed chaos and $\beta =3/5$.

Two types of spontaneous breaking of the space translational symmetry in distributed chaos have been considered for turbulent thermal convection at large values of Rayleigh number. First type is related to boundaries and second type is related to appearance of inertial range of scales. The first type is dominated by vorticity correlation integral: $\int_{V} \langle {\boldsymbol \omega} ({\bf x},t) \cdot {\boldsymbol \omega} ({\bf x} + {\bf r},t) \rangle_{V} d{\bf r}$ and is characterized by stretched exponential spectrum $\exp-(k/k_{\beta})^{\beta }$ with $\beta =1/2$. The second type is dominated by energy correlation integral: $\int_{V} \langle {\bf u}^2 ({\bf x},t) ~ {\bf u}^2({\bf x} + {\bf r},t) \rangle_{V} d{\bf r}$ and is characterized by $\beta =3/5$. Good agreement has been established with laboratory experimental data obtained at large values of Rayleigh number $Ra \sim 10^{11}-10^{14}$ (the range relevant to solar photosphere) in upright cylinder cells. Taylor hypothesis transforms the wavenumber spectrum into frequency spectrum $\exp-(f/f_{\beta})^{1/2}$. It is shown that turnover frequency of large-scale circulation (wind): $f_w = f_{\beta}/2$. Results of an experiment in horizontal cylinder are also briefly discussed. The analysis suggests that in this case the large-scale circulation can be considered as a natural (harmonic) part of the distributed chaos.

It is shown that in turbulent flows the distributed chaos with spontaneously broken translational space symmetry (homogeneity) has a stretched exponential spectrum $\exp-(k/k_{\beta})^{\beta }$ with $\beta =1/2$. Good agreement has been established between the theory and the data of direct numerical simulations of isotropic homogeneous turbulence (energy dissipation rate field), of a channel flow (velocity field), of a fully developed boundary layer flow (velocity field), and the experimental data at the plasma edges of different fusion devices (stellarators and tokamaks). An astrophysical application to the large-scale galaxies distribution has been briefly discussed and good agreement with the data of recent Sloan Digital Sky Survey SDSS-III has been established.

It is shown, that the both angular CMB Doppler spectrum: $C_l$ (Planck space telescope - cosmic microwave background [1]) and the 3D galaxy-galaxy power spectrum: $P(k)$ (Sloan Digital Sky Survey SDSS-II [2]), exhibit a considerable range with an exponential decay: $370 < l < 2500$ and $0.05 < k < 0.27~~(h/Mpc)$, respectively. The rates of the exponential decay are $l_c \simeq 300$ for $C_l \sim \exp-(l/l_c)$ and $k_c \simeq 0.09~~(h/Mpc)$ for $P(k) \sim \exp-(k/k_c)$. A waviness is observed along a straight line representing the exponential decay in the log-linear scale graphs of these spectra. In both cases the waviness has period (distance between peaks) equal to the same $l_c$ and $k_c$ as for the exponential decay. It means that the waviness is generated by the same, presumably chaotic, mechanism that generates the exponential decay. A more complex, distributed, chaos is observed in the Baryon Oscillation Spectroscopic Survey (the largest component of the SDSS-III containing nearly one million galaxies) [3]. In this case the $P(k)$ spectrum is a weighted superposition of the exponentials. At assumption that dynamics of the dispersive waves, driving the distributed chaos, is dominated by the effects of a surface tension this weighted superposition of the exponentials is converged to a compact form of a stretched exponential $P(k) \sim \exp-(k/k_b)^{1/2}$, in good agreement with the data.

Universal role of the nonlinear one-third subharmonic resonance mechanism in generation of the strong fluctuations in such complex natural dynamical systems as global climate and global solar activity is discussed using wavelet regression detrended data. Role of the oceanic Rossby waves in the year-scale global temperature fluctuations and the nonlinear resonance contribution to the El Nino phenomenon have been discussed in detail. The large fluctuations of the reconstructed temperature on the millennial time-scales (Antarctic ice cores data for the past 400,000 years) are also shown to be dominated by the one-third subharmonic resonance, presumably related to Earth precession effect on the energy that the intertropical regions receive from the Sun. Effects of Galactic turbulence on the temperature fluctuations are discussed in this content. It is also shown that the one-third subharmonic resonance can be considered as a background for the 11-years solar cycle, and again the global (solar) rotation and chaotic propagating waves play significant role in this phenomenon. Finally, a multidecadal chaotic coherence between the detrended solar activity and global temperature has been briefly discussed.

The wavelet regression detrended fluctuations of the reconstructed temperature for the past three ice ages: approximately 340000 years (Antarctic ice cores isotopic data), exhibit clear evidences of the galactic turbulence modulation up to 2500 years time-scales. The observed strictly Kolmogorov turbulence features indicates the Kolmogorov nature of galactic turbulence, and provide explanation to random-like fluctuations of the global temperature on the millennial time scales.

It is shown that, the wavelet regression detrended fluctuations of the monthly sunspot number for 1749-2009 years exhibit strong periodicity with a period approximately equal to 3.7 years. The wavelet regression method detrends the data from the approximately 11-years period. Therefore, it is suggested that the one-third subharmonic resonance can be considered as a background for the 11-years solar cycle. It is also shown that the broad-band part of the wavelet regression detrended fluctuations spectrum exhibits an exponential decay that, together with the positive largest Lyapunov exponent, are the hallmarks of chaos. Using a complex-time analytic approach the rate of the exponential decay of the broad-band part of the spectrum has been theoretically related to the Carrington solar rotation period. Relation of the driving period of the subharmonic resonance (3.7-years) to the active longitude flip-flop phenomenon, in which the dominant part of the sunspot activity changes the longitude every 3.7 years on average, has been briefly discussed.

It is shown that correlation function of the mean wind velocity in a turbulent thermal convection (Rayleigh number $Ra \sim 10^{11}$) exhibits exponential decay with a very long correlation time, while corresponding largest Lyapunov exponent is certainly positive. These results together with the reconstructed phase portrait indicate presence of a chaotic component in the examined mean wind. Telegraph approximation is also used to study relative contribution of the chaotic and stochastic components to the mean wind fluctuations and an equilibrium between these components has been studied. Since solar activity is based on the thermal convection processes, it is reasoned that the observed solar activity long-term correlations can be an imprint of the mean wind chaotic properties. In particular, correlation function of the daily sunspots number exhibits exponential decay with a very long correlation time and corresponding largest Lyapunov exponent is certainly positive, also relative contribution of the chaotic and stochastic components follows the same pattern as for the convection mean wind.

Solar activity is studied using a cluster analysis of the time-fluctuations of sunspot number. In the historic period (1850-1932) the cluster exponent was approximately equal to 0.37 (strong clustering) for the high activity components of the solar cycles. In the modern period (last seven solar cycles: 1933-2007) the cluster exponent was approximately equal to 0.50 (random, white noise-like). Comparing these results with the corresponding data from laboratory experiments on convection it is shown, that in the historic period emergence of sunspots in the solar photosphere was dominated by turbulent photospheric convection. In the modern period, this domination was broken by a new more active dynamics of the inner layers of the convection zone. Cluster properties of the solar wind magnetic field and the aa-geomagnetic-index also support this result. Long-range chaotic dynamics in the solar activity is briefly discussed.

Solar activity is studied using a cluster analysis of the time-fluctuations of the sunspot number. It is shown that in an Historic period the high activity components of the solar cycles exhibit strong clustering, whereas in a Modern period (last seven solar cycles: 1933-2007) they exhibit a white-noise (non-)clustering behavior. Using this observation it is shown that in the Historic period, emergence of the sunspots in the solar photosphere was strongly dominated by turbulent photospheric convection. In the Modern period, this domination was broken by a new more active dynamics of the inner layers of the convection zone. Then, it is shown that the dramatic change of the sun dynamics at the transitional period (between the Historic and Modern periods, solar cycle 1933-1944yy) had a clear detectable impact on Earth climate. A scenario of a chain of transitions in the solar convective zone is suggested in order to explain the observations, and a forecast for the global warming is suggested on the basis of this scenario. A relation between the recent transitions and solar long-period chaotic dynamics has been found. Contribution of the galactic turbulence (due to galactic cosmic rays) has been discussed. These results are also considered in a content of chaotic climate dynamics at millennial timescales.

Coherent properties of the baryon-photon fluid decoupling are considered in the terms of an effective nonlinear Schr\"{o}dinger equation for a macroscopic wave function that specifies the index of the coherent state. Generation of a transitional acoustic turbulence preceding formation of large-scale condensate in the plasma and its influence on the CMB power spectrum has been studied. A scaling $k^{-1}$ law is derived for the CMB Doppler spectrum $E(k)$ (angle-averaged) in the {\it wavenumber} space, for sufficiently large wavenumber $k$ and for the weak nonlinear and completely disordered initial conditions. Using the recent WMAP data it is shown that the so-called first acoustic peak represents (in a compensated spectral form) a pre-condensate fraction of the spectrum $E(k)$ at a rather advance stage of the condensate formation process.

It is shown that nonlocal interactions determine energy spectrum in isotropic turbulence at small Reynolds numbers. It is also shown that for moderate Reynolds numbers the bottleneck effect is determined by the same nonlocal interactions. Role of the large and small scales covariance at the nonlocal interactions and in energy balance has been investigated. A possible hydrodynamic mechanism of the nonlocal solution instability at large scales has been briefly discussed. A quantitative relationship between effective strain of the nonlocal interactions and viscosity has been found. All results are supported by comparison with the data of experiments and numerical simulations.

An analytic perturbation theory is suggested in order to find finite-size corrections to the scaling power laws. In the frame of this theory it is shown that the first order finite-size correction to the scaling power laws has following form $S(r) \cong cr^{\alpha_0}[\ln(r/\eta)]^{\alpha_1}$, where $\eta$ is a finite-size scale (in particular for turbulence, it can be the Kolmogorov dissipation scale). Using data of laboratory experiments and numerical simulations it is shown shown that a degenerate case with $\alpha_0 =0$ can describe turbulence statistics in the near-dissipation range $r > \eta$, where the ordinary (power-law) scaling does not apply. For moderate Reynolds numbers the degenerate scaling range covers almost the entire range of scales of velocity structure functions (the log-corrections apply to finite Reynolds number). Interplay between local and non-local regimes has been considered as a possible hydrodynamic mechanism providing the basis for the degenerate scaling of structure functions and extended self-similarity. These results have been also expanded on passive scalar mixing in turbulence. Overlapping phenomenon between local and non-local regimes and a relation between position of maximum of the generalized energy input rate and the actual crossover scale between these regimes are briefly discussed.

Since the strong clustering of luminous matter in the observable universe is a consequence of an initial non-uniformity of the baryon-photon fluid in the last scattering surface, an investigation of clustering of the isotherms in the cosmic microwave background has been performed. The isotherms clustering has been related to the baryon-photon fluid dynamics and the Taylor-microscale Reynolds number of this motion is estimated to be $10^2$. 3-year WMAP cosmic microwave background map has been used in this investigation.

Using the WMAP cosmic microwave background data it is shown that collisions between Alfven wave packets and the cascades generated by these collisions (the Iroshnikov model) can determine the photon temperature fluctuations for arcminute scales on the last scattering surface.

Broad theoretical arguments are proposed to show, formally, that the magnitude G of the temperature gradients in turbulent thermal convection at high Rayleigh numbers obeys the same advection-diffusion equation that governs the temperature fluctuation T, except that the velocity field in the new equation is substantially smoothed. This smoothed field leads to a -1 scaling of the spectrum of G in the same range of scales for which the spectral exponent of T lies between -7/5 and -5/3. This result is confirmed by measurements in a confined container with cryogenic helium gas as the working fluid for Rayleigh number Ra=1.5x10^{11}. Also confirmed is the logarithmic form of the autocorrelation function of G. The anomalous scaling of dissipation-like quantities of T and G are identical in the inertial range, showing that the analogy between the two fields is quite deep.

The Extended Self-Similarity (ESS) of cosmic microwave background (CMB) radiation has been studied using recent data obtained by the space-craft based Wilkinson Microwave Anisotropy Probe. Using the ESS and the high angular scale resolution (arcminutes) of the data it is shown that the CMB temperature space {\it increments} exhibit considerable and systematic declination from Gaussianity for high order moments at the small angular scales. Moreover, the CMB space increment ESS exponents have remarkably close values to the ESS exponents observed in turbulence (in magnetohydrodynamic turbulence).

A combination of spectral and moments analysis of the continuous X-ray flux data is used to show consistency of statistical properties of the powerful solar flares with 2D BTW prototype model of self-organized criticality.

Multiscaling analysis of differential flux dissipation rate of galactic cosmic rays (Carbon nuclei) is performed in the energy ranges: 56.3-73.4 Mev/nucleon and 183.1-198.7 MeV/nucleon, using the data collected by ACE/CRIS spacecraft instrument for 2000 year. The analysis reveals strong (turbulence-like) intermittency of the flux dissipation rate for the short-term intervals: 1-30 hours. It is also found that type of the intermittency can be different in different energy ranges.

Multiscaling properties of the Anderson localization of the cosmic electromagnetic fields before the recombination time are studied and results of a numerical simulation for a random banded matrix ensemble are found to be in good agreement with the large-scale cosmic data.

Scaling properties of the cosmic microwave background (CMB) radiation are studied using satellite (COBE-DMR maps), balloon-borne and ground-based (combined QMASK map) data. Quantitative consistency is found between the multiscaling properties of the COBE-DMR and QMASK CMB maps. Surprisingly, it is found that the observed CMB temperature multiscaling resembles quantitatively the multiscaling properties of fluid turbulence, that indicates primordial plasma turbulence as an origin of the CMB temperature space anisotropy.