Some tests and improvements to the VFISV: Very Fast Inversion of the Stokes Vector for the Helioseismic and Magnetic Imager

Some experimental tests and improvements to the Very Fast Inversion of the Stokes Vector program, which is designed for the inversion calculation used by the Helioseismic and Magnetic Imager instrument on the Solar Dynamics Observatory,are given. On one hand, the interpolation for calculating the Voigt function is not smooth, which may occasionally cause the iteration process to converge to different minima although they are very close to initial values. This problem can be solved by a smoother interpolation. On the other hand, in order to improve the performance of this program, we have tried to abandon the randomly-jump-out strategy and set the initial value properly to avoid non-global minima. The resulting method costs only1/4 of the computational time, and will be very competitive when the users are only interested in the vectorial magnetic fields and the velocities along the line of sight.     

Magnetic Field Vector Retrieval With the Helioseismic and Magnetic Imager

We investigate the accuracy to which we can retrieve the solar photospheric magnetic field vector using the Helioseismic and Magnetic Imager (HMI) that will fly onboard of the Solar Dynamics Observatory by inverting simulated HMI profiles. The simulated profiles realistically take into account the effects of the photon noise, limited spectral resolution, instrumental polarization modulation, solar p modes, and temporal averaging. The accuracy of the determination of the magnetic field vector is studied by considering the different operational modes of the instrument.     

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: Optimization of the Spectral Line Inversion Code

The Very Fast Inversion of the Stokes Vector (VFISV) is a Milne–Eddington spectral line inversion code used to determine the magnetic and thermodynamic parameters of the solar photosphere from observations of the Stokes vector in the 6173 Å Fe i line by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We report on the modifications made to the original VFISV inversion code in order to optimize its operation within the HMI data pipeline and provide the smoothest solution in active regions. The changes either sped up the computation or reduced the frequency with which the algorithm failed to converge to a satisfactory solution. Additionally, coding bugs which were detected and fixed in the original VFISV release are reported here.     

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: Overview and Performance

The Helioseismic and Magnetic Imager (HMI) began near-continuous full-disk solar measurements on 1 May 2010 from the Solar Dynamics Observatory (SDO). An automated processing pipeline keeps pace with observations to produce observable quantities, including the photospheric vector magnetic field, from sequences of filtergrams. The basic vector-field frame list cadence is 135 seconds, but to reduce noise the filtergrams are combined to derive data products every 720 seconds. The primary 720 s observables were released in mid-2010, including Stokes polarization parameters measured at six wavelengths, as well as intensity, Doppler velocity, and the line-of-sight magnetic field. More advanced products, including the full vector magnetic field, are now available. Automatically identified HMI Active Region Patches (HARPs) track the location and shape of magnetic regions throughout their lifetime. The vector field is computed using the Very Fast Inversion of the Stokes Vector (VFISV) code optimized for the HMI pipeline; the remaining 180^° azimuth ambiguity is resolved with the Minimum Energy (ME0) code. The Milne–Eddington inversion is performed on all full-disk HMI observations. The disambiguation, until recently run only on HARP regions, is now implemented for the full disk. Vector and scalar quantities in the patches are used to derive active region indices potentially useful for forecasting; the data maps and indices are collected in the SHARP data series, hmi.sharp_720s. Definitive SHARP processing is completed only after the region rotates off the visible disk; quick-look products are produced in near real time. Patches are provided in both CCD and heliographic coordinates. HMI provides continuous coverage of the vector field, but has modest spatial, spectral, and temporal resolution. Coupled with limitations of the analysis and interpretation techniques, effects of the orbital velocity, and instrument performance, the resulting measurements have a certain dynamic range and sensitivity and are subject to systematic errors and uncertainties that are characterized in this report.     

Lossless Compression Method for the Magnetic and Helioseismic Imager(MHI) Payload

The Solar Polar-orbit Observatory (SPO), proposed by Chinese scientists, is designed to observe the solar polar regions in an unprecedented way with a spacecraft traveling in a large solar inclination angle and a small ellipticity.However, one of the most significant challenges lies in ultra-long-distance data transmission, particularly for the Magnetic and Helioseismic Imager (MHI), which is the most important payload and generates the largest volume of data in SPO. In this paper, we propose a ...     

Interpreting the Helioseismic and Magnetic Imager (HMI) Multi-Height Velocity Measurements

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe i 6173.3 Å line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multi-height velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection-simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams. We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30?–?40 km above (line center) and 30?–?40 km below (average wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic-cutoff frequency.     

Initial Results from Fitting $$-Modes Using Intensity Observations from the Helioseismic and Magnetic Imager

The Helioseismic and Magnetic Imager project recently started processing the continuum-intensity images following global helioseismology procedures similar to those used to process the velocity images. The spatial decomposition of these images has produced time series of spherical harmonic coefficients for degrees up to \(\ell=300\), using a different apodization than the one used for velocity observations. The first 360 days of observations were processed and are made available. I present initial results from fitting these time series using my fitting method and compare the derived mode characteristics to those estimated using coeval velocity observations.     

黑子磁场的Stokes光谱反演技术

对太阳大气磁场的可靠测量有助于人们更好地理解太阳活动区内外的许多活动现象,如耀斑的触发和能量释放过程、黑子的形态和黑子大气的平衡、日珥的形成等.由于原子在磁场中的一些能级会产生分裂(Zeeman效应),使对应这些能级的谱线分裂成若干个具有不同偏振特性的分量,因此目前对黑子磁场的测量主要是通过偏振光,即Stokes参量I、Q、U、V的观测来实现的.该文主要介绍近30年来太阳黑子光谱反演的方法以及所取得的成就;同时也对光谱反演和滤光器型的望远镜矢量磁场的测量进行了简单的比较.     

Transient Artifacts in a Flare Observed by the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory

The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) provides a new tool for the systematic observation of white-light flares, including Doppler and magnetic information as well as continuum. In our initial analysis of the highly impulsive $\mathrm{\gamma}$ -ray flare SOL2010-06-12T00:57 (Martínez Oliveros et al., Solar Phys. 269, 269, 2011), we reported the signature of a strong blueshift in the two footpoint sources. Concerned that this might be an artifact due to aliasing peculiar to the HMI instrument, we undertook a comparative analysis of Global Oscillation Network Group (GONG++) observations of the same flare, using the PArametric Smearing Correction ALgorithm (PASCAL) algorithm to correct for artifacts caused by variations in atmospheric smearing. This analysis confirms the artifactual nature of the apparent blueshift in the HMI observations, finding weak redshifts at the footpoints instead. We describe the use of PASCAL with GONG++ observations as a complement to the SDO observations and discuss constraints imposed by the use of HMI far from its design conditions. With proper precautions, these data provide rich information on flares and transients.     

Comparison of Line-of-Sight Magnetograms Taken by the Solar Dynamics Observatory/Helioseismic and Magnetic Imager and Solar and Heliospheric Observatory/Michelson Doppler Imager

We compare line-of-sight magnetograms from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) and the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO). The line-of-sight magnetic signal inferred from the calibrated MDI data is greater than that derived from the HMI data by a factor of 1.40. This factor varies somewhat with center-to-limb distance. An upper bound to the random noise for the 1′′ resolution HMI 720-second magnetograms is 6.3 Mx?cm^?2, and 10.2 Mx?cm^?2 for the 45-second magnetograms. Virtually no p-mode leakage is seen in the HMI magnetograms, but it is significant in the MDI magnetograms. 12-hour and 24-hour periodicities are detected in strong fields in the HMI magnetograms. The newly calibrated MDI full-disk magnetograms have been corrected for the zero-point offset and underestimation of the flux density. The noise is 26.4 Mx?cm^?2 for the MDI one-minute full-disk magnetograms and 16.2 Mx?cm^?2 for the five-minute full-disk magnetograms observed with four-arcsecond resolution. The variation of the noise over the Sun’s disk found in MDI magnetograms is likely due to the different optical distortions in the left- and right-circular analyzers, which allows the granulation and p-mode to leak in as noise. Saturation sometimes seen in sunspot umbrae in MDI magnetograms is caused by the low intensity and the limitation of the onboard computation. The noise in the HMI and MDI line-of-sight magnetic-field synoptic charts appears to be fairly uniform over the entire map. The noise is 2.3 Mx?cm^?2 for HMI charts and 5.0 Mx?cm^?2 for MDI charts. No evident periodicity is found in the HMI synoptic charts.     

Polarization Calibration of the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO)

As part of the overall ground-based calibration of the Helioseismic and Magnetic Imager (HMI) instrument an extensive set of polarimetric calibrations were performed. This paper describes the polarimetric design of the instrument, the test setup, the polarimetric model, the tests performed, and some results. It is demonstrated that HMI achieves an accuracy of 1% or better on the crosstalks between Q, U, and V and that our model can reproduce the intensities in our calibration sequences to about 0.4%. The amount of depolarization is negligible when the instrument is operated as intended which, combined with the flexibility of the polarimeter design, means that the polarimetric efficiency is excellent.     

Image Quality of the Helioseismic and Magnetic Imager (HMI) Onboard the Solar Dynamics Observatory (SDO)

We describe the imaging quality of the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) as measured during the ground calibration of the instrument. We describe the calibration techniques and report our results for the final configuration of HMI. We present the distortion, modulation transfer function, stray light, image shifts introduced by moving parts of the instrument, best focus, field curvature, and the relative alignment of the two cameras. We investigate the gain and linearity of the cameras, and present the measured flat field.     

Very Fast and Accurate Azimuth Disambiguation of Vector Magnetograms

We present a method for fast and accurate azimuth disambiguation of vector magnetogram data regardless of the location of the analyzed region on the solar disk. The direction of the transverse field is determined with the principle of minimum deviation of the field from the reference (potential) field. The new disambiguation (NDA) code is examined on the well-known models of Metcalf et al. (Solar Phys. 237, 267, 2006) and Leka et al. (Solar Phys. 260, 83, 2009), and on an artificial model based on the observed magnetic field of AR 10930 (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We compare Hinode/SOT-SP vector magnetograms of AR 10930 disambiguated with three codes: the NDA code, the nonpotential magnetic-field calculation (NPFC: Georgoulis, Astrophys. J. Lett. 629, L69, 2005), and the spherical minimum-energy method (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We then illustrate the performance of NDA on SDO/HMI full-disk magnetic-field observations. We show that our new algorithm is more than four times faster than the fastest algorithm that provides the disambiguation with a satisfactory accuracy (NPFC). At the same time, its accuracy is similar to that of the minimum-energy method (a very slow algorithm). In contrast to other codes, the NDA code maintains high accuracy when the region to be analyzed is very close to the limb.     

Wavelength Dependence of the Helioseismic and Magnetic Imager (HMI) Instrument onboard the Solar Dynamics Observatory (SDO)

The Helioseismic and Magnetic Imager (HMI) instrument will produce Doppler-velocity and vector-magnetic-field maps of the solar surface, whose accuracy is dependent on a thorough knowledge of the transmission profiles of the components of the HMI optical-filter system. Here we present a series of wavelength-dependence calibration tests, performed on the instrument from 2005 onwards, to obtain these profiles. We obtained the transmittances as a function of wavelength for the tunable and non-tunable filter elements, as well as the variation of these transmittances with temperature and the angle of incidence of rays of light. We also established the presence of fringe patterns produced by interferences inside the blocking filter and the front window, as well as a change in transmitted intensity with the tuning position. This thorough characterization of the HMI-filter system confirmed the very high quality of the instrument, and showed that its properties are well within the required specifications to produce superior data with high spatial and temporal resolution.     

Helioseismic Constraints on Rotation and Magnetic Fields in the Solar Core

In recent years, more and more precise measurements have been made of solar oscillation frequencies and line widths. From space, the Solar and Heliospheric Observatory/Michelson Doppler Imager (MDI) data has led to much progress. From the ground, networks, like Global Oscillation Network Group (GONG), Taiwanese Oscillation Network (TON), and Birmingham Solar Oscillations Network (BiSON) have also led to much progress. The sharpened and enriched oscillation spectrum of data have been critically complemented by advances in the treatments of the opacities and the equation of state. All of this has led to a significantly more precise probing of the solar core. Here we discuss the progress made and suggest how the core may be better probed with seismic data on-hand. In particular, we review our knowledge of the rotation and structure of the core. We further argue that much may be learned about the core by exploiting the line width data from the aforementioned sources. Line-width data can be used to place sharper constraints on core properties, like the degree to which the Sun rotates on a single axis and the upper limit on magnetic fields that may be buried in the core.     

Analysis of Supergranule Sizes and Velocities Using Solar Dynamics Observatory (SDO)/Helioseismic Magnetic Imager (HMI) and Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) Dopplergrams

Co-temporal Doppler images from Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) and Solar Dynamics Observatory (SDO)/Helioseismic Magnetic Imager (HMI) have been analyzed to extract quantitative information about global properties of the spatial and temporal characteristics of solar supergranulation. Preliminary comparisons show that supergranules appear to be smaller and have stronger horizontal velocity flows within HMI data than was measured with MDI. There appears to be no difference in their evolutionary timescales. Supergranule sizes and velocities were analyzed over a ten-day time period at a 15-minute cadence. While the averages of the time-series retain the aforementioned differences, fluctuations of these parameters first observed in MDI data were seen in both MDI and HMI time-series, exhibiting a strong cross-correlation. This verifies that these fluctuations are not instrumental, but are solar in origin. The observed discrepancies between the averaged values from the two sets of data are a consequence of instrument resolution. The lower spatial resolution of MDI results in larger observed structures with lower velocities than is seen in HMI. While these results offer a further constraint on the physical nature of supergranules, they also provide a level of calibration between the two instruments.     

Vector Magnetic Field Synoptic Charts from the <Emphasis Type="Italic">Helioseismic and Magnetic Imager</Emphasis> (HMI)

Vector magnetic field synoptic charts from the Helioseismic and Magnetic Imager (HMI) are now available for each Carrington Rotation (CR) starting from CR 2097 in May 2010. Synoptic charts are produced using 720-second cadence full-disk vector magnetograms remapped to Carrington coordinates. The vector field is derived from the Stokes parameters (\(I, Q, U, V\)) using a Milne–Eddington-based inversion model. The \(180^{\circ}\) azimuth ambiguity is resolved using the minimum energy algorithm for pixels in active regions and for strong-field pixels (the field is greater than about 150 G) in quiet-Sun regions. Three other methods are used for the rest of the pixels: the potential-field method, the radial acute-angle method, and the random method. The vector field synoptic charts computed using these three disambiguation methods are evaluated. The noise in the three components of the vector magnetic field is generally much higher in the potential-field method charts. The component noise levels are significantly different in the radial-acute charts. However, the noise levels in the random-method charts are lower and comparable. The assumptions used in the potential-field and radial-acute methods to disambiguate the weak transverse field introduce bias that propagates differently into the three vector-field components, leading to unreasonable pattern and artifacts, whereas the random method appears not to introduce any systematic bias. The current sheet on the source surface, computed using the potential-field source-surface model applied to random-method charts, agrees with the best solution (the result computed from the synoptic charts with the minimum energy algorithm applied to each and every pixel in the vector magnetograms) much better than the other two. Differences in the synoptic charts determined with the best method and the random method are much smaller than those from the best method and the other two. This comparison indicates that the random method is better for vector field synoptic maps computed from near-central meridian data. The vector field synoptic charts provided by the Joint Science Operations Center (JSOC) are therefore produced with the random method.     

Experimental Set-Up for Detecting Very Fast and Dispersed Millisecond Pulsars

We have developed observational and data processingtechniques for detecting millisecond pulsars. Themethod of data processing consists in correlating timeseries of data by folding them in two steps accordingto a trial period. While very time consuming,especially for very short periods, this methods allowsus in principle to detect very fast and dispersedmillisecond pulsars. The experimental set-up has beenintegrated into a data acquisition system developed atCentre national d'études spatiales (CNES).Observations of known millisecond pulsars atNançay radio observatory are used to validate thesystem and assess the sensitivity of the survey, whichis 3.5 mJy for 240 sec of integration. Futureimprovements in the radio-telescope and the systemshould bring it down to below 1 mJy. The applicationof the method to a millisecond pulsar survey is underway at Nançay radio observatory in collaborationbetween the Bureau international des poids et mesure(BIPM), Observatoire Midi-Pyrénées (OMP), and CNES.