Here we present datasets provided by a SCINDA GNSS receiver installed in the Lisbon airport area from November of 2014 to July of 2019. The installed equipment is a NovAtel EURO4 with a JAVAD Choke-Ring antenna. The data are in an archived format and include the general messages on quality of records (*.msg), RANGE files (*.rng), raw observables as the signal-to-noise (S/N) ratios, pseudoranges and phases (*.obs), receiver position information (*.psn), ionosphere scintillations monitor (ISMRB; *.ism) and ionospheric parameters: total electron content (TEC), rate of change of TEC index (ROTI), and the scintillation index S4 (*.scn). The presented data cover the full 2015 year. The raw data are of 1-minute resolution and available for each of the receiver-satellite pairs. The processing and the analysis of the ionosphere scintillation datasets can be done using a specific "SCINDA-Iono" toolbox for the MATLAB developed by T. Barlyaeva (2019) and available online via MathWorks File Exchange system. The toolbox calculates 1-hour means for ionospheric parameters for each of the available receiver-satellite pairs and averaged over all available satellites during the analyzed hour. Here we present the processed data for the following months in 2015: March, June, October, and December. The months were selected as containing most significant geomagnetic events of 2015. The 1-hour means for other months can be obtained from the raw data using the aforementioned toolbox. The provided datasets are interesting for the GNSS and ionosphere based scientific communities.

Here we present a "SCINDA-Iono" toolbox for the MATLAB. This is a software to analyze ionosphere scintillation indices provided by a SCINDA GNSS receiver. The toolbox is developed in the MATLAB R2018b. This software allows to preprocess the original data and analyze ionosphere scintillations on the 1-minute and 1-hour time scales both for averaged over all available satellites values and separately for each receiver-satellite pair.

Recently the atmospheric muon spectra at high energies were reconstructed for two ranges of zenith angles, basing on the events collected with the IceCube detector. These measurements reach high energies at which the contribution to atmospheric muon fluxes from decays of short-lived hadrons is expected. Latest IceCube measurements of the high-energy atmospheric muon spectrum indicate the presence of prompt muon component at energies above 500 TeV. In this work, the atmospheric conventional muon flux in the energy range 10 GeV - 10 PeV is calculated using a set of hadronic models in combination with known parameterizations of the cosmic ray spectrum by Zatsepin $\&$ Sokolskaya and by Hillas $\&$ Gaisser. The calculation of the prompt muons with use of the quark-gluon string model (QGSM) reproduces the muon data of the IceCube experiment. Nevertheless, an additional contribution to the prompt muon component is required to describe the IceCube muon spectra in case if a charm production model predicts the appreciably lower prompt lepton flux as compared with QGSM. This addition, apparently originating from rare decay modes of the short-lived unflavored mesons $\eta, \eta^\prime, \rho, \omega, \phi$, might ensure the competing contribution to the high-energy atmospheric muon flux.

Cosmic rays originating from extraterrestrial sources are permanently arriving at Earth atmosphere, where they produce up to billions of secondary particles. The analysis of the secondary particles reaching to the surface of the Earth may provide a very valuable information about the Sun activity, changes in the geomagnetic field and the atmosphere, among others. In this article, we present the first preliminary results of the analysis of the cosmic rays measured with a high resolution tracking detector, TRAGALDABAS, located at the Univ. of Santiago de Compostela, in Spain.

We calculate the atmospheric neutrino fluxes in the energy range $100$ GeV -- $10$ PeV with the use of several known hadronic models and a few parametrizations of the cosmic ray spectra which take into account the knee. The calculations are compared with the atmospheric neutrino measurements by Frejus, AMANDA, IceCube and ANTARES. An analytic description is presented for the conventional ($\nu_\mu+\bar\nu_\mu$) and ($\nu_e+\bar\nu_e$) energy spectra, averaged over zenith angles, which can be used to obtain test data of the neutrino event reconstruction in neutrino telescopes. The sum of the calculated atmospheric $\nu_\mu$ flux and the IceCube best-fit astrophysical flux gives the evidently higher flux as compared to the IceCube59 data, giving rise the question concerning the hypothesis of the equal flavor composition of the high-energy astrophysical neutrino flux. Calculations show that the transition from the atmospheric electron neutrino flux to the predominance of the astrophysical neutrinos occurs at $30-100$ TeV if the prompt neutrino component is taken into consideration. The neutrino flavor ratio, extracted from the IceCube data, does not reveal the trend to increase with the energy as is expected for the conventional neutrino flux in the energy range $100$ GeV - $30$ TeV. A depression of the ratio $R_{\nu_\mu/\nu_e}$ possibly indicates that the atmospheric electron neutrino flux obtained in the IceCube experiment contains an admixture of the astrophysical neutrinos in the range $10-50$ TeV.

We present the analysis of temperature variations in Portugal for 140 years (from 1865 to 2005). The two stations with the longest time series of temperature measurements (Lisbon and Coimbra) have been used to study the dependence of the portuguese climate variations on the changes of both global circulation and solar activity. Monthly averaged temperature series have been analyzed together with monthly North-Atlantic Oscillation index data, sunspot numbers and cosmic ray flux intensity. Different statistical methods (the correlation analysis and the multiple regression analysis) were used. Our results show that the temperature in Portugal depends not only on the atmospheric variations in the studied region but also on the variations of the solar parameters such as sunspot numbers and the cosmic rays flux intensity. Also, the dependence of temperature on solar parameters is strong during the cold season (November to February), while much weaker during the warm season. For some months, strong correlations between the temperature variations of the current month and the North-Atlantic Oscillation index values of the previous month have been found. The correlation between climatic and solar parameters shows up best on the decadal and decadal-to-centennial timescale. It is found that the temperature correlates positively with the sunspot numbers and negatively with the CR flux intensity throughout the year. Besides, the absolute values of the correlation coefficients between the temperature and the CR are higher than those between the temperature and the sunspot numbers. Our results are consistent with some of the proposed mechanisms that relate solar activity to Earth climate and could be explained through the effect of the solar UV radiation and stratosphere-troposphere coupling or/and through the effect of the CR particles on clouds and stratospheric and tropospheric conditions.