No Arabic abstract
Holm (ASR, 2018) claims that Scafetta (ASR 57, 2121-2135, 2016) is irreproducible because I would have left undocumented the values of two parameters (a reduced-rank index p and a regularization term) that he claimed to be requested in the Magnitude Squared Coherence Canonical Correlation Analysis (MSC-CCA). Yet, my analysis did not require such two parameters. In fact: 1) using the MSC-CCA reduced-rank option neither changes the result nor was needed since Scafetta (2016) statistically evaluated the significance of the coherence spectral peaks; 2) the analysis algorithm neither contains nor needed the regularization term. Herein, I show that Holm could not replicate Scafetta (2016) because he used different analysis algorithms. In fact, although Holm claimed to be using MSC-CCA, for his figures 2-4 he used a MatLab code labeled gcs_cca_1D.m (see paragraph 2 of his Section 3), which Holm also modified, that implements a different methodology known as the Generalized Coherence Spectrum using the Canonical Correlation Analysis (GCS-CCA). This code is herein demonstrated to be unreliable under specific statistical circumstances such as those required to replicate Scafetta (2016). On the contrary, the MSC-CCA method is stable and reliable. Moreover, Holm could not replicate my result also in his figure 5 because there he used the basic Welch MSC algorithm by erroneously equating it to MSC-CCA. Herein I clarify step-by-step how to proceed with the correct analysis, and I fully confirm the 95% significance of my results. I add data and codes to easily replicate my results.
During the last few years a number of works have proposed that planetary harmonics regulate solar oscillations and the Earth climate. Herein I address some critiques. Detailed analysis of the data do support the planetary theory of solar and climate variation. In particular, I show that: (1) high-resolution cosmogenic 10Be and 14C solar activity proxy records both during the Holocene and during the Marine Interglacial Stage 9.3 (MIS 9.3), 325-336 kyr ago, present four common spectral peaks at about 103, 115, 130 and 150 yrs (this is the frequency band that generates Maunder and Dalton like grand solar minima) that can be deduced from a simple solar model based on a generic non-linear coupling between planetary and solar harmonics; (2) time-frequency analysis and advanced minimum variance distortion-less response (MVDR) magnitude squared coherence analysis confirm the existence of persistent astronomical harmonics in the climate records at the decadal and multidecadal scales when used with an appropriate window length (110 years) to guarantee a sufficient spectral resolution. However, the best coherence test can be currently made only by comparing directly the temperature and astronomical spectra as done in Scafetta (J. Atmos. Sol. Terr. Phys. 72(13), 951-970, 2010). The spectral coherence between planetary, solar and climatic oscillations is confirmed at the following periods: 5.2 yr, 5.93 yr, 6.62 yr, 7.42 yr, 9.1 yr (main lunar tidal cycle), 10.4 yr (related to the 9.93-10.87-11.86 yr solar cycle harmonics), 13.8-15.0 yr, 20 yr, 30 yr and 61 yr, 103 yr, 115 yr, 130 yr, 150 yr and about 1000 year. This work responds to the critiques of Cauquoin et al. (Astron. Astrophys. 561, A132, 2014) who ignored alternative planetary theories of solar variations, and of Holm (J. Atmos. Sol. Terr. Phys. 110-111, 23-27, 2014) who used inadequate physical and time frequency analysis of the data.
Here we respond briefly to a comment on our work in arXiv:cond-mat/0405501.
Here we use a very simple conceptual model in an attempt to reduce essential parts of the complex nonlinearity of abrupt glacial climate changes (the so-called Dansgaard-Oeschger events) to a few simple principles, namely (i) a threshold process, (ii) an overshooting in the stability of the system and (iii) a millennial-scale relaxation. By comparison with a so-called Earth system model of intermediate complexity (CLIMBER-2), in which the events represent oscillations between two climate states corresponding to two fundamentally different modes of deep-water formation in the North Atlantic, we demonstrate that the conceptual model captures fundamental aspects of the nonlinearity of the events in that model. We use the conceptual model in order to reproduce and reanalyse nonlinear resonance mechanisms that were already suggested in order to explain the characteristic time scale of Dansgaard-Oeschger events. In doing so we identify a new form of stochastic resonance (i.e. an overshooting stochastic resonance) and provide the first explicitly reported manifestation of ghost resonance in a geosystem, i.e. of a mechanism which could be relevant for other systems with thresholds and with multiple states of operation. Our work enables us to explicitly simulate realistic probability measures of Dansgaard-Oeschger events (e.g. waiting time distributions, which are a prerequisite for statistical analyses on the regularity of the events by means of Monte-Carlo simulations). We thus think that our study is an important advance in order to develop more adequate methods to test the statistical significance and the origin of the proposed glacial 1470-year climate cycle.
Seismic attributes calculated by conventional methods are susceptible to noise. Conventional filtering reduces the noise in the cost of losing the spectral bandwidth. The challenge of having a high-resolution and robust signal processing tool motivated us to propose a sparse time-frequency decomposition while is stabilized for random noise. The procedure initiates by using Sparsity-based adaptive S-transform to regularize abrupt variations in frequency content of the nonstationary signals. Then, considering the fact that a higher amplitude of a frequency component results in a higher signal to noise ratio, an adaptive filter is applied to the time-frequency spectrum which is sparcified previously. The proposed zero adaptive filter enhances the high amplitude frequency components while suppresses the lower ones. The performance of the proposed method is compared to the sparse S-transform and the robust window Hilbert transform in estimation of instantaneous attributes by applying on synthetic and real data sets. Seismic attributes estimated by the proposed method is superior to the conventional ones in terms of its robustness and high resolution image. The proposed approach has a vast application in interpretation and identification of geological structures.
In April and May 2019, as a part of the National Geographic and Roxel Perpetual Planet Everest Expedition, the most interdisciplinary scientific ever was launched. This research identified changing dynamics, including emergent risks resulting from natural and anthropogenic change to the natural system. We have identified compounded risks to ecosystem and human health, geologic hazards, and changing climate conditions that impact the local community, climbers, and trekkeers in the future. This review brings together perspectives from across the biological, geological, and health sciences to better understand emergent risks on Mt. Everest and in the Khumbu region. Understanding and mitigating these risks is critical for the ~10,000 people living in the Khumbu region, as well as the thousands of visiting trekkers and the hundreds of climbers who attempt to summit each year.