Atmospheric concentrations of greenhouse gases are rising, leading to a positive radiative forcing of climate and an expected warming of surface temperatures. In general, fractal properties may be observed in the time series of the dynamics of complex systems. To study the relation between the atmospheric CO 2 concentration and the climate indices, we investigated the change of fractal behavior of the CO 2, the carbon isotope ratio (δ 13C) of atmospheric CO 2, the El Niño-Southern Oscillation (ENSO), the Pacific Decadal Oscillation (PDO), and the North Atlantic Oscillation (NAO) indices using the multifractal analysis. When the atmospheric CO 2 growth rate was large, the multifractality of CO 2, δ 13C in CO 2, ENSO, and NAO was large and the changes were large from the change of fractality. The changes of CO 2 and ENSO were closely related and the influence of the CO 2 on the ENSO was strong from the change in fractality and wavelet coherence. When the El Niño occurred, the CO 2 growth rate was large. The CO 2 related to PDO, NAO, and global temperature from the change in fractality and wavelet coherence. Especially, the changes of CO 2 and global temperature were closely related. When the global warming hiatus occurred, the multifractality of the global temperature was weaker than that of CO 2 and the change of the global temperature was stable. These findings will contribute to the research of the relation between the atmospheric CO 2 and climate change.
Atmospheric concentrations of greenhouse gases are rising, leading to a positive radiative forcing of climate and an expected warming of surface temperatures.
Many of the El Niño events appear to be associated with a net transfer of CO2 from the biosphere to the atmosphere. ENSO events produce droughts, particularly in tropical regions, which can cause drought and forest fires and increase CO2 (Keeling et al., 1995; Langenfelds et al., 2002) .
The North Atlantic Oscillation (NAO) and the Pacific-North America (PNA) indices are correlated with the observed atmospheric CO2 growth rate (Murayama et al., 2004) .
Various objects in nature show the so-called self-similarity or fractal property. Monofractal shows a roughly similar pattern at different scales and is characterized through a fractal dimension. Multifractal is a non-uniform, more complex fractal and is separated into many sub-sets characterized through different fractal dimensions. Fractal property can be observed in the time series representing dynamics of complex systems as well. A change of fractality occurs with a phase transition and changes of state. The multifractal properties of daily rain were studied in two contrasting climates: an East Asian monsoon climate with drastic rain variability and a mild climate with moderate rain variability (Svensson et al., 1996) . In both the climates, the frontal rain shows monofractality and the convective-type rain shows multifractality.
Hence, climate change can be interpreted from the view of fractals. A change of fractality may be observed when the climate changes. We use the wavelet transform to analyze the multifractal behavior of the climate index. Wavelet methods are useful for the analysis of complex non-stationary time series. The wavelet transform allows good multifractal analysis to be performed (Muzy et al., 1991) . We used the wavelet transform to analyze the multifractal behavior of the climate index. We concluded that a climatic regime shift corresponds to a change from multifractality to monofractality of the Pacific Decadal Oscillation (PDO) index (Maruyama et al., 2015) .
To study the relation between the atmospheric CO2 concentration and the climate indices, we investigated the change of multifractal behavior of the CO2, the carbon isotope ratio (δ13C) of atmospheric CO2, the Southern Oscillation Index (SOI), and the NAO index using the multifractal analysis.
The SOI, Niño3.4, PDO, NAO indices, and global mean surface air temperature anomalies provided by NOAA’s Climate Prediction Center, USA (CPC) were used. The SOI is a standardized index based on the observed sea level pressure differences between Tahiti and Darwin, Australia. The SOI is one measure of the large-scale fluctuations in air pressure occurring between the western and eastern tropical Pacific (i.e., the state of the Southern Oscillation) during El Niño and La Niña episodes. Prolonged periods of negative (positive) SOI values coincide with abnormally warm (cold) ocean waters across the eastern tropical Pacific typical of El Niño (La Niña) episodes.
Atmospheric CO2 concentrations (ppm) derived from in situ air measurements at Mauna Loa, Observatory, Hawaii was used. Monthly atmospheric 13C concentrations (per mil) in CO2 derived from flask air samples at Mauna Loa Observatory was used. Annual mean growth rate of CO2 for Mauna Loa obtained from NOAA was used.
For the examination, we used the Daubechies wavelet, which is widely used in solving a broad range of problems, e.g., self-similarity properties of a signal and signal discontinuities. We made use of a discrete signal that was fitted the Daubechies mother wavelet with the capacity of correct inverse transformation. Thus, we can precisely calculate the following best τ(q), which can be regarded as a characteristic function of the fractal behavior. We can define the τ(q) from the power-law behavior of the partition function, as shown in equation (2). We then computed the scaling of the partition function Zq(a), which is defined as the sum of the q-th powers of the modulus of the wavelet transform coefficients at scale a, where q is the q-th moment. In our computation, the wavelet-transform coefficients did not grow zero. Thus, for a correct calculation, the summation was considered for the whole set. Muzy et al. (1991) defined Zq(a) as the sum of the q-th powers of the local maxima of the modulus to avoid dividing by zero. We got the following partition function Zq(a):
Z q ( a ) = ∑ | W φ [ f ] ( a , b ) | q (1)
where W φ [ f ] ( a , b ) , a and b are the wavelet coefficient of function f, a scale parameter, and a space parameter, respectively. W φ [ f ] ( a , b ) is defined as below.
W φ [ f ] ( a , b ) = 1 | a | ∫ − ∞ + ∞ f ( t ) ϕ * ( t − b a ) d t (2)
where f ( t ) is data and φ is wavelet function. For small scales, we expect
Z q ( a ) ~ a τ ( q ) (3)
First, we examined the changes in Zq(a) in the time series at a different scale a for each moment q. We plotted the logarithm of Zq(a) against that of time scale a. Here τ(q) is the slope of the fitted straight line for each q. Next, we plotted τ(q) versus q. The time window was advanced by one year, which was repeated. The time window was fixed to 6 years, when a moderate change in fractality was observed. Monofractal and multifractal signals were defined as follows. For τ(q), a monofractal signal corresponds to a straight line, while a multifractal signal is nonlinear (Frish & Parisi, 1985) . We calculated the R2 value, which is the coefficient of determination, for the fitted straight line. If R2 ≥ 0.98, the time series is monofractal and if 0.98 > R2, it is multifractal.
We calculated the τ(q) of different moments q for individual records in the Niño3.4 index. In
was calculated from n80 up to n89. A monofractal signal would correspond to a straight line for τ(q), while a multifractal signal would be nonlinear. Nearly all the multifractality observed is due to the negative value of q, i.e., small fluctuations are more inhomogeneous than large fluctuations. In
We plotted the value of the τ(−6) in each index. The negative large values of the τ(−6) show large multifractality. For the τ(q), q = −6 is the appropriate number to show the change of τ.
The time series of atmospheric CO2 and δ13C in CO2 are shown in
We investigated the relation between the CO2 and δ13C. The τ(−6) of the CO2, and δ13C are shown in
We studied the relationship between the CO2 and δ13C by means of wavelet coherence. We show the wavelet coherence and phase using the Morlet wavelet between the CO2 and δ13C in
We investigated the relation between the annual mean CO2 growth rate and the multifractality of CO2. The annual mean CO2 growth rate for Mauna Loa and the τ(−6) of the CO2 are shown in
We investigated the relation between the CO2 and SOI. The τ(−6) of the CO2, and SOI are shown in
We investigated the relation between the CO2 and Niño3.4. The τ(−6) of the CO2, and Niño3.4 are shown in
We investigated the relation between the CO2 and PDO. The τ(−6) of the CO2,
and PDO are shown in
We studied the relationship between the CO2 and PDO by means of wavelet coherence. We show the wavelet coherence and phase using the Morlet wavelet between the CO2 and PDO in
We investigated the relation between the CO2 and NAO. The τ(−6) of the CO2, and NAO are shown in
We studied the relationship between the CO2 and NAO by means of wavelet coherence. We show the wavelet coherence and phase using the Morlet wavelet between the CO2 and NAO in
We investigated the relation between the CO2 and global temperature. The global mean surface temperature anomalies is shown in
We studied the relationship between the CO2 and global temperature by means of wavelet coherence. We show the wavelet coherence and phase using the Morlet wavelet between the CO2 and global temperature in
(middle) and (bottom), respectively. The coherence between the CO2 and global temperature in two year scale was very strong, when the rate of global temperature increased largely. The lead of the global temperature was observed.
The CO2 increased and δ13C decreased with time, so there was an inverse relationship between CO2 and δ13C, which indicated that these changes were mainly generated by activities of terrestrial biosphere. The δ13C decreased due to combustion of fossil fuel with isotopically light CO2 (Nakazawa et al., 1993) .
When the CO2 growth rate was large, the τ(−6) of the CO2 was small and multifractality was large and the change was large. The changes in fractality of the CO2 and δ13C were very similar for 1985-2010. When the CO2 growth rate was large, the τ(−6) of δ13C was small and the multifractality was large. The changes of CO2 and δ13C in CO2 were closely related.
The changes in fractality of the CO2 and SOI were very similar for 1965-2010. The coherence between the CO2 and SOI in 2 - 4 year scale was strong for 1970-2000 and the lead of the CO2 was observed. The CO2 differences between the central and western Pacific Ocean correlate well with the SOI. There is more (less) midtropospheric CO2 in the central Pacific and less (more) midtropospheric CO2 in the western Pacific during El Niño (La Nina) events (Jiang et al., 2013) . When the CO2 growth rate was large, the ENSO event occurred and the τ(−6) of SOI was small and multifractal was strong. The changes of CO2 and ENSO were closely related. The influence of the CO2 on the ENSO (SOI and Niño3.4) was strong from the change of fractality and wavelet coherence. When the El Niño (La Niña) occurred, the CO2 growth rate was large (small). The relationship between the CO2 and the SOI values is shown (Matsueda et al., 2015) .
The changes in fractality of the CO2 and PDO were very similar for 1965-2010 especially for the 1970s and 1990s and the changes of CO2 and PDO were related to each other. The changes in fractality of the CO2 and NAO were very similar for 1965-2000. The coherence between the CO2 and NAO in one year scale was strong and the leads of the CO2 and NAO were observed. When the CO2 growth rate was large, the τ(−6) of NAO was small and multifractality was strong. The changes of CO2 and NAO were related to each other. When the rate of global temperature increased largely, the changes in fractality of the CO2 and global temperature were very similar and the coherence between the CO2 and global temperature was very strong. For 1998-2012 (global warming hiatus), the multifractality of the global temperature was weaker than that of CO2. The change of the global temperature was more stable than that of CO2. The changes of CO2 and global temperature were closely related. Hence, the CO2 related to PDO, NAO, and global temperature from the change in fractality and wavelet coherence.
To study the relation between the atmospheric CO2 concentration and the climate indices, we investigated the change of fractal behavior of the CO2, the carbon isotope ratio (δ13C) of atmospheric CO2, SOI, Niño3.4, PDO, and NAO indices using the multifractal analysis. We showed the change of factuality by plotting the τ-function and used the wavelet coherence. The main findings are summarized below.
1) When the atmospheric CO2 growth rate was large, the multifractality of CO2, δ13C in CO2, ENSO, and NAO was large and the changes were large from the change of fractality.
2) The changes of CO2 and ENSO were closely related and the influence of the CO2 on the ENSO was strong from the change in fractality and wavelet coherence. When the El Niño occurred, the CO2 growth rate was large.
3) The CO2 related to PDO, NAO, and global temperature from the change in fractality and wavelet coherence. Especially, the changes of CO2 and global temperature were closely related. When the global warming hiatus occurred, the multifractality of the global temperature was weaker than that of CO2 and the change of the global temperature was stable.
These findings will contribute to the research of the relation between the atmospheric CO2 and climate change.
The author declares no conflicts of interest regarding the publication of this paper.
Maruyama, F. (2019). Relationship between the Atmospheric CO2 and Climate Indices by Wavelet-Based Multifractal Analysis. Journal of Geoscience and Environment Protection, 7, 38-51. https://doi.org/10.4236/gep.2019.71004