Next generation radio interferometric telescopes are entering an era of big data with extremely large data sets. While these telescopes can observe the sky in higher sensitivity and resolution than before, computational challenges in image reconstruction need to be overcome to realize the potential of forthcoming telescopes. New methods in sparse image reconstruction and convex optimization techniques (cf. compressive sensing) have shown to produce higher fidelity reconstructions of simulations and real observations than traditional methods. This article presents distributed and parallel algorithms and implementations to perform sparse image reconstruction, with significant practical considerations that are important for implementing these algorithms for Big Data. We benchmark the algorithms presented, showing that they are considerably faster than their serial equivalents. We then pre-sample gridding kernels to scale the distributed algorithms to larger data sizes, showing application times for 1 Gb to 2.4 Tb data sets over 25 to 100 nodes for up to 50 billion visibilities, and find that the run-times for the distributed algorithms range from 100 milliseconds to 3 minutes per iteration. This work presents an important step in working towards computationally scalable and efficient algorithms and implementations that are needed to image observations of both extended and compact sources from next generation radio interferometers such as the SKA. The algorithms are implemented in the latest versions of the SOPT (https://github.com/astro-informatics/sopt) and PURIFY (https://github.com/astro-informatics/purify) software packages {(Versions 3.1.0)}, which have been released alongside of this article.

Comparison of appropriate theoretical derived line ratios with observational data can yield estimates of a plasma's physical parameters, such as electron density or temperature. The usual practice in the calculation of the line ratio is the assumption of excitation by electrons/protons followed by radiative decay. Furthermore, it is normal to use the so-called coronal approximation, i.e. one only considers ionization and recombination to and from the ground state. A more accurate treatment is to include the ionization/recombination to and from meta-stable levels. Here, we apply this to two lines from adjacent ionization stages; Mg IX 368A and Mg X 625A, which has been shown to be a very useful temperature diagnostic. At densities typical of coronal hole conditions, the difference between the electron temperature derived assuming the zero density limit compared with the electron density dependent ionization/recombination is small. This however is not the case for flares where the electron density is orders of magnitude larger. The derived temperature for the coronal hole at solar maximum is around 1.04 MK compared to just below 0.82 MK at solar minimum.

Context. We diagnose the properties of the plume and interplume regions in a polar coronal hole and the role of waves in the acceleration of the solar wind. Aims. We attempt to detect whether Alfven waves are present in the polar coronal holes through variations in EUV line widths. Methods. Using spectral observations performed over a polar coronal hole region with the EIS spectrometer on Hinode, we study the variation in the line width and electron density as a function of height. We use the density sensitive line pairs of Fe xii 186.88 A & 195.119 A and Fe xiii 203.82 A & 202.04 A . Results. For the polar region, the line width data show that the nonthermal line-of-sight velocity increases from 26 km/s at 1000 above the limb to 42 km/s some 15000 (i.e. 110,000 km) above the limb. The electron density shows a decrease from 3:3 10^9 cm^-3 to 1:9 10^8 cm^-3 over the same distance. Conclusions. These results imply that the nonthermal velocity is inversely proportional to the quadratic root of the electron density, in excellent agreement with what is predicted for undamped radially propagating linear Alfven waves. Our data provide signatures of Alfven waves in the polar coronal hole regions, which could be important for the acceleration of the solar wind.