Open Access
Issue
A&A
Volume 683, March 2024
Article Number A127
Number of page(s) 10
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202348443
Published online 13 March 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Numerical simulations where the heating is concentrated chiefly towards the footpoints of coronal loops can lead to a phenomenon known as coronal rain (e.g. Antiochos & Klimchuk 1991; Müller et al. 2003, 2004; Mikić et al. 2013). Footpoint heating leads the chromospheric plasma to “evaporate” and rise along the loop, ultimately reaching coronal temperatures and therefore filling it with coronal-temperature material. As the heating is predominantly concentrated at the footpoints, the material near the coronal loop apex is unable to maintain its thermodynamic state due to the radiative losses. These increase as the density of the plasma increases and its temperature diminishes, a dynamics that can sometimes lead to a runaway cooling process if the temperature falls below a certain threshold, usually at ∼106 K. In this situation, the plasma will continue to cool down to transition region and chromospheric temperatures, while at the same time increasing its density. Such dense condensation is out of the mechanical equilibrium and flows down along the loop in the form of coronal rain.

The coronal rain material consists of a denser and colder core of chromospheric temperature surrounded by a warm envelope, whose temperature transitions from that of the clump’s core to coronal values (Antolin et al. 2022). This temperature inhomogeneity allows for the plasma to be visible in a myriad of spectral lines that sample coronal and chromospheric temperatures. Coronal rain appears in these spectral lines as clump-like concentrations of material falling along the magnetic field of coronal loops with velocities of up to 100 − 150 km s−1 (Antolin et al. 2010, 2015; Antolin & Rouppe van der Voort 2012; Froment et al. 2020; Kriginsky et al. 2021). Spectroscopic observations of coronal rain in chromospheric lines, such as the Ca II k and Hα lines, with high-resolution instruments such as the Solar Optical Telescope on board of Hinode (Tsuneta et al. 2008) and the CRISP imaging spectropolarimeter (Scharmer et al. 2008) at the Swedish 1-m solar telescope (Scharmer et al. 2003) have allowed for the determination of multiple properties of coronal rain clumps, such as their temperature and downfall velocities. The use of off-limb observations in chromospheric lines forming in the optically thin regime allowed for several estimates of the average temperature of coronal rain clumps (Antolin & Rouppe van der Voort 2012; Antolin et al. 2015; Froment et al. 2020). These studies have found temperatures of coronal rain plasma in the range between 5000 K up to values beyond 30 000 K.

The focus of this work is to use another diagnostic technique that still has not been exploited in coronal rain observations: spectroscopic inversions. When the coronal rain material reaches chromospheric temperatures and can be detected on the disc with enough signal in a spectral line whose formation properties are strongly correlated with those of the plasma constituting the coronal rain clumps, the use of inversion codes might serve as an alternative diagnostic tool. The diagnostic of the properties of the coronal rain plasma is most easily done in off-limb observations. However, if the line under study has enough opacity in coronal rain material for it to be formed under optically thick conditions, it could also potentially be used in on-disc observations. The Mg II h & k lines, which form in the upper chromosphere and possess stronger opacity than the more extensively used Ca II h & k lines are possibly among the most promising candidates.

The launch of the Interface Region Imaging Spectrograph (IRIS, De Pontieu et al. 2014) mission has made systematic observations of the chromosphere and transition region broadly available, including spectra from lines that are inaccessible from ground-based solar observatories. The near-ultraviolet (near-UV) spectral range of IRIS includes spectroscopic data from both the Mg II h & k lines and the secondary triplet of Mg II located in their vicinity. With IRIS completing a decade of observations, the online database now includes a multitude of coronal rain observations, most easily seen off-limb.

Before and after the launch of IRIS, a series of publications have been devoted to the study of the different spectral regions and lines accessible to its spectrograph. In particular, the works of Leenaarts et al. (2013a,b) focused on the formation properties of the Mg II h & k lines and the connection of such properties to the physical state of the plasma. This was done with the use of synthesis of the spectral profiles computed numerically using as the input atmosphere a simulation snapshot from the Bifrost code (Gudiksen et al. 2011). The aforementioned works found a strong relationship between the temperature at the formation height of the emission peaks of the Mg II h & k lines and their intensity; therefore, their inclusion in inversions could lead to a better constraint of the temperature at their formation height. This temperature sensitivity has been exploited in a number of works (de la Cruz Rodríguez et al. 2016; da Silva Santos et al. 2018; Vissers et al. 2019). In addition to constraining the temperature of the plasma, Leenaarts et al. (2013b) showed that, at least for an atmosphere arising from a simulation snapshot, there are some constraints about the vertical velocity of the plasma that can be obtained by combining different spectral features of the Mg II h & k lines.

The goal of the present work is to exploit the diagnostic capabilities of the Mg II h & k lines in order to infer properties about the coronal rain clumps forming in active regions observed with the IRIS spectrograph. This work is organised as follows. Section 2 details the different observations used and their calibration, with the inversion scheme and data analysis of spectral features detailed in Sect. 3. The results are discussed in Sect. 4 and our conclusions are presented in Sect. 5.

2. Observations

Coronal rain is nearly always present in active regions (Şahin et al. 2023). However, its detection on the disc is not straightforward. The presence of a myriad of different phenomena in active regions hinders the possibility of detecting coronal rain events automatically. In this work, a comprehensive manual inspection of the IRIS public data archive1 was undertaken in order to detect on-disc coronal rain events that crossed the slit position at any given point in their evolution.

The study of the IRIS observations in coordination with data in the 30.4 nm passband from the Atmospheric Imaging Assembly (AIA, Boerner et al. 2012; Lemen et al. 2012), on board the Solar Dynamics Observatory (SDO, Pesnell et al. 2012), demonstrated that the easiest way to detect coronal rain clumps in active regions is when the slit is positioned in the neighbourhood of the bright plage areas and the coronal loop hosting the coronal rain event is crossing the field of view of this neighbourhood. This situation allows for the easiest detection of the presence of coronal rain in the slitjaw images of IRIS as well as in the spectra of the Mg II h & k lines. The presence of coronal rain in the field of view corresponds to a diminution of over 40% of the intensity compared to the background. This diminution is much smaller or nonexistent in the case of coronal rain clumps crossing a quiet Sun area. Therefore, it becomes nearly impossible to unambiguously detect the motion and position of coronal rain in the slitjaw in such a scenario.

The data used in this work consisted of nine observations undertaken with the IRIS spacecraft between 2014 and 2015. All observations targeted active regions. In each one, coronal rain is seen crossing the position of the slit, with the slit either in a sit-and-stare mode or performing raster scans. A detailed description of each observation is provided in Table 1. The calibration of the data was done in line with the reduction methods detailed in Wülser et al. (2018). Radiometric calibration of the raster data was carried out using version four of the calibration files employed by the iris_get_response routine in SolarSoft (SSW; Freeland & Handy 1998).

Table 1.

Description of the datasets used in this work.

Referring to the labels in Table 1, observations 1–6 correspond to a tracking campaign of active region 12007. The remainder of the observations correspond to single observations of different active regions where coronal rain was visibly crossing the slit. In observations 1–3 and observation 7, the scan mode was sit-and-stare; whereas for observations 4–6 a sparse (1″ steps) raster scan was undertaken and for observations 8–9 dense (0.33″ steps) and coarse (2″ steps) scans were performed, respectively. Figure 1 provides a context of the observations for each of the different active regions listed in Table 1. For brevity and clarity, only one observation corresponding to active region number 12007 is shown, as these observations are largely similar to each other, differing only with respect to the particular coronal loop that hosts the coronal rain event during the observation. Coronal rain material in the shape of long strands can be clearly seen in the AIA 30.4 nm images (first and second rows of Fig. 1). Detecting this material in the 279.6 nm passband is much harder (third row of Fig. 1).

thumbnail Fig. 1.

Observations used in this study. Each column corresponds to one of four of the observations listed in Table 1, namely: observations 1, 7, 8, and 9. The first row displays the field of view of the IRIS slitjaw imager as a coloured rectangle within a co-aligned image in the 30.4 nm passband obtainedw with SDO/AIA. The second row panels display a zoomed view into the field of view inside the rectangles of the first row panels. The corresponding slitjaw images in the 279.6 nm passband are displayed in the third row panels. The thick vertical lines in each panel of the second and third rows outline the spatial portion of the slit used in the fourth row panels. Arrows indicate the location of coronal rain as it crosses the slit. The fourth row panels display the intensity in the core of the Mg II h line as a function of time. The vertical coloured lines in each panel of the fourth row provides the time instant used for the images in the top three corresponding panels. Finally, the last row shows a mask of the fourth row panels for each observation, with the pixels used for the inversions in each case shown in black. The areas coloured white were not used in the inversions. Two accompanying movies display the coronal rain formation and its crossing of the slit position for observations 7 and 8. The two movies have been trimmed to show the time when coronal rain is present.

3. Methods

3.1. Inversions

The inversion of the IRIS data was made with the MPI-parallel Stockholm inversion code (STiC; de la Cruz Rodríguez et al. 2016, 2019), which is a regularised Levenberg-Marquardt code. It uses an optimised version of the RH code (Uitenbroek 2001) to solve the atom population densities assuming statistical equililibrium and plane-parallel geometry for multi-atom inversions of multiple spectral lines assuming non-local thermodynamic equilibrium (non-LTE). It also allows for the inclusion of partial redistribution effects of scattered photons (Leenaarts et al. 2012). The radiative transport equation is solved using cubic Bezier solvers (de la Cruz Rodríguez & Piskunov 2013). The inversion engine of STiC includes an equation of state extracted from the Spectroscopy Made Easy code (SME; Piskunov & Valenti 2017). The electron densities were derived with the assumption of non-LTE hydrogen ionisation by solving the statistical equilibrium equations imposing charge conservation (Leenaarts et al. 2007). As a 1.5D plane parallel code, STiC cannot account for some 3D radiative transfer effects that could be important for the inversion of the Mg II h & k lines (Leenaarts et al. 2013a,b; Pereira et al. 2015).

A model atom of Mg II consisting of 11 atomic levels including the ground level of Mg III was used. In addition, a photospheric Ni I line at 281.4 nm was used in order to have some constraint about the state of the photosphere, which was treated under the assumption of LTE. In order to save computation time given the large number of inversions performed, only the region around the cores of the Mg II h & k lines was used for the inversions, using the same spectral regions as those used by Kriginsky et al. (2023, in particular, their Fig. 4).

A regular depth grid of 36 points covering the range log τ500 = [ − 7, 0] (where τ500 is the optical depth of the 500 nm radiation) was used in the inversions. The initial temperature depth profile was interpolated from a FAL-C (Avrett 1985; Fontenla et al. 1993) atmosphere and used as an initial guess for the inversions. Given the different view points in the observations, the initial guess for the velocity along the line of sight (LoS) was modified in order to provide a suitable guess for each observation. The parameters that were actively inverted were the temperature, T, the LoS component of the velocity, vLOS, and the microturbulent velocity, vturb. The node distribution was not the same for every observation. However, the same philosophy was followed, namely an increasing number of nodes for each cycle. Most of the inversions were performed in no more than three inversion cycles, except for observation 9, which required up to five inversion cycles in order to reach a reasonable fit to the observations. This is probably due to to the point of view of the observation being much closer to the limb than the other observations.

Given the large amount of pixels in the data from all the observations, only the pixels containing coronal rain were inverted. They were selected taking advantage of the large intensity difference between pixels containing coronal rain and background plage pixels through an intensity masking procedure (see bottom row of Fig. 1). Additionally, in order to reduce the number of wavelength points used in the inversions, a region centred at the rest wavelengths of the Mg II h & k lines with a total width of 0.2 nm was used. For the two Mg II triplet lines between the h & k lines, an additional region with a width of 0.1 nm was selected, centred at the 279.9 nm wavelength.

3.2. Velocity diagnostics

Macroscopic velocity fields can be inferred from the Doppler displacement of a spectral line. This procedure has been extensively used before with spectra obtained in chromospheric lines, such as the Ca II H line (Antolin et al. 2010) or the Hα line, to determine the typical velocities of coronal rain clumps (Antolin et al. 2010, 2015; Antolin & Rouppe van der Voort 2012; Froment et al. 2020; Kriginsky et al. 2021). The results of these works are in great agreement with each other, with clump total velocities of up to 200 km s−1 obtained.

The spectra in the Mg II h & k lines also allows for a similar study. In fact, there are several additional diagnostics that are possible by using both lines individually and by combining them. Using spectra computed with the Multi3D radiative transfer code (Leenaarts & Carlsson 2009) with data from a simulation snapshot performed with the Bifrost code (Gudiksen et al. 2011), Leenaarts et al. (2013b) performed a series of tests using several spectral features of these lines in order to assess what information can be extracted from them. The diagnostics used in this work are described below (see Leenaarts et al. 2013b for the complete discussion and for the definition of the nomenclature used below):

  • With respect to the h3 and k3 Doppler shifts: The Doppler displacements of the line core position of the Mg II h & k lines are found to be strongly correlated with the velocity along the LoS at the height of formation of the line cores. This value is hereafter referred to as either v3 when grouping together results for both h & k lines or as vh3 and vk3 when specifying the line becomes relevant.

  • Difference between h3 and k3: The core of the k line was found to be formed a few hundred km below the transition region, with the core of the h line forming at slightly lower heights. For small differences (smaller than 2 km s−1) between the velocities at the formation heights of h3 and k3, there is a strong correlation between such differences and the difference between the velocities measured from Doppler shifts at h3 and k3.

  • Average peak Doppler shift: A correlation between the average Doppler shift of the red and blue emission peaks of both lines was found to be tightly related to the average vertical velocity at the height of formation of the peaks. This velocity value is referred to here as v2 as a general term for either the h or k lines and explicitly as vh2 or vk2 when the value for each line is referenced.

Access to these diagnostics makes it possible to estimate the vertical velocities and their differences inside the coronal clump plasma.

3.3. Identification of spectral features

In order to make use of the diagnostics given above, the different spectral features (emission peaks and central reversal) had to be determined for each pixel. This was done by means of fitting the region around the line cores of the Mg II h & k lines with a sum of splines. There are mainly four types of Mg II h & k line core profiles, namely: the typically double-peaked profile with a central reversal (top panel of Fig. 2), a single-peaked profile with a secondary lobe (bottom panel of Fig. 2), a single-peaked profile without a secondary lobe, and profiles with more than two emission peaks. In order to minimise the uncertainties when determining the position of the h2, h3, k2, and k3 spectral features, only the first two types of profiles were used, and the rest were discarded. This represented a loss of about 20% of all coronal rain pixels.

thumbnail Fig. 2.

Spectral features, featuring a typical Mg II h line profile with two emission peaks and a line core minimum (shown in the top panel) and a single-peaked profile (shown in the bottom panel). The positions of the extremes of the fitted profile are marked with a red star, with the vertical red lines denoting the retrieved h3 and h2 features. For the single-peaked profile, the position of the only extreme is also marked with a red star, while the extremes of the first and second derivatives are marked with green and blue stars, respectively. The vertical lines also identify the retrieved h3 and h2 features. The black dashed curves in both panels represent the fitted curve to the spectral profile, with the black dots representing observed spectral points.

For the regular, double-peaked profiles, the detection algorithm is straightforward: the extremes of the first derivative of the fitted curve are determined, with the two maxima and the minimum between them extracted. Profiles with one maximum and a secondary lobe required the use of the second and third derivatives. The position of either the red or blue emission peak was determined as the zero crossing position of the first derivative. These profiles are characterised by four points where the second derivative is zero. If the profile maximum is in the red (blue) side, then the second point with null second derivative to its right (left) was taken as the position of the remaining emission peak. The line centre was taken to be the position of the zero crossing of the third derivative in between the two emission peaks (see Fig. 2).

4. Results

4.1. Response functions

The optical depth of formation of each wavelength was approximated by the depth of maximum response to temperature variations as computed numerically with STiC. The response function to temperature was calculated for every inverted pixel. The distributions of the optical depth of peak response to temperature of the h2r and k2r features are shown in the third panel of Fig. 3 for observations 7 and 9. For both observations there are clearly two groups of pixels; one for which the formation depth of the h2r and k2r features is in the upper atmosphere and another for which it is in the middle chromosphere. This was an expected result, as there is the possibility that the emission peaks of the Mg II h & k lines form deeper in the atmosphere (Leenaarts et al. 2013b). Pixels pertaining to the second aforementioned group were discarded as they formed deeper in the atmosphere, arguably below the height of the coronal rain material, with the division between both groups set at log τ500 = −4.5, where the distribution of formation depths of observation 9 has a maximum (see bottom right panel of Fig. 3). The coronal rain material in observation 7 is much more abundant than the small rain event of observation 9; therefore, the coronal rain pixels in observation 7 mostly belong to the first group mentioned above, with the opposite being true for observation 9. With the very small amount of pixels in observation 9 forming at upper layers of the atmosphere, this observation was not used. In observation 7, the discarded pixels belong to the edges of the coronal rain material (yellow-coloured pixels in the upper middle panel of Fig. 3), where (arguably) the coronal rain material exhibits a lower opacity; therefore, the formation height of the h2r and k2r features is in the middle chromosphere, not in the coronal rain itself.

thumbnail Fig. 3.

Response functions. Left columns display the line core intensity in the Mg II k line as a function of time, with the top one referring to observation 7, and the bottom to observation 9. The middle panels show the optical depth of maximum response to temperature in these observations for the h2r feature. The right columns show the distributions of maximum response to temperature for each observation, in each case, for both h2r and k2r features.

4.2. Velocities

The methods described in Sects. 3.2 and 3.3 were applied to all the observations, with the results displayed in Figs. 4 and 5. Positive values of the velocity in the first two panels correspond to movements away from the observer. In the first two panels of each row, the velocities at the formation height of the emission peaks (first panels) and of the line centre (second panels) are mostly positive, indicating that the material is falling towards the chromosphere. The derived velocity distributions for observations 4–6 do have a sizeable negative proportion of values. As they are all close together in time and location (see Table 1), this could be caused by the orientation of the coronal loops respect to the observer. The measured velocities with respect to the observer are rather small in comparison to the values reported by studies in the Ca II H line and in the Hα line (Antolin et al. 2010, 2015; Antolin & Rouppe van der Voort 2012; Froment et al. 2020; Kriginsky et al. 2021). The reason for this discrepancy is that in all observations reported here the motion of the coronal rain clump is constrained mostly to the plane of the sky in a perpendicular direction to the observations.

thumbnail Fig. 4.

Velocity diagnostics. Each row represents a different observation, from observation 1 to 4 in ascending order from the top to the bottom. The first two columns represent the distribution of the LoS velocity as measured in the position of the emission peaks (first column) and in the line core (second column). The third column shows the distribution of the relative values of the line core velocities between the h & k lines, and the fourth column shows the relative velocity as measured in the line core and the emission peaks.

thumbnail Fig. 5.

Same as Fig. 4, for observations 5 to 8.

The larger opacity of the Mg II k line leads to a slightly higher formation height than the h line (Leenaarts et al. 2013b). Additionally, the emission peaks of each line form deeper than the line cores. Assuming that this is also true in the case of coronal rain, it is (in principle) possible to estimate velocity differences between different parts of coronal rain clumps. Those are shown in the third and fourth columns of Figs. 4 and 5. The relative velocities between different heights are measured with the expression:

V r = | v h | | v l | × sign ( v h ) × sign ( v l ) , $$ \begin{aligned} V_{\rm r} = |v_{\rm h}|-|v_{\rm l}|\times \mathrm{sign} (v_{\rm h})\times \mathrm{sign} (v_{\rm l}), \end{aligned} $$(1)

where vh is the velocity assumed to be measured higher in the atmosphere and vl is the velocity assumed to be measured at a lower height. With Eq. (1), a positive relative velocity means that the higher layer is moving faster towards or away from the observer than the velocity at the lower layer and a negative velocity means the opposite is true.

The relative velocity between the core of the h & k lines is relatively small, as the cores of the lines tend to form very close to each other (Leenaarts et al. 2013b). However, the distributions seem to all be slightly shifted towards positive values, indicating that the plasma at the outer layers of the coronal rain clump is moving faster than the plasma below it. The right-most panels of Figs. 4 and 5 show the relative velocity between the core of each line and its emission peaks, again estimated using Eq. (1). The relative velocities are slightly larger than the ones measured between the two line cores, with the shift towards positive values is even larger than that between the line cores, indicating that the layer where the line core is formed is moving with larger speeds than the layer where the emission peaks are formed. Such result should be confirmed and justified in the future thanks to further observations and coronal rain simulations.

4.3. Temperatures and densities

The temperature inversions were performed with the methodology described in Sect. 3. For each pixel in the different observations inversions (for which both the h2 and k2 are formed in the upper atmosphere), Fig. 6 shows the temperature distribution at the height of formation of these features for each of the observations. Observation 9 was discarded because for its pixels, both the h2 and k2 features are mostly formed in the middle chromosphere.

thumbnail Fig. 6.

Temperature diagnostics. From left to right and top to bottom: histogram of the average temperature at the depth of maximum response to temperature of the h2 and k2 features for observations 1 through 8.

The inverted temperatures are consistent between the different observations and all remain within the narrow temperature range between 5200 K up to 7500 K. This could indicate that in the different condensation events, the rain becomes visible in the Mg II h & k lines in a very specific moment during its cooling phase and that different condensation events pertaining to independent coronal loops reach similar thermal conditions. The inverted temperature values are consistent with the temperature at the denser core of coronal rain found in simulations (Antolin et al. 2022), which could mean that that is where the emission peaks of the Mg II h & k lines are being formed. The inverted temperature of the k2 feature seems to be slightly larger than that inferred at the h2 feature (with the exception of observation 6), indicating that the h2 forms closer to the core of the coronal rain clump. Some of the temperature distributions show apparent peaks centred at multiple temperature values. For example, the temperature distribution at the formation height of the k2 feature for observation 3 has a peak around 6500 K and another apparent peak at 7100 K. The temperature distributions for observation 7 have multiple peaks for both the h2 and k2 features. This multiplicity in distribution peaks can be caused by the variation of the geometrical height of formation of the spectral features even for a single observation, where they might have their formation at locations of different temperatures. Additionally, observation 7 stems from a flare event which led to the formation of at least three condensation events that crossed the position of the slit. Small differences in the properties of the plasma constituting these events could also serve as another cause for the variety of peaks in the corresponding temperature distribution.

Figure 7 shows the derived densities, which are also similar to those at the denser cores of the coronal rain plasma found in simulations (Antolin et al. 2022) and with values obtained from estimates of the density based on the correlation between the emission measure and intensity in the Hα line (Froment et al. 2020). Similarly to the temperature values, there seems to be some shift between the densities inferred in the h2 and k2 features, with denser values at the formation height of the h2 feature being inferred.

thumbnail Fig. 7.

Density diagnostics. Same as Fig. 6, but for densities.

5. Conclusions

The strongly stratified heating that coronal loops are subject to can give rise to the ubiquitously-observed phenomenon of coronal rain. The multi-thermal nature of coronal rain allows for it to be studied with lines formed in the chromosphere, the transition region and the corona. This study has made use of successful observations of coronal rain formation in the Mg II h & k lines as observed with the spectrograph mounted in the IRIS satellite.

As opposed to the off-limb observations, which rely in the inference of the temperature assigning a single temperature value to the emitting plasma, inversions can in theory provide a depth-dependent temperature stratification. The non-LTE code STiC was used for the inversion of coronal rain spectra in nine different active region observations obtained on the solar disc. The sensitivity of the Mg II h & k lines allowed for the inference of the temperature of coronal rain clumps using this alternative method. The results largely agree with inferences made using thermal width estimates from the emission line profiles. The results also agree with those obtained in numerical simulations of coronal loops undergoing catastrophic cooling. The results are specially aligned with the cool, dense cores of such studies, which might suggest that this is where the Mg II h & k lines are formed. This view also agrees with the density estimates obtained from the inversions, which agree with the values obtained at the core of coronal rain clumps from simulations. From the distributions of density and temperature obtained from the spectral inversions, it seems that there is a noticeable change in the temperature and density of the plasma between the formation height of the h2 and k2 features, with the h2 appearing to form in a slightly colder and denser region.

The formation properties of the different spectral features of the Mg II h & k lines have also allowed for a study of the velocity in the LoS direction. Several relationships obtained from the study of Leenaarts et al. (2013b) based in the comparison of synthetic line profiles obtained from 3D radiation transfer calculations in a snapshot of a 3D simulation were used in this study. If these relationships between spectral features and vertical velocities still hold for the case of the Sun in the presence of coronal rain clumps, the results of this work are suggestive of an interesting scenario at work. The plasma formed in the outer part of the coronal rain clump, where the line cores of the Mg II h & k lines are formed, falls towards the chromosphere with a greater speed than that of the plasma constituting the denser part of the clumps, where the emission peaks of the lines form. This situation seems to take place in all the observations, suggesting that it is a common phenomenon in the dynamics of coronal rain. It might not be present in all stages during the fall of the material and it might only occur in the initial part of the descent along magnetic field lines. The existence, cause, and possible significance of this velocity gradient should be analysed in future studies that would include numerical simulations and additional observations of on-disc coronal rain. Some numerical simulations (e.g. Antolin et al. 2022; Martínez-Gómez et al. 2020) show that the falling velocity of the outer plasma layers is smaller than that of the inner parts. The interpretation of this fact is that less dense material falls more slowly than denser one (see Oliver et al. 2014). The results presented here, however, show that the material at the outer part of coronal rain falls faster than that closer to its core, where we can safely assume that the highest density is found. This discrepancy between numerical simulations and the present results warrants further attention in subsequent studies.

Movies

Movie 1 associated with Fig. 1 (2014-05-13) Access here

Movie 2 associated with Fig. 1 (2015-07-10) Access here


Acknowledgments

This publication is part of the R+D+i project PID2020-112791GB-I00, financed by MCIN/AEI/10.13039/501100011033. M.K. acknowledges the support from the Vicepresidéncia i Conselleria d’Innovació, Recerca i Turisme del Govern de les Illes Balears and the Fons Social Europeu 2014-2020 de les Illes Balears. IRIS is a NASA small explorer mission developed and operated by LMSAL with mission operations executed at NASA Ames Research center and major contributions to downlink communications funded by ESA and the Norwegian Space Centre. This research was supported by the International Space Science Institute (ISSI) in Bern, through ISSI International Team project #545 (“Observe Local Think Global: What Solar Observations can Teach us about Multiphase Plasmas across Physical Scales”). The authors thank the anonymous referee for some comments and suggestions that helped improve the manuscript.

References

  1. Antiochos, S. K., & Klimchuk, J. A. 1991, ApJ, 378, 372 [Google Scholar]
  2. Antolin, P., & Rouppe van der Voort, L. 2012, ApJ, 745, 152 [Google Scholar]
  3. Antolin, P., Shibata, K., & Vissers, G. 2010, ApJ, 716, 154 [NASA ADS] [CrossRef] [Google Scholar]
  4. Antolin, P., Vissers, G., Pereira, T. M. D., Rouppe van der Voort, L., & Scullion, E. 2015, ApJ, 806, 81 [Google Scholar]
  5. Antolin, P., Martínez-Sykora, J., & Şahin, S. 2022, ApJ, 926, L29 [NASA ADS] [CrossRef] [Google Scholar]
  6. Avrett, E. H. 1985, in Chromospheric Diagnostics and Modelling, ed. B. W. Lites, 67 [Google Scholar]
  7. Boerner, P., Edwards, C., Lemen, J., et al. 2012, Sol. Phys., 275, 41 [Google Scholar]
  8. da Silva Santos, J. M., de la Cruz Rodríguez, J., & Leenaarts, J. 2018, A&A, 620, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. de la Cruz Rodríguez, J., & Piskunov, N. 2013, ApJ, 764, 33 [Google Scholar]
  10. de la Cruz Rodríguez, J., Leenaarts, J., & Asensio Ramos, A. 2016, ApJ, 830, L30 [Google Scholar]
  11. de la Cruz Rodríguez, J., Leenaarts, J., Danilovic, S., & Uitenbroek, H. 2019, A&A, 623, A74 [Google Scholar]
  12. De Pontieu, B., Title, A. M., Lemen, J. R., et al. 2014, Sol. Phys., 289, 2733 [Google Scholar]
  13. Fontenla, J. M., Avrett, E. H., & Loeser, R. 1993, ApJ, 406, 319 [Google Scholar]
  14. Freeland, S. L., & Handy, B. N. 1998, Sol. Phys., 182, 497 [Google Scholar]
  15. Froment, C., Antolin, P., Henriques, V. M. J., Kohutova, P., & Rouppe van der Voort, L. H. M. 2020, A&A, 633, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Gudiksen, B. V., Carlsson, M., Hansteen, V. H., et al. 2011, A&A, 531, A154 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Kriginsky, M., Oliver, R., Antolin, P., Kuridze, D., & Freij, N. 2021, A&A, 650, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Kriginsky, M., Oliver, R., & Kuridze, D. 2023, A&A, 672, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Leenaarts, J., & Carlsson, M. 2009, in The Second Hinode Science Meeting: Beyond Discovery-Toward Understanding, eds. B. Lites, M. Cheung, T. Magara, J. Mariska, & K. Reeves, ASP Conf. Ser., 415, 87 [Google Scholar]
  20. Leenaarts, J., Carlsson, M., Hansteen, V., & Rutten, R. J. 2007, A&A, 473, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Leenaarts, J., Pereira, T., & Uitenbroek, H. 2012, A&A, 543, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Leenaarts, J., Pereira, T. M. D., Carlsson, M., Uitenbroek, H., & De Pontieu, B. 2013a, ApJ, 772, 89 [NASA ADS] [CrossRef] [Google Scholar]
  23. Leenaarts, J., Pereira, T. M. D., Carlsson, M., Uitenbroek, H., & De Pontieu, B. 2013b, ApJ, 772, 90 [NASA ADS] [CrossRef] [Google Scholar]
  24. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [Google Scholar]
  25. Martínez-Gómez, D., Oliver, R., Khomenko, E., & Collados, M. 2020, A&A, 634, A36 [EDP Sciences] [Google Scholar]
  26. Mikić, Z., Lionello, R., Mok, Y., Linker, J. A., & Winebarger, A. R. 2013, ApJ, 773, 94 [Google Scholar]
  27. Müller, D. A. N., Hansteen, V. H., & Peter, H. 2003, A&A, 411, 605 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Müller, D. A. N., Peter, H., & Hansteen, V. H. 2004, A&A, 424, 289 [CrossRef] [EDP Sciences] [Google Scholar]
  29. Oliver, R., Soler, R., Terradas, J., Zaqarashvili, T. V., & Khodachenko, M. L. 2014, ApJ, 784, 21 [Google Scholar]
  30. Pereira, T. M. D., Carlsson, M., De Pontieu, B., & Hansteen, V. 2015, ApJ, 806, 14 [NASA ADS] [CrossRef] [Google Scholar]
  31. Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3 [Google Scholar]
  32. Piskunov, N., & Valenti, J. A. 2017, A&A, 597, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Şahin, S., Antolin, P., Froment, C., & Schad, T. A. 2023, ApJ, 950, 171 [CrossRef] [Google Scholar]
  34. Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., & Petterson, B. 2003, Proc. SPIE, 4853, 341 [Google Scholar]
  35. Scharmer, G. B., Narayan, G., Hillberg, T., et al. 2008, ApJ, 689, L69 [Google Scholar]
  36. Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167 [Google Scholar]
  37. Uitenbroek, H. 2001, ApJ, 557, 389 [Google Scholar]
  38. Vissers, G. J. M., de la Cruz Rodríguez, J., Libbrecht, T., et al. 2019, A&A, 627, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Wülser, J. P., Jaeggli, S., De Pontieu, B., et al. 2018, Sol. Phys., 293, 149 [Google Scholar]

All Tables

Table 1.

Description of the datasets used in this work.

All Figures

thumbnail Fig. 1.

Observations used in this study. Each column corresponds to one of four of the observations listed in Table 1, namely: observations 1, 7, 8, and 9. The first row displays the field of view of the IRIS slitjaw imager as a coloured rectangle within a co-aligned image in the 30.4 nm passband obtainedw with SDO/AIA. The second row panels display a zoomed view into the field of view inside the rectangles of the first row panels. The corresponding slitjaw images in the 279.6 nm passband are displayed in the third row panels. The thick vertical lines in each panel of the second and third rows outline the spatial portion of the slit used in the fourth row panels. Arrows indicate the location of coronal rain as it crosses the slit. The fourth row panels display the intensity in the core of the Mg II h line as a function of time. The vertical coloured lines in each panel of the fourth row provides the time instant used for the images in the top three corresponding panels. Finally, the last row shows a mask of the fourth row panels for each observation, with the pixels used for the inversions in each case shown in black. The areas coloured white were not used in the inversions. Two accompanying movies display the coronal rain formation and its crossing of the slit position for observations 7 and 8. The two movies have been trimmed to show the time when coronal rain is present.

In the text
thumbnail Fig. 2.

Spectral features, featuring a typical Mg II h line profile with two emission peaks and a line core minimum (shown in the top panel) and a single-peaked profile (shown in the bottom panel). The positions of the extremes of the fitted profile are marked with a red star, with the vertical red lines denoting the retrieved h3 and h2 features. For the single-peaked profile, the position of the only extreme is also marked with a red star, while the extremes of the first and second derivatives are marked with green and blue stars, respectively. The vertical lines also identify the retrieved h3 and h2 features. The black dashed curves in both panels represent the fitted curve to the spectral profile, with the black dots representing observed spectral points.

In the text
thumbnail Fig. 3.

Response functions. Left columns display the line core intensity in the Mg II k line as a function of time, with the top one referring to observation 7, and the bottom to observation 9. The middle panels show the optical depth of maximum response to temperature in these observations for the h2r feature. The right columns show the distributions of maximum response to temperature for each observation, in each case, for both h2r and k2r features.

In the text
thumbnail Fig. 4.

Velocity diagnostics. Each row represents a different observation, from observation 1 to 4 in ascending order from the top to the bottom. The first two columns represent the distribution of the LoS velocity as measured in the position of the emission peaks (first column) and in the line core (second column). The third column shows the distribution of the relative values of the line core velocities between the h & k lines, and the fourth column shows the relative velocity as measured in the line core and the emission peaks.

In the text
thumbnail Fig. 5.

Same as Fig. 4, for observations 5 to 8.

In the text
thumbnail Fig. 6.

Temperature diagnostics. From left to right and top to bottom: histogram of the average temperature at the depth of maximum response to temperature of the h2 and k2 features for observations 1 through 8.

In the text
thumbnail Fig. 7.

Density diagnostics. Same as Fig. 6, but for densities.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.