RHESSI Data Analysis Software: Rationale and Methods

The Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) performs imaging spectroscopy of the Sun with high spatial and spectral resolution from 3 keV to 17 MeV using indirect Fourier-transform techniques. We review the rationale behind the RHESSI data analysis software, and explain the underlying structure of the software tools. Our goal was to make the large data set available within weeks after the RHESSI launch, and to make it possible for any member of the scientific community to analyze it easily. This paper describes the requirements for the software and explores our decisions to use the SolarSoftWare and Interactive Data Language programming packages, to support both Windows and Unix platforms, and to use object-oriented programming. We also describe how the data are rapidly disseminated and how ancillary data sets are used to enhance the RHESSI science. Finally, we give a schematic overview of some of the data flow through the high-level analysis tools. More information on the data and analysis procedures can be found at the RHESSI Data Center website, http://hesperia.gsfc.nasa.gov/rhessidatacenter. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022444531435     

Reconstruction of RHESSI Solar Flare Images with a Forward Fitting Method

We describe a forward-fitting method that has been developed to reconstruct hard X-ray images of solar flares from the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI), a Fourier imager with rotation-modulated collimators that was launched on 5 February 2002. The forward-fitting method is based on geometric models that represent a spatial map by a superposition of multiple source structures, which are quantified by circular gaussians (4 parameters per source), elliptical gaussians (6 parameters), or curved ellipticals (7 parameters), designed to characterize real solar flare hard X-ray maps with a minimum number of geometric elements. We describe and demonstrate the use of the forward-fitting algorithm. We perform some 500 simulations of rotation-modulated time profiles of the 9 RHESSI detectors, based on single and multiple source structures, and perform their image reconstruction. We quantify the fidelity of the image reconstruction, as function of photon statistics, and the accuracy of retrieved source positions, widths, and fluxes. We outline applications for which the forward-fitting code is most suitable, such as measurements of the energy-dependent altitude of energy loss near the limb, or footpoint separation during flares. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022469811115     

Spectral and Spatial Variations of Flare Hard X-ray Footpoints

In a sample of strong RHESSI M-class flares we have made a study of the relationship between the `hardness' of the HXR spectrum and the intensity in the 30–50 keV energy range. In all events we find clear evidence for a `soft–hard–soft' pattern of correlation between hardness and flux, on time scales as short as 10 s. We investigate whether or not this pattern is intrinsic to the acceleration mechanism. The RHESSI images in this energy range are dominated by footpoint brightenings, and we have searched for a correlation between footpoint separation velocity and spectral hardness, to be compared qualitatively with theoretical flare models. We find quite systematic footpoint motions, and also note that episodes in which footpoint separation varies rapidly often correspond with episodes of significant change in the flare spectral index, though not as the simplest flare models would predict. We report also on one of our events, on 14 March 2002, which exhibits highly sheared HXR footpoint ribbons extending over a scale of 100 arc sec. For this flare we find a correlation between footpoint motion and hard X-ray flux. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022479610710     

Energy budget and imaging spectroscopy of a compact flare

We present the analysis of a compact flare that occurred on 26 February 2002 at 10:26 UT, seen by both RHESSI and TRACE. The size of the nearly circular hard X-ray source is determined to be 5.6 (±0.8)′′, using different methods. The power-law distribution of non-thermal photons is observed to extend down to 10 keV without flattening, and to soften with increasing distance from the flare kernel. The former indicates that the energy of the precipitating flare electron population is larger than previously estimated: it amounts to 2.6 (±0.8)×10³⁰ erg above 10 keV, assuming thick-target emission. The thermal energy content of the soft X-ray source (isothermal temperature of 20.8 (±0.9) MK) and its radiated power were derived from the thermal emission at low energies. TRACE has observed a low-temperature ejection in the form of a constricted bubble, which is interpreted as a reconnection jet. Its initial energy of motion is estimated. Using data from both satellites, an energy budget for this flare is derived. The kinetic energy of the jet bulk motion and the thermal and radiated energies of the flare kernel were more than an order of magnitude smaller than the derived electron beam energy. A movie is available on the CD-ROM accompanying this volume. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022478300679     

The Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI)

RHESSI is the sixth in the NASA line of Small Explorer (SMEX) missions and the first managed in the Principal Investigator mode, where the PI is responsible for all aspects of the mission except the launch vehicle. RHESSI is designed to investigate particle acceleration and energy release in solar flares, through imaging and spectroscopy of hard X-ray/gamma-ray continua emitted by energetic electrons, and of gamma-ray lines produced by energetic ions. The single instrument consists of an imager, made up of nine bi-grid rotating modulation collimators (RMCs), in front of a spectrometer with nine cryogenically-cooled germanium detectors (GeDs), one behind each RMC. It provides the first high-resolution hard X-ray imaging spectroscopy, the first high-resolution gamma-ray line spectroscopy, and the first imaging above 100 keV including the first imaging of gamma-ray lines. The spatial resolution is as fine as ∼ 2.3 arc sec with a full-Sun (≳ 1°) field of view, and the spectral resolution is ∼ 1–10 keV FWHM over the energy range from soft X-rays (3 keV) to gamma-rays (17 MeV). An automated shutter system allows a wide dynamic range (>10⁷) of flare intensities to be handled without instrument saturation. Data for every photon is stored in a solid-state memory and telemetered to the ground, thus allowing for versatile data analysis keyed to specific science objectives. The spin-stabilized (∼ 15 rpm) spacecraft is Sun-pointing to within ∼ 0.2° and operates autonomously. RHESSI was launched on 5 February 2002, into a nearly circular, 38° inclination, 600-km altitude orbit and began observations a week later. The mission is operated from Berkeley using a dedicated 11-m antenna for telemetry reception and command uplinks. All data and analysis software are made freely and immediately available to the scientific community. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022428818870     

Relative Timing and Spectra of Solar Flare Hard X-ray Sources

Hard X-ray lightcurves, spectrograms, images, and spectra of three medium-sized flares observed by the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) are presented. Imaging spectroscopy of the 20 February 2002, 11:06 UT flare at 10′′ spatial resolution, comparable to the best previous hard X-ray imaging from Yohkoh, shows two footpoints with an ∼ 8 s delay of peak emission between footpoints. Subsequent imaging at le4′′ shows three sources consistent with two separate loops and simultaneous brightening in connected footpoints. Imaging for the simple two footpoint flare of 2 June 2002 also shows simultaneous footpoint brightening. The more complex 17 March 2002 flare shows at least four different sources during the main peak of the event, and it is difficult to clearly demonstrate simultaneous brightening of connected footpoints. Non-thermal power laws are observed down to ∼ 12–13 keV without flattening in all these events, indicating the energy content in energetic electrons may be significantly greater than previously estimated from assumed 25 keV low energy cutoff. Simultaneously brightening footpoints show similar spectra, at least in the three flares investigated. Double-power-law spectra with a relatively sharp break are often observed. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022469902940     

RHESSI Observations of the Size Scales of Solar Hard X-ray Sources

The Reuven Ramaty High-Energy Solar Spectroscopic Imager RHESSI telescope produces hard X-ray images by Fourier imaging techniques that are capable of determining the sizes and shapes of sources with spatial scales in the range ∼ 2′′–180′′. Applying the method of Unpixelized Forward Fitting to RHESSI modulation profiles from simple flares, we have identified the presence of `halo' sources whose size scale (∼ 40′′) greatly exceeds the `core' sizes (≤ 6′′–14′′). Although such `core-halo' structures have been observed at radio wavelengths using a similar technique, the radio and hard X-ray phenomena may be different. These observations raise questions about the nature of these `halos'. Among the possibilities are that they are albedo sources, thin-target loops, or unidentified diffuse structures. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022484822851     

On the Photometric Accuracy of RHESSI Imaging and Spectrosocopy

We compare the photometric accuracy of spectra and images in flares observed with the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) spacecraft. We test the accuracy of the photometry by comparing the photon fluxes obtained in different energy ranges from the spectral-fitting software SPEX with those fluxes contained in the images reconstructed with the Clean, MEM, MEM-Vis, Pixon, and Forward-fit algorithms. We quantify also the background fluxes, the fidelity of source geometries, and spatial spectra reconstructed with the five image reconstruction algorithms. We investigate the effects of grid selection, pixel size, field of view, and time intervals on the quality of image reconstruction. The detailed parameters and statistics are provided in an accompanying CD-ROM and web page. We find that Forward-fit, Pixon, and Clean have a robust convergence behavior and a photometric accuracy in the order of a few percent, while MEM does not converge optimally for large degrees of freedom (for large field of view and/or small pixel sizes), and MEM-Vis suffers in the case of time-variable sources. This comparative study documents the current status of the RHESSI spectral and imaging software, one year after launch. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/B:SOLA.0000021801.83038.aa     

Observations of Rotating Sunspots from TRACE

Recent observations from TRACE in the photospheric white-light channel have shown sunspots that rotate up to 200° about their umbral centre over a period of 3–5 days. The corresponding loops in the coronal fan are often seen to twist and can erupt as flares. In an ongoing study, seven cases of rotating sunspots have been identified, two of which can be associated with sigmoid structures appearing in Yohkoh/SXT and six with events seen by GOES. This paper analyzes the rotation rates of the sunspots using TRACE white-light data. Observations from AR 9114 are presented in detail in the main text and a summary of the results for the remaining six sunspots is presented in Appendixes A–F. Discussion of the key results, particularly common features, are presented, as well as possible mechanisms for sunspot rotation. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1026138413791     

A Multiscale Vision Model applied to analyze EIT images of the solar corona

The large dynamic range provided by the SOHO/EIT CCD (1 : 5000) is needed to observe the large EUV zoom of coronal structures from coronal homes up to flares. Histograms show that often a wide dynamic range is present in each image. Extracting hidden structures in the background level requires specific techniques such as the use of the Multiscale Vision Model (MVM, Bijaoui et al., 1998). This method, based on wavelet transformations optimizes detection of various size objects, however complex they may be. Bijaoui et al. built the Multiscale Vision Model to extract small dynamical structures from noise, mainly for studying galaxies. In this paper, we describe requirements for the use of this method with SOHO/EIT images (calibration, size of the image, dynamics of the subimage, etc.). Two different areas were studied revealing hidden structures: (1) classical coronal mass ejection (CME) formation and (2) a complex group of active regions with its evolution. The aim of this paper is to define carefully the constraints for this new method of imaging the solar corona with SOHO/EIT. Physical analysis derived from multi-wavelength observations will later complete these first results. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1017579625208 with J.P. Delaboudinière P.I.     

Global Simulation of Magnetospheric Space Weather Effects of the Bastille Day Storm

We present results from a global simulation of the interaction of the solar wind with Earth's magnetosphere, ionosphere, and thermosphere for the Bastille Day geomagnetic storm and compare the results with data. We find that during this event the magnetosphere becomes extremely compressed and eroded, causing 3 geosynchronous GOES satellites to enter the magnetosheath for an extended time period. At its extreme, the magnetopause moves at local noon as close as 4.9 R _E to Earth which is interpreted as the consequence of the combined action of enhanced dynamic pressure and strong dayside reconnection due to the strong southward interplanetary magnetic field component B _z, which at one time reaches a value of −60 nT. The lobes bulge sunward and shield the dayside reconnection region, thereby limiting the reconnection rate and thus the cross polar cap potential. Modeled ground magnetic perturbations are compared with data from 37 sub-auroral, auroral, and polar cap magnetometer stations. While the model can not yet predict the perturbations and fluctuations at individual ground stations, its predictions of the fluctuation spectrum in the 0–3 mHz range for the sub-auroral and high-latitude regions are remarkably good. However, at auroral latitudes (63° to 70° magnetic latitude) the predicted fluctuations are slightly too high. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1014228230714     

High-resolution Observations of Plasma Jets in the Solar Corona

We present recent observations of coronal jets, made by TRACE and Yohkoh/SXT on 28 May and 19 August 1998. The high spatial resolution of TRACE enables us to see in detail the process of material ejection; in the line of Fe ix (one million degrees) we see both bright emitting material and dark absorbing/scattering material being ejected, i.e., both hot and cold material, highly collimated and apparently ejected along the direction of the overlying field lines. Bright ejecta are seen simultaneously in Lyman α for one event and Yohkoh/SXT in the other. The jets on the two days are different in that the 19 August jet displays the morphology typical of a one-sided anemone jet while the 28 May jet exhibits a two-sided jet morphology. The 19 August jet shows evidence for rotation and an interesting bifurcation at large distances from the energy release site. We study the physical properties and energetics of these jetting events, and conclude that existing theoretical models capture the essential physics of the jet phenomena. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005213826793     

High-Speed x-ray Jets Associated with the 18 June 1999 Limb Flares

Peculiar high-speed X-ray jets associated with the 18 June 1999 limb flares have directly been observed with the soft X-ray Telescope (SXT) aboard Yohkoh. The jets have a much shorter lifetime (within 200 s) and a much larger velocity (∼ 1700 km s⁻¹) than the previously reported jets. Judged from their large velocity, the kinetic energy of the jets is approximately one order of magnitude larger than the thermal energy content, which is far different from the jets previously reported by other researchers. Here we present the preliminary results and discuss the particular features of the jets. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1013324126229     

Microflares and hot component in solar active regions

Open-shutter RHESSI observations of 3–15 keV X-rays are found to exhibit active-region transient brightenings and microflares at a rate of a least 10 per hour occurring even during the periods of lowest solar activity so far in the mission. A thermal component fitted by temperatures of 6–14 MK dominates from 3 keV to about 9 keV, but can be traced up to 14 keV in some cases, and has an average duration of 131(±103) s at 7–8 keV. The duration increases with decreasing photon energy. The peak count rate defined by cross-correlation is delayed at low energies. The temperature peaks early in the event and then decreases, whereas the emission measure increases throughout the event. The properties are consistent with thermal conduction dominating the evolution. In some of the bigger events, a second component was found in the 11–14 keV range extending down to 8 keV in some cases. The duration is typically 3 times shorter and ends near the peak time of the thermal component consistent with the Neupert effect of regular flares. Therefore the second component is suggested to be of non-thermal origin, presumably causing the beam-driven evaporation of the first component. The two components can be separated and analyzed in detail for the first time. Low-keV measurements allow a reliable estimate of the energy input by microflares necessary to assess their relevance for coronal heating. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022496515506     

Forward modeling of the coronal response to reconnection in an X-ray bright point

We use MDI magnetic field observations and the theory of reconnection through a separator to constrain a numerical simulation of an X-ray bright point observed in EUV by TRACE. A gasdynamic model is employed to describe the corona and transition region in the bright point loop. Nonlocal effects such as opacity and ambipolar diffusion are important to the transition region; these effects are approximated locally by modification of the radiative loss and thermal conduction. A straightforward comparison of measured light curves versus those generated by the simulation shows that the reconnection model is unable to account for the observations. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005205807984     

Spectral Energy Distribution fitting: Application to Lyα-emitting galaxies

Spectral Energy Distribution (SED) fitting is a well-developed astrophysical tool that has recently been applied to high-redshift Lyα-emitting galaxies. If rest-frame ultraviolet through near-infrared photometry is available, it allows the simultaneous determination of the star formation history and dust extinction of a galaxy. Lyα-emitter SED fitting results from the literature find star formation rates 3 M yr⁻¹, stellar masses 10⁹ M for the general population but 10¹⁰ M for the subset detected by IRAC, and very low dust extinction, A_V  0.3, although a couple of outlying analyses prefer significantly more dust and higher intrinsic star formation rates. A checklist of 14 critical choices that must be made when performing SED fitting is discussed.     

Solar Microflares Observed by SphinX and RHESSI

In 2009, the Russian Complex Orbital Observations Near-Earth of Activity of the Sun (CORONAS-Photon) spacecraft was launched, carrying the Polish Solar PHotometer In X-rays (SphinX). The SphinX was most sensitive in the spectral range 1.2?–?15 keV, thus an excellent opportunity appeared for comparison with the low-energy end of Ramaty High Energy Solar Spectroscopic Imager (RHESSI) spectra. Common spectral measurements with these instruments cover the range where most of the flare energy is accumulated. We have chosen four consecutive small solar events observed on 4 July 2009 at 13:43 UT, 13:48 UT, 13:52 UT, and 13:55 UT (RHESSI flare peak times) and used them to compare the data and results from the two instruments. Moreover, we included Geostationary Operational Environmental Satellite (GOES) records in our analysis. In practice, the range of comparison performed for SphinX and RHESSI is limited roughly to 3?–?6 keV. RHESSI fluxes measured with a use of one, four, and nine detectors in the 3?–?4 keV energy band agree with SphinX measurements. However, we observed that SphinX spectral irradiances are three times higher than those of RHESSI in the 4?–?6 keV energy band. This effect contributes to the difference in obtained emission measures, but the derived temperatures of plasma components are similar. RHESSI spectra were fitted using a model with two thermal components. We have found that the RHESSI hot component is in agreement with GOES, and the RHESSI hotter component fits the SphinX flaring component well. Moreover, we calculated the so-called thermodynamic measure and the total thermal energy content in the four microflares that we studied. The results obtained show that SphinX is a very sensitive complementary observatory for RHESSI and GOES.     

Faint UBVRI CCD sequences for wide-field surveys – II. UBVR sequences at δ=−30°

We present results from a continuing campaign to secure deep multi-colour CCD sequences for photometric calibration in UK Schmidt fields with galactic latitudes | b |>50°. In this paper we present UBVR photometry in 12 fields and BR photometry in a further 14 fields observed within UK Schmidt survey fields centred at δ =−30°. Photometric errors are at the 0.05 level at 20.2, 21, 20.5 and 20 mag for UBVR sequences respectively. Our observations were carried out with the 0.9-m telescope at the Cerro-Tololo Inter-American Observatory. These data are not intended for use as highly accurate individual photometric standards, but rather for use as sequences, using a large number of stars to calibrate wide-area data such as photographic plates. The data are available electronically at http://icstar5.ph.ic.ac.uk/~scroom/phot/photom.html .     

The Astro-WISE approach to quality control for astronomical data

We present a novel approach to quality control during the processing of astronomical data. Quality control in the Astro-WISE Information System is integral to all aspects of data handing and provides transparent access to quality estimators for all stages of data reduction from the raw image to the final catalog. The implementation of quality control mechanisms relies on the core features in this Astro-WISE Environment (AWE): an object-oriented framework, full data lineage, and both forward and backward chaining. Quality control information can be accessed via the command-line awe-prompt and the web-based Quality-WISE service. The quality control system is described and qualified using archive data from the 8-CCD Wide Field Imager (WFI) instrument (http://www.eso.org/lasilla/instruments/wfi/) on the 2.2-m MPG/ESO telescope at La Silla and (pre-)survey data from the 32-CCD OmegaCAM instrument (http://www.astro-wise.org/~omegacam/) on the VST telescope at Paranal.     

TomograPy: A Fast, Instrument-Independent, Solar Tomography Software

Solar tomography has progressed rapidly in recent years thanks to the development of robust algorithms and the availability of more powerful computers. It can today provide crucial insights in solving issues related to the line-of-sight integration present in the data of solar imagers and coronagraphs. However, there remain challenges such as the increase of the available volume of data, the handling of the temporal evolution of the observed structures, and the heterogeneity of the data in multi-spacecraft studies. We present a generic software package that can perform fast tomographic inversions that scales linearly with the number of measurements, linearly with the length of the reconstruction cube (and not the number of voxels), and linearly with the number of cores and can use data from different sources and with a variety of physical models: TomograPy ( http://nbarbey.github.com/TomograPy/ ), an open-source software freely available on the Python Package Index. For performance, TomograPy uses a parallelized-projection algorithm. It relies on the World Coordinate System standard to manage various data sources. A variety of inversion algorithms are provided to perform the tomographic-map estimation. A test suite is provided along with the code to ensure software quality. Since it makes use of the Siddon algorithm it is restricted to rectangular parallelepiped voxels but the spherical geometry of the corona can be handled through proper use of priors. We describe the main features of the code and show three practical examples of multi-spacecraft tomographic inversions using STEREO/EUVI and STEREO/COR1 data. Static and smoothly varying temporal evolution models are presented.