For Science / Math Lovers: Measurement of a Solar Filament Temperature.

I LOVE finding out about different ways to appreciate the Sun and light in general. Use this forum to post your info or questions about various outside the mainstream ways to appreciate our life giving star!
Post Reply
User avatar
p_zetner
Way More Fun to Share It!!
Way More Fun to Share It!!
Posts: 1644
Joined: Sun Mar 18, 2012 4:59 pm
Location: Winnipeg, Manitoba, Canada
Has thanked: 166 times
Been thanked: 719 times

For Science / Math Lovers: Measurement of a Solar Filament Temperature.

Post by p_zetner »

Hi Everyone.

In my last imaging session with the spectroheliograph, I observed a strong filament in both the H alpha and Ca K images. Here are two images with the filament at, roughly, the solar equator towards the E limb chosen to be the one for study.
Montage_Finals_Ha-CaK lvls siz.png
Montage_Finals_Ha-CaK lvls siz.png (501.82 KiB) Viewed 4834 times

An interesting aspect of the H alpha filament was its persistence over a large detuning from line centre. This can be appreciated in the spectral sweep animation presented here:

viewtopic.php?f=8&t=24946

I thought the broad absorption (in wavelength span) from the H alpha filament might be due to a high temperature and corresponding large thermal (Doppler) broadening. Researching the subject a little, I came across the paper by Park et al.: “Temperature of Solar Prominences Obtained with the Fast Imaging Solar Spectrograph on the 1.6m New Solar Telescope at the Big Bear Solar Observatory” pgs105-116 in “Initial Results from the Fast Imaging Solar Spectrograph (FISS)” (2015) ed. Jongchul Chae, Springer International Publishing. These authors measured prominence temperatures and I decided to apply a similar analysis to the filament.

https://web.njit.edu/~sirenko/PapersNJIT/Cao_5_2013.pdf

The analysis requires some elementary radiation transport theory and a measurement of the spectral profile for absorption by the filament from two atomic species, in the present case H and Ca(II). The starting point makes use of equation (2-19) in Foukal’s “Solar Astrophysics”, and reduces to the equation shown in the following figure.
filament spectral analysis fig1 new siz.png
filament spectral analysis fig1 new siz.png (23.82 KiB) Viewed 4712 times

The equation relates the transmitted intensity to the incident intensity on an absorbing medium (the filament) in terms of the optical depth / thickness and source function, at a given wavelength lambda. Note added 08Nov2018: I've since learned that this approach was introduced by J.M. Beckers (1964, PhD thesis , Univ. Utrecht ) as the "cloud model" and is discussed in detail in the article by Jongchul Chae, Astrophysical Journal, v780, p109 (11pp), 2014. The critical underlying assumption of validity in the cloud model, used in the present analysis, is that the source function is constant in position and independent of wavelength.

Foukal and others define and explain the significance of the optical depth / thickness and the source function. In fact, here are a couple of instructive videos:

https://www.youtube.com/watch?v=hGqKOwBRKwA
https://www.youtube.com/watch?v=I_90kOpj7J4

Park et al. give the lambda dependence of the optical depth in terms of a Doppler broadened profile (Gaussian) as presented in the next figure.
filament spectral analysis fig2 siz.png
filament spectral analysis fig2 siz.png (28.44 KiB) Viewed 4770 times

In addition to the temperature,T, there is a “non thermal” velocity component contributing to the Doppler width which, according to Park et al. “arises from all kinds of unresolved motion along the line of sight, including random motions of fine-structure threads”. This is something I'm not too familiar with so, time for more research. Note added 08Nov2018: This non thermal velocity component is also referred to as "microturbulence". It is described, for example, in the lecture notes of J.B.Tatum (University of Victoria, Canada):

http://www.astro.uvic.ca/~tatum/stellatm.html

The measurement of filament spectra arising from two different types of absorbing atoms (H and Ca) allows one to determine both the temperature and the non thermal velocity. The exercise involved fitting the spectral profile of a ratio of intensities: filament spectral intensity to underlying chromospheric spectral intensity. The model equation used is shown in the figure below has, as fitting parameters, the optical depth, the absorbing linewidth, the line centre wavelength and the source function. The fitting parameters were determined for both the H and Ca absorption spectra. As shown in the figure, a correction was applied to account for the instrumental spectral profile of the spectroheliograph.
filament spectral analysis fig4 new siz.png
filament spectral analysis fig4 new siz.png (45.8 KiB) Viewed 4712 times

Once the spectral linewidths are determined, by the fitting, for both H and Ca spectra, the two unknowns which determine the linewidth (temperature and non-thermal velocity) can be calculated. The assumption is that the H and Ca filaments are in thermal equilibrium (same temperature) and have the same non thermal velocities. The figure below shows the steps used to calculate the temperature, T. Once T is known, the non thermal velocity can be determined from either the H or Ca linewidth.
filament spectral analysis fig3 siz.png
filament spectral analysis fig3 siz.png (24.44 KiB) Viewed 4770 times

The next figure shows the intensity ratio spectra and the location in the filament and adjacent portion of the chromosphere where the measurements were made. The H filament is shown, but, of course, the same measurement locations were used for both H and Ca filaments. The yellow squares in the surface image indicate the area over which an average spectrum was measured. For the spectrum plots, squares represent the actual measured intensities while curves represent the results of applying the fitting model described above. Plotting the spectra in terms of wavelength detuning from (chromospheric) line centre allowed me to superimpose the H and Ca results on to the same plot. Note that there’s an interesting “hump” near the centre of the H alpha filament absorption line. Park et al. noticed a similar result in some of their prominence emission lines (a dip near emission centre was observed) for large optical depth. This type of feature is called a “central reversal” (for both absorption and emission) and its absence in the model spectrum is due to the breakdown of a simplifying assumption used, namely that “the source function is constant along the line of sight”.
Montage_Ha_CaK_Fil-spec_comparison_v3_siz.png
Montage_Ha_CaK_Fil-spec_comparison_v3_siz.png (292.22 KiB) Viewed 4834 times

Fits to the spectral line profiles returned the values of four parameters: optical depth, line centre wavelength, linewidth and source function. In both H and Ca cases the line centre wavelength of the filament absorption was close to zero, indicating the filament was relatively stationary with respect to the underlying chromosphere. The results for optical depth, line centre wavelength and linewidth are shown in the next figure. The numerical value of the source function isn’t particularly relevant to the discussion here.
filament spectral analysis fig5 siz.png
filament spectral analysis fig5 siz.png (32.71 KiB) Viewed 4770 times

Values obtained for the temperature (8930 K) and the non thermal velocity (10.2 km/s) are consistent with those measured by Park et al. in prominences. Park et al. were able to work with very high angular resolution and, consequently, could probe many separate regions in a given prominence. They showed a distribution of temperatures and non thermal velocities. My result is averaged over a fairly large region of the filament but falls well within their distribution.

So, what have I learned from this exercise?

First, the absorption line profile can have a relatively complicated shape but it can be well modelled by some simple radiation transport theory applied to the filament as an optically thick absorber. An optical depth of 1.0 usually serves as the dividing line between optically thick (> 1.0) and optically thin (< 1.0) absorbing media. The filament is clearly very optically thick in H and much less so in Ca. Hence the filament shows up more easily (with greater contrast) in H alpha images than in Ca images.

Second, the persistence of the filament in H alpha over a wide range of wavelengths around line centre (ie. the broad filament absorption lineshape) has little to do with temperature but a lot to do with optical depth. The filament is so highly absorbing near H alpha that light incident upon it (from below) is well attenuated even for wavelengths relatively far from line centre.

Third, an amateur spectroheliograph is capable of producing interesting solar images as well as some real scientific analysis of solar features!

(Note added on October16: Fourth, feedback in a forum is sometimes invaluable! The analysis above has been edited in response to some of the discussion here: https://www.cloudynights.com/topic/6366 ... -filament/
Thank you to those who posted their ideas!)

Cheers.
Peter
Last edited by p_zetner on Thu Nov 08, 2018 2:31 pm, edited 5 times in total.


User avatar
marktownley
Librarian
Librarian
Posts: 42269
Joined: Tue Oct 18, 2011 5:27 pm
Location: Brierley Hills, UK
Has thanked: 20424 times
Been thanked: 10242 times
Contact:

Re: For Science / Math Lovers: Measurement of a Solar Filament Temperature.

Post by marktownley »

Very interesting Peter, some excellent science here!


Image
http://brierleyhillsolar.blogspot.co.uk/
Solar images, a collection of all the most up to date live solar data on the web, imaging & processing tutorials - please take a look!
LTHB
Im an EXPERT!
Im an EXPERT!
Posts: 433
Joined: Thu Mar 03, 2016 7:31 pm
Been thanked: 581 times

Re: For Science / Math Lovers: Measurement of a Solar Filament Temperature.

Post by LTHB »

Marvelous! Reading this almost makes me want to seek early retirement so I can pursue small science projects like this, too! Thanks also for the concise summary of the theoretical background, very interesting!

Regards,

Frank


User avatar
p_zetner
Way More Fun to Share It!!
Way More Fun to Share It!!
Posts: 1644
Joined: Sun Mar 18, 2012 4:59 pm
Location: Winnipeg, Manitoba, Canada
Has thanked: 166 times
Been thanked: 719 times

Re: For Science / Math Lovers: Measurement of a Solar Filament Temperature.

Post by p_zetner »

Thanks Mark and Frank.

I’m still kind of amazed at the kind of science results I can get from the spectroheliograph.
I’ve measured solar disk rotation (https://www.cloudynights.com/topic/5745 ... try7837091),
magnetic fields (https://www.cloudynights.com/topic/5723 ... try7821145), velocity fields of prominences and the chromosphere, and temperatures! I have to revisit these studies to follow up a bit more thoroughly. I view my stuff so far as sort of proof of principle.

Frank: As a fairly new retiree I have to admit that I love doing these studies. I was always a science nerd, building questionable devices and frightening my parents when I was a kid. It’s nice to return to science as a hobby in my old age!

Cheers.
Peter


christian viladrich
Way More Fun to Share It!!
Way More Fun to Share It!!
Posts: 2150
Joined: Sun Jun 14, 2015 4:46 pm
Location: France
Has thanked: 1 time
Been thanked: 2712 times
Contact:

Re: For Science / Math Lovers: Measurement of a Solar Filament Temperature.

Post by christian viladrich »

Marvellous !!


Christian Viladrich
Co-author of "Planetary Astronomy"
http://planetary-astronomy.com/
Editor of "Solar Astronomy"
http://www.astronomiesolaire.com/
User avatar
p_zetner
Way More Fun to Share It!!
Way More Fun to Share It!!
Posts: 1644
Joined: Sun Mar 18, 2012 4:59 pm
Location: Winnipeg, Manitoba, Canada
Has thanked: 166 times
Been thanked: 719 times

Re: For Science / Math Lovers: Measurement of a Solar Filament Temperature.

Post by p_zetner »

Thank you, Christian. Glad you enjoyed it.

Cheers.
Peter


User avatar
marktownley
Librarian
Librarian
Posts: 42269
Joined: Tue Oct 18, 2011 5:27 pm
Location: Brierley Hills, UK
Has thanked: 20424 times
Been thanked: 10242 times
Contact:

Re: For Science / Math Lovers: Measurement of a Solar Filament Temperature.

Post by marktownley »

Maths is fun!


Image
http://brierleyhillsolar.blogspot.co.uk/
Solar images, a collection of all the most up to date live solar data on the web, imaging & processing tutorials - please take a look!
Post Reply