Source author record

S. K. Solanki

S. K. Solanki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

148works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

148 published item(s)

preprint2025arXiv

Formation of chromospheric Fe I lines in the near ultraviolet in 1D atmospheres

In the near ultraviolet (NUV) part of the solar spectrum, there are several Fe I lines with very broad profiles, typical of chromospheric lines. These lines are largely unexplored due to the lack of high-resolution data in this region. This changed with the successful Sunrise III flight in 2024, when spectro-polarimetric data were recorded with high spatial, spectral, and temporal resolution covering a large variety of solar targets. The aim of this work is to investigate the formation of the lines and lay the groundwork for further studies. We compute the spectrum of the lines at 358.12 nm, 371.99 nm, 406.36 nm, and 407.17 nm emerging from the standard 1D FAL-atmospheres using the non-local thermodynamic equilibrium (NLTE) radiative transfer code RH. We find that the lines are affected by overionization in the wings, but have line cores strongly affected by scattering. The line cores form well into the chromosphere in the tested atmosphere models except the colder FALX model where the line core forms in the temperature minimum (which lies at traditional chromospheric heights). In the presence of a vertical magnetic field, the Stokes $V$ signal is multi-lobed. The profile can be decomposed into two broad photospheric lobes and two sharper lobes forming in the flanks of the chromospheric line core. We have investigated the properties of four lines of Fe I that form in the lower chromosphere. The results provide a basic understanding of the formation of the lines, which will be useful for later analysis of formation in 3D magnetohydrodynamic simulations and an eventual investigation into their diagnostic potential.

preprint2023arXiv

Focused Space Weather Strategy for Securing Earth, and Human Exploration of the Moon and Mars

This white paper recognizes gaps in observations that will, when addressed, much improve solar radiation hazard and geomagnetic storm forecasting. Radiation forecasting depends on observations of the entire "Solar Radiation Hemisphere" that we will define. Mars exploration needs strategic placement of radiation-relevant observations. We also suggest an orbital solution that will improve geomagnetic storm forecasting through improved in situ and solar/heliospheric remote sensing.

preprint2022arXiv

Chromospheric Extension of the MURaM Code

Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines.

preprint2022arXiv

Faculae cancel out on the surfaces of active Suns

Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres.

preprint2022arXiv

Magnetized supersonic downflows in the chromosphere. A statistical study using the He I 10830 Å lines

The chromosphere above active regions (ARs) on the Sun hosts magnetized supersonic downflows. Studies of these supersonic downflows help to decipher the magnetic fine structure and dynamics of the chromosphere. We perform a statistical analysis of the magnetized supersonic downflows in a number of ARs and survey their characteristics. We analyze spectro-polarimetric scans of parts of 13 ARs obtained in the infrared He I 10830 Å triplet formed in the upper chromosphere recorded with the GREGOR Infrared Spectrograph (GRIS) mounted at the GREGOR solar telescope. We retrieve the line-of-sight velocities and the magnetic field vector using the HeLIx+ inversion code that assumes Milne-Eddington atmospheres. We find magnetized supersonic downflows in all the ARs, with larger area coverage by such flows in ARs observed during their emerging phase. The fact that supersonic downflows were detected in all scans, though they cover only a small fraction, 0.2--6.4%, of the observed field-of-view, suggests that they are a common phenomenon in the upper chromospheres of ARs. The supersonic downflows are found to be associated with many AR features such as pores, sunspot umbrae, sunspot penumbrae, light bridges, plages, He I loops as part of arch filament systems characteristic of emerging fields, and filaments. Although several mechanisms are identified to be causing the supersonic downflows, by far the most common one appears to be the draining of plasma along the legs of rising magnetic loops. The loops mainly drain into forming pores. The line-of-sight velocities of the supersonic downflows reach up to 49 kms-1 and the velocity distribution shows multiple populations. Almost 92% of these supersonic downflows coexist with a subsonic flow component. The weaker, more horizontal fields associated with the supersonic component suggests that it is formed above the subsonic component.

preprint2022arXiv

Predictions of Astrometric Jitter for Sun-like Stars. III. Fast Rotators

A breakthrough in exoplanet detections is foreseen with the unprecedented astrometric measurement capabilities offered by instrumentation aboard Gaia space observatory. Besides, astrometric discoveries of exoplanets are expected from the planned space mission, Small-JASMINE. In this setting, the present series of papers focuses on estimating the effect of magnetic activity of G2V-type host stars on the astrometric signal. This effect interferes with the astrometric detections of Earth-mass planets. While the first two papers considered stars rotating at the solar rotation rate, this paper focuses on stars having solar effective temperature and metallicity but rotating faster than the Sun, and consequently more active. By simulating the distribution of active regions on such stars using the Flux Emergence And Transport model, we show that the contribution of magnetic activity to the astrometric measurements becomes increasingly significant with increasing rotation rates. We further show that the jitter for the most variable periodic Kepler stars is high enough to be detected by Gaia. Furthermore, due to a decrease in the facula-to-spot area ratio for more active stars, the magnetic jitter is found to be spot-dominated for rapid rotators. Our simulations of the astrometric jitter has the potential to aid the interpretation of data from Gaia and upcoming space astrometry missions.

preprint2022arXiv

The magnetic drivers of campfires seen by the Polarimetric and Helioseismic Imager (PHI) on Solar Orbiter

The Extreme Ultraviolet Imager (EUI) on board the Solar Orbiter (SO) spacecraft observed small extreme ultraviolet (EUV) bursts, termed campfires, that have been proposed to be brightenings near the apexes of low-lying loops in the quiet-Sun atmosphere. The underlying magnetic processes driving these campfires are not understood. During the cruise phase of SO and at a distance of 0.523\,AU from the Sun, the Polarimetric and Helioseismic Imager on Solar Orbiter (SO/PHI) observed a quiet-Sun region jointly with SO/EUI, offering the possibility to investigate the surface magnetic field dynamics underlying campfires at a spatial resolution of about 380~km. In 71\% of the 38 isolated events, campfires are confined between bipolar magnetic features, which seem to exhibit signatures of magnetic flux cancellation. The flux cancellation occurs either between the two main footpoints, or between one of the footpoints of the loop housing the campfire and a nearby opposite polarity patch. In one particularly clear-cut case, we detected the emergence of a small-scale magnetic loop in the internetwork followed soon afterwards by a campfire brightening adjacent to the location of the linear polarisation signal in the photosphere, that is to say near where the apex of the emerging loop lays. The rest of the events were observed over small scattered magnetic features, which could not be identified as magnetic footpoints of the campfire hosting loops. The majority of campfires could be driven by magnetic reconnection triggered at the footpoints, similar to the physical processes occurring in the burst-like EUV events discussed in the literature. About a quarter of all analysed campfires, however, are not associated to such magnetic activity in the photosphere, which implies that other heating mechanisms are energising these small-scale EUV brightenings.

preprint2022arXiv

The on-ground data reduction and calibration pipeline for SO/PHI-HRT

The ESA/NASA Solar Orbiter space mission has been successfully launched in February 2020. Onboard is the Polarimetric and Helioseismic Imager (SO/PHI), which has two telescopes, a High Resolution Telescope (HRT) and the Full Disc Telescope (FDT). The instrument is designed to infer the photospheric magnetic field and line-of-sight velocity through differential imaging of the polarised light emitted by the Sun. It calculates the full Stokes vector at 6 wavelength positions at the Fe I 617.3 nm absorption line. Due to telemetry constraints, the instrument nominally processes these Stokes profiles onboard, however when telemetry is available, the raw images are downlinked and reduced on ground. Here the architecture of the on-ground pipeline for HRT is presented, which also offers additional corrections not currently available on board the instrument. The pipeline can reduce raw images to the full Stokes vector with a polarimetric sensitivity of $10^{-3}\cdot I_{c}$ or better.

preprint2021arXiv

A solar coronal loop in a box: Energy generation and heating

Coronal loops are the basic building block of the upper solar atmosphere. Comprehending how these are energized, structured, and evolve is key to understanding stellar coronae. Here we investigate how the energy to heat the loop is generated by photospheric magneto-convection, transported into the upper atmosphere, and how the internal structure of a coronal loop forms. In a 3D magnetohydrodynamics (MHD) model, we study an isolated coronal loop rooted with both footpoints in a shallow layer within the convection zone using the MURaM code. To resolve its internal structure, we limited the computational domain to a rectangular box containing a single coronal loop as a straightened magnetic flux tube. Field-aligned heat conduction, gray radiative transfer in the photosphere and chromosphere, and optically thin radiative losses in the corona were taken into account. The footpoints were allowed to interact self-consistently with the granulation surrounding them. The loop is heated by a Poynting flux that is self-consistently generated through small-scale motions within individual magnetic concentrations in the photosphere. Turbulence develops in the upper layers of the atmosphere as a response to the footpoint motions. We see little sign of heating by large-scale braiding of magnetic flux tubes from different photospheric concentrations at a given footpoint. The synthesized emission, as it would be observed by the Atmospheric Imaging Assembly or the X-ray Telescope, reveals transient bright strands that form in response to the heating events. Overall, our model roughly reproduces the properties and evolution of the plasma as observed within coronal loops. With this model we can build a coherent picture of how the energy flux to heat the upper atmosphere is generated near the solar surface and how this process drives and governs the heating and dynamics of a coronal loop.

preprint2021arXiv

Observing and modelling the young solar analogue EK Draconis: starspot distribution, elemental abundances, and evolutionary status

Observations and modelling of stars with near-solar masses in their early phases of evolution is critical for a better understanding of how dynamos of solar-type stars evolve. We examine the chemical composition and the spot distribution of the pre-main-sequence solar analogue EK Dra. Using spectra from the HERMES Spectrograph (La Palma), we obtain the abundances of 23 elements with respect to the solar ones, which lead to a $[{\rm Fe/H}]=0.03$, with significant overabundance of Li and Ba. The s-process elements Sr, Y, and Ce are marginally overabundant, while Co, Ni, Cu, Zn are marginally deficient compared to solar abundances. The overabundance of Ba is most likely due to the assumption of depth-independent microturbulent velocity. Li abundance is consistent with the age and the other abundances may indicate distinct initial conditions of the pre-stellar nebula. We estimate a mass of 1.04 $M_\odot$ and an age of $27^{+11}_{-8}$\,Myr using various spectroscopic and photometric indicators. We study the surface distribution of dark spots, using 17 spectra collected during 15 nights using the CAFE Spectrograph (Calar Alto). We also conduct flux emergence and transport (FEAT) simulations for EK Dra's parameters and produce 15-day-averaged synoptic maps of the likely starspot distributions. Using Doppler imaging, we reconstruct the surface brightness distributions for the observed spectra and FEAT simulations, which show overall agreement for polar and mid-latitude spots, while in the simulations there is a lack of low-latitude spots compared to the observed image. We find indications that cross-equatorial extensions of mid-latitude spots can be artefacts of the less visible southern-hemisphere activity.

preprint2021arXiv

The Dimmest State of the Sun

How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm-2 to a rise of 6.3 Wm-2. Consequently, the exact contribution by solar forcing to the rise in global temperatures over the past centuries remains inconclusive. Adopting a novel approach based on state-of-the-art solar imagery and numerical simulations, we establish the TSI level of the Sun when it is in its least-active state to be 2.0 +/- 0.7 Wm-2 below the 2019 level. This means TSI could not have risen since the Maunder minimum by more than this amount, thus restricting the possible role of solar forcing in global warming.

preprint2020arXiv

Connecting measurements of solar and stellar brightness variations

Comparing solar and stellar brightness variations is hampered by the difference in spectral passbands used in observations as well as by the possible difference in the inclination of their rotation axes from the line of sight. We calculate the rotational variability of the Sun as it would be measured in passbands used for stellar observations. In particular, we consider the filter systems used by the CoRoT, $Kepler$, TESS, and $Gaia$ space missions. We also quantify the effect of the inclination of the rotation axis on the solar rotational variability. We employ the Spectral And Total Irradiance REconstructions (SATIRE) model to calculate solar brightness variations in different filter systems as observed from the ecliptic plane. We then combine the simulations of the surface distribution of the magnetic features at different inclinations using a surface flux transport model (SFTM) with the SATIRE calculations to compute the dependence of the variability on the inclination. For an ecliptic-bound observer, the amplitude of the solar rotational variability, as observed in the total solar irradiance (TSI) is 0.68 mmag (averaged over solar cycles 21-24). We obtained corresponding amplitudes in the $Kepler$ (0.74 mmag), CoRoT (0.73 mmag), TESS (0.62 mmag), $Gaia~ $ (0.74 mmag), $Gaia~ G_{RP}$ (0.62 mmag), and ), $Gaia~ G_{BP}$ (0.86 mmag) passbands. Decreasing the inclination of the rotation axis decreases the rotational variability. For a sample of randomly inclined stars, the variability is on average 15% lower in all filter systems considered in this work. This almost compensates for the difference in the amplitudes of the variability in TSI and $Kepler$ passbands, making the amplitudes derived from the TSI records an ideal representation of the solar rotational variability for comparison to $Kepler$ stars with unknown inclinations.

preprint2020arXiv

Connecting the Wilson depression to the magnetic field of sunspots

Context: In sunspots, the geometric height of continuum optical depth unity is depressed compared to the quiet Sun. This so-called Wilson depression is caused by the Lorentz force of the strong magnetic field inside the spots. However, it is not understood in detail yet, how the Wilson depression is related to the strength and geometry of the magnetic field or to other properties of the sunspot. Aims: We aim to study the dependence of the Wilson depression on the properties of the magnetic field of the sunspots and how exactly the magnetic field contributes to balancing the Wilson depression with respect to the gas pressure of the surroundings of the spots. Methods: Our study is based on 24 spectropolarimetric scans of 12 individual sunspots performed with Hinode. We derived the Wilson depression for each spot using both, a recently developed method that is based on minimizing the divergence of the magnetic field, and an approach developed earlier that enforces an equilibrium between the gas pressure and the magnetic pressure inside the spot and the gas pressure in the quiet Sun, thus neglecting the influence of the curvature force. We then performed a statistical analysis by comparing the Wilson depression resulting from the two techniques with each other and by relating them to various parameters of the sunspots, such as their size or the strength of the magnetic field. Results: We find that the Wilson depression becomes larger for spots with a stronger magnetic field, but not as much as one would expect from the increased magnetic pressure. This suggests that the curvature integral provides an important contribution to the Wilson depression, particularly for spots with a weak magnetic field. Our results indicate that the geometry of the magnetic field in the penumbra is different between spots with different strengths of the average umbral magnetic field.

preprint2020arXiv

Detection of the strongest magnetic field in a sunspot light bridge

Traditionally, the strongest magnetic fields on the Sun have been measured in sunspot umbrae. More recently, however, much stronger fields have been measured at the ends of penumbral filaments carrying the Evershed and counter-Evershed flows. Superstrong fields have also been reported within a light bridge separating two umbrae of opposite polarities. We aim to accurately determine the strengths of the strongest fields in a light bridge using an advanced inversion technique and to investigate their detailed structure. We analyze observations from the spectropolarimeter on board the Hinode spacecraft of the active region AR 11967. The thermodynamic and magnetic configurations are obtained by inverting the Stokes profiles using an inversion scheme that allows multiple height nodes. Both the traditional 1D inversion technique and the so-called 2D coupled inversions, which take into account the point spread function of the Hinode telescope, are used. We find a compact structure with an area of 32.7 arcsec$^2$ within a bipolar light bridge with field strengths exceeding 5 kG, confirming the strong fields in this light bridge reported in the literature. Two regions associated with downflows of $\sim$5 km s$^{-1}$ harbor field strengths larger than 6.5 kG, covering a total area of 2.97 arcsec$^2$. The maximum field strength found is 8.2 kG, which is the largest ever observed field in a bipolar light bridge up to now.

preprint2020arXiv

Effect of metallicity on the detectability of rotational periods in solar-like stars

Accurate determination of stellar rotation periods is important for estimating stellar ages as well as for understanding stellar activity and evolution. While for about thirty thousand stars in the Kepler field rotation periods can be determined, there are over hundred thousand stars, especially with low photometric variability and irregular pattern of variations, for which rotational periods are unknown. Here, we investigate the effect of metallicity on the detectability of rotation periods. This is done by synthesising light curves of hypothetical stars, which are identical to our Sun, with the exception of the metallicity. These light curves are then used as an input to the period determination algorithms. We find that the success rate for recovering the rotation signal has a minimum close to the solar metallicity value. This can be explained by the compensation effect of facular and spot contributions. In addition, selecting solar-like stars with near-solar effective temperature, near solar photometric variability, and with metallicity between M/H = -0.35 and M/H = 0.35 from the Kepler sample, we analyse the fraction of stars for which rotational periods have been detected as a function of metallicity. In agreement with our theoretical estimate we found a local minimum for the detection fraction close to the solar metallicity. We further report rotation periods of 87~solar-like Kepler stars for the first time.

preprint2020arXiv

No universal connection between the vertical magnetic field and the umbra-penumbra boundary in sunspots

Context. It has been reported that the boundary between the umbra and the penumbra of sunspots occurs at a canonical value of the strength of the vertical magnetic field, independently of the size of the spot. This critical field strength is interpreted as to be the threshold for the onset of magnetoconvection. Aims. Here we investigate the reasons why this criterion, also called the Jurčák criterion in the literature, does not always identify the boundary between umbra and penumbra. Methods. We perform a statistical analysis of 23 sunspots observed with Hinode/SOT. We compare the properties of the continuum intensity and the vertical magnetic field between filaments and spines and how they vary between spots of different sizes. Results. We find that the inner boundary of the penumbra is not related to a universal value of the vertical magnetic field. The properties of spines and filaments vary between spots of different sizes. Both components are darker in larger spots and the spines exhibit stronger vertical magnetic field. These variations of the properties of filaments and spines with spot size are also the reason for the reported invariance of the averaged vertical magnetic field at 50% of the mean continuum intensity. Conclusions. The formation of filaments and the onset of magnetoconvection are not related to a canonical value of the strength of the vertical magnetic field. Such a seemingly unique magnetic field strength is rather an effect of the filling factor of spines and penumbral filaments.

preprint2020arXiv

On the Magnetic Nature of an Exploding Granule as Revealed by Sunrise/IMaX

We study the photospheric evolution of an exploding granule observed in the quiet Sun at high spatial ($0.3^{\prime\prime}$) and temporal (31.5 s) resolution by the imaging magnetograph Sunrise/IMaX in June 2009. These observations show that the exploding granule is cospatial to a magnetic flux emergence event occurring at mesogranular scale (up to 12 Mm$^{2}$ area). Using a modified version of the SIR code for inverting the IMaX spectropolarimetric measurements, we obtain information about the magnetic configuration of this photospheric feature. In particular, we find evidence of highly inclined emerging fields in the structure, carrying a magnetic flux content up to $4 \times 10^{18}$ Mx. The balance between gas and magnetic pressure in the region of flux emergence, compared with a very quiet region of the Sun, indicates that the additional pressure carried by the emerging flux increases by about 5% the total pressure and appears to allow the granulation to be modified, as predicted by numerical simulations. The overall characteristics suggest that a multi-polar structure emerges into the photosphere, resembling an almost horizontal flux sheet. This seems to be associated with exploding granules. Finally, we discuss the origin of such flux emergence events.

preprint2020arXiv

Power spectra of solar brightness variations at different inclinations

Magnetic features on the surfaces of cool stars cause variations of their brightness. Such variations have been extensively studied for the Sun. Recent planet-hunting space telescopes allowed measuring brightness variations in hundred thousands of other stars. The new data posed the question of how typical is the Sun as a variable star. Putting solar variability into the stellar context suffers, however, from the bias of solar observations being made from its near-equatorial plane, whereas stars are observed at all possible inclinations. We model solar brightness variations at timescales from days to years as they would be observed at different inclinations. In particular, we consider the effect of the inclination on the power spectrum of solar brightness variations. The variations are calculated in several passbands routinely used for stellar measurements. We employ the Surface Flux Transport Model (SFTM) to simulate the time-dependent spatial distribution of magnetic features on both near- and far-sides of the Sun. This distribution is then used to calculate solar brightness variations following the SATIRE (Spectral And Total Irradiance REconstruction) approach. We have quantified the effect of the inclination on solar brightness variability at timescales down to a day. Thus, our results allow making solar brightness records directly comparable to those obtained by the planet-hunting space telescopes. Furthermore, we decompose solar brightness variations into the components originating from the solar rotation and from the evolution of magnetic features.

preprint2020arXiv

Power spectrum of turbulent convection in the solar photosphere

The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law.

preprint2020arXiv

Reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars"

This is our reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars" by T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gomez. We hope that both the comment and our reply lead to fruitful discussions which of the two presented scenarios is more likely.

preprint2020arXiv

Rotation periods from the inflection point in the power spectrum of stellar brightness variations: II. The Sun

Young and active stars generally have regular, almost sinusoidal, patterns of variability attributed to their rotation, while the majority of older and less active stars, including the Sun, have more complex and non-regular light-curves which do not have clear rotational-modulation signals. Consequently, the rotation periods have been successfully determined only for a small fraction of the Sun-like stars observed by transit-based planet-hunting missions, such as CoRoT, Kepler, and TESS. This suggests that only a small fraction of such systems have been properly identified as solar-like analogs. We apply a new method for determining rotation periods of low-activity stars, like the Sun. The method is based on calculating the gradient of the power spectrum (GPS) of stellar brightness variations and identifying a tell-tale inflection point in the spectrum. The rotation frequency is then proportional to the frequency of that inflection point. In this paper test this GPS method against available photometric records of the Sun. We apply GPS, autocorrelation functions, Lomb-Scargle periodograms, and wavelet analyses to the total solar irradiance (TSI) time series obtained from the Total Irradiance Monitor (TIM) on the Solar Radiation and Climate Experiment (SORCE) and the Variability of solar IRradiance and Gravity Oscillations (VIRGO) experiment on the SOlar and Heliospheric Observatory (SoHO) missions. We analyse the performance of all methods at various levels of solar activity. We show that the GPS method returns accurate values of solar rotation independently of the level of solar activity. In particular, it performs well during periods of high solar activity, when TSI variability displays an irregular pattern and other methods fail. Our results suggest that the GPS method can successfully determine the rotational periods of stars with both regular and non-regular light-curves.

preprint2020arXiv

Simulations Show that Vortex Flows could Heat the Chromosphere in Solar Plage

The relationship between vortex flows at different spatial scales and their contribution to the energy balance in the chromosphere is not yet fully understood. We perform three-dimensional (3D) radiation-magnetohydrodynamic (MHD) simulations of a unipolar solar plage region at a spatial resolution of 10 km using the MURaM code. We use the swirling-strength criterion that mainly detects the smallest vortices present in the simulation data. We additionally degrade our simulation data to smooth-out the smaller vortices, so that also the vortices at larger spatial scales can be detected. Vortex flows at various spatial scales are found in our simulation data for different effective spatial resolutions. We conclude that the observed large vortices are likely clusters of much smaller ones that are not yet resolved by observations. We show that the vertical Poynting flux decreases rapidly with reduced effective spatial resolutions and is predominantly carried by the horizontal plasma motions rather than vertical flows. Since the small-scale horizontal motions or the smaller vortices carry most of the energy, the energy transported by vortices deduced from low resolution data is grossly underestimated. In full resolution simulation data, the Poynting flux contribution due to vortices is more than adequate to compensate for the radiative losses in plage, indicating their importance for chromospheric heating.

preprint2020arXiv

The influence of NLTE effects in Fe I lines on an inverted atmosphere I. 6301 A and 6302 A lines formed in 1D NLTE

Ultraviolet over-ionisation of iron atoms in the solar atmosphere leads to deviations in their level populations from the Saha-Boltzmann statistics. This causes their line profiles to form in Non-Local Thermodynamic Equilibrium (NLTE) conditions. While inverting such profiles to determine atmospheric parameters, the NLTE effects are often neglected and deviations from LTE are compensated for by tweaking other quantities. We investigate how the routinely employed LTE inversion of iron lines formed in NLTE under- or over-estimates atmospheric quantities such as temperature (T), line-of-sight velocity (v_LOS), magnetic field strength (B) and inclination (gamma) while the previous papers have focused mainly on T. We synthesize the Stokes profiles of Fe I 6301.5 A and 6302.5 A lines in both LTE and NLTE using a snapshot of a 3D MHD simulation. The profiles are then inverted in LTE. By considering the atmosphere inferred from inversion of LTE profiles to be the fiducial model, we compare atmosphere from the inversion of NLTE profiles with it. Any differences observed are attributed to NLTE effects. Neglecting the NLTE effects introduces errors in the inverted atmosphere. While the errors in T can go up to 13%, in v_LOS and B the errors can be as high as 50% or more. We find these errors to be present at all three inversion nodes. Importantly, they survive degradation from spatial averaging of the profiles. We give an overview of how the neglect of NLTE effects influences the values of T, v_LOS, B and gamma determined by inverting Fe I 6300 A line pair, as observed, e.g., by Hinode. Errors are found at the sites of granules, intergranular lanes, magnetic elements, basically in every region susceptible to NLTE effects. For an accurate determination of atmospheric quantities and their stratification, it is therefore important to take account of NLTE effects.

preprint2019arXiv

Inflection point in the power spectrum of stellar brightness variations. I. The model

Context. Considerable effort has been put into using light curves observed by space telescopes such as CoRoT, Kepler and TESS for determining stellar rotation periods. While rotation periods of active stars can be reliably determined, the light curves of many older and less active stars (e.g. stars similar to the Sun) are quite irregular, which hampers determination of their periods. Aims. We examine the factors causing the irregularities in stellar brightness variations and develop a method for determining rotation periods of low activity stars with irregular light curves. Methods. We extend the Spectral And Total Irradiance Reconstruction (SATIRE) approach for modelling solar brightness variations to Sun-like stars. We calculate the power spectra of stellar brightness variations for various combinations of parameters defining the surface configuration and evolution of stellar magnetic features. Results. The short lifetime of spots in comparison to the stellar rotation period as well as the interplay between spot and facular contributions to brightness variations of stars with near solar activity cause irregularities in their light curves. The power spectra of such stars often lack a peak associated with the rotation period. Nevertheless, the rotation period can still be determined by measuring the period where the concavity of the power spectrum plotted in the log-log scale changes sign, i.e. by identifying the position of the inflection point. Conclusions. The inflection point of the (log-log) power spectrum is found to be a new diagnostic for stellar rotation periods that is shown to work even in cases where the power spectrum shows no peak at the rotation rate.

preprint2016arXiv

Active region fine structure observed at 0.08 arcsec resolution

The various mechanisms of magneto-convective energy transport determines the structure of sunspots and active regions. We characterise the appearance of light bridges and other fine-structure details and elaborate on their magneto-convective nature. We present speckle-reconstructed images taken with the broad-band imager at the 1.5 m GREGOR telescope in the 486nm and 589nm bands. We estimate the spatial resolution from the noise characteristics of the image bursts and obtain 0.08" at 589nm. We describe structure details in individual best images as well as the temporal evolution of selected features. We find branched dark lanes extending along thin (~1") light bridges in sunspots at various heliocentric angles. In thick (~2") light bridges the branches are disconnected from the central lane and have a Y shape with a bright grain toward the umbra. The images reveal that light bridges exist on varying intensity levels and that their small-scale features evolve on timescales of minutes. Faint light bridges show dark lanes outlined by the surrounding bright features. Dark lanes are very common and are also found in the boundary of pores. They have a characteristic width of 0.1" or smaller. Intergranular dark lanes of that width are seen in active region granulation. While central dark lanes in thin light bridges are elevated and associated with an density increase above upflows, the dark lane branches correspond to locations of downflows and are depressed relative to the adjacent bright plasma. Thick light bridges with central dark lanes show no projection effect. They have a flat elevated plateau that falls off steeply at the umbral boundary. There, Y-shaped filaments form as they do in the inner penumbra. This indicates the presence of inclined magnetic fields, meaning that the umbral magnetic field is wrapped around the convective light bridge.

preprint2016arXiv

Are solar brightness variations faculae- or spot-dominated?

Regular spaceborne measurements have revealed that solar brightness varies on multiple timescales, variations on timescales greater than a day being attributed to surface magnetic field. Independently, ground-based and spaceborne measurements suggest that Sun-like stars show a similar, but significantly broader pattern of photometric variability. To understand whether the broader pattern of stellar variations is consistent with the solar paradigm we assess relative contributions of faculae and spots to solar magnetically-driven brightness variability. We investigate how the solar brightness variability as well as its facular and spot contributions depend on the wavelength, timescale of variability, and position of the observer relative to the ecliptic plane. We perform calculations with the SATIRE model, which returns solar brightness with daily cadence from solar disc area coverages of various magnetic features. Moving the observer away from the ecliptic plane increases the amplitude of 11-year variability as it would be seen in Strömgren $(b+y)/2$ photometry, but decreases the amplitude of the rotational brightness variations as it would appear in Kepler and CoRoT passbands. The spot and facular contributions to the 11-year solar variability in the Strömgren $(b+y)/2$ photometry almost fully compensate each other so that the Sun appears anomalously quiet with respect to its stellar cohort. Such a compensation does not occur on the rotational timescale. The rotational solar brightness variability as it would appear in Kepler and CoRoT passband from the ecliptic plane is spot-dominated but the relative contribution of faculae increases for out-of-ecliptic viewing so that the apparent brightness variations are faculae-dominated for inclinations less than about $i=45^{\circ}$.

preprint2016arXiv

Deep probing of the photospheric sunspot penumbra: no evidence for magnetic field-free gaps

Some models for the topology of the magnetic field in sunspot penumbrae predict the existence of field-free or dynamically weak-field regions in the deep Photosphere. To confirm or rule out the existence of weak-field regions in the deepest photospheric layers of the penumbra. The magnetic field at $\logτ_5=0$ is investigated by means of inversions of spectropolarimetric data of two different sunspots located very close to disk center with a spatial resolution of approximately 0.4-0.45 arcsec. The data have been recorded using the GRIS instrument attached to the 1.5-meters GREGOR solar telescope at El Teide observatory. It includes three Fe I lines around 1565 nm, whose sensitivity to the magnetic field peaks at half a pressure-scale-height deeper than the sensitivity of the widely used Fe I spectral line pair at 630 nm. Prior to the inversion, the data is corrected for the effects of scattered light using a deconvolution method with several point spread functions. At $\logτ_5=0$ we find no evidence for the existence of regions with dynamically weak ($B<500$ Gauss) magnetic fields in sunspot penumbrae. This result is much more reliable than previous investigations done with Fe I lines at 630 nm. Moreover, the result is independent of the number of nodes employed in the inversion, and also independent of the point spread function used to deconvolve the data, and does not depend on the amount of straylight (i.e. wide-angle scattered light) considered.

preprint2016arXiv

Diagnostics of Coronal Magnetic Fields Through the Hanle Effect in UV and IR Lines

The plasma thermodynamics in the solar upper atmosphere, particularly in the corona, are dominated by the magnetic field, which controls the flow and dissipation of energy. The relative lack of knowledge of the coronal vector magnetic field is a major handicap for progress in coronal physics. This makes the development of measurement methods of coronal magnetic fields a high priority in solar physics. The Hanle effect in the UV and IR spectral lines is a largely unexplored diagnostic. We use magnetohydrodynamic (MHD) simulations to study the magnitude of the signal to be expected for typical coronal magnetic fields for selected spectral lines in the UV and IR wavelength ranges, namely the H I Ly-$α$ and the He I 10830 Å lines. We show that the selected lines are useful for reliable diagnosis of coronal magnetic fields. The results show that the combination of polarization measurements of spectral lines with different sensitivities to the Hanle effect may be most appropriate for deducing coronal magnetic properties from future observations.

preprint2016arXiv

Fitting peculiar spectral profiles in He I 10830 Å absorption features

The new generation of solar instruments provides better spectral, spatial, and temporal resolution for a better understanding of the physical processes that take place on the Sun. Multiple-component profiles are more commonly observed with these instruments. Particularly, the He I 10830 Å triplet presents such peculiar spectral profiles, which give information on the velocity and magnetic fine structure of the upper chromosphere. The purpose of this investigation is to describe a technique to efficiently fit the two blended components of the He I 10830 Å triplet, which are commonly observed when two atmospheric components are located within the same resolution element. The observations used in this study were taken on 2015 April 17 with the very fast spectroscopic mode of the GREGOR Infrared Spectrograph (GRIS) attached to the 1.5-meter GREGOR solar telescope, located at the Observatorio del Teide, Tenerife, Spain. We apply a double-Lorentzian fitting technique using Levenberg-Marquardt least-squares minimization. This technique is very simple and much faster than inversion codes. Line-of-sight Doppler velocities can be inferred for a whole map of pixels within just a few minutes. Our results show sub- and supersonic downflow velocities of up to 32 km/s for the fast component in the vicinity of footpoints of filamentary structures. The slow component presents velocities close to rest.

preprint2016arXiv

Flows in and around active region NOAA12118 observed with the GREGOR solar telescope and SDO/HMI

Accurate measurements of magnetic and velocity fields in and around solar active regions are key to unlocking the mysteries of the formation and the decay of sunspots. High spatial resolution image and spectral sequences with a high cadence obtained with the GREGOR solar telescope give us an opportunity to scrutinize 3-D flow fields with local correlation tracking and imaging spectroscopy. We present GREGOR early science data acquired in 2014 July - August with the GREGOR Fabry-Pérot Interferometer and the Blue Imaging Channel. Time-series of blue continuum (? 450.6 nm) images of the small active region NOAA 12118 were restored with the speckle masking technique to derive horizontal proper motions and to track the evolution of morphological changes. In addition, high-resolution observations are discussed in the context of synoptic data from the Solar Dynamics Observatory.

preprint2016arXiv

Formation of a solar Ha filament from orphan penumbrae

The formation of an Ha filament in active region (AR) 10953 is described. Observations from the Solar Optical Telescope (SOT) aboard the Hinode satellite starting on 27th April 2007 until 1st May 2007 were analysed. 20 scans of the 6302A Fe I line pair recorded by SOT/SP were inverted using the SPINOR code. The inversions were analysed together with SOT/BFI G-band and Ca II H and SOT/NFI Ha observations. Following the disappearance of an initial Ha filament aligned along the polarity inversion line (PIL) of the AR, a new Ha filament formed in its place some 20 hours later, which remained stable for at least 1.5 days. The creation of the new Ha filament was driven by the ascent of horizontal magnetic fields from the photosphere into the chromosphere at three separate locations along the PIL. The magnetic fields at two of these locations were situated directly underneath the initial Ha filament and formed orphan penumbrae already aligned along the Ha filament channel. The 700 G orphan penumbrae were stable and trapped in the photosphere until the disappearance of the overlying initial Ha filament, after which they started to ascend into the chromosphere at 10pm5 m/s. Each ascent was associated with a simultaneous magnetic flux reduction of up to 50% in the photosphere. The ascended orphan penumbrae formed dark 'seed' structures in Ha in parallel with the PIL, which elongated and merged to form an Ha filament. The filament channel featured horizontal magnetic fields of on average 260 G at log(tau)=-2 suspended above the nearly field-free lower photosphere. The fields took on an 'inverse' configuration at log(tau)=-2 suggesting a flux rope for the new Ha filament. The orphan penumbral fields ascend into the chromosphere 9-24 hours before the Ha filament is fully formed. The destruction of the initial Ha filament was likely caused by the flux emergence at the third location along the PIL.

preprint2016arXiv

Formation of the penumbra and start of the Evershed flow

We studied the variations of line-of-sight photospheric plasma flows during the formation phase of the penumbra around a pore in Active Region NOAA 11490. We used a high spatial, spectral, and temporal resolution data set acquired by the Interferometric BIdimensional Spectrometer (IBIS) operating at the NSO/Dunn Solar Telescope as well as data taken by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory satellite (SDO/HMI). Before the penumbra formed we observed a redshift of the spectral line in the inner part of the annular zone surrounding the pore as well as a blueshift of material associated with opposite magnetic polarity further away from the pore. We found that the onset of the classical Evershed flow occurs in a very short time scale -- 1-3 hours -- while the penumbra is forming. During the same time interval we found changes in the magnetic field inclination in the penumbra, with the vertical field actually changing sign near the penumbral edge, while the total magnetic field showed a significant increase, about 400 G. To explain these and other observations related to the formation of the penumbra and the onset of the Evershed flow we propose a scenario in which the penumbra is formed by magnetic flux dragged down from the canopy surrounding the initial pore. The Evershed flow starts when the sinking magnetic field dips below the solar surface and magnetoconvection sets in.

preprint2016arXiv

High-resolution, high-sensitivity, ground-based solar spectropolarimetry with a new fast imaging polarimeter

Context. Remote sensing of weak and small-scale solar magnetic fields is of utmost relevance for a number of important open questions in solar physics. This requires the acquisition of spectropolarimetric data with high spatial resolution (0.1 arcsec) and low noise (1e-3 to 1e-5 of the continuum intensity). The main limitations to obtain these measurements from the ground, are the degradation of the image resolution produced by atmospheric seeing and the seeing-induced crosstalk (SIC). Aims. We introduce the prototype of the Fast Solar Polarimeter (FSP), a new ground-based, high-cadence polarimeter that tackles the above-mentioned limitations by producing data that are optimally suited for the application of post-facto image restoration, and by operating at a modulation frequency of 100 Hz to reduce SIC. Results. The pnCCD camera reaches 400 fps while keeping a high duty cycle (98.6 %) and very low noise (4.94 erms). The modulator is optimized to have high (> 80%) total polarimetric efficiency in the visible spectral range. This allows FSP to acquire 100 photon-noise-limited, full-Stokes measurements per second. We found that the seeing induced signals present in narrow-band, non-modulated, quiet-sun measurements are (a) lower than the noise (7e-5) after integrating 7.66 min, (b) lower than the noise (2.3e-4) after integrating 1.16 min and (c) slightly above the noise (4e-3) after restoring case (b) by means of a multi-object multi-frame blind deconvolution. In addition, we demonstrate that by using only narrow-band images (with low SNR of 13.9) of an active region, we can obtain one complete set of high-quality restored measurements about every 2 s.

preprint2016arXiv

Horizontal flow fields in and around a small active region-- The transition period between flux emergence and decay

Aims. Combining high-resolution and synoptic observations aims to provide a comprehensive description of flux emergence at photospheric level and of the growth process that eventually leads to a mature active region. Methods. Small active region NOAA 12118 was observed on 2014 July 18 with the 1.5-meter GREGOR solar telescope on 2014 July 18. High-resolution time-series of blue continuum and G-band images acquired in the blue imaging channel (BIC) of the GREGOR Fabry-Pérot Interferometer (GFPI) were complemented by LOS magnetograms and continuum images obtained with the HMI onboard the SDO. Horizontal proper motions and horizontal plasma velocities were computed with local correlation tracking (LCT) and the differential affine velocity estimator, respectively. Morphological image processing was employed to measure the photometric/magnetic area, magnetic flux, and the separation profile of the EFR during its evolution. Results. The computed growth rates for photometric area, magnetic area, and magnetic flux are about twice as high as the respective decay rates. The space-time diagram using HMI magnetograms of five days traces a leaf-like structure, which is determined by the initial separation of the two polarities, a rapid expansion phase, a time when the spread stalls, and a period when the region slowly shrinks again. The separation rate of 0.26 km\s is highest in the initial stage, and it decreases when the separation comes to a halt. Horizontal plasma velocities computed at four evolutionary stages indicate a changing pattern of inflows. In LCT maps we find persistent flow patterns such as outward motions in the outer part of the two major pores, a diverging feature near the trailing pore marking the site of upwelling plasma and flux emergence, and low velocities in the interior of pores. We detected many elongated rapidly expanding granules between the two major polarities.

preprint2016arXiv

Magnetic fields of opposite polarity in sunspot penumbrae

Context. A significant part of the penumbral magnetic field returns below the surface in the very deep photosphere. For lines in the visible, a large portion of this return field can only be detected indirectly by studying its imprints on strongly asymmetric and three-lobed Stokes V profiles. Infrared lines probe a narrow layer in the very deep photosphere, providing the possibility of directly measuring the orientation of magnetic fields close to the solar surface. Aims. We study the topology of the penumbral magnetic field in the lower photosphere, focusing on regions where it returns below the surface. Methods. We analyzed 71 spectropolarimetric datasets from Hinode and from the GREGOR infrared spectrograph. We inferred the quality and polarimetric accuracy of the infrared data after applying several reduction steps. Techniques of spectral inversion and forward synthesis were used to test the detection algorithm. We compared the morphology and the fractional penumbral area covered by reversed-polarity and three-lobed Stokes V profiles for sunspots at disk center. We determined the amount of reversed-polarity and three-lobed Stokes V profiles in visible and infrared data of sunspots at various heliocentric angles. From the results, we computed center-to-limb variation curves, which were interpreted in the context of existing penumbral models. Results. Observations in visible and near-infrared spectral lines yield a significant difference in the penumbral area covered by magnetic fields of opposite polarity. In the infrared, the number of reversed-polarity Stokes V profiles is smaller by a factor of two than in the visible. For three-lobed Stokes V profiles the numbers differ by up to an order of magnitude.

preprint2016arXiv

Probing deep photospheric layers of the quiet Sun with high magnetic sensitivity

Context. Investigations of the magnetism of the quiet Sun are hindered by extremely weak polarization signals in Fraunhofer spectral lines. Photon noise, straylight, and the systematically different sensitivity of the Zeeman effect to longitudinal and transversal magnetic fields result in controversial results in terms of the strength and angular distribution of the magnetic field vector. Aims. The information content of Stokes measurements close to the diffraction limit of the 1.5 m GREGOR telescope is analyzed. We took the effects of spatial straylight and photon noise into account. Methods. Highly sensitive full Stokes measurements of a quiet-Sun region at disk center in the deep photospheric Fe I lines in the 1.56 μm region were obtained with the infrared spectropolarimeter GRIS at the GREGOR telescope. Noise statistics and Stokes V asymmetries were analyzed and compared to a similar data set of the Hinode spectropolarimeter (SOT/SP). Simple diagnostics based directly on the shape and strength of the profiles were applied to the GRIS data. We made use of the magnetic line ratio technique, which was tested against MHD simulations. Results. About 80% of the GRIS spectra of a very quiet solar region show polarimetric signals above a 3σ level. Area and amplitude asymmetries agree well with small-scale surface dynamo MHD simulations. The magnetic line ratio analysis reveals ubiquitous magnetic regions in the ten to hundred Gauss range with some concentrations of kilo-Gauss fields. Conclusions. The GRIS spectropolarimetric data at a spatial resolution of 0.40" are so far unique in the combination of high spatial resolution scans and high magnetic field sensitivity. Nevertheless, the unavoidable effect of spatial straylight and the resulting dilution of the weak Stokes profiles means that inversion techniques still bear a high risk of misinterpretating the data.

preprint2016arXiv

Reconstruction of spectral solar irradiance since 1700 from simulated magnetograms

We present a reconstruction of the spectral solar irradiance since 1700 using the SATIRE-T2 (Spectral And Total Irradiance REconstructions for the Telescope era version 2) model. This model uses as input magnetograms simulated with a surface flux transport model fed with semi-synthetic records of emerging sunspot groups. We used statistical relationships between the properties of sunspot group emergence, such as the latitude, area, and tilt angle, and the sunspot cycle strength and phase to produce semi-synthetic sunspot group records starting in the year 1700. The semisynthetic records are fed into a surface flux transport model to obtain daily simulated magnetograms that map the distribution of the magnetic flux in active regions (sunspots and faculae) and their decay products on the solar surface. The magnetic flux emerging in ephemeral regions is accounted for separately based on the concept of extended cycles whose length and amplitude are linked to those of the sunspot cycles through the sunspot number. The magnetic flux in each surface component (sunspots, faculae and network, and ephemeral regions) was used to compute the spectral and total solar irradiance between the years 1700 and 2009. This reconstruction is aimed at timescales of months or longer although the model returns daily values. We found that SATIRE-T2, besides reproducing other relevant observations such as the total magnetic flux, reconstructs the total solar irradiance (TSI) on timescales of months or longer in good agreement with the PMOD composite of observations, as well as with the reconstruction starting in 1878 based on the RGO-SOON data. The model predicts an increase in the TSI of 1.2[+0.2, -0.3] Wm-2 between 1700 and the present. The spectral irradiance reconstruction is in good agreement with the UARS/SUSIM measurements as well as the Lyman-alpha composite.

preprint2016arXiv

Search for chameleons with CAST

In this work we present a search for (solar) chameleons with the CERN Axion Solar Telescope (CAST). This novel experimental technique, in the field of dark energy research, exploits both the chameleon coupling to matter ($β_{\rm m}$) and to photons ($β_γ$) via the Primakoff effect. By reducing the X-ray detection energy threshold used for axions from 1$\,$keV to 400$\,$eV CAST became sensitive to the converted solar chameleon spectrum which peaks around 600$\,$eV. Even though we have not observed any excess above background, we can provide a 95% C.L. limit for the coupling strength of chameleons to photons of $β_γ\!\lesssim\!10^{11}$ for $1<β_{\rm m}<10^6$.

preprint2016arXiv

Slipping reconnection in a solar flare observed in high resolution with the GREGOR solar telescope

A small flare ribbon above a sunspot umbra in active region 12205 was observed on November 7, 2014, at 12:00 UT in the blue imaging channel of the 1.5 m GREGOR telescope, using a 1 A Ca II H interference filter. Context observations from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO), the Solar Optical Telescope (SOT) onboard Hinode, and the Interface Region Imaging Spectrograph (IRIS) show that this ribbon is part of a larger one that extends through the neighboring positive polarities and also participates in several other flares within the active region. We reconstructed a time series of 140 seconds of Ca II H images by means of the multiframe blind deconvolution method, which resulted in spatial and temporal resolutions of 0.1 arcsec and 1 s. Light curves and horizontal velocities of small-scale bright knots in the observed flare ribbon were measured. Some knots are stationary, but three move along the ribbon with speeds of 7-11 km/s. Two of them move in the opposite direction and exhibit highly correlated intensity changes, which provides evidence of a slipping reconnection at small spatial scales.

preprint2016arXiv

Solar science with the Atacama Large Millimeter/submillimeter Array - A new view of our Sun

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

preprint2016arXiv

Spectropolarimetric observations of an arch filament system with the GREGOR solar telescope

Arch filament systems occur in active sunspot groups, where a fibril structure connects areas of opposite magnetic polarity, in contrast to active region filaments that follow the polarity inversion line. We used the GREGOR Infrared Spectrograph (GRIS) to obtain the full Stokes vector in the spectral lines Si I 1082.7 nm, He I 1083.0 nm, and Ca I 1083.9 nm. We focus on the near-infrared calcium line to investigate the photospheric magnetic field and velocities, and use the line core intensities and velocities of the helium line to study the chromospheric plasma. The individual fibrils of the arch filament system connect the sunspot with patches of magnetic polarity opposite to that of the spot. These patches do not necessarily coincide with pores, where the magnetic field is strongest. Instead, areas are preferred not far from the polarity inversion line. These areas exhibit photospheric downflows of moderate velocity, but significantly higher downflows of up to 30 km/s in the chromospheric helium line. Our findings can be explained with new emerging flux where the matter flows downward along the fieldlines of rising flux tubes, in agreement with earlier results.

preprint2016arXiv

Three-dimensional structure of a sunspot light bridge

Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 Å and Ca I 10839 Å lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$β$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution.

preprint2016arXiv

Upper Chromospheric Magnetic Field of a Sunspot Penumbra: Observations of Fine Structure

The fine-structure of magnetic field of a sunspot penumbra in the upper chromosphere is to be explored and compared to that in the photosphere. High spatial resolution spectropolarimetric observations were recorded with the 1.5-meter GREGOR telescope using the GREGOR Infrared Spectrograph (GRIS). The observed spectral domain includes the upper chromospheric He I triplet at 1083.0 nm and the photospheric Si I 1082.7 nm and Ca I 1083.3 nm spectral lines. The upper chromospheric magnetic field is obtained by inverting the He I triplet assuming a Milne-Eddington type model atmosphere. A height dependent inversion was applied to the Si I 1082.7 nm and Ca I 1083.3 nm lines to obtain the photospheric magnetic field. We find that the inclination of the magnetic field shows variations in the azimuthal direction both in the photosphere, but also in the upper chromosphere. The chromospheric variations remarkably well coincide with the variations in the inclination of the photospheric field and resemble the well-known spine and inter-spine structure in the photospheric layers of penumbrae. The typical peak-to-peak variations in the inclination of the magnetic field in the upper chromosphere is found to be 10-15 degree, i.e., roughly half the variation in the photosphere. In contrast, the magnetic field strength of the observed penumbra does not show variations on small spatial scales in the upper chromosphere. Thanks to the high spatial resolution observations possible with the GREGOR telescope at 1.08 microns, we find that the prominent small-scale fluctuations in the magnetic field inclination, which are a salient part of the property of sunspot penumbral photospheres, also persist in the chromosphere, although at somewhat reduced amplitudes. Such a complex magnetic configuration may facilitate penumbral chromospheric dynamic phenomena, such as penumbral micro-jets or transient bright dots.

preprint2015arXiv

A New Calibrated Sunspot Group Series Since 1749: Statistics of Active Day Fractions

Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a novel method to calibrate the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900-1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds [S_S] defined such that the observer is assumed to miss all of the groups with an area smaller than S_S and report all the groups larger than S_S. Next, using a Monte-Carlo method we construct, from the reference data set, a correction matrix for each observer. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot group records for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series displays secular minima around 1800 (Dalton minimum) and 1900 (Gleissberg minimum), as well as the Modern grand maximum of activity in the second half of the 20th century. The uniqueness of the grand maximum is confirmed for the last 250 years. It is shown that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the 18th and 19th centuries in some reconstructions.

preprint2015arXiv

Flarelike brightenings of active region loops observed with SUMER

Coronal loops on the east limb of the Sun were observed by SUMER on SOHO for several days. Small flare-like brightenings are detected very frequently in the hot flare line Fe~{\small XIX}. We find that the relatively intense events are in good coincidence with the transient brightenings seen by Yohkoh/SXT. A statistical analysis shows that these brightenings have durations of 5-84 min and extensions along the slit of 2-67 Mm. The integrated energy observed in Fe~{\small XIX} for each event is in the range of $3\times10^{18}-5\times10^{23}$ ergs, and the estimated thermal energy ranges from $10^{26}-10^{29}$ ergs. Application of the statistical method proposed by Parnell \& Jupp (2000) yields a value of 1.5 to 1.8 for the index of a power law relation between the frequency of the events and the radiated energy in Fe~{\small XIX}, and a value of 1.7 to 1.8 for the index of the frequency distribution of the thermal energy in the energy range $>10^{27}$ ergs. We examine the possibility that these small brightenings give a big contribution to heating of the active region corona.

preprint2015arXiv

Magnetic reconnection as a source of jets from a penumbral intrusion into a sunspot umbra

We present the results of high resolution co-temporal and co-spatial photospheric and chromospheric observations of sunspot penumbral intrusions. The data was taken with the Swedish Solar Telescope (SST) on the Canary Islands. Time series of Ca\,II H images show a series of transient jets extending roughly 3000 km above a penumbral intrusion into the umbra. For most of the time series jets were seen along the whole length of the intruding bright filament. Some of these jets develop a clear $λ$-shaped structure, with a small loop appearing at their footpoint and lasting for around a minute. In the framework of earlier studies, the observed transient $λ$ shape of these jets strongly suggests that they are caused by magnetic reconnection between a curved arcade-like or flux-rope like field in the lower part of the penumbral intrusion and the more vertical umbral magnetic field forming a cusp-shaped structure above the penumbral intrusion.

preprint2015arXiv

Magneto-static modelling of the mixed plasma Beta solar atmosphere based on SUNRISE/IMaX data

Our aim is to model the 3D magnetic field structure of the upper solar atmosphere, including regions of non-negligible plasma beta. We use high-resolution photospheric magnetic field measurements from SUNRISE/IMaX as boundary condition for a magneto-static magnetic field model. The high resolution of IMaX allows us to resolve the interface region between photosphere and corona, but modelling this region is challenging for the following reasons. While the coronal magnetic field is thought to be force-free (the Lorentz-force vanishes), this is not the case in the mixed plasma $β$ environment in the photosphere and lower chromosphere. In our model, pressure gradients and gravity forces are taken self-consistently into account and compensate the non-vanishing Lorentz-force. Above a certain height (about 2 Mm) the non-magnetic forces become very weak and consequently the magnetic field becomes almost force-free. Here we apply a linear approach, where the electric current density consists of a superposition of a field-line parallel current and a current perpendicular to the Sun's gravity field. We illustrate the prospects and limitations of this approach and give an outlook for an extension towards a non-linear model.

preprint2015arXiv

New solar axion search in CAST with $^4$He filling

The CERN Axion Solar Telescope (CAST) searches for $a\toγ$ conversion in the 9 T magnetic field of a refurbished LHC test magnet that can be directed toward the Sun. Two parallel magnet bores can be filled with helium of adjustable pressure to match the X-ray refractive mass $m_γ$ to the axion search mass $m_a$. After the vacuum phase (2003--2004), which is optimal for $m_a\lesssim0.02$ eV, we used $^4$He in 2005--2007 to cover the mass range of 0.02--0.39 eV and $^3$He in 2009--2011 to scan from 0.39--1.17 eV. After improving the detectors and shielding, we returned to $^4$He in 2012 to investigate a narrow $m_a$ range around 0.2 eV ("candidate setting" of our earlier search) and 0.39--0.42 eV, the upper axion mass range reachable with $^4$He, to "cross the axion line" for the KSVZ model. We have improved the limit on the axion-photon coupling to $g_{aγ}< 1.47\times10^{-10} {\rm GeV}^{-1}$ (95% C.L.), depending on the pressure settings. Since 2013, we have returned to vacuum and aim for a significant increase in sensitivity.

preprint2015arXiv

Properties of solar plage from a spatially coupled inversion of Hinode SP data

The properties of magnetic fields forming an extended plage region in AR 10953 were investigated. Stokes spectra of the Fe I line pair at 6302 Årecorded by the spectropolarimeter aboard the Hinode satellite were inverted using the SPINOR code. The code performed a 2D spatially coupled inversion on the Stokes spectra, allowing the retrieval of gradients in optical depth within the atmosphere of each pixel, whilst accounting for the effects of the instrument's PSF. Consequently, no magnetic filling factor was needed. The inversion results reveal that plage is composed of magnetic flux concentrations (MFCs) with typical field strengths of 1520 G at log(τ)=-0.9 and inclinations of 10-15 degrees. The MFCs expand by forming magnetic canopies composed of weaker and more inclined magnetic fields. The expansion and average temperature stratification of isolated MFCs can be approximated well with an empirical plage thin flux-tube model. The highest temperatures of MFCs are located at their edges in all log(τ) layers. Whilst the plasma inside MFCs is nearly at rest, each is surrounded by a ring of downflows of on average 2.4 km/s at log(τ)=0 and peak velocities of up to 10 km/s, which are supersonic. The downflow ring of an MFC weakens and shifts outwards with height, tracing the MFC's expansion. Such downflow rings often harbour magnetic patches of opposite polarity to that of the main MFC with typical field strengths below 300 G at log(τ)=0. These opposite polarity patches are situated beneath the canopy of their main MFC. We found evidence of a strong broadening of the Stokes profiles in MFCs and particularly in the downflow rings surrounding MFCs (expressed by a microturbulence in the inversion). This indicates the presence of strong unresolved velocities. Larger magnetic structures such as sunspots cause the field of nearby MFCs to be more inclined.

preprint2015arXiv

Reconstruction of total and spectral solar irradiance from 1974 to 2013 based on KPVT, SoHO/MDI and SDO/HMI observations

Total and spectral solar irradiance are key parameters in the assessment of solar influence on changes in the Earth's climate. We present a reconstruction of daily solar irradiance obtained using the SATIRE-S model spanning 1974 to 2013 based on full-disc observations from the KPVT, SoHO/MDI and SDO/HMI. SATIRE-S ascribes variation in solar irradiance on timescales greater than a day to photospheric magnetism. The solar spectrum is reconstructed from the apparent surface coverage of bright magnetic features and sunspots in the daily data using the modelled intensity spectra of these magnetic structures. We cross-calibrated the various data sets, harmonizing the model input so as to yield a single consistent time series as the output. The model replicates 92% of the variability in the PMOD TSI composite including the secular decline between the 1996 and 2008 solar cycle minima. The model also reproduces most of the variability in observed Lyman-alpha irradiance and the Mg II index. The UV solar irradiance measurements from the UARS and SORCE missions are mutually consistent up to about 180 nm before they start to exhibit discrepant rotational and cyclical variability, indicative of unresolved instrumental effects. As a result, the agreement between model and measurement starts to deteriorate above this wavelength. As with earlier similar investigations, the reconstruction cannot reproduce the overall trends in SORCE/SIM SSI. We argue, from the lack of clear solar cycle modulation in the SIM record and the inconsistency between the total flux recorded by the instrument and TSI, that unaccounted instrumental trends are present. The reconstruction is consistent with observations from multiple sources, demonstrating its validity. It also provides further evidence that photospheric magnetism is the prime driver of variation in solar irradiance on timescales greater than a day.

preprint2015arXiv

Solar extreme ultraviolet variability of the quiet Sun

The last solar minimum has been unusually quiet compared to the previous minima (since space-based radiometric measurements are available). The Sun's magnetic flux was substantially lower during this minimum. Some studies also show that the total solar irradiance during the minimum after cycle 23 may have dropped below the values known from the two minima prior to that. For chromospheric and coronal radiation, the situation is less clear-cut. The Sun's 10.7\,cm flux shows a decrease of $\sim4\%$ during the solar minimum in 2008 compared to the previous minimum, but \ion{Ca}{II} K does not. Here we consider additional wavelengths in the extreme ultraviolet (EUV), specifically transitions of \ion{He}{I} at 584.3\,Å and \ion{O}{V} at 629.7\,Å, of which the CDS spectrometer aboard SOHO has been taking regular scans along the solar central meridian since 1996. We analysed this unique dataset to verify if and how the radiance distribution undergoes measurable variations between cycle minima. To achieve this aim we determined the radiance distribution of quiet areas around the Sun centre. Concentrating on the last two solar minima, we found out that there is very little variation in the radiance distribution of the chromospheric spectral line \ion{He}{I} between these minima. The same analysis shows a modest, although significant, 4\% variation in the radiance distribution of the transition region spectral line \ion{O}{V}. These results are comparable to those obtained by earlier studies employing other spectral features, and they confirm that chromospheric indices display a small variation, whereas in the TR a more significant reduction of the brighter features is visible.

preprint2015arXiv

Solar magnetoconvection and small-scale dynamo: Recent developments in observation and simulation

A number of observational and theoretical aspects of solar magnetoconvection are considered in this review. We discuss recent developments in our understanding of the small-scale structure of the magnetic field on the solar surface and its interaction with convective flows, which is at the centre of current research. Topics range from plage areas in active regions over the magnetic network shaped by supergranulation to the ubiquituous `turbulent' internetwork fields. On the theoretical side, we focus upon magnetic field generation by small-scale dynamo action.

preprint2015arXiv

The role of the Fraunhofer lines in solar brightness variability

The solar brightness varies on timescales from minutes to decades. A clear identification of the physical processes behind such variations is needed for developing and improving physics-based models of solar brightness variability and reconstructing solar brightness in the past. This is, in turn, important for better understanding the solar-terrestrial and solar-stellar connections. We estimate the relative contributions of the continuum, molecular, and atomic lines to the solar brightness variations on different timescales. Our approach is based on the assumption that variability of the solar brightness on timescales greater than a day is driven by the evolution of the solar surface magnetic field. We calculated the solar brightness variations employing the solar disc area coverage of magnetic features deduced from the MDI/SOHO observations. The brightness contrasts of magnetic features relative to the quiet Sun were calculated with a non-LTE radiative transfer code as functions of disc position and wavelength. By consecutive elimination of molecular and atomic lines from the radiative transfer calculations, we assessed the role of these lines in producing solar brightness variability. We show that the variations in Fraunhofer lines define the amplitude of the solar brightness variability on timescales greater than a day and even the phase of the total solar irradiance variability over the 11-year cycle. We also demonstrate that molecular lines make substantial contribution to solar brightness variability on the 11-year activity cycle and centennial timescales. In particular, our model indicates that roughly a quarter of the total solar irradiance variability over the 11-year cycle originates in molecular lines. The maximum of the absolute spectral brightness variability on timescales greater than a day is associated with the CN violet system between 380 and 390 nm.

preprint2015arXiv

UV solar irradiance in observations and the NRLSSI and SATIRE-S models

Total solar irradiance and UV spectral solar irradiance have been monitored since 1978 through a succession of space missions. This is accompanied by the development of models aimed at replicating solar irradiance by relating the variability to solar magnetic activity. The NRLSSI and SATIRE-S models provide the most comprehensive reconstructions of total and spectral solar irradiance over the period of satellite observation currently available. There is persistent controversy between the various measurements and models in terms of the wavelength dependence of the variation over the solar cycle, with repercussions on our understanding of the influence of UV solar irradiance variability on the stratosphere. We review the measurement and modelling of UV solar irradiance variability over the period of satellite observation. The SATIRE-S reconstruction is consistent with spectral solar irradiance observations where they are reliable. It is also supported by an independent, empirical reconstruction of UV spectral solar irradiance based on UARS/SUSIM measurements from an earlier study. The weaker solar cycle variability produced by NRLSSI between 300 and 400 nm is not evident in any available record. We show that although the method employed to construct NRLSSI is principally sound, reconstructed solar cycle variability is detrimentally affected by the uncertainty in the SSI observations it draws upon in the derivation. Based on our findings, we recommend, when choosing between the two models, the use of SATIRE-S for climate studies.

preprint2015arXiv

Variation of the Mn I 539.4 nm line with the solar cycle

As a part of the long-term program at Kitt Peak National Observatory (KPNO), the Mn I 539.4 nm line has been observed for nearly three solar cycles using the McMath telescope and the 13.5 m spectrograph in double-pass mode. These full-disk spectrophotometric observations revealed an unusually strong change of this line's parameters over the solar cycle. Optical pumping by the Mg II k line was originally proposed to explain these variations. More recent studies have proposed that this is not required and that the magnetic variability might explain it. Magnetic variability is also the mechanism that drives the changes in total solar irradiance variations (TSI). With this work we investigate this proposition quantitatively by using using the model SATIRE-S. We applied exactly the same model atmospheres and value of the free parameter as were used in previous solar irradiance reconstructions to now model the variation in the Mn I 539.4 nm line profile and in neighboring Fe I lines. We compared the results of the theoretical model with KPNO observations. Our result confirms that optical pumping of the Mn I 539.4 nm line by Mg II k is not the main cause of its solar cycle change. It also provides independent confirmation of solar irradiance models which are based on the assumption that irradiance variations are caused by the evolution of the solar surface magnetic flux. The result obtained here also supports the spectral irradiance variations computed by these models.

preprint2014arXiv

CAST solar axion search with 3^He buffer gas: Closing the hot dark matter gap

The CERN Axion Solar Telescope (CAST) has finished its search for solar axions with 3^He buffer gas, covering the search range 0.64 eV < m_a <1.17 eV. This closes the gap to the cosmological hot dark matter limit and actually overlaps with it. From the absence of excess X-rays when the magnet was pointing to the Sun we set a typical upper limit on the axion-photon coupling of g_ag < 3.3 x 10^{-10} GeV^{-1} at 95% CL, with the exact value depending on the pressure setting. Future direct solar axion searches will focus on increasing the sensitivity to smaller values of g_a, for example by the currently discussed next generation helioscope IAXO.

preprint2014arXiv

Comparison of solar photospheric bright points between SUNRISE observations and MHD simulations

Bright points (BPs) in the solar photosphere are radiative signatures of magnetic elements described by slender flux tubes located in the darker intergranular lanes. They contribute to the ultraviolet (UV) flux variations over the solar cycle and hence may influence the Earth's climate. Here we combine high-resolution UV and spectro-polarimetric observations of BPs by the SUNRISE observatory with 3D radiation MHD simulations. Full spectral line syntheses are performed with the MHD data and a careful degradation is applied to take into account all relevant instrumental effects of the observations. It is demonstrated that the MHD simulations reproduce the measured distributions of intensity at multiple wavelengths, line-of-sight velocity, spectral line width, and polarization degree rather well. Furthermore, the properties of observed BPs are compared with synthetic ones. These match also relatively well, except that the observations display a tail of large and strongly polarized BPs not found in the simulations. The higher spatial resolution of the simulations has a significant effect, leading to smaller and more numerous BPs. The observation that most BPs are weakly polarized is explained mainly by the spatial degradation, the stray light contamination, and the temperature sensitivity of the Fe I line at 5250.2 Å. The Stokes $V$ asymmetries of the BPs increase with the distance to their center in both observations and simulations, consistent with the classical picture of a production of the asymmetry in the canopy. This is the first time that this has been found also in the internetwork. Almost vertical kilo-Gauss fields are found for 98 % of the synthetic BPs. At the continuum formation height, the simulated BPs are on average 190 K hotter than the mean quiet Sun, their mean BP field strength is 1750 G, supporting the flux-tube paradigm to describe BPs.

preprint2014arXiv

Discriminant analysis of solar bright points and faculae II. Contrast and morphology analysis

Taken at a high spatial resolution of 0.1 arcsec, Bright Points (BPs) are found to coexist with faculae in images and the latter are often resolved as adjacent striations. Understanding the properties of these different features is fundamental to carrying out proxy magnetometry. To shed light on the relationship between BPs and faculae, we studied them separately after the application of a classification method, developed and described in a previous paper) on active region images at various heliocentric angles. In this Paper, we explore different aspects of the photometric properties of BPs and faculae, namely their G-band contrast profiles, their peak contrast in G-band and continuum, as well as morphological parameters. We find that: (1) the width of the contrast profiles of the classified BPs and faculae are consistent with studies of disk center BPs at and limb faculae, which indirectly confirms the validity of our classification, (2) the profiles of limb faculae are limbward skewed on average, while near disk center they exhibit both centerward and limbward skewnesses due to the distribution of orientations of the faculae, (3) the relation between the peak contrasts of BPs and faculae and their apparent area discloses a trend reminiscent of magnetogram studies. The skewness of facular profiles provides a novel constraint for 3D MHD models of faculae. As suggested by the asymmetry and orientation of their contrast profiles, faculae near disk center could be induced by inclined fields, while apparent BPs near the limb seem to be in fact small faculae misidentified. The apparent area of BPs and faculae could be possibly exploited for proxy magnetometry.

preprint2014arXiv

Inclinations of small quiet-Sun magnetic features based on a new geometric approach

High levels of horizontal magnetic flux have been reported in the quiet-Sun internetwork, often based on Stokes profile inversions. Here we introduce a new method for deducing the inclination of magnetic elements and use it to test magnetic field inclinations from inversions. We determine accurate positions of a set of small, bright magnetic elements in high spatial resolution images sampling different photospheric heights obtained by the Sunrise balloon-borne solar observatory. Together with estimates of the formation heights of the employed spectral bands, these provide us with the inclinations of the magnetic features. We also compute the magnetic inclination angle of the same magnetic features from the inversion of simultaneously recorded Stokes parameters. Our new, geometric method returns nearly vertical fields (average inclination of around 14 deg with a relatively narrow distribution having a standard deviation of 6 deg). In strong contrast to this, the traditionally used inversions give almost horizontal fields (average inclination of 75+-8 deg) for the same small magnetic features, whose linearly polarised Stokes profiles are adversely affected by noise. The almost vertical field of bright magnetic features from our geometric method is clearly incompatible with the nearly horizontal magnetic fields obtained from the inversions. This indicates that the amount of magnetic flux in horizontal fields deduced from inversions is overestimated in the presence of weak Stokes signals, in particular if Stokes Q and U are close to or under the noise level. By combining the proposed method with inversions we are not just improving the inclination, but also the field strength. This technique allows us to analyse features that are not reliably treated by inversions, thus greatly extending our capability to study the complete magnetic field of the quiet Sun.

preprint2014arXiv

Magnetic Flux Transport at the Solar Surface

After emerging to the solar surface, the Sun's magnetic field displays a complex and intricate evolution. The evolution of the surface field is important for several reasons. One is that the surface field, and its dynamics, sets the boundary condition for the coronal and heliospheric magnetic fields. Another is that the surface evolution gives us insight into the dynamo process. In particular, it plays an essential role in the Babcock-Leighton model of the solar dynamo. Describing this evolution is the aim of the surface flux transport model. The model starts from the emergence of magnetic bipoles. Thereafter, the model is based on the induction equation and the fact that after emergence the magnetic field is observed to evolve as if it were purely radial. The induction equation then describes how the surface flows -- differential rotation, meridional circulation, granular, supergranular flows, and active region inflows -- determine the evolution of the field (now taken to be purely radial). In this paper, we review the modeling of the various processes that determine the evolution of the surface field. We restrict our attention to their role in the surface flux transport model. We also discuss the success of the model and some of the results that have been obtained using this model.

preprint2014arXiv

Migration of Ca II H bright points in the internetwork

The migration of magnetic bright point-like features (MBP) in the lower solar atmosphere reflects the dispersal of magnetic flux as well as the horizontal flows of the atmospheric layer they are embedded in. We analyse trajectories of the proper motion of intrinsically magnetic, isolated internetwork Ca II H MBPs (mean lifetime 461 +- 9 s) to obtain their diffusivity behaviour. We use seeing-free high spatial and temporal resolution image sequences of quiet-Sun, disc-centre observations obtained in the Ca II H 3968 Å passband of the Sunrise Filter Imager (SuFI) onboard the Sunrise balloon-borne solar observatory. Small MBPs in the internetwork are automatically tracked. The trajectory of each MBP is then calculated and described by a diffusion index (γ) and a diffusion coefficient (D). We further explore the distribution of the diffusion indices with the help of a Monte Carlo simulation. We find γ = 1.69 +- 0.08 and D = 257 +- 32 km^2/s averaged over all MBPs. Trajectories of most MBPs are classified as super-diffusive, i.e., γ > 1, with the determined γ being to our knowledge the largest obtained so far. A direct correlation between D and time-scale (τ) determined from trajectories of all MBPs is also obtained. We discuss a simple scenario to explain the diffusivity of the observed, relatively short-lived MBPs while they migrate within a small area in a supergranule (i.e., an internetwork area). We show that the scatter in the γ values obtained for individual MBPs is due to their limited lifetimes. The super-diffusive MBPs can be well-described as random walkers (due to granular evolution and intergranular turbu- lence) superposed on a large systematic (background) velocity, caused by granular, mesogranular and supergranular flows.

preprint2014arXiv

Modelling total solar irradiance since 1878 from simulated magnetograms

We present a new model of total solar irradiance (TSI) based on magnetograms simulated with a surface flux transport model (SFTM) and the SATIRE (Spectral And Total Irradiance REconstructions) model. Our model provides daily maps of the distribution of the photospheric field and the TSI starting from 1878. We first calculate the magnetic flux on the solar surface emerging in active and ephemeral regions. The evolution of the magnetic flux in active regions is computed using a surface flux transport model fed with the observed record of sunspot group areas and positions. The magnetic flux in ephemeral regions is treated separately using the concept of overlapping cycles. To model the ephemeral region cycles, we assume that their length and amplitude are related to that of the sunspot cycles. We then use a version of the SATIRE model to compute the TSI. The area coverage and the distribution of different magnetic features as a function of time, which are required by SATIRE, are extracted from the simulated magnetograms and the modelled ephemeral region magnetic flux. Previously computed intensity spectra of the various types of magnetic features are employed. Our model reproduces the PMOD composite of TSI measurements starting from 1978 at daily and rotational timescales more accurately than the previous version of the SATIRE model computing TSI over this period of time. The simulated magnetograms provide a more realistic representation of the evolution of the magnetic field on the photosphere and also allow us to make use of information on the spatial distribution of the magnetic fields before the times when observed magnetograms were available. We find that the secular increase in TSI since 1878 is fairly stable to modifications of the treatment of the ephemeral region magnetic flux.

preprint2014arXiv

Nonlinear Force-Free Field Modeling of the Solar Magnetic Carpet and Comparison with SDO/HMI and Sunrise/IMaX Observations

In the quiet solar photosphere, the mixed polarity fields form a magnetic carpet, which continuously evolves due to dynamical interaction between the convective motions and magnetic field. This interplay is a viable source to heat the solar atmosphere. In this work, we used the line-of-sight (LOS) magnetograms obtained from the Helioseismic and Magnetic Imager (HMI) on the \textit{Solar Dynamics Observatory} (\textit{SDO}), and the Imaging Magnetograph eXperiment (IMaX) instrument on the \textit{Sunrise} balloon-borne observatory, as time dependent lower boundary conditions, to study the evolution of the coronal magnetic field. We use a magneto-frictional relaxation method, including hyperdiffusion, to produce time series of three-dimensional (3D) nonlinear force-free fields from a sequence of photospheric LOS magnetograms. Vertical flows are added up to a height of 0.7 Mm in the modeling to simulate the non-force-freeness at the photosphere-chromosphere layers. Among the derived quantities, we study the spatial and temporal variations of the energy dissipation rate, and energy flux. Our results show that the energy deposited in the solar atmosphere is concentrated within 2 Mm of the photosphere and there is not sufficient energy flux at the base of the corona to cover radiative and conductive losses. Possible reasons and implications are discussed. Better observational constraints of the magnetic field in the chromosphere are crucial to understand the role of the magnetic carpet in coronal heating.

preprint2014arXiv

Point spread function of SDO/HMI and the effects of stray light correction on the apparent properties of solar surface phenomena

We present a PSF for SDO/HMI and discuss the effects of its removal on the apparent properties of solar surface phenomena in HMI data. The PSF was retrieved from observations of Venus in transit by matching it to the convolution of a model of the venusian disc and solar background with a guess PSF. Observations recorded near in time to the transit of Venus were corrected for instrumental scattered light by the deconvolution with the PSF. Granulation contrast in restored HMI data is greatly enhanced relative to the original data and exhibit reasonable agreement with numerical simulations. Image restoration enhanced the apparent intensity and pixel averaged magnetic field strength of photospheric magnetic features significantly. For small-scale magnetic features, restoration enhanced intensity contrast in the continuum and core of the Fe I 617.3 nm line by a factor of 1.3, and the magnetogram signal by a factor of 1.7. For sunspots and pores, the enhancement varied strongly within and between features, being more acute for smaller features. Magnetic features are also rendered smaller, as signal smeared onto the surrounding quiet Sun is recovered. Image restoration increased the apparent amount of magnetic flux above the noise floor by a factor of about 1.2, most of the gain coming from the quiet Sun. Line-of-sight velocity due to granulation and supergranulation is enhanced by a factor of 1.4 to 2.1, depending on position on the solar disc. There are indications that the PSF varies across the FOV, with time and between the two CCDs. However, all these variations were found to be relatively small, such that the derived PSF can be applied to HMI data from both CCDs, over the period examined without introducing significant error. Image restoration affects the measured radiant, magnetic and dynamic properties of solar surface phenomena sufficiently to significantly impact interpretation.

preprint2014arXiv

Scattered Lyman-alpha radiation of comet 2012/S1 (ISON) observed by SUMER/SOHO

During its sungrazing perihelion passage, comet ISON appeared in the field of view of the SUMER spectrometer and allowed unique observations at far-ultraviolet wavelengths with high spatial and temporal resolution. We report results of these observations completed on November 28, 2013, when the comet was only 2.82 R_Sun away from the Sun. Our data show the arrow-shaped dust tail in Ly-$α$ emission trailing behind the predicted position of the nucleus, but offset from the trajectory. We interpret the emission as sunlight that is scattered at micron-sized dust particles. We modeled the dust emission and dynamics to reproduce the appearance of the tail. We were unable to detect any signature of cometary gas or plasma around the expected position of the nucleus and conclude that the outgassing processes must have stopped before the observation started. Moreover, the model we used to reproduce the observed dust tail needs a sharp fall-off of the dust production hours before perihelion transit. We compare the radiances of the disk and the dust tail for an estimate of the dust column density and tail mass.

preprint2014arXiv

Simulated Magnetic Flows in the Solar Photosphere

Recent Sunrise/IMaX observations have revealed the existence of supersonic magnetic flows. Our aim is to determine the origin of such flows by using realistic MHD simulations. We simulate cancellation and emergence of magnetic flux through the solar photosphere. Our first numerical experiment starts with magnetic field of both polarities. To simulate emergence into a region with pre-existing field, we introduce a large-scale horizontally uniform sheet of horizontal field. We follow the subsequent evolution, creating synthetic polarimetric observations, including known instrumental effects of the Sunrise/IMaX and Hinode/SP instruments. We compare the simulated and observed spectropolarimetric signals. Strongly blue- and redshifted Stokes V signals are produced in locations where strong line-of-sight velocities coincide with the strong line-of-sight component of magnetic field. The size and strength of simulated events is smaller than observed and they are mostly associated with downflows, contrary to observations. In a few cases where they appear above a granule, single blue lobed Stokes V are produced due to strong gradients in magnetic field and velocity. No change of magnetic field sign is detected along the line of sight in these instances. More high-speed magnetized flows occur in the case where emergence is simulated than in the case where no horizontal field was added. The simulations indicate that the observed events result from magnetic flux emergence, where reconnection may take place but does not seem to be necessary.

preprint2014arXiv

Solar cycle variation in solar irradiance

The correlation between solar irradiance and the 11-year solar activity cycle is evident in the body of measurements made from space, which extend over the past four decades. Models relating variation in solar irradiance to photospheric magnetism have made significant progress in explaining most of the apparent trends in these observations. There are, however, persistent discrepancies between different measurements and models in terms of the absolute radiometry, secular variation and the spectral dependence of the solar cycle variability. We present an overview of solar irradiance measurements and models, and discuss the key challenges in reconciling the divergence between the two.

preprint2014arXiv

The continuum intensity as a function of magnetic field II. Local magnetic flux and convective flows

To deepen our understanding of the role of small-scale magnetic fields in active regions (ARs) and in the quiet Sun (QS) on the solar irradiance, it is fundamental to investigate the physical processes underlying their continuum brightness. Previous results showed that magnetic elements in the QS reach larger continuum intensities than in ARs at disk center, but left this difference unexplained. We use Hinode/SP disk center data to study the influence of the local amount of magnetic flux on the vigour of the convective flows and the continuum intensity contrasts. The apparent (i.e. averaged over a pixel) longitudinal field strength and line-of-sight (LOS) plasma velocity were retrieved by means of Milne-Eddington inversions (VFISV code). We analyzed a series of boxes taken over AR plages and the QS, to determine how the continuum intensity contrast of magnetic elements, the amplitude of the vertical flows and the box-averaged contrast were affected by the mean longitudinal field strength in the box (which scales with the total unsigned flux in the box). Our results suggest that as the magnetic flux increases locally (e.g. from weak network to strong plage), the heating of the magnetic elements is reduced by the intermediate of a more suppressed convective energy transport within the larger and stronger magnetic patches. This, together with the known presence of larger magnetic features, could explain the previously found lower contrasts of the brightest magnetic elements in ARs compared to the QS. The inhibition of convection also affects the average continuum brightness of a photospheric region, so that at disk center, an area of photosphere in strong network or plage appears darker than a purely quiet one. This is qualitatively consistent with the predictions of 3D MHD simulations.

preprint2014arXiv

The formation and disintegration of magnetic bright points observed by Sunrise/IMaX

The evolution of the physical parameters of magnetic bright points (MBPs) located in the quiet Sun (mainly in the interwork) during their lifetime is studied. First we concentrate on the detailed description of the magnetic field evolution of three MBPs. This reveals that individual features follow different, generally complex, and rather dynamic scenarios of evolution. Next we apply statistical methods on roughly 200 observed MBP evolutionary tracks. MBPs are found to be formed by the strengthening of an equipartition field patch, which initially exhibits a moderate downflow. During the evolution, strong downdrafts with an average velocity of 2.4 km/s set in. These flows, taken together with the concurrent strengthening of the field, suggest that we are witnessing the occurrence of convective collapses in these features, although only 30% of them reach kG field strengths. This fraction might turn out to be larger when the new 4 m class solar telescopes are operational as observations of MBPs with current state of the art instrumentation could still be suffering from resolution limitations. Finally, when the bright point disappears (although the magnetic field often continues to exist) the magnetic field strength has dropped to the equipartition level and is generally somewhat weaker than at the beginning of the MBP's evolution. Noteworthy is that in about 10% of the cases we observe in the vicinity of the downflows small-scale strong (exceeding 2 km/s) intergranular upflows related spatially and temporally to these downflows.

preprint2014arXiv

The variability of Sun-like stars: reproducing observed photometric trends

The Sun and stars with low magnetic activity levels, become photometrically brighter when their activity increases. Magnetically more active stars display the opposite behaviour and get fainter when their activity increases. We reproduce the observed photometric trends in stellar variations with a model that treats stars as hypothetical Suns with coverage by magnetic features different from that of the Sun. The presented model attributes the variability of stellar spectra to the imbalance between the contributions from different components of the solar atmosphere, such as dark starspots and bright faculae. A stellar spectrum is calculated from spectra of the individual components, by weighting them with corresponding disc area coverages. The latter are obtained by extrapolating the solar dependences of spot and facular disc area coverages on chromospheric activity to stars with different levels of mean chromospheric activity. We have found that the contribution by starspots to the variability increases faster with chromospheric activity than the facular contribution. This causes the transition from faculae-dominated variability and direct activity--brightness correlation to spot-dominated variability and inverse activity--brightness correlation with increasing chromospheric activity level. We have shown that the regime of the variability also depends on the angle between the stellar rotation axis and the line-of-sight and on the latitudinal distribution of active regions on the stellar surface. Our model can be used as a tool to extrapolate the observed photometric variability of the Sun to Sun-like stars at different activity levels, which makes possible the direct comparison between solar and stellar irradiance data.

preprint2013arXiv

3-D non-LTE radiative transfer effects in Fe I lines: II. Line formation in 3-D radiation hydrodynamic simulations

We investigate the effects of horizontal radiative transfer (RT) a NLTE on important diagnostic iron lines in a realistic 3-D HD simulation. Using a multi-level atom we compute and compare widely used Fe I line profiles at 3 different levels of approximation (3-D NLTE, 1-D NLTE, LTE). We find that the influence of horizontal RT is of the same order of magnitude as that of NLTE, although spatially more localized. Also, depending on the temperature of the surroundings, horizontal RT is found to weaken or strengthen spectral lines. Line depths and equivalent width may differ by up to 20% against the corresponding LTE value if 3-D RT is applied. Residual intensity contrasts in LTE are found to be larger than those in 3-D NLTE by up to a factor of two. When compared to 1-D NLTE, we find that horizontal RT weakens the contrast by up to 30% almost independently of the angle of line of sight. While the CLV of the 1-D and 3-D NLTE contrasts are of similar form, the LTE contrast CLV shows a different run. The determination of temperatures by 1-D NLTE inversions of spatially resolved observations may produce errors of up to 200 K if one neglects 3-D RT. We find a linear correlation between the intensity difference of 1-D and 3-D NLTE and a simple estimate of the temperature in the horizontal environment of the line formation region. This correlation could be used to coarsely correct for the effects of horizontal RT in inversions done in 1-D NLTE. Horizontal RT is less important if one considers spatially averaged line profiles because local line strengthening and weakening occur with similar frequency in our HD atmosphere. Thus, the iron abundance is underestimated by 0.012 dex if calculated using 1-D NLTE RT. Since effects of horizontal RT are largest for spatially resolved quantities, the use of 3-D RT is particularly important for the interpretation of high spatial resolution observations.

preprint2013arXiv

Fine structures in the atmosphere above a sunspot umbra

We present simultaneous photospheric and chromospheric observations of the trailing sunspot in NOAA 10904, obtained with the Swedish Solar Telescope (SST) La Palma, Canary Islands. Time series of high resolution \ion{Ca}{ii}\,$H$ images show transient jet-like structures in sunspot umbrae are elongated, which we call umbral microjets. These jets are directed roughly parallel to nearby penumbral microjets, suggesting that they are aligned with the background magnetic field. In general, first a bright dot-like structure appears, from which a jet later emerges, although some jets appear without an associated chromospheric dot. Bright photospheric umbral dots are associated with umbral microjets arising in the outer umbra. Nevertheless, a one-to-one correspondence between jet-like events and underlying umbral dots is not seen. They are typically less than 1\arcsec ~long and less than 0\farcs3 wide. The typical lifetime of umbral microjets is around one minute. The brightness of these structures increases from the center of the umbra towards the umbra-penumbra boundary along with the brightness of the local background.

preprint2013arXiv

First evidence of interaction between longitudinal and transverse waves in solar magnetic elements

Small-scale magnetic fields are thought to play an important role in the heating of the outer solar atmosphere. By taking advantage of the unprecedented high-spatial and temporal cadence of IMaX, the filter vector polarimeter on board the Sunrise balloon-borne observatory, we study the transversal and longitudinal velocity oscillations in small magnetic elements. The results of this analysis are then compared to MHD simulations, showing excellent agreement. We found buffeting-induced transverse oscillations with velocity amplitudes of the order of 1-2 km/s, to be common along with longitudinal oscillations with amplitudes 0.4 km/s. Moreover, we also found an interaction between transverse oscillations and longitudinal velocity oscillations, showing a +-90 degrees phase lag at the frequency at which they exhibit the maximum coherence in the power spectrum. Our results are consistent with the theoretical picture in which MHD longitudinal waves are excited inside small magnetic elements as a response of the flux tube to the forcing action of the granular flows.

preprint2013arXiv

First high-resolution images of the Sun in the 2796 Å Mg II k line

We present the first high-resolution solar images in the Mg II k 2796 Å line. The images, taken through a 4.8 Å broad interference filter, were obtained during the second science flight of SUNRISE in June 2013 by the SuFI instrument. The Mg II k images display structures that look qualitatively very similar to images taken in the core of Ca II H. The Mg II images exhibit reversed granulation (or shock waves) in the internetwork regions of the quiet Sun, at intensity contrasts that are similar to those found in Ca II H. Very prominent in Mg II are bright points, both in the quiet Sun and in plage regions, particularly near disk center. These are much brighter than at other wavelengths sampled at similar resolution. Furthermore, Mg II k images also show fibril structures associated with plage regions. Again, the fibrils are similar to those seen in Ca II H images, but tend to be more pronounced, particularly in weak plage.

preprint2013arXiv

Intensity contrast of solar network and faculae

This study aims at setting observational constraints on the continuum and line core intensity contrast of network and faculae, specifically, their relationship with magnetic field and disc position. Full-disc magnetograms and intensity images by the HMI instrument onboard SDO were employed. Bright magnetic features, representing network and faculae, were identified and the relationship between their intensity contrast at continuum and line core with magnetogram signal and heliocentric angle examined. Care was taken to minimize the inclusion of the magnetic canopy and straylight from sunspots and pores as network and faculae. In line with earlier studies, network features, on a per unit magnetic flux basis, appeared brighter than facular features. Intensity contrasts in the continuum and line core differ considerably, most notably, they exhibit opposite centre-to-limb variations. We found this difference in behaviour to likely be due to the different mechanisms of the formation of the two spectral components. From a simple model based on bivariate polynomial fits to the measured contrasts we confirmed spectral line changes to be a significant driver of facular contribution to variation in solar irradiance. The discrepancy between the continuum contrast reported here and in the literature was shown to arise mainly from differences in spatial resolution and treatment of magnetic signals adjacent to sunspots and pores. HMI is a source of accurate contrasts and low-noise magnetograms covering the full solar disc. For irradiance studies it is important to consider not just the contribution from the continuum but also from the spectral lines. In order not to underestimate long-term variations in solar irradiance, irradiance models should take the greater contrast per unit magnetic flux associated with magnetic features with low magnetic flux into account.

preprint2013arXiv

Magnetic structure of an activated filament in a flaring active region

While the magnetic field in quiescent prominences has been widely investigated, less is known about the field in activated prominences. We introduce observational results on the magnetic field structure of an activated filament in a flaring active region. We study, in particular, its magnetic structure and line-of-sight flows during its early activated phase, shortly before it displays signs of rotation. We invert the Stokes profiles of the chromospheric He I 10830 A triplet and the photospheric Si I 10827 A line observed in this filament by the VTT on Tenerife. Using these inversion results we present and interpret the first maps of velocity and magnetic field obtained in an activated filament, both in the photosphere and the chromosphere. Up to 5 different magnetic components are found in the chromospheric layers of the filament, while outside the filament a single component is sufficient to reproduce the observations. Magnetic components displaying an upflow are preferentially located towards the centre of the filament, while the downflows are concentrated along its periphery. Also, the upflowing gas is associated with an opposite-polarity magnetic configuration with respect to the photosphere, while the downflowing gas is associated with a same-polarity configuration. The activated filament has a rather complex structure. Nonetheless, it is compatible with a flux rope, although with a distorted one, in the normal configuration. The observations are best explained by a rising flux rope in which a part of the filament material is still stably stored (upflowing material, rising with the field), while a part is no longer stably stored and flows down along the field lines.

preprint2013arXiv

Quiet Sun magnetic fields observed by Hinode: Support for a local dynamo

The Hinode mission has revealed copious amounts of horizontal flux covering the quiet Sun. Local dynamo action has been proposed to explain the presence of this flux. We sought to test whether the quiet Sun flux detected by Hinode is due to a local or the global dynamo by studying long-term variations in the polarisation signals detectable at the disc centre of the quiet Sun between November 2006 and May 2012, with particular emphasis on weak signals in the internetwork. The investigation focusses on line-integrated circular polarisation V_tot and linear polarisation LP_tot profiles obtained from the Fe I 6302.5 Åabsorption line in Hinode SOT/SP. Both circular and linear polarisation signals show no overall variation in the fraction of selected pixels from 2006 until 2012. There is also no variation in the magnetic flux in this interval of time. The probability density functions (PDF) of the line-of-sight magnetic flux can be fitted with a power law from 1.17 x 10^17 Mx to 8.53 x 10^18 Mx with index α=-1.82 \pm 0.02 in 2007. The variation of α's across all years does not exceed a significance of 1σ. Linearly polarised features are also fitted with a power law, with index α=-2.60 \pm 0.06 in 2007. Indices derived from linear polarisation PDFs of other years also show no significant variation. Our results show that the ubiquitous horizontal polarisation on the edges of bright granules seen by Hinode are invariant during the minimum of cycle 23. This supports the notion that the weak circular and linear polarisation is primarily caused by an independent local dynamo.

preprint2013arXiv

Recent variability of the solar spectral irradiance and its impact on climate modelling

The lack of long and reliable time series of solar spectral irradiance (SSI) measurements makes an accurate quantification of solar contributions to recent climate change difficult. Whereas earlier SSI observations and models provided a qualitatively consistent picture of the SSI variability, recent measurements by the SORCE satellite suggest a significantly stronger variability in the ultraviolet (UV) spectral range and changes in the visible and near-infrared (NIR) bands in anti-phase with the solar cycle. A number of recent chemistry-climate model (CCM) simulations have shown that this might have significant implications on the Earth's atmosphere. Motivated by these results, we summarize here our current knowledge of SSI variability and its impact on Earth's climate. We present a detailed overview of existing SSI measurements and provide thorough comparison of models available to date. SSI changes influence the Earth's atmosphere, both directly, through changes in shortwave (SW) heating and therefore, temperature and ozone distributions in the stratosphere, and indirectly, through dynamical feedbacks. We investigate these direct and indirect effects using several state-of-the art CCM simulations forced with measured and modeled SSI changes. A unique asset of this study is the use of a common comprehensive approach for an issue that is usually addressed separately by different communities. Omissis. Finally, we discuss the reliability of the available data and we propose additional coordinated work, first to build composite SSI datasets out of scattered observations and to refine current SSI models, and second, to run coordinated CCM experiments.

preprint2013arXiv

Solar Irradiance Variability and Climate

The brightness of the Sun varies on all time scales on which it has been observed, and there is increasing evidence that it has an influence on climate. The amplitudes of such variations depend on the wavelength and possibly on the time scale. Although many aspects of this variability are well established, the exact magnitude of secular variations (going beyond a solar cycle) and the spectral dependence of variations are under discussion. The main drivers of solar variability are thought to be magnetic features at the solar surface. The climate reponse can be, on a global scale, largely accounted for by simple energetic considerations, but understanding the regional climate effects is more difficult. Promising mechanisms for such a driving have been identified, including through the influence of UV irradiance on the stratosphere and dynamical coupling to the surface. Here we provide an overview of the current state of our knowledge, as well as of the main open questions.

preprint2013arXiv

Temporal relation between quiet-Sun transverse fields and the strong flows detected by IMaX/SUNRISE

Localized strongly Doppler-shifted Stokes V signals were detected by IMaX/SUNRISE. These signals are related to newly emerged magnetic loops that are observed as linear polarization features. We aim to set constraints on the physical nature and causes of these highly Doppler-shifted signals. In particular, the temporal relation between the appearance of transverse fields and the strong Doppler shifts is analyzed in some detail. We calculated the time difference between the appearance of the strong flows and the linear polarization. We also obtained the distances from the center of various features to the nearest neutral lines and whether they overlap or not. These distances were compared with those obtained from randomly distributed points on observed magnetograms. Various cases of strong flows are described in some detail. The linear polarization signals precede the appearance of the strong flows by on average 84+-11 seconds. The strongly Doppler-shifted signals are closer (0.19") to magnetic neutral lines than randomly distributed points (0.5"). Eighty percent of the strongly Doppler-shifted signals are close to a neutral line that is located between the emerging field and pre-existing fields. That the remaining 20% do not show a close-by pre-existing field could be explained by a lack of sensitivity or an unfavorable geometry of the pre-existing field, for instance, a canopy-like structure. Transverse fields occurred before the observation of the strong Doppler shifts. The process is most naturally explained as the emergence of a granular-scale loop that first gives rise to the linear polarization signals, interacts with pre-existing fields (generating new neutral line configurations), and produces the observed strong flows. This explanation is indicative of frequent small-scale reconnection events in the quiet Sun.

preprint2013arXiv

The AD775 cosmic event revisited: the Sun is to blame

Miyake et al. (henceforth M12) recently reported, based on 14C data, an extreme cosmic event ca. AD775. Using a simple model, M12 claimed that the event was too strong to be caused by a solar flare within the standard theory. This implied a new paradigm of either an impossibly strong solar flare or a very strong cosmic ray event of unknown origin occurred ca. AD775. We show that the strength of the event was significantly overestimated by M12. Several subsequent works have attempted to find a possible exotic source for such an event, but they are all based on incorrect estimates by M12. We revisit this event with analysis of new datasets and consistent theoretical modelling. We verified the experimental result for the AD775 event using independent datasets including 10Be series and newly measured 14C annual data. We surveyed available historical chronicles for astronomical observations for the AD770s to identify potential sightings of aurorae or supernovae. We interpreted the 14C measurements using an appropriate carbon cycle model. We show that: (1) The reality of the AD775 event is confirmed by new measurements of 14C; (2) by using an inappropriate carbon cycle model, M12 strongly overestimated the event's strength; (3) The revised magnitude of the event is consistent with different independent datasets (14C, 10Be, 36Cl) and can be associated with a strong, but not inexplicably strong, SEP event (or a sequence of events), and provides the first evidence for an event of this magnitude (the fluence >30 MeV was about 4.5*10^{10} /cm2) in multiple datasets; (4) This is in agreement with increased auroral activity identified in historical chronicles. This point to the likely solar origin of the event, which is the greatest solar event on a multi-millennial time scale, placing a strong observational constraint on the theory of explosive energy releases on the Sun and cool stars.

preprint2013arXiv

Vertical flows and mass flux balance of sunspot umbral dots

A new Stokes inversion technique that greatly reduces the effect of the spatial point spread function of the telescope is used to constrain the physical properties of umbral dots (UDs). The depth-dependent inversion of the Stokes parameters from a sunspot umbra recorded with Hinode SOT/SP revealed significant temperature enhancements and magnetic field weakenings in the core of the UDs in deep photospheric layers. Additionally, we found upflows of around 960 m/s in peripheral UDs (i.e., UDs close to the penumbra) and $\approx$ 600 m/s in central UDs. For the first time, we also detected systematic downflows for distances larger than 200 km from the UD center that balance the upflowing mass flux. In the upper photosphere, we found almost no difference between the UDs and their diffuse umbral background.

preprint2012arXiv

A retrospective of the GREGOR solar telescope in scientific literature

In this review, we look back upon the literature, which had the GREGOR solar telescope project as its subject including science cases, telescope subsystems, and post-focus instruments. The articles date back to the year 2000, when the initial concepts for a new solar telescope on Tenerife were first presented at scientific meetings. This comprehensive bibliography contains literature until the year 2012, i.e., the final stages of commissioning and science verification. Taking stock of the various publications in peer-reviewed journals and conference proceedings also provides the "historical" context for the reference articles in this special issue of Astronomische Nachrichten/Astronomical Notes.

preprint2012arXiv

Break up of returning plasma after the 7 June 2011 filament eruption by Rayleigh-Taylor instabilities

A prominence eruption on 7 June 2011 produced spectacular curtains of plasma falling through the lower corona. At the solar surface they created an incredible display of extreme ultraviolet brightenings. The aim is to identify and analyze some of the local instabilities which produce structure in the falling plasma. The structures were investigated using SDO/AIA 171A and 193A images in which the falling plasma appeared dark against the bright coronal emission. Several instances of the Rayleigh-Taylor instability were investigated. In two cases the Alfven velocity associated with the dense plasma could be estimated from the separation of the Rayleigh-Taylor fingers. A second type of feature, which has the appearance of self-similar branching horns, is also discussed.

preprint2012arXiv

CAST search for sub-eV mass solar axions with 3He buffer gas

The CERN Axion Solar Telescope (CAST) has extended its search for solar axions by using 3He as a buffer gas. At T=1.8 K this allows for larger pressure settings and hence sensitivity to higher axion masses than our previous measurements with 4He. With about 1 h of data taking at each of 252 different pressure settings we have scanned the axion mass range 0.39 eV < m_a < 0.64 eV. From the absence of excess X-rays when the magnet was pointing to the Sun we set a typical upper limit on the axion-photon coupling of g_ag < 2.3 x 10^{-10} GeV^{-1} at 95% CL, the exact value depending on the pressure setting. KSVZ axions are excluded at the upper end of our mass range, the first time ever for any solar axion search. In future we will extend our search to m_a < 1.15 eV, comfortably overlapping with cosmological hot dark matter bounds.

preprint2012arXiv

Long-term magnetic field monitoring of the Sun-like star ξBootis A

Aims. We aim to investigate the long-term temporal evolution of the magnetic field of the solar-type star ξBootis A, both from direct magnetic field measurements and from the simultaneous estimate of indirect activity indicators. Methods. We obtained seven epochs of high-resolution, circularly-polarized spectra from the NARVAL spectropolarimeter between 2007 and 2011, for a total of 76 spectra. Using approximately 6,100 photospheric spectral lines covering the visible domain, we employed a cross-correlation procedure to compute a mean polarized line profile from each spectrum. The large-scale photospheric magnetic field of the star was then modelled by means of Zeeman-Doppler Imaging, allowing us to follow the year-to-year evolution of the reconstructed magnetic topology. Simultaneously, we monitored the width of several magnetically sensitive spectral lines, the radial velocity, the line asymmetry of intensity line profiles, and the chromospheric emission in the cores of the Ca II H and Halpha lines. Results. During the highest observed activity states, in 2007 and 2011, the large-scale field of ξBoo A is almost completely axisymmetric and is dominated by its toroidal component. The magnetic topologies reconstructed for these activity maxima are very similar, suggesting a form of short cyclicity in the large-scale field distribution. Correlated temporal evolution, due to both rotational modulation and seasonal variability, is observed between the Ca II emission, the Halpha emission and the width of magnetically sensitive lines. When measurable, the differential rotation reveals a strong latitudinal shear in excess of 0.2 rad/d.

preprint2012arXiv

Resolving the internal magnetic structure of the solar network

We analyze the spectral asymmetry of Stokes V (circularly polarized) profiles of an individual network patch in the quiet Sun observed by Sunrise/IMaX. At a spatial resolution of 0.15"-0.18", the network elements contain substructure which is revealed by the spatial distribution of Stokes V asymmetries. The area asymmetry between the red and blue lobes of Stokes V increases from nearly zero at the core of the structure to values close to unity at its edges (one-lobed profiles). Such a distribution of the area asymmetry is consistent with magnetic fields expanding with height, i.e., an expanding magnetic canopy (which is required to fulfill pressure balance and flux conservation in the solar atmosphere). Inversion of the Stokes I and V profiles of the patch confirms this picture, revealing a decreasing field strength and increasing height of the canopy base from the core to the periphery of the network patch. However, the non-roundish shape of the structure and the presence of negative area and amplitude asymmetries reveal that the scenario is more complex than a canonical flux tube expanding with height surrounded by downflows.

preprint2012arXiv

Results and perspectives of the solar axion search with the CAST experiment

The status of the solar axion search with the CERN Axion Solar Telescope (CAST) will be presented. Recent results obtained by the use of $^3$He as a buffer gas has allowed us to extend our sensitivity to higher axion masses than our previous measurements with $^4$He. With about 1 h of data taking at each of 252 different pressure settings we have scanned the axion mass range 0.39 eV$ \le m_{a} \le $ 0.64 eV. From the absence of an excess of x rays when the magnet was pointing to the Sun we set a typical upper limit on the axion-photon coupling of g$_{aγ} \le 2.3\times 10^{-10}$ GeV$^{-1}$ at 95% C.L., the exact value depending on the pressure setting. CAST published results represent the best experimental limit on the photon couplings to axions and other similar exotic particles dubbed WISPs (Weakly Interacting Slim Particles) in the considered mass range and for the first time the limit enters the region favored by QCD axion models. Preliminary sensitivities for axion masses up to 1.16 eV will also be shown reaching mean upper limits on the axion-photon coupling of g$_{aγ} \le 3.5\times 10^{-10}$ GeV$^{-1}$ at 95% C.L. Expected sensibilities for the extension of the CAST program up to 2014 will be presented. Moreover long term options for a new helioscope experiment will be evoked.

preprint2012arXiv

Spectroscopic Observations of Propagating Disturbances in a Polar Coronal Hole: Evidence of Slow Magneto-acoustic Waves

We focus on detecting and studying quasi-periodic propagating features that have been interpreted both in terms of slow magneto-acoustic waves and of high speed upflows. We analyze long duration spectroscopic observations of the on-disk part of the south polar coronal hole taken on 1997 February 25 by the SUMER spectrometer aboard SOHO. We calibrated the velocity with respect to the off-limb region and obtain time--distance maps in intensity, Doppler velocity and line width. We also perform a cross correlation analysis on different time series curves at different latitudes. We study average spectral line profiles at the roots of propagating disturbances and along the propagating ridges, and perform a red-blue asymmetry analysis. We find the clear presence of propagating disturbances in intensity and Doppler velocity with a projected propagation speed of about $60\pm 4.8$ km s$^{-1}$ and a periodicity of $\approx$14.5 min. To our knowledge, this is the first simultaneous detection of propagating disturbances in intensity as well as in Doppler velocity in a coronal hole. During the propagation, an intensity enhancement is associated with a blue-shifted Doppler velocity. These disturbances are clearly seen in intensity also at higher latitudes (i.e. closer to the limb), while disturbances in Doppler velocity becomes faint there. The spectral line profiles averaged along the propagating ridges are found to be symmetric, to be well fitted by a single Gaussian, and have no noticeable red-blue asymmetry. Based on our analysis, we interpret these disturbances in terms of propagating slow magneto-acoustic waves.

preprint2012arXiv

Structure and Dynamics of Isolated Internetwork Ca II H Bright Points Observed by Sunrise

We aim to improve our picture of the low chromosphere in the quiet-Sun internetwork by investigating the intensity, horizontal velocity, size and lifetime variations of small bright points (BPs; diameter smaller than 0.3 arcsec) observed in the Ca II H 3968 Å passband along with their magnetic field parameters, derived from photospheric magnetograms. Several high-quality time series of disc-centre, quiet-Sun observations from the Sunrise balloon-borne solar telescope, with spatial resolution of around 100 km on the solar surface, have been analysed to study the dynamics of BPs observed in the Ca II H passband and their dependence on the photospheric vector magnetogram signal. Parameters such as horizontal velocity, diameter, intensity and lifetime histograms of the isolated internetwork and magnetic Ca II H BPs were determined. Mean values were found to be 2.2 km/s, 0.2 arcsec (150 km), 1.48 average Ca II H quiet-Sun and 673 sec, respectively. Interestingly, the brightness and the horizontal velocity of BPs are anti-correlated. Large excursions (pulses) in horizontal velocity, up to 15 km/s, are present in the trajectories of most BPs. These could excite kink waves travelling into the chromosphere and possibly the corona, which we estimate to carry an energy flux of 310 W/m^2, sufficient to heat the upper layers, although only marginally. The stable observing conditions of Sunrise and our technique for identifying and tracking BPs have allowed us to determine reliable parameters of these features in the internetwork. Thus we find, e.g., that they are considerably longer lived than previously thought. The large velocities are also reliable, and may excite kink waves. Although these wave are (marginally) energetic enough to heat the quiet corona, we expect a large additional contribution from larger magnetic elements populating the network and partly also the internetwork.

preprint2012arXiv

Supersonic Magnetic Flows in the Quiet Sun

In this contribution we describe some recent observations of high-speed magnetized flows in the quiet Sun granulation. These observations were carried out with the Imaging Magnetograph eXperiment (IMaX) onboard the stratospheric balloon {\sc Sunrise}, and possess an unprecedented spatial resolution and temporal cadence. These flows were identified as highly shifted circular polarization (Stokes $V$) signals. We estimate the LOS velocity responsible for these shifts to be larger than 6 km s$^{-1}$, and therefore we refer to them as {\it supersonic magnetic flows}. The average lifetime of the detected events is 81.3 s and they occupy an average area of about 23\,000 km$^2$. Most of the events occur within granular cells and correspond therefore to upflows. However some others occur in intergranular lanes or bear no clear relation to the convective velocity pattern. We analyze a number of representative examples and discuss them in terms of magnetic loops, reconnection events, and convective collapse.

preprint2012arXiv

The Height of Chromospheric Loops in an Emerging Flux Region

Context. The chromospheric layer observable with the He I 10830 Å triplet is strongly warped. The analysis of the magnetic morphology of this layer therefore requires a reliable technique to determine the height at which the He I absorption takes place. Aims. The He I absorption signature connecting two pores of opposite polarity in an emerging flux region is investigated. This signature is suggestive of a loop system connecting the two pores. We aim to show that limits can be set on the height of this chromospheric loop system. Methods. The increasing anisotropy in the illumination of a thin, magnetic structure intensifies the linear polarization signal observed in the He I triplet with height. This signal is altered by the Hanle effect. We apply an inversion technique incorporating the joint action of the Hanle and Zeeman effects, with the absorption layer height being one of the free parameters. Results. The observed linear polarization signal can be explained only if the loop apex is higher than \approx5 Mm. Best agreement with the observations is achieved for a height of 6.3 Mm. Conclusions. The strength of the linear polarization signal in the loop apex is inconsistent with the assumption of a He I absorption layer at a constant height level. The determined height supports the earlier conclusion that dark He 10830 Å filaments in emerging flux regions trace emerging loops.

preprint2012arXiv

Three-dimensional non-LTE radiative transfer effects in Fe I lines I. Flux sheet and flux tube geometries

In network and active region plages, the magnetic field is concentrated into structures often described as flux tubes (FTs) and sheets (FSs). 3-D radiative transfer (RT) is important for energy transport in these concentrations. It is also expected to be important for diagnostic purposes but has rarely been applied for that purpose. Using true 3-D, non-LTE (NLTE) RT in FT/FS models, we compute Fe line profiles commonly used to diagnose the Sun's magnetic field by comparing the results with those obtained from LTE/1-D (1.5-D) NLTE calculations. Employing a multilevel iron atom, we study the influence of basic parameters such as Wilson depression, wall thickness, radius/width, thermal stratification or magnetic field strength on all Stokes $I$ parameters in the thin-tube approximation. The use of different levels of approximations of RT may lead to considerable differences in profile shapes, intensity contrasts, equivalent widths, and the determination of magnetic field strengths. In particular, LTE, which often provides a good approach in planar 1-D atmospheres, is a poor approximation in our flux sheet model for some of the most important diagnostic Fe I lines (524.7nm, 525.0nm, 630.1nm, and 630.2nm). The observed effects depend on parameters such as the height of line formation, field strength, and internal temperature stratification. Differences between the profile shapes may lead to errors in the determination of magnetic fields on the order of 10 to 20%, while errors in the determined temperature can reach 300-400K. The empirical FT models NET and PLA turn out to minimize the effects of 3D RT, so that results obtained with these models by applying LTE may also remain valid for 3-D NLTE calculations. Finally, horizontal RT is found to only insignificantly smear out structures such as the optically thick walls of FTs and FSs, allowing features as narrow as 10km to remain visible.

preprint2012arXiv

Waves as the source of apparent twisting motions in sunspot penumbrae

The motion of dark striations across bright filaments in a sunspot penumbra has become an important new diagnostic of convective gas flows in penumbral filaments. The nature of these striations has, however, remained unclear. Here we present an analysis of small scale motions in penumbral filaments in both simulations and observations. The simulations, when viewed from above, show fine structure with dark lanes running outwards from the dark core of the penumbral filaments. The dark lanes either occur preferentially on one side or alternate between both sides of the filament. We identify this fine structure with transverse (kink) oscillations of the filament, corresponding to a sideways swaying of the filament. These oscillations have periods in the range of 5-7 min and propagate outward and downward along the filament. Similar features are found in observed G-band intensity time series of penumbral filaments in a sunspot located near disk center obtained by the Broadband Filter Imager (BFI) on board {\it Hinode}. We also find that some filaments show dark striations moving to both sides of the filaments. Based on the agreement between simulations and observations we conclude that the motions of these striations are caused by transverse oscillations of the underlying bright filaments.

preprint2011arXiv

Continuous upflows and sporadic downflows observed in active regions

We present a study of the temporal evolution of coronal loops in active regions and its implications for the dynamics in coronal loops. We analyzed images of the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) at multiple temperatures to detect apparent motions in the coronal loops. Quasi-periodic brightness fluctuations propagate upwards from the loop footpoint in hot emission at 1MK, while sporadic downflows are seen in cool emission below 1MK. The upward motion in hot emission increases just after the cool downflows. The apparent propagating pattern suggests a hot upflow from the loop footpoints, and is considered to supply hot plasma into the coronal loop, but a wavelike phenomenon cannot be ruled out. Coronal condensation occasionally happens in the coronal loop, and the cool material flows down to the footpoint. Emission from cool plasma could have a significant contribution to hot AIA channels in the event of coronal condensation.

preprint2011arXiv

Convective Nature of Sunspot Penumbral Filaments: Discovery of Downflows in the Deep Photosphere

We study the velocity structure of penumbral filaments in the deep photosphere to obtain direct evidence for the convective nature of sunspot penumbrae. A sunspot was observed at high spatial resolution with the 1-m Swedish Solar Telescope in the deep photospheric C I 5380 Å absorption line. The Multi-Object Multi-Frame Blind Deconvolution (MOMFBD) method is used for image restoration and straylight is filtered out. We report here the discovery of clear redshifts in the C I 5380 Å line at multiple locations in sunspot penumbral filaments. For example, bright head of filaments show larger concentrated blueshift and are surrounded by darker, redshifted regions, suggestive of overturning convection. Elongated downflow lanes are also located beside bright penumbral fibrils. Our results provide the strongest evidence yet for the presence of overturning convection in penumbral filaments and highlight the need to observe the deepest layers of the penumbra in order to uncover the energy transport processes taking place there.

preprint2011arXiv

EUV jets, type III radio bursts and sunspot waves investigated using SDO/AIA observations

Images from the Solar Dynamics Observatory (SDO) at 211A are used to identify the solar source of the type III radio bursts seen in WIND/WAVES dynamic spectra. We analyse a 2.5 hour period during which six strong bursts are seen. The radio bursts correlate very well with the EUV jets coming from the western side of a sunspot in AR11092. The EUV jet emission also correlates well with brightening at what looks like their footpoint at the edge of the umbra. For 10-15 min after strong EUV jets are ejected, the footpoint brightens at roughly 3 min intervals. In both the EUV images and the extracted light curves, it looks as though the brightening is related to the 3-min sunspot oscillations, although the correlation coefficient is rather low. The only open field near the jets is rooted in the sunspot. We conclude that active region EUV/X-ray jets and interplanetary electron streams originate on the edge of the sunspot umbra. They form along a current sheet between the sunspot open field and closed field connecting to underlying satellite flux. Sunspot running penumbral waves cause roughly 3-min jet footpoint brightening. The relationship between the waves and jets is less clear.

preprint2011arXiv

Kink waves in an active region dynamic fibril

We present high spatial and temporal resolution Ca II 8542 observations of a kink wave in an on-disk chromospheric active region fibril. The properties of the wave are similar to those observed in off-limb spicules. From the observed phase and period of the wave we determine a lower limit for the field strength in the chromospheric active region fibril located at the edge of a sunspot to be a few hundred Gauss. We find indications that the event was triggered by a small-scale reconnection event higher up in the atmosphere.

preprint2011arXiv

Magnetic field emergence in mesogranular-sized exploding granules observed with SUNRISE/IMaX data

We report on magnetic field emergences covering significant areas of exploding granules. The balloon-borne mission SUNRISE provided high spatial and temporal resolution images of the solar photosphere. Continuum images, longitudinal and transverse magnetic field maps and Dopplergrams obtained by IMaX onboard SUNRISE are analyzed by Local Correlation Traking (LCT), divergence calculation and time slices, Stokes inversions and numerical simulations are also employed. We characterize two mesogranular-scale exploding granules where $\sim$ 10$^{18}$ Mx of magnetic flux emerges. The emergence of weak unipolar longitudinal fields ($\sim$100 G) start with a single visible magnetic polarity, occupying their respective granules' top and following the granular splitting. After a while, mixed polarities start appearing, concentrated in downflow lanes. The events last around 20 min. LCT analyses confirm mesogranular scale expansion, displaying a similar pattern for all the physical properties, and divergence centers match between all of them. We found a similar behaviour with the emergence events in a numerical MHD simulation. Granule expansion velocities are around 1 \kms while magnetic patches expand at 0.65 \kms. One of the analyzed events evidences the emergence of a loop-like structure. Advection of the emerging magnetic flux features is dominated by convective motion resulting from the exploding granule due to the magnetic field frozen in the granular plasma. Intensification of the magnetic field occurs in the intergranular lanes, probably because of being directed by the downflowing plasma.

preprint2011arXiv

Milne-Eddington inversions of the He I 10830 Å Stokes profiles: Influence of the Paschen-Back effect

The Paschen-Back effect influences the Zeeman sublevels of the He I multiplet at 10830 Å, leading to changes in strength and in position of the Zeeman components of these lines. We illustrate the relevance of this effect using synthetic Stokes profiles of the He I 10830 Å multiplet lines and investigate its influence on the inversion of polarimetric data. We invert data obtained with the Tenerife Infrared Polarimeter (TIP) at the German Vacuum Tower Telescope (VTT). We compare the results of inversions based on synthetic profiles calculated with and without the Paschen-Back effect being included. We find that when taking into account the incomplete Paschen-Back effect, on average 16% higher field strength values are obtained. We also show that this effect is not the main cause for the area asymmetry exhibited by many He I 10830 Stokes V-profiles. This points to the importance of velocity and magnetic field gradients over the formation height range of these lines.

preprint2011arXiv

Multicomponent He I 10830 Å profiles in an active filament

We present new spectropolarimetric observations of the chromospheric He I 10830 Å multiplet observed in a filament during its phase of activity. The data were recorded with the new Tenerife Infrared Polarimeter (TIP-II) at the German Vacuum Tower Telescope (VTT) on 2005 May 18. We inverted the He Stokes profiles using multiple atmospheric components. The observed He Stokes profiles display a remarkably wide variety of shapes. Most of the profiles show very broad Stokes I absorptions and complex and spatially variable Stokes V signatures. The inversion of the profiles shows evidence of different atmospheric blue- and redshifted components of the He I lines within the resolution element (1 arcsec), with supersonic velocities of up to 100 km/s. Up to five different atmospheric components are found in the same profile. We show that even these complex profiles can be reliably inverted.

preprint2011arXiv

Solar total irradiance in cycle 23

The apparently unusual behaviour of the TSI during the most recent minimum of solar activity has been interpreted as evidence against solar surface magnetism as the main driver of the secular change in the TSI. We test claims that the evolution of the solar surface magnetic field does not reproduce the observed TSI in cycle 23. We use sensitive, 60-minute averaged MDI magnetograms and quasi-simultaneous continuum images as an input to our SATIRE-S model and calculate the TSI variation over cycle 23, sampled roughly twice-monthly. The computed TSI is then compared to the PMOD composite of TSI measurements and to the data from two individual instruments, SORCE/TIM and UARS/ACRIM II, that monitored the TSI during the declining phase of cycle 23 and over the previous minimum in 1996, respectively. Excellent agreement is found between the trends shown by the model and almost all sets of measurements. The only exception is the early, i.e. 1996 to 1998, PMOD data. Whereas the agreement between the model and the PMOD composite over the period 1999-2009 is almost perfect, the modelled TSI shows a steeper increase between 1996 and 1999 than implied by the PMOD composite. On the other hand, the steeper trend in the model agrees remarkably well with the ACRIM II data. A closer look at the VIRGO data, that make the basis of the PMOD composite after 1996, reveals that only one of the two VIRGO instruments, the PMO6V, shows the shallower trend present in the composite, whereas the DIARAD measurements indicate a steeper trend. We conclude that (1) the sensitivity changes of the PMO6V radiometers within VIRGO during the first two years have very likely not been correctly evaluated, and that (2) the TSI variations over cycle 23 and the change in the TSI levels between the minima in 1996 and 2008 are consistent with the solar surface magnetism mechanism.

preprint2011arXiv

The continuum intensity as a function of magnetic field I. Active region and quiet Sun magnetic elements

Small-scale magnetic fields are major contributors to the solar irradiance variations. Hence, the continuum intensity contrast of magnetic elements in the quiet Sun (QS) network and in active region (AR) plage is an essential quantity that needs to be measured reliably. By using Hinode/SP disk center data at a constant, high spatial resolution, we aim at updating results of earlier ground-based studies of contrast vs. magnetogram signal, and to look for systematic differences between AR plages and QS network. The field strength, filling factor and inclination of the field was retrieved by Milne-Eddington inversion (VFISV). As in earlier studies, we performed a pixel-by-pixel study of 630.2 nm continuum contrast vs. apparent (i.e. averaged over a pixel) longitudinal magnetic field over large fields of view in ARs and in the QS. The contrast of magnetic elements reaches larger values in the QS (on average 3.7%) than in ARs (1.3%). This could not be attributed to any systematic difference in the chosen contrast references. At Hinode's spatial resolution, the relationship between contrast and apparent longitudinal field strength exhibits a peak at around 700 G in both the QS and ARs, whereas earlier lower resolution studies only found a peak in the QS and a monotonous decrease in ARs. We attribute this discrepancy both to our careful removal of the pores and their close surroundings affected by the telescope diffraction, as well as to the enhanced spatial resolution and very low scattered light of the Hinode Solar Optical Telescope. According to our inversions, the magnetic elements producing the peak of the contrast curves have similar properties (field strength, inclination, filling factor) in ARs and in the QS, so that the larger brightness of magnetic elements in the QS remains unexplained.

preprint2011arXiv

The frontier between Small-scale bipoles and Ephemeral Regions in the solar photosphere: Emergence and Decay of an Intermediate-scale bipole observed with IMaX/SUNRISE

We report on the photospheric evolution of an intermediate-scale (~4 Mm footpoint separation) magnetic bipole, from emergence to decay, observed in the quiet Sun at high spatial 0".3 and temporal (33 s) resolution. The observations were acquired by the IMaX imaging magnetograph during the first science flight of the Sunrise balloon-borne solar observatory. The bipole flux content is 6 x 10^17 Mx, representing a structure bridging the gap between granular scale bipoles and the smaller ephemeral regions. Footpoints separate at a speed of 3.5 km s-1 and reach a maximum distance of 4.5 Mm before the field dissolves. The evolution of the bipole is revealed to be very dynamic: we found a proper motion of the bipole axis and detected a change of the azimuth angle of 90° in 300 seconds. The overall morphology and behaviour are in agreement with previous analyses of bipolar structures emerging at granular scale, but we also found several similarities with emerging flux structures at larger scale. The flux growth rate is 2.6 x 15 Mx s-1, while the mean decay rate is one order of magnitude smaller. We describe in some detail the decay phase of the bipole footpoints which includes break up into smaller structures, interaction with pre-existing fields leading to cancellation but appears to be dominated by an as-yet unidentified diffusive process that removes most of the flux with an exponential flux decay curve. The diffusion constant (8 x 10^2 km^2 s-1) associated with this decay is similar to the values used to describe the large scale diffusion in flux transport models.

preprint2011arXiv

Transport of magnetic flux from the canopy to the internetwork

Recent observations have revealed that 8% of linear polarization patches in the internetwork quiet Sun are fully embedded in downflows. These are not easily explained with the typical scenarios for the source of internetwork fields which rely on flux emergence from below. We explore using radiative MHD simulations a scenario where magnetic flux is transported from the magnetic canopy overlying the internetwork into the photosphere by means of downward plumes associated with convective overshoot. We find that if a canopy-like magnetic field is present in the simulation, the transport of flux from the canopy is an important process for seeding the photospheric layers of the internetwork with magnetic field. We propose that this mechanism is relevant for the Sun as well, and it could naturally explain the observed internetwork linear polarization patches entirely embedded in downflows.

preprint2011arXiv

Unnoticed magnetic field oscillations in the very quiet Sun revealed by Sunrise/IMaX

We present observational evidence for oscillations of magnetic flux density in the quiet areas of the Sun. The majority of magnetic fields on the solar surface have strengths of the order of or lower than the equipartition field (300-500 G). This results in a myriad of magnetic fields whose evolution is largely determined by the turbulent plasma motions. When granules evolve they squash the magnetic field lines together or pull them apart. Here we report on the periodic deformation of the shapes of features in circular polarization observed at high resolution with Sunrise. In particular, we note that the area of patches with constant magnetic flux oscillates with time, which implies that the apparent magnetic field intensity oscillates in antiphase. The periods associated to this oscillatory pattern is compatible with the granular life-time and change abruptly, which suggests that these oscillations might not correspond to characteristic oscillatory modes of magnetic structures, but to the forcing by granular motions. In one particular case, we find three patches around the same granule oscillating in phase, which means that the spatial coherence of these oscillations can reach 1600 km. Interestingly, the same kind of oscillatory phenomenon is found also in the upper photosphere.

preprint2010arXiv

A homogeneous database of sunspot areas covering more than 130 years

The historical record of sunspot areas is a valuable and widely used proxy of solar activity and variability. The Royal Greenwich Observatory (RGO) regularly measured this and other parameters between 1874 and 1976. After that time records from a number of different observatories are available. These, however, show systematic differences and often have significants gaps. Our goal is to obtain a uniform and complete sunspot area time series by combining different data sets. A homogeneus composite of sunspot areas is essential for different applications in solar physics, among others for irradiance reconstructions. Data recorded simultaneously at different observatories are statistically compared in order to determine the intercalibration factors. Using these data we compile a complete and cross-calibrated time series. The Greenwich data set is used as a basis until 1976, the Russian data (a compilation of observations made at stations in the former USSR) between 1977 and 1985 and data compiled by the USAF network since 1986. Other data sets (Rome, Yunnan, Catania) are used to fill up the remaining gaps. Using the final sunspot areas record the Photometric Sunspot Index is calculated. We also show that the use of uncalibrated sunspot areas data sets can seriously affect the estimate of irradiance variations. Our analysis implies that there is no basis for the claim that UV irradiance variations have a much smaller influence on climate than total solar irradiance variations.

preprint2010arXiv

Bright points in the quiet Sun as observed in the visible and near-UV by the balloon-borne observatory Sunrise

Bright points (BPs) are manifestations of small magnetic elements in the solar photosphere. Their brightness contrast not only gives insight into the thermal state of the photosphere (and chromosphere) in magnetic elements, but also plays an important role in modulating the solar total and spectral irradiance. Here we report on simultaneous high-resolution imaging and spectropolarimetric observations of BPs using Sunrise balloon-borne observatory data of the quiet Sun at disk center. BP contrasts have been measured between 214 nm and 525 nm, including the first measurements at wavelengths below 388 nm. The histograms of the BP peak brightness show a clear trend toward broader contrast distributions and higher mean contrasts at shorter wavelengths. At 214 nm we observe a peak brightness of up to five times the mean quiet-Sun value, the highest BP contrast so far observed. All BPs are associated with a magnetic signal, although in a number of cases it is surprisingly weak. Most of the BPs show only weak downflows, the mean value being 240 m/s, but some display strong down- or upflows reaching a few km/s.

preprint2010arXiv

Detection of large acoustic energy flux in the solar atmosphere

We study the energy flux carried by acoustic waves excited by convective motions at sub-photospheric levels. The analysis of high-resolution spectropolarimetric data taken with IMaX/Sunrise provides a total energy flux of ~ 6400--7700 Wm$^{-2}$ at a height of ~ 250 km in the 5.2-10 mHz range, i.e. at least twice the largest energy flux found in previous works. Our estimate lies within a factor of 2 of the energy flux needed to balance radiative losses from the chromosphere according to Anderson & Athay (1989) and revives interest in acoustic waves for transporting energy to the chromosphere. The acoustic flux is mainly found in the intergranular lanes but also in small rapidly-evolving granules and at the bright borders, forming dark dots and lanes of splitting granules.

preprint2010arXiv

Detection of vortex tubes in solar granulation from observations with Sunrise

We have investigated a time series of continuum intensity maps and corresponding Dopplergrams of granulation in a very quiet solar region at the disk center, recorded with the Imaging Magnetograph eXperiment (IMaX) on board the balloon-borne solar observatory Sunrise. We find that granules frequently show substructure in the form of lanes composed of a leading bright rim and a trailing dark edge, which move together from the boundary of a granule into the granule itself. We find strikingly similar events in synthesized intensity maps from an ab initio numerical simulation of solar surface convection. From cross sections through the computational domain of the simulation, we conclude that these `granular lanes' are the visible signature of (horizontally oriented) vortex tubes. The characteristic optical appearance of vortex tubes at the solar surface is explained. We propose that the observed vortex tubes may represent only the large-scale end of a hierarchy of vortex tubes existing near the solar surface.

preprint2010arXiv

Discriminant analysis of solar bright points and faculae I. Classification method and center-to-limb distribution

While photospheric magnetic elements appear mainly as Bright Points (BPs) at the disk center and as faculae near the limb, high-resolution images reveal the coexistence of BPs and faculae over a range of heliocentric angles. This is not explained by a "hot wall" effect through vertical flux tubes, and suggests that the transition from BPs to faculae needs to be quantitatively investigated. To achieve this, we made the first recorded attempt to discriminate BPs and faculae, using a statistical classification approach based on Linear Discriminant Analysis(LDA). This paper gives a detailed description of our method, and shows its application on high-resolution images of active regions to retrieve a center-to-limb distribution of BPs and faculae. Bright "magnetic" features were detected at various disk positions by a segmentation algorithm using simultaneous G-band and continuum information. By using a selected sample of those features to represent BPs and faculae, suitable photometric parameters were identified in order to carry out LDA. We thus obtained a Center-to-Limb Variation (CLV) of the relative number of BPs and faculae, revealing the predominance of faculae at all disk positions except close to disk center (mu > 0.9). Although the present dataset suffers from limited statistics, our results are consistent with other observations of BPs and faculae at various disk positions. The retrieved CLV indicates that at high resolution, faculae are an essential constituent of active regions all across the solar disk. We speculate that the faculae near disk center as well as the BPs away from disk center are associated with inclined fields.

preprint2010arXiv

Expansion of magnetic flux concentrations: a comparison of Hinode SOT d ata and models

Context: The expansion of network magnetic fields with height is a fundamental property of flux tube models. A rapid expansion is required to form a magnetic canopy. Aims: We characterize the observed expansion properties of magnetic network elements and compare them with the thin flux tube and sheet approximations, as well as with magnetoconvection simulations. Methods: We used data from the Hinode SOT NFI NaD1 channel and spectropolarimeter to study the appearance of magnetic flux concentrations seen in circular polarization as a function of position on the solar disk. We compared the observations with synthetic observables from models based on the thin flux tube approximation and magnetoconvection simulations with two different upper boundary conditions for the magnetic field (potential and vertical). Results: The observed circular polarization signal of magnetic flux concentrations changes from unipolar at disk center to bipolar near the limb, which implies an expanding magnetic field. The observed expansion agrees with expansion properties derived from the thin flux sheet and tube approximations. Magnetoconvection simulations with a potential field as the upper boundary condition for the magnetic field also produce bipolar features near the limb while a simulation with a vertical field boundary condition does not. Conclusions: The near-limb apparent bipolar magnetic features seen in high-resolution Hinode observations can be interpreted using a simple flux sheet or tube model. This lends further support to the idea that magnetic features with vastly varying sizes have similar relative expansion rates. The numerical simulations presented here are less useful in interpreting the expansion since the diagnostics we are interested in are strongly influenced by the choice of the upper boundary condition for the magnetic field in the purely photospheric simulations.

preprint2010arXiv

Exploring the magnetic topologies of cool stars

Magnetic fields of cool stars can be directly investigated through the study of the Zeeman effect on photospheric spectral lines using several approaches. With spectroscopic measurement in unpolarised light, the total magnetic flux averaged over the stellar disc can be derived but very little information on the field geometry is available. Spectropolarimetry provides a complementary information on the large-scale component of the magnetic topology. With Zeeman-Doppler Imaging (ZDI), this information can be retrieved to produce a map of the vector magnetic field at the surface of the star, and in particular to assess the relative importance of the poloidal and toroidal components as well as the degree of axisymmetry of the field distribution. The development of high-performance spectropolarimeters associated with multi-lines techniques and ZDI allows us to explore magnetic topologies throughout the Hertzsprung-Russel diagram, on stars spanning a wide range of mass, age and rotation period. These observations bring novel constraints on magnetic field generation by dynamo effect in cool stars. In particular, the study of solar twins brings new insight on the impact of rotation on the solar dynamo, whereas the detection of strong and stable dipolar magnetic fields on fully convective stars questions the precise role of the tachocline in this process.

preprint2010arXiv

Fully resolved quiet-Sun magnetic flux tube observed with the Sunrise IMaX instrument

Until today, the small size of magnetic elements in quiet Sun areas has required the application of indirect methods, such as the line-ratio technique or multi-component inversions, to infer their physical properties. A consistent match to the observed Stokes profiles could only be obtained by introducing a magnetic filling factor that specifies the fraction of the observed pixel filled with magnetic field. Here, we investigate the properties of a small magnetic patch in the quiet Sun observed with the IMaX magnetograph on board the balloon-borne telescope Sunrise with unprecedented spatial resolution and low instrumental stray light. We apply an inversion technique based on the numerical solution of the radiative transfer equation to retrieve the temperature stratification and the field strength in the magnetic patch. The observations can be well reproduced with a one-component, fully magnetized atmosphere with a field strength exceeding 1 kG and a significantly enhanced temperature in the mid- to upper photosphere with respect to its surroundings, consistent with semi-empirical flux tube models for plage regions. We therefore conclude that, within the framework of a simple atmospheric model, the IMaX measurements resolve the observed quiet-Sun flux tube.

preprint2010arXiv

Intensity contrast from MHD simulations and from HINODE observations

Changes in the solar surface area covered by small-scale magnetic elements are thought to cause long-term changes in the solar spectral irradiance, which are important for determining the impact on Earth's climate. To study the effect of small-scale magnetic elements on total and spectral irradiance, we derive their contrasts from 3-D MHD simulations of the solar atmosphere. Such calculations are necessary since measurements of small-scale flux tube contrasts are confined to a few wavelengths and suffer from scattered light and instrument defocus, even for space observations. To test the contrast calculations, we compare rms contrasts from simulations with those obtained with the broad-band filter imager mounted on the Solar Optical Telescope (SOT) onboard the Hinode satellite and also analyse centre-to-limb variations (CLV). The 3-D MHD simulations include the interaction between convection and magnetic flux tubes. They have been run with non-grey radiative transfer using the MURaM code. Simulations have an average vertical magnetic field of 0G, 50G, and 200G. Emergent intensities are calculated with the spectral synthesis code ATLAS9 and are convolved with a theoretical point-spread function to account for the properties of the observations' optical system. We find reasonable agreement for simulated and observed intensity distributions in the visible continuum bands. Agreement is poorer for the CN and G-Bands. The analysis of the simulations exhibits a potentially more realistic centre-to-limb behaviour than calculations based on 1-D model atmospheres. We conclude that starting from 3-D MHD simulations represents a powerful approach to obtaining intensity contrasts for a wide wavelength coverage and for different positions on the solar disk. This also paves the way for future calculations of facular and network contrast as a function of magnetic fluxes.

preprint2010arXiv

Long-term magnetic field monitoring of the Sun-like star Ksi Bootis A

Phase-resolved observations of the solar-type star Ksi Bootis A were obtained using the NARVAL spectropolarimeter at the Telescope Bernard Lyot (Pic du Midi, France) during years 2007, 2008, 2009 and 2010. The data sets enable us to study both the rotational modulation and the long-term evolution of various magnetic and activity tracers. Here, we focus on the large-scale photospheric magnetic field (reconstructed by Zeeman-Doppler Imaging), the Zeeman broadening of the FeI 846.84 nm magnetic line, and the chromospheric CaII H and H alpha emission.

preprint2010arXiv

Mesogranulation and the solar surface magnetic field distribution

The relation of the solar surface magnetic field with mesogranular cells is studied using high spatial (~ 100 km) and temporal (~ 30 sec) resolution data obtained with the IMaX instrument aboard SUNRISE. First, mesogranular cells are identified using Lagrange tracers (corks) based on horizontal velocity fields obtained through Local Correlation Tracking. After ~ 20 min of integration, the tracers delineate a sharp mesogranular network with lanes of width below about 280 km. The preferential location of magnetic elements in mesogranular cells is tested quantitatively. Roughly 85% of pixels with magnetic field higher than 100 G are located in the near neighborhood of mesogranular lanes. Magnetic flux is therefore concentrated in mesogranular lanes rather than intergranular ones. Secondly, magnetic field extrapolations are performed to obtain field lines anchored in the observed flux elements. This analysis, therefore, is independent of the horizontal flows determined in the first part. A probability density function (PDF) is calculated for the distribution of distances between the footpoints of individual magnetic field lines. The PDF has an exponential shape at scales between 1 and 10 Mm, with a constant characteristic decay distance, indicating the absence of preferred convection scales in the mesogranular range. Our results support the view that mesogranulation is not an intrinsic convective scale (in the sense that it is not a primary energy-injection scale of solar convection), but also give quantitative confirmation that, nevertheless, the magnetic elements are preferentially found along mesogranular lanes.

preprint2010arXiv

Observations of a rotating macrospicule associated with an X-ray jet

We attempt to understand the driving mechanism of a macrospicule and its relationship with a coronal jet. We study the dynamics of a macrospicule and an associated coronal jet captured by multi-spacecraft observations. Doppler velocities both in the macrospicule and the coronal jet are determined by EIS and SUMER spectra. Their temporal evolution is studied using X-ray and He II 304 images. A blueshift of -120+/-15 km/s is detected on one side of the macrospicule, while a redshift of 50+/-6 km/s is found at the base of the other side. The inclination angle of the macrospicule inferred from a stereoscopic analysis with STEREO suggests that the measured Doppler velocities can be attributed to a rotating motion of the macrospicule rather than a radial flow or an expansion. The macrospicule is driven by the unfolding motion of a twisted magnetic flux rope, while the associated X-ray jet is a radial outflow.

preprint2010arXiv

Probing quiet Sun magnetism using MURaM simulations and Hinode/SP results: support for a local dynamo

We obtain information about the magnetic flux present in the quiet Sun by comparing radiative MHD simulations with Hinode/SP observations, with particular emphasis on the role of surface dynamo action. Simulation runs with different magnetic Reynolds numbers (Rm) are used together with observations at different heliocentric angles with different levels of noise. The results show that simulations with an imposed mixed-polarity field and Rm below the threshold for dynamo action reproduce the observed vertical flux density, but do not display a sufficiently high horizontal flux density. Surface dynamo simulations at the highest Rm feasible at the moment yield a ratio of the horizontal and vertical flux density consistent with observational results, but the overall amplitudes are too low. Based on the properties of the local dynamo simulations, a tentative scaling of the magnetic field strength by a factor 2 - 3 reproduces the signal observed in the internetwork regions. We find an agreement with observations at different heliocentric angles. The mean field strength in internetwork, implied by our analysis, is roughly 170 G at the optical depth unity. Our study shows that surface dynamo could be responsible for most of the magnetic flux in the quiet Sun outside the network given that the extrapolation to higher Rm is valid.

preprint2010arXiv

Radiative emission of solar features in the Ca II K line: comparison of measurements and models

We study the radiative emission of various types of solar features, such as quiet Sun, enhanced network, plage, and bright plage regions, identified on filtergrams taken in the Ca II K line. We analysed fulldisk images obtained with the PSPT, by using three interference filters that sample the Ca II K line with different bandpasses. We studied the dependence of the radiative emission of disk features on the filter bandpass. We also performed a NLTE spectral synthesis of the Ca II K line integrated over the bandpass of PSPT filters. The synthesis was carried out by utilizing both the PRD and CRD with the most recent set of semi empirical atmosphere models in the literature and some earlier atmosphere models. We measured the CLV of intensity values for various solar features identified on PSPT images and compared the results obtained with those derived from the synthesis. We find that CRD calculations derived using the most recent quiet Sun model, on average, reproduce the measured values of the quiet Sun regions slightly more accurately than PRD computations with the same model. This may reflect that the utilized atmospheric model was computed assuming CRD. Calculations with PRD on earlier quiet Sun model atmospheres reproduce measured quantities with a similar accuracy as to that achieved here by applying CRD to the recent model. We also find that the median contrast values measured for most of the identified bright features, disk positions, and filter widths are, on average, a factor 1.9 lower than those derived from PRD simulations performed using the recent bright feature models. The discrepancy between measured and modeled values decreases by 12% after taking into account straylight effects on PSPT images. PRD computations on either the most recent or the earlier atmosphere models of bright features reproduce measurements from plage and bright plage regions with a similar accuracy.

preprint2010arXiv

Sideways displacement of penumbral fibrils by the solar flare on 2006 December 13

Flares are known to restructure the magnetic field in the corona and to accelerate the gas between the field lines, but their effect on the photosphere is less well studied. New data of the Solar Optical Telescope (SOT) onboard Hinode provide unprecedented opportunity to uncover the photospheric effect of a solar flare, which associates with an active region NOAA AR 10930 on 2006 December 13. We find a clear lateral displacement of sunspot penumbral regions scanned by two flare ribbons. In the impulsive phase of the flare, the flare ribbons scan the sunspot at a speed of around 18 km/s, derived from Ca II and G-band images. We find instantaneous horizontal shear of penumbral fibrils, with initial velocities of about 1.6 km/s, produced when a flare ribbon passes over them. This velocity decreases rapidly at first, then gradually decays, so that about one hour later, the fibrils return to a new equilibrium. During the one hour interval, the total displacement of these fibrils is around 2.0 Mm, with an average shear velocity of 0.55 km/s. This lateral motion of the penumbral fibrils indicates that the magnetic footpoints of these field lines being rearranged in the corona also move.

preprint2010arXiv

Sunrise: instrument, mission, data and first results

The Sunrise balloon-borne solar observatory consists of a 1m aperture Gregory telescope, a UV filter imager, an imaging vector polarimeter, an image stabilization system and further infrastructure. The first science flight of Sunrise yielded high-quality data that reveal the structure, dynamics and evolution of solar convection, oscillations and magnetic fields at a resolution of around 100 km in the quiet Sun. After a brief description of instruments and data, first qualitative results are presented. In contrast to earlier observations, we clearly see granulation at 214 nm. Images in Ca II H display narrow, short-lived dark intergranular lanes between the bright edges of granules. The very small-scale, mixed-polarity internetwork fields are found to be highly dynamic. A significant increase in detectable magnetic flux is found after phase-diversity-related reconstruction of polarization maps, indicating that the polarities are mixed right down to the spatial resolution limit, and probably beyond.

preprint2010arXiv

SUNRISE/IMaX observations of convectively driven vortex flows in the Sun

We characterize the observational properties of the convectively driven vortex flows recently discovered on the quiet Sun, using magnetograms, Dopplergrams and images obtained with the 1-m balloon-borne Sunrise telescope. By visual inspection of time series, we find some 3.1e-3 vortices/(Mm^2 min), which is a factor of 1.7 larger than previous estimates. The mean duration of the individual events turns out to be 7.9 min, with a standard deviation of 3.2 min. In addition, we find several events appearing at the same locations along the duration of the time series (31.6 min). Such recurrent vortices show up in the proper motion flow field map averaged over the time series. The typical vertical vorticities are <= 6e-3 1/sec, which corresponds to a period of rotation of some 35 min. The vortices show a preferred counterclockwise sense of rotation, which we conjecture may have to do with the preferred vorticity impinged by the solar differential rotation.

preprint2010arXiv

Supersonic Magnetic Upflows in Granular Cells Observed with Sunrise/IMaX

Using the IMaX instrument on-board the Sunrise stratospheric balloon-telescope we have detected extremely shifted polarization signals around the Fe I 5250.217 Å spectral line within granules in the solar photosphere. We interpret the velocities associated with these events as corresponding to supersonic and magnetic upflows. In addition, they are also related to the appearance of opposite polarities and highly inclined magnetic fields. This suggests that they are produced by the reconnection of emerging magnetic loops through granular upflows. The events occupy an average area of 0.046 arcsec$^2$ and last for about 80 seconds, with larger events having longer lifetimes. These supersonic events occur at a rate of $1.3\times10^{-5}$ occurrences per second per arcsec$^{2}$.

preprint2010arXiv

Surface waves in solar granulation observed with {\sc Sunrise}

Solar oscillations are expected to be excited by turbulent flows in the intergranular lanes near the solar surface. Time series recorded by the IMaX instrument aboard the {\sc Sunrise} observatory reveal solar oscillations at high resolution, which allow studying the properties of oscillations with short wavelengths. We analyze two times series with synchronous recordings of Doppler velocity and continuum intensity images with durations of 32\thinspace min and 23\thinspace min, resp., recorded close to the disk center of the Sun to study the propagation and excitation of solar acoustic oscillations. In the Doppler velocity data, both the standing acoustic waves and the short-lived, high-degree running waves are visible. The standing waves are visible as temporary enhancements of the amplitudes of the large-scale velocity field due to the stochastic superposition of the acoustic waves. We focus on the high-degree small-scale waves by suitable filtering in the Fourier domain. Investigating the propagation and excitation of $f$- and $p_1$-modes with wave numbers $k > 1.4$\thinspace 1/Mm we find that also exploding granules contribute to the excitation of solar $p$-modes in addition to the contribution of intergranular lanes.

preprint2010arXiv

The Filter Imager SuFI and the Image Stabilization and Light Distribution System ISLiD of the Sunrise Balloon-Borne Observatory: Instrument Description

We describe the design of the Sunrise Filter Imager (SuFI) and the Image Stabilization and Light Distribution (ISLiD) unit onboard the Sunrise balloon borne solar observatory. This contribution provides the necessary information which is relevant to understand the instruments working principles, the relevant technical data, and the necessary information about calibration issues directly related to the science data.

preprint2010arXiv

The Imaging Magnetograph eXperiment (IMaX) for the Sunrise balloon-borne solar observatory

The Imaging Magnetograph eXperiment (IMaX) is a spectropolarimeter built by four institutions in Spain that flew on board the Sunrise balloon-borne telesocope in June 2009 for almost six days over the Arctic Circle. As a polarimeter IMaX uses fast polarization modulation (based on the use of two liquid crystal retarders), real-time image accumulation, and dual beam polarimetry to reach polarization sensitivities of 0.1%. As a spectrograph, the instrument uses a LiNbO3 etalon in double pass and a narrow band pre-filter to achieve a spectral resolution of 85 mAA. IMaX uses the high Zeeman sensitive line of Fe I at 5250.2 AA and observes all four Stokes parameters at various points inside the spectral line. This allows vector magnetograms, Dopplergrams, and intensity frames to be produced that, after reconstruction, reach spatial resolutions in the 0.15-0.18 arcsec range over a 50x50 arcsec FOV. Time cadences vary between ten and 33 seconds, although the shortest one only includes longitudinal polarimetry. The spectral line is sampled in various ways depending on the applied observing mode, from just two points inside the line to 11 of them. All observing modes include one extra wavelength point in the nearby continuum. Gauss equivalent sensitivities are four Gauss for longitudinal fields and 80 Gauss for transverse fields per wavelength sample. The LOS velocities are estimated with statistical errors of the order of 5-40 m/s. The design, calibration and integration phases of the instrument, together with the implemented data reduction scheme are described in some detail.

preprint2010arXiv

The Wave-Front Correction System for the Sunrise Balloon-Borne Solar Observatory

This paper describes the wave-front correction system developed for the Sunrise balloon telescope, and provides information about its in-flight performance. For the correction of low-order aberrations, a Correlating Wave-Front Sensor (CWS) was used. It consisted of a six-element Shack-Hartmann wave-front sensor (WFS), a fast tip-tilt mirror for the compensation of image motion, and an active telescope secondary mirror for focus correction. The CWS delivered a stabilized image with a precision of 0.04 arcsec (rms), whenever the coarse pointing was better than 90 arcsec peak-to-peak. The automatic focus adjustment maintained a focus stability of 0.01 waves in the focal plane of the CWS. During the 5.5 day flight, good image quality and stability was achieved during 33 hours, containing 45 sequences that lasted between 10 and 45 minutes.

preprint2010arXiv

Transverse component of the magnetic field in the solar photosphere observed by Sunrise

We present the first observations of the transverse component of photospheric magnetic field acquired by the imaging magnetograph Sunrise/IMaX. Using an automated detection method, we obtain statistical properties of 4536 features with significant linear polarization signal. Their rate of occurrence is 1-2 orders of magnitude larger than values reported by previous studies. We show that these features have no characteristic size or lifetime. They appear preferentially at granule boundaries with most of them being caught in downflow lanes at some point in their evolution. Only a small percentage are entirely and constantly embedded in upflows (16%) or downflows (8%).

preprint2010arXiv

Where the granular flows bend

Based on IMaX/Sunrise data, we report on a previously undetected phenomenon in solar granulation. We show that in a very narrow region separating granules and intergranular lanes the spectral line width of the Fe I 5250.2 A line becomes extremely small. We offer an explanation of this observation with the help of magneto-convection simulations. These regions with extremely small line widths correspond to the places where the granular flows bend from mainly upflow in granules to downflow in intergranular lanes. We show that the resolution and image stability achieved by IMaX/Sunrise are important requisites to detect this interesting phenomenon.

preprint2009arXiv

A comparison of measured and simulated solar network contrast

Long-term trends in the solar spectral irradiance are important to determine the impact on Earth's climate. These long-term changes are thought to be caused mainly by changes in the surface area covered by small-scale magnetic elements. The direct measurement of the contrast to determine the impact of these small-scale magnetic elements is, however, limited to a few wavelengths, and is, even for space instruments, affected by scattered light and instrument defocus. In this work we calculate emergent intensities from 3-D simulations of solar magneto-convection and validate the outcome by comparing with observations from Hinode/SOT. In this manner we aim to construct the contrast at wavelengths ranging from the NUV to the FIR.

preprint2009arXiv

ACRIM-gap and total solar irradiance revisited: Is there a secular trend between 1986 and 1996?

A gap in the total solar irradiance (TSI) measurements between ACRIM-1 and ACRIM-2 led to the ongoing debate on the presence or not of a secular trend between the minima preceding cycles 22 (in 1986) and 23 (1996). It was recently proposed to use the SATIRE model of solar irradiance variations to bridge this gap. When doing this, it is important to use the appropriate SATIRE-based reconstruction, which we do here, employing a reconstruction based on magnetograms. The accuracy of this model on months to years timescales is significantly higher than that of a model developed for long-term reconstructions used by the ACRIM team for such an analysis. The constructed `mixed' ACRIM - SATIRE composite shows no increase in the TSI from 1986 to 1996, in contrast to the ACRIM TSI composite.

preprint2009arXiv

Bright fibrils in Ca II K

Context: Except for the Ca II resonance lines, fibrils are ubiquitously present in most high-resolution observations of chromospheric lines. Aims: We show that fibrils are also a prevailing feature in Ca II K, provided the spatial-resolution is sufficiently high. Methods: We present high spatial resolution observations of an active region in the Ca I} K line from the Swedish Solar Telescope. Through a comparison between photospheric intensity and magnetic field data, we study the connection between bright chromospheric fibrils and photospheric structures. Additionally, using Fourier analysis we study how the fibrils are linked to the observed dynamics. Results: We find that very narrow, bright fibrils are a prevailing feature over large portions of the observed field. We also find a clear connection between the fibril footpoints and photospheric magnetic features. We show that the fibrils play two distinct roles in the observed dynamics: depending on their location they can act as a canopy suppressing oscillations or they can channel low-frequency oscillations into the chromosphere. Conclusions: The Ca II K fibrils share many characteristics with fibrils observed in other chromospheric lines, but some features, such as the very small widths, are unique to these observations.

preprint2009arXiv

Comparison of the thin flux tube approximation with 3D MHD simulations

The structure and dynamics of small vertical photospheric magnetic flux concentrations has been often treated in the framework of an approximation based upon a low-order truncation of the Taylor expansions of all quantities in the horizontal direction, together with the assumption of instantaneous total pressure balance at the boundary to the non-magnetic external medium. Formally, such an approximation is justified if the diameter of the structure (a flux tube or a flux sheet) is small compared to all other relevant length scales (scale height, radius of curvature, wavelength, etc.). The advent of realistic 3D radiative MHD simulations opens the possibility of checking the consistency of the approximation with the properties of the flux concentrations that form in the course of a simulation. We carry out a comparative analysis between the thin flux tube/sheet models and flux concentrations formed in a 3D radiation-MHD simulation. We compare the distribution of the vertical and horizontal components of the magnetic field in a 3D MHD simulation with the field distribution in the case of the thin flux tube/sheet approximation. We also consider the total (gas plus magnetic) pressure in the MHD simulation box. Flux concentrations with super-equipartition fields are reasonably well reproduced by the second-order thin flux tube/sheet approximation. The differences between approximation and simulation are due to the asymmetry and the dynamics of the simulated structures.

preprint2009arXiv

Convective motions and net circular polarization in sunspot penumbrae

We have employed a penumbral model, that includes the Evershed flow and convective motions inside penumbral filaments, to reproduce the azimuthal variation of the net circular polarization (NCP) in sunspot penumbrae at different heliocentric angles for two different spectral lines. The theoretical net circular polarization fits the observations as satisfactorily as penumbral models based on flux-tubes. The reason for this is that the effect of convective motions on the NCP is very small compared to the effect of the Evershed flow. In addition, the NCP generated by convective upflows cancels out the NCP generated by the downflows. We have also found that, in order to fit the observed NCP, the strength of the magnetic field inside penumbral filaments must be very close to 1000 G. In particular, field-free or weak-field filaments fail to reproduce both the correct sign of the net circular polarization, as well as its dependence on the azimuthal and heliocentric angles.

preprint2009arXiv

Magnetic field intensification: comparison of 3D MHD simulations with Hinode/SP results

Recent spectro-polarimetric observations have provided detailed measurements of magnetic field, velocity and intensity during events of magnetic field intensification in the solar photosphere. We consider the temporal evolution of the relevant physical quantities for three cases of magnetic field intensification in a numerical simulation. We determine the evolution of the intensity, magnetic flux density and zero-crossing velocity derived from the synthetic Stokes parameters by taking into account the spectral and spatial resolution of the spectropolarimeter (SP) on board Hinode. The three events considered show a similar evolution: advection of magnetic flux to a granular vertex, development of a strong downflow, evacuation of the magnetic feature, increase of the field strength and the appearance of the bright point. We find that synthetic and real observations are qualitatively consistent and, for one of the cases considered, agree very well also quantitatively. The effect of finite resolution (spatial smearing) is most pronounced in the case of small features, for which the synthetic Hinode/SP observations miss the bright point formation and also the high-velocity downflows during the formation of the smaller magnetic features.

preprint2009arXiv

Nonlinear force-free modelling: influence of inaccuracies in the measured magnetic vector

Context: Solar magnetic fields are regularly extrapolated into the corona starting from photospheric magnetic measurements that can suffer from significant uncertainties. Aims: Here we study how inaccuracies introduced into the maps of the photospheric magnetic vector from the inversion of ideal and noisy Stokes parameters influence the extrapolation of nonlinear force-free magnetic fields. Methods: We compute nonlinear force-free magnetic fields based on simulated vector magnetograms, which have been produced by the inversion of Stokes profiles, computed froma 3-D radiation MHD simulation snapshot. These extrapolations are compared with extrapolations starting directly from the field in the MHD simulations, which is our reference. We investigate how line formation and instrumental effects such as noise, limited spatial resolution and the effect of employing a filter instrument influence the resulting magnetic field structure. The comparison is done qualitatively by visual inspection of the magnetic field distribution and quantitatively by different metrics. Results: The reconstructed field is most accurate if ideal Stokes data are inverted and becomes less accurate if instrumental effects and noise are included. The results demonstrate that the non-linear force-free field extrapolation method tested here is relatively insensitive to the effects of noise in measured polarization spectra at levels consistent with present-day instruments. Conclusions heading: Our results show that we can reconstruct the coronal magnetic field as a nonlinear force-free field from realistic photospheric measurements with an accuracy of a few percent, at least in the absence of sunspots.

preprint2009arXiv

Reconstructed and measured total solar irradiance: Is there a secular trend between 1978 and 2003?

Total solar irradiance reconstructed between 1978 and 2003 using solar surface magnetic field distributions is compared with three composites of total solar irradiance measurements. A good correspondence is found with the total solar irradiance composite from PMOD/WRC, with no bias between the three cycles. The agreement with the other composites (the ACRIM composite, mainly based on the Active Cavity Radiometer Irradiance Monitors I, II & III, and the IRMB composite from the Institut Royal Meteorologique Belgique) is significantly poorer. In particular, a secular increase in the irradiance exhibited by these composites is not present in the reconstructions. Hence any secular trend in total solar irradiance between 1978 and 2003 is not due to magnetic fields at the solar surface.

preprint2009arXiv

Simulation of a flux emergence event and comparison with observations by Hinode

We study the observational signature of flux emergence in the photosphere using synthetic data from a 3D MHD simulation of the emergence of a twisted flux tube. Several stages in the emergence process are considered. At every stage we compute synthetic Stokes spectra of the two iron lines Fe I 6301.5 Å and Fe I 6302.5 Å and degrade the data to the spatial and spectral resolution of Hinode's SOT/SP. Then, following observational practice, we apply Milne-Eddington-type inversions to the synthetic spectra in order to retrieve various atmospheric parameters and compare the results with recent Hinode observations. During the emergence sequence, the spectral lines sample different parts of the rising flux tube, revealing its twisted structure. The horizontal component of the magnetic field retrieved from the simulations is close to the observed values. The flattening of the flux tube in the photosphere is caused by radiative cooling, which slows down the ascent of the tube to the upper solar atmosphere. Consistent with the observations, the rising magnetized plasma produces a blue shift of the spectral lines during a large part of the emergence sequence.

preprint2009arXiv

Spectropolarimetric Diagnostics at the Solar Photosphere near the Limb

In the present work, we investigate the formation of Stokes profiles and spectro-polarimetric diagnostics in an active region plage near the limb. We use 3-D radiation-MHD simulations with unipolar fields of an average strength of 400 G, which is largely concentrated in flux tubes in which the field reaches typical kilo-Gauss values. We generate synthetic Stokes spectra by radiative transfer calculations, then we degrade the simulated Stokes signal to account for observational conditions. The synthetic data treated in this manner are compared with and found to roughly reproduce spectro-polarimetric high-resolution observations at Mu=0.39 obtained by the SOUP instrument with the Swedish 1-m Solar Telescope at the beginning of 2006.

preprint2009arXiv

Stray light correction and contrast analysis of Hinode broad band images

The contrasts of features in the quiet Sun are studied using filtergrams recorded by the Broad-band Filter Imager on the Hinode/Solar Optical Telescope. In a first step, the scattered light originating in the instrument is modeled using Mercury transit data. Combinations of four two-dimensional Gaussians with different widths and weights were employed to retrieve the point-spread functions (PSF) of the instrument at different wavelengths, which also describe instrumental scattered light. The parameters of PSFs at different wavelengths are tabulated. The observed images were then deconvolved using the PSFs. The corrected images were used to obtain contrasts of features such as bright points and granulation in different wavelength bands. After correction, rms contrasts of the granulation of between 0.11 (at 668 nm) and 0.22 (at 388 nm) are obtained. Similarly, bright point contrasts ranging from 0.07 (at 668 nm) to 0.78 (at 388 nm) are found, which are a factor of 1.8 to 2.8 higher than those obtained before PSF deconvolution. The mean contrast of the bright points is found to be somewhat higher in the CN-band than in the G-band, which confirms theoretical predictions.

preprint2009arXiv

Towards a long-term record of solar total and spectral irradiance

The variation of total solar irradiance (TSI) has been measured since 1978 and that of the spectral irradiance for an even shorter amount of time. Semi-empirical models are now available that reproduce over 80% of the measured irradiance variations. An extension of these models into the more distant past is needed in order to serve as input to climate simulations. Here we review our most recent efforts to model solar total and spectral irradiance on time scales from days to centuries and even longer. Solar spectral irradiance has been reconstructed since 1947. Reconstruction of solar total irradiance goes back to 1610 and suggests a value of about 1-1.5 Wm$^{-2}$ for the increase in the cycle-averaged TSI since the end of the Maunder minimum, which is significantly lower than previously assumed but agrees with other modern models. First steps have also been made towards reconstructions of solar total and spectral irradiance on time scales of millennia.

preprint2008arXiv

Identification of different types of kink modes in coronal loops: principles and application to TRACE results

We explore the possible observational signatures of different types of kink modes (horizontal and vertical oscillations in their fundamental mode and second harmonic) that may arise in coronal loops, with the aim of determining how well the individual modes can be uniquely identified from time series of images. A simple, purely geometrical model is constructed to describe the different types of kink-mode oscillations. These are then `observed' from a given direction. In particular, we employ the 3D geometrical parameters of 14 TRACE loops of transverse oscillations to try to identify the correct observed wave mode. We find that for many combinations of viewing and loop geometry it is not straightforward to distinguish between at least two types of kink modes just using time series of images. We also considered Doppler signatures and find that these can help obtain unique identifications of the oscillation modes when employed in combination with imaging. We then compare the modeled spatial signatures with the observations of 14 TRACE loops. We find that out of three oscillations previously identified as fundamental horizontal mode oscillations, two cases appear to be fundamental vertical mode oscillations (but possibly combined with the fundamental horizontal mode), and one case appears to be a combination of the fundamental vertical and horizontal modes, while in three cases it is not possible to clearly distinguish between the fundamental mode and the second-harmonic of the horizontal oscillation. In five other cases it is not possible to clearly distinguish between a fundamental horizontal mode and the second-harmonic of a vertical mode.

preprint2008arXiv

Solar axion search with the CAST experiment

The CAST (CERN Axion Solar Telescope) experiment is searching for solar axions by their conversion into photons inside the magnet pipe of an LHC dipole. The analysis of the data recorded during the first phase of the experiment with vacuum in the magnet pipes has resulted in the most restrictive experimental limit on the coupling constant of axions to photons. In the second phase, CAST is operating with a buffer gas inside the magnet pipes in order to extent the sensitivity of the experiment to higher axion masses. We will present the first results on the $^{4}{\rm He}$ data taking as well as the system upgrades that have been operated in the last year in order to adapt the experiment for the $^{3}{\rm He}$ data taking. Expected sensitivities on the coupling constant of axions to photons will be given for the recent $^{3}{\rm He}$ run just started in March 2008.

preprint2008arXiv

Stratification of sunspot umbral dots from inversion of Stokes profiles recorded by Hinode

This work aims to constrain the physical nature of umbral dots (UDs) using high-resolution spectropolarimetry. Full Stokes spectra recorded by the spectropolarimeter on Hinode of 51 UDs in a sunspot close to the disk center are analyzed. The height dependence of the temperature, magnetic field vector, and line-of-sight velocity across each UD is obtained from an inversion of the Stokes vectors of the two FeI lines at 630 nm. No difference is found at higher altitudes (-3 <= log(tau) <= -2) between the UDs and the diffuse umbral background. Below that level the difference rapidly increases, so that at the continuum formation level (log(tau) = 0) we find on average a temperature enhancement of 570 K, a magnetic field weakening of 510 G, and upflows of 800 m/s for peripheral UDs, whereas central UDs display an excess temperature of on average 550 K, a field weakening of 480 G, and no significant upflows. The results for, in particular, the peripheral UDs, including cuts of magnetic vector and velocity through them, look remarkably similar to the output of recent radiation MHD simulations. They strongly suggest that UDs are produced by convective upwellings.

preprint2007arXiv

Grand minima and maxima of solar activity: New observational constraints

Using a reconstruction of sunspot numbers stretching over multiple millennia, we analyze the statistics of the occurrence of grand minima and maxima and set new observational constraints on long-term solar and stellar dynamo models. We present an updated reconstruction of sunspot number over multiple millennia, from $^{14}$C data by means of a physics-based model, using an updated model of the evolution of the solar open magnetic flux. A list of grand minima and maxima of solar activity is presented for the Holocene (since 9500 BC) and the statistics of both the length of individual events as well as the waiting time between them are analyzed. The occurrence of grand minima/maxima is driven not by long-term cyclic variability, but by a stochastic/chaotic process. The waiting time distribution of the occurrence of grand minima/maxima deviates from an exponential distribution, implying that these events tend to cluster together with long event-free periods between the clusters. Two different types of grand minima are observed: short (30--90 years) minima of Maunder type and long ($>$110 years) minima of Spörer type, implying that a deterministic behaviour of the dynamo during a grand minimum defines its length. The duration of grand maxima follows an exponential distribution, suggesting that the duration of a grand maximum is determined by a random process. These results set new observational constraints upon the long-term behaviour of the solar dynamo.

preprint2007arXiv

Properties of solar polar coronal plumes constrained by Ultraviolet Coronagraph Spectrometer data

We investigate the plasma dynamics (outflow speed and turbulence) inside polar plumes. We compare line profiles (mainly of \ion{O}{6}) observed by the UVCS instrument on SOHO at the minimum of solar cycle 22-23 with model calculations. We consider Maxwellian velocity distributions with different widths in plume and inter-plume regions. Electron densities are assumed to be enhanced in plumes and to approach inter-plume values with increasing height. Different combinations of the outflow and turbulence velocity in the plume regions are considered. We compute line profiles and total intensities of the \ion{H}{1} Ly$α$ and the \ion{O}{6} doublets. The observed profile shapes and intensities are reproduced best by a small solar wind speed at low altitudes in plumes that increases with height to reach ambient inter-plume values above roughly 3-4 $R_\sun$ combined with a similar variation of the width of the velocity distribution of the scattering atoms/ions. We also find that plumes very close to the pole give narrow profiles at heights above 2.5 $R_\sun$, which are not observed. This suggests a tendency for plumes to be located away from the pole. We find that the inclusion of plumes in the model computations provides an improved correspondence with the observations and confirms previous results showing that published UVCS observations in polar coronal holes can be roughly reproduced without the need for large temperature anisotropy. The latitude distributions of plumes and magnetic flux distributions are studied by analyzing data from different instruments on SOHO and with SOLIS.

preprint2004arXiv

Discovery of kilogauss magnetic fields in three DA white dwarfs

We have detected longitudinal magnetic fields between 2 and 4 kG in three (WD 0446$-$790, WD 1105$-$048, WD 2359$-$434) out of a sample of 12 normal DA white dwarfs by using optical spectropolarimetry done with the VLT Antu 8 m telescope equipped with FORS1. With the exception of 40 Eri B (4 kG) these are the first positive detections of magnetic fields in white dwarfs below 30 kG. Although suspected, it was not clear whether a significant fraction of white dwarfs contain magnetic fields at this level. These fields may be explained as fossil relics from magnetic fields in the main-sequence progenitors considerably enhanced by magnetic flux conservation during the shrinkage of the core. A detection rate of 25 % (3/12) may indicate now for the first time that a substantial fraction of white dwarfs have a weak magnetic field. This result, if confirmed by future observations, would form a cornerstone for our understanding on the evolution of stellar magnetic fields. Keywords: stars: white dwarfs - stars: magnetic fields - stars: individual: WD0446-790, WD1105-048, WD2359-434