Supernova remnants (SNRs) are believed to be the source of Galactic cosmic rays (CRs). SNR shocks accelerate CR protons and electrons which reveal key insights into the non-thermal physics by means of their synchrotron and $\gamma$-ray emission. The remnant SN 1006 is an ideal particle acceleration laboratory because it is observed across all electromagnetic wavelengths from radio to $\gamma$-rays. We perform three-dimensional (3D) magnetohydrodynamics (MHD) simulations where we include CR protons and follow the CR electron spectrum. By matching the observed morphology and non-thermal spectrum of SN 1006 in radio, X-rays and $\gamma$-rays, we gain new insight into CR electron acceleration and magnetic field amplification. 1. We show that a mixed leptonic-hadronic model is responsible for the $\gamma$-ray radiation: while leptonic inverse-Compton emission and hadronic pion-decay emission contribute equally at GeV energies observed by Fermi, TeV energies observed by imaging air Cherenkov telescopes are hadronically dominated. 2. We show that quasi-parallel acceleration (i.e., when the shock propagates at a narrow angle to the upstream magnetic field) is preferred for CR electrons and that the electron acceleration efficiency of radio-emitting GeV electrons at quasi-perpendicular shocks is suppressed at least by a factor ten. This poses a challenge for the interpretation of plasma particle-in-cell simulations of particle acceleration at shocks. 3. To match the radial emission profiles and the $\gamma$-ray spectrum, we require a volume-filling, turbulently amplified magnetic field and that the Bell-amplified magnetic field is damped in the immediate post-shock region. Our work connects micro-scale plasma physics simulations to the scale of SNRs.

Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ray morphology of SNRs such as Vela Jr. and SN1006 solely by varying the magnetic morphology. This implies that gamma-ray bright regions result from quasi-parallel shocks (i.e., when the shock propagates at a narrow angle to the upstream magnetic field), which are known to efficiently accelerate CR protons, and that gamma-ray dark regions point to quasi-perpendicular shock configurations. Comparison of the simulated gamma-ray morphology to observations allows us to constrain the magnetic coherence scale $\lambda_B$ around Vela Jr. and SN1006 to $\lambda_B \simeq 13_{-4.3}^{+13}$ pc and $\lambda_B >200_{-10}^{+80}$ pc, respectively, where the ambient magnetic field of SN1006 is consistent with being largely homogeneous. We find consistent pure hadronic and mixed hadronic-leptonic models that both reproduce the multi-frequency spectra from the radio to TeV gamma rays and match the observed gamma-ray morphology. Finally, to capture the propagation of a SNR shock in a clumpy interstellar medium, we study the interaction of a shock with a dense cloud with numerical simulations and analytics. We construct an analytical gamma-ray model for a core collapse SNR propagating through a structured interstellar medium, and show that the gamma-ray luminosity is only biased by 30% for realistic parameters.

Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code {\sc arepo} for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.

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.

The increased amount of information provided by ongoing missions such as the Solar Dynamics Observatory (SDO) represents a great challenge for the understanding of basic questions such as the internal structure of sunspots and how they evolve with time. Here, we contribute with the exploitation of new data, to provide a better understanding of the separate growth and decay of sunspots, umbra and penumbra. Using fuzzy sets to compute separately the areas of sunspot umbra and penumbra, the growth and decay rates for active regions NOAA 11117, NOAA 11428, NOAA 11429, and NOAA 11430 are computed, from the analysis of intensitygrams obtained by the Helioseismic and Magnetic Imager onboard SDO. A simplified numerical model is proposed for the decay phase, whereby an empirical irrotational and uniformly convergent horizontal velocity field interacting with an axially symmetric and height invariant magnetic field, reproduces the large-scale features of the much more complex convection observed inside sunspots.

The emergence of the magnetic field through the photosphere has multiple manifestations and sunspots are the most prominent examples of this. One of the most relevant sunspot properties, to study both its structure and evolution, is the sunspot area: either total, umbra or penumbra area. Recently Schlichenmaier et al. (2010) studied the evolution of the active region (AR) NOAA 11024 concluding that during the penumbra formation the umbra area remains constant and that the increase of the total sunspot area is caused exclusively by the penumbra growth. In this presentation the Schlichenmaier's conclusion is firstly tested, investigating the evolution of four different ARs. Hundreds of Intensitygram images from the Helioseismic and Magnetic Imager (HMI) images are used, obtained by the Solar Dynamics Observatory, in order to describe the area evolution of the above ARs and estimate the increase and decrease rates for umbra and penumbra areas, separately. A simple magnetohydrodynamic model is then tentatively used in a first approximation to explain the observed results.

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.