J. de la Cruz Rodriguez

CRISPRED: A data pipeline for the CRISP imaging spectropolarimeter

The production of science-ready data from major solar telescopes requires expertise beyond that of the typical observer. This is a consequence of the increasing complexity of instruments and observing sequences, which require calibrations and corrections for instrumental and seeing effects that are not only difficult to measure, but are also coupled in ways that require careful analysis in the design of the correction procedures. Modern space-based telescopes have data-processing pipelines capable of routinely producing well-characterized data products. High resolution imaging spectropolarimeters at ground-based telescopes need similar data pipelines. We present new methods for flat-fielding spectropolarimetric data acquired with telecentric Fabry-Perot instruments and a new approach for accurate camera co-alignment for image restoration. We document a procedure that forms the basis of current state-of- the-art processing of data from the CRISP imaging spectropolarimeter at the Swedish 1 m Solar Telescope (SST). By collecting, implementing, and testing a suite of computer programs, we have defined a data reduction pipeline for this instrument. This pipeline, CRISPRED, streamlines the process of making science-ready data. It is implemented and operated in IDL, with time-consuming steps delegated to C. CRISPRED will also be the basis for the data pipeline of the forthcoming CHROMIS instrument.

Non-local thermodynamic equilibrium inversions from a 3D magnetohydrodynamic chromospheric model

The structure of the solar chromosphere is believed to be governed by magnetic fields, even in quiet-Sun regions that have a relatively weak photospheric field. During the past decade inversion methods have emerged as powerful tools for analyzing the chromosphere of active regions. The applicability of inversions to infer the stratification of the physical conditions in a dynamic 3D solar chromosphere has not yet been studied in detail. This study aims to establish the diagnostic capabilities of non-local thermodynamical equilibrium (NLTE) inversion techniques of Stokes profiles induced by the Zeeman effect in the Ca II 8542 line. We computed the Ca II atomic level populations in a snapshot from a 3D radiation-MHD simulation of the quiet solar atmosphere in non-LTE using the 3D radiative transfer code Multi3d. These populations were used to compute synthetic full-Stokes profiles in the Ca II 8542 line using 1.5D radiative transfer and the inversion code Nicole. The profiles were then spectrally degraded to account for finite filter width and Gaussian noise was added to account for finite photon flux. These profiles were inverted using Nicole and the results were compared with the original model atmosphere. Our NLTE inversions applied to quiet-Sun synthetic observations provide reasonably good estimates of the chromospheric magnetic field, line-of-sight velocities and somewhat less accurate, but still very useful, estimates of the temperature. Three dimensional scattering of photons cause cool pockets in the chromosphere to be invisible in the line profile and consequently they are also not recovered by the inversions. To successfully detect Stokes linear polarization in this quiet snapshot, a noise level below 10^{-3.5} is necessary.

An Open Source, Massively Parallel Code for Non-LTE Synthesis and Inversion of Spectral Lines and Zeeman-induced Stokes Profiles

With the advent of a new generation of solar telescopes and instrumentation, interpreting chromospheric observations (in particular, spectropolarimetry) requires new, suitable diagnostic tools. This paper describes a new code, NICOLE, that has been designed for Stokes non-LTE radiative transfer, for synthesis and inversion of spectral lines and Zeeman-induced polarization profiles, spanning a wide range of atmospheric heights from the photosphere to the chromosphere. The code features a number of unique features and capabilities and has been built from scratch with a powerful parallelization scheme that makes it suitable for application on massive datasets using large supercomputers. The source code is written entirely in Fortran 90/2003 and complies strictly with the ANSI standards to ensure maximum compatibility and portability. It is being publicly released, with the idea of facilitating future branching by other groups to augment its capabilities.

Opposite polarity field with convective downflow and its relation to magnetic spines in a sunspot penumbra

We discuss NICOLE inversions of Fe I 630.15 nm and 630.25 nm Stokes spectra from a sunspot penumbra recorded with the CRISP imaging spectropolarimeter on the Swedish 1-m Solar Telescope at a spatia ...

Are solar chromospheric fibrils tracing the magnetic field

Fibrils are thin elongated features visible in the solar chromosphere in and around magnetized regions. Because of their visual appearance, they have been traditionally considered a tracer of the m ...

Solar velocity references from 3D HD photospheric models

Context. The measurement of Doppler velocities in spectroscopic solar observations requires a reference for the local frame of rest. The rotational and radial velocities of the Earth and the rotation of the Sun introduce velocity offsets in the observations. Normally, good references for velocities are missing (e.g. telluric lines), especially in filter-based spectropolarimetric observations. Aims. We determine an absolute reference for line-of-sight velocities measured from solar observations for any heliocentric angle, calibrating the convective line shift of spatially-averaged profiles on quiet sun from a 3D hydrodynamical simulation. This method works whenever there is quiet sun in the field-of-view, and it has the advantage of being relatively insensitive to uncertainties in the atomic data. Methods. We carry out radiative transfer computations in LTE for selected C i and Fe i lines, whereas the Ca ii infrared lines are synthesized in non-LTE. Radiative transfer calculations are done with a modified version of Multi, using the snapshots of a nonmagnetic 3D hydrodynamical simulation of the photosphere. Results. The resulting synthetic profiles show the expected C-shaped bisector at disk center. The degree of asymmetry and the line shifts, however, show a clear dependence on the heliocentric angle and the properties of the lines. The profiles at μ = 1 are compared with observed profiles to prove their reliability, and they are tested against errors induced by the LTE calculations, inaccuracies in the atomic data and the 3D simulation. Conclusions. Theoretical quiet-sun profiles of lines commonly used by solar observers are provided to the community. Those can be used as absolute references for line-of-sight velocities. The limb effect is produced by the projection of the 3D atmosphere along the line of sight. Non-LTE effects on Fe i lines are found to have a small impact on the convective shifts of the lines, reinforcing the usability of the LTE approximation in this case. We estimate the precision of the disk-center line shifts to be approximately 50 m s −1 , but the off-center profiles remain to be tested against observations.

SHORT DYNAMIC FIBRILS IN SUNSPOT CHROMOSPHERES

Sunspot chromospheres display vigorous oscillatory signature when observed in chromospheric diagnostics like the strong Ca II lines and H . New high-resolution sunspot observations from the Swedish 1-m Solar Telescope show the ubiquitous presence of small-scale periodic jet-like features that move up and down. This phenomenon has not been described before. Their typical width is about 0: 3 and they display clear parabolic trajectories in space-time diagrams. The maximum extension of the top of the jets is lowest in the umbra, a few 100 km, and progressively longer further away from the umbra in the penumbra, with the longest more than 1000 km. These jets resemble dynamic fibrils found in plage regions but at smaller extensions. LTE inversion of spectro-polarimetric Ca II 8542 observations enabled for a comparison of the magnetic field inclination and the properties of these short jets. We find that the most extended of these jets also have longer periods and tend to be located in regions with more horizontal magnetic fields. This is a direct observational confirmation of the mechanism of long-period waves propagating along inclined magnetic fields into the solar chromosphere. This mechanism was identified earlier as the driver of dynamic fibrils in plage, part of the mottles in quiet Sun, and type I spicules at the limb. The sunspot dynamic fibrils that we report here represent a new class of manifestation of this mechanism. They are not the same as the transient penumbral and umbral micro-jets reported earlier.

High-order aberration compensation with multi-frame blind deconvolution and phase diversity image restoration techniques

Context. For accurately measuring intensities and determining magnetic field strengths of small-scale solar (magnetic) structure, knowledge of and compensation for the point spread function is crucial. For images recorded with the Swedish 1-meter Solar Telescope (SST), restoration with multi-frame blind deconvolution (MFBD) and joint phase diverse speckle (JPDS) methods lead to remarkable improvements in image quality but granulation contrasts that are too low, indicating additional stray light. Aims. We propose a method to compensate for stray light from high-order atmospheric aberrations not included in MFBD and JPDS processing. Methods. To compensate for uncorrected aberrations, a reformulation of the image restoration process is proposed that allows the average effect of hundreds of high-order modes to be compensated for by relying on Kolmogorov statistics for these modes. The applicability of the method requires simultaneous measurements of Frieds parameter r(0). The method is tested with simulations as well as real data and extended to include compensation for conventional stray light. Results. We find that only part of the reduction of granulation contrast in SST images is due to uncompensated high-order aberrations. The remainder is still unaccounted for and attributed to stray light from the atmosphere, the telescope with its re-imaging system and to various high-altitude seeing effects. Conclusions. We conclude that statistical compensation of high-order modes is a viable method to reduce the loss of contrast occurring when a limited number of aberrations is explicitly compensated for with MFBD and JPDS processing. We show that good such compensation is possible with only 10 recorded frames. The main limitation of the method is that already MFBD and JPDS processing introduces high-order compensation that, if not taken into account, can lead to over-compensation.

DELO-Bezier Formal Solutions of the Polarized Radiative Transfer Equation

We present two new accurate and efficient methods to compute the formal solution of the polarized radiative transfer equation. In this work, the source function and the absorption matrix are approximated using quadratic and cubic Bezier spline interpolants. These schemes provide second- and third-order approximations, respectively, and do not suffer from erratic behavior of the polynomial approximation (overshooting). The accuracy and the convergence of the new method are studied along with other popular solutions of the radiative transfer equation, using stellar atmospheres with strong gradients in the line-of-sight velocity and in the magnetic-field vector.

Spectropolarimetric capabilities of Ca ii 8542 Å line

The next generation of space-and ground-based solar missions aim to study the magnetic properties of the solar chromosphere using the infrared Ca II lines and the He I 10830 angstrom line. The former seem to be the best candidates to study the stratification of magnetic fields in the solar chromosphere and their relation to the other thermodynamical properties underlying the chromospheric plasma. The purpose of this work is to provide a detailed analysis of the diagnostic capabilities of the Ca II 8542 angstrom line, anticipating forthcoming observational facilities. We study the sensitivity of the Ca II 8542 angstrom line to perturbations applied to the physical parameters of reference semi-empirical 1D model atmospheres using response functions and we make use of 3D magnetohydrodynamics simulations to examine the expected polarization signals for moderate magnetic field strengths. Our results indicate that the Ca II 8542 angstrom line is mostly sensitive to the layers enclosed in the range log tau = [0, -5.5], under the physical conditions that are present in our model atmospheres. In addition, the simulated magnetic flux tube generates strong longitudinal signals in its centre and moderate transversal signals, due to the vertical expansion of magnetic field lines, in its edge. Thus, observing the Ca II 8542 angstrom line we will be able to infer the 3D geometry of moderate magnetic field regions.