# Spatially Decomposed Low Cloud Amounts and Wavelet Analysis

^{*}

**Corresponding Author:**Neïla Zarrouk, Laboratoire de Physique de la Matière Condensée, Faculté des Sciences de Tunis, Université Tunis El Manar, Tunisia, Tel: 0021623066470, Email: neila.zarrouk@yahoo.fr

*
Received: 20-Sep-2012 /
Accepted Date: 19-Nov-2012 /
Published Date: 21-Nov-2012 DOI: 10.4172/2157-7617.1000125 *

### Abstract

In this work we have used monthly data series of low cloud amounts (LCA) describing the repartition and evolution of these cloud amounts in space and time. Thus we have decomposed values of low cloud amounts into spherical harmonics. For each month we have decomposed in spherical harmonics the values of LCA varying in space. Limiting our development to the fourth spherical wave we have obtained 35 components for each month covering the period 01/1990-12/1992. Each correspondent spherical component describing monthly variations of LCA is then submitted to Morlet wavelet analysis. Generally we can notice that the most LCA are obtained for winter and spring. Indeed Sun angle, solar intensity and clouds are important controls of the weather. At the equinox the Sun passes directly overhead at solar noon at the equator. Although the noon solar Sun angle on the equinox is greater than that on the first day of winter, the first remains low inducing also high amounts of clouds.

**Keywords:** Low cloud amounts; Spherical harmonics; Wavelets

#### Introduction

Spherical harmonics have been used to solve potential problems in physics. One important property of SH is how projections interact with rotations. The SH are closed under rotation. The correspondent rotation matrices are in a block structure, where each band is rotationally independent [1-4]. The rotational invariance for SH is similar to the translational invariance in the Fourier transform [5,6].

In association with SH we have investigated wavelet analysis already used in our previous works [7-11] to show the hidden characteristics in time evolution.

The Continuous Wavelet Transform (CWT) is an ideal tool for mapping the changing properties of non-stationary signals and also to determine whether or not a signal is stationary in a global sense. CWT is then used to build a time-frequency representation of a signal that offers very good time and frequency localization [12-15]. An adjustable parameter controlling the properties of the localized oscillation is generally included in the mother wavelet which can be complex or real.

Morlet wavelets are more powerful in revealing hidden detailed structures. This type of wavelet has been used in our previous work [10] and in other works to examine the processes, models and structures of cosmic rays, solar activity and modulations in low cloud amounts [16-19].

Thus instead of analyzing with Morlet wavelets the apparent time series of LCA we increased our analysis quality by studying the decomposition of LCA in spherical harmonics (eq.1). The obtained components are then submitted to Morlet analysis [12] for the 36 months covering the period between 1990-1992.

**Spherical harmonics components of LCA submitted to Morlet analysis**

**Low cloud amounts expanded into spherical harmonics:** Just as the Fourier basis represents an important tool for evaluation of convolutions in a one or two dimensional space, the spherical harmonics basis is a similar tool but defined on the surface of a sphere [20].

Spherical waves propagate isotropically, that is, the wave amplitude is spherically symmetric. We have decomposed values of low cloud amounts (LCA) into spherical harmonics. Then instead of analyzing LCA values directly we have proceeded by decomposing these values in spherical harmonics. Therefore, they are the correspondent components which were analyzed with Morlet wavelets.

Since, spherical harmonics form an orthonormal basis, projecting a spherical function into spherical harmonic coefficients is simple and a straight-forward application of definition (eq.1)

A continuous function can be expanded as a linear combination of spherical harmonics

(1)

Where* l *,*m* are integer indices, is a constant coefficient known in mathematics under the name of generalized Fourier coefficient relative to this base.

The functions where *l *= *0*, are called zonal harmonics (ZH), they have rotational symmetry around the z axis and the zeros (locations where the function is zero) are contours on the sphere parallel to the XY plane.

The functions where are called sectorial harmonics [21] and the zeros define regions like apple slices.

Everything needed to approximate a given function *f *with arbitrary accuracy is to compute the coefficients describing how much each basis function is like *f* .

Generalized and normalized spherical harmonics are given ^{(1-4)} in the final form:

(2)

is the polynomial of Legendre defined by :

(3)

We have examined monthly averaged low cloud cover data [22] from the International Satellite Cloud Climatology Project (ISCCP). For each month the corresponding data table gives the spatial repartition of LCA in latitude and longitude. We have decomposed into spherical harmonics the values of LCA for 36 months covering the period 1990- 1992. We have limited our spherical harmonics decomposition to the first 5 spherical harmonic bands *l* = 0,.......4 , on the same hand the different* m* vary : −*l*..... + *l* resulting in a total of 25 spherical harmonic components for each month.

**Morlet analysis of spherical components for LCA:** At this step we have investigated Morlet wavelet already used in previous works [7-11] Morlet wavelet is defined as a complex sine wave, localized with a Gaussian. The frequency domain representation is a single symmetric Gaussian peak. This wavelet incorporates a wave of a certain period, it is finite in extent and is given by:

(4)

Already wavelet analysis is a useful tool both to find the dominant mode of variation and also to study how it varies with time [15].

In this paper we have directly decomposed in Morlet wavelets the low cloud amounts relative to some countries for the 36 months covering the period 1990-1992 (**Figure 1**).

The yearly modulation is clear for all reconstructed variations of LCA. Likewise, a structure of 6 months appears also for Morlet reconstructed LCA variations relative particularly to Rome, Alger and Rabat.

On the other side and looking further for our analysis, we have used Morlet wavelets to study the time variation of the resulting 25 spherical harmonic components (t) of LCA during 36 months between 1990-1992.

For analysis of a discrete signal (t) and particularly for our case of study we have investigated limited and discrete time series. Therefore the discretized expressions of real parts for Morlet coefficients [10] are then given as follows:

, (5)

The phase constant was taken ω_{0} = 24*π* . The period was fixed to *T* = 1 *month*

After decomposing the monthly time series relative to each spherical component of LCA, we have reconstructed into Morlet wavelets the 25 spherical harmonics components of LCA. We have obtained 35 reconstructed components considering both parts for the complex components the module and the angle.

(6)

It’s clear that a yearly modulation occurs in all figures (**Figures 2-4**) for reconstructed spherical harmonics components of LCA and for the whole period of 36 months. The relative maxima of LCA for this period are generally localised around months 11-12 or months 1-3, this is clearer for module components of spherical harmonics. For all sectorial harmonics (**Figure 2**) and particularly for module components, the shapes of LCA variations are the same, the maxima appear around months 2-3. The angle components present maxima around months 8-9. However the modulation is not very pronounced and is clearer for non-sectorial harmonics. The highest LCA are obtained for the module component *f*_{44} , *f*_{4−4} . Thus LCA are more condensed in winter, autumn and the beginning of spring relative to spherical harmonics components covering the 3 years 1990-1992.

A secondary structure around 6 months occurs clearly in the module part of the fourth spherical harmonic component of LCA (**Figure 4**). Although, the spherical harmonic *f*_{3−1} is a maximum in month 6, the minima and the secondary minima relative to these harmonics are generally around months 6-7, thus this means that the lowest amounts of clouds are found in middle of summer which is expected.

Generally we can notice that the most LCA are obtained for winter and spring.

Indeed Sun angle, solar intensity and clouds are important controls of the weather. Clouds are great reflectors of solar energy. This reflected energy is energy that cannot be used to warm the earth’s surface.

The solar intensity is much lower when the sun is closer to the horizon. When the sun is low on the horizon there is a tendency for more clouds in moist climates. This is because a weaker solar intensity produces a weaker evaporative potential. When the solar intensity is strong (sun high in the sky) and the atmosphere is stable, the warming tends to produce less clouds since moisture is evaporated away into the dry air. At the equinox the Sun passes directly overhead at solar noon at the equator. Although the noon solar Sun angle on the equinox is greater than that on the first day of winter, the first remains low inducing also high amounts of clouds [22].

#### Conclusion

From reconstructed spherical harmonics components of LCA we have detected the yearly modulation. The relative maxima of LCA for this period are generally localised around months 11-12 or months 1-3. The shapes of LCA variations are the same for module components, the maxima appear around months 2-3. The highest LCA are obtained for the module component *f*_{44} , *f*_{4−4} . Thus LCA are more condensed in winter, autumn and the beginning of spring relative to spherical harmonics.

A secondary structure of 6 months occurs also. The minima and the secondary minima are generally around months 6-7, thus the lowest amounts of clouds are found in middle of summer.

The most LCA are then obtained for winter and spring. Indeed Sun angle, solar intensity and clouds are important controls of the weather. Although the noon solar Sun angle on the equinox is greater than that on the first day of winter, the first remains low inducing also high amounts of clouds.

Spherical harmonics decomposition, already known as “Swiss army knife of mathematical physics” associated with Morlet wavelets expansion greatly increase our analysing tool power. Indeed it seems more rigorous to analyse in time the fine structure of LCA represented by spherical components than to analyse an averaged value of LCA.

We plan to study in another work the correlation between low cloud amounts and solar activity through spherical harmonics decomposition.

#### References

- MacRobert TM (1948) An elementary treatise on spherical harmonics: harmonic functions with applications. Dover Publications, NewYork, USA
- Ruffini G, Marco J, Grau C (2002) Spherical Harmonics Interpolation, Computation of Laplacians and Gauge Theory. Starlab Technical Notes
- Axler S, Bourdon P, Ramey W (2000) Harmonic Function Theory. Springer-Verlag
- Anen T, Kautz J, Durand F, Seidel H (2004) Spherical Harmonic Gradients for Mid-Range Illumination. Eurographics Symposium on Rendering
- Green R (2003) Spherical Harmonic Lighting: The Gritty Details. Archives of the Game Developers Conference, 2003
- Wang J, Xu K, Zhou K, Lin S, Hu S, et al. (2006) Spherical Harmonic Scaling. Visual Comput 22: 713-720
- Zarrouk N, Bennaceur R (2009) A wavelet based analysis of cosmic rays modulation. Acta Astronaut 65: 262-272
- Zarrouk N, Bennaceur R (2010) Neural network and wavelets in prediction of cosmic ray variability: the North Africa as study case. Acta Astronaut 66: 1008–1016
- Zarrouk N, Bennaceur R (2009) Extrapolating cosmic ray variations and impacts on life : Morlet wavelet analysis. Int J Astrobio 8: 169-174
- Zarrouk N, Bennaceur R (2010) Link nature between low cloud amounts and cosmic rays through wavelet analysis. Acta astronaut 66: 1311–1319
- Zarrouk N, Bennaceur R (2011) Cosmic rays transport equation solution and wavelets analysis: A case of pulsar source. Acta Astronaut 68: 1650–1659
- Morlet J, Arehs G, Fourgeau I, Giard D (1982) Wave propagation and sampling theory. Geophys. 47: 203-221
- Kozlov VI, Markov VV (2007) Wavelet image of the fine structure of the 11-year cycle based on studying cosmic ray fluctuations during cycles 20-23. Geomagn Aeron+ 47: 47-55
- Kozlov VI, Markov VV (2007) Wavelet image of a Heliospheric Storm in Cosmic rays. Geomagn Aeron+ 47: 56–65
- Prokoph A, Patterson RT (2004) Application of Wavelet and regression Analysis in assessing temporal and geographic climate variability: Eastern, Ontario, Canada as a case study. Atmos Ocean 42: 201-212
- Frick P, Galyagin D, Hoyt DV, Nesme-ribes E, Schatten KH, et al. (1997) Wavelet analysis of solar activity recorded by sunspot groups. Astron Astrophys 328: 670-681
- Kudela K, Storini M, Antalova A, Rybak J (2001) On the wavelet approach to cosmic ray variability. 27th International Cosmic Ray Conference (ICRC) Hamburg, Germany
- Boberg F, Lundstedt H, Hoeksema JT, Scherrer PH, Lui W (2002) Solar mean magnetic field variability: A wavelet approach to Wilcox Solar Observatory and SOHO/Michelson Doppler Image observations. J Geophys Res 107: 1318
- Sternberg WJ, Smith TL (1946) The Theory of Potential and Spherical Harmonics. University of Toronto Press, Toronto, Canada

Citation: Zarrouk N, Bennaceur R (2012) Spatially Decomposed Low Cloud Amounts and Wavelet Analysis. J Earth Sci Climate Change 3: 125. Doi: 10.4172/2157-7617.1000125

Copyright: ©2012 Zarrouk N, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.