May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area...

36
Magnetic structure over plage region May 5 2020 Tetsu Anan Photosphere Chromosphere HAZEL result Really?

Transcript of May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area...

Page 1: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Magnetic structure over plage region May 5 2020 Tetsu Anan

Photosphere

Chromosphere

HAZEL result Really?

Page 2: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Contents

• Observed plage region• IRIS inversion• GREGOR/GRIS data– Stray light– IQUV crosstalks– Expected choromospheric B field from Hanle diagram

• HAZEL inversion– Photosphere– Chromosphere

DefinitionΘ! : Inclination of B in LOS frame)! : Inclination of B in solar frameΦ! : Azimuth of B in LOS frame+! : Azimuth of B in solar frame

Page 3: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Where is the observed active region• Heliocentric coordinates (693”, −237”)• Solar radius 958”

Page 4: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Where is the observed active region

Page 5: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Plage has single polarityReply to a comment

Page 6: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

IRIS inversion ’IRIS2’Reply to a comment

Page 7: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Stray light of GRIS �4 %Reply to a comment

Stray light derived by comparing intensitiesat the line center is 13%Wing is not consistent with the atlas

Stray light derived by comparing equivalent widths in a λ range between dashed lines is 4%

Atlas: Sac Peak Solar Flux Atlas by Kurucz et al.(1984)

GRIS Disk centerATLASCalibrated

Page 8: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Definition of Stokes parametersComparison with HMI SHARP

• GRIS (email from Manolo)– +Q is the solar N-S direction– +U rotate counter-clockwise from +Q on the solar surface– +V is right circular polarization on the solar surface

Not plage data

Reply to a comment

HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates

Page 9: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Cross talk IóQ,U,V �0.01-0.02%Plage data

Polarization degree in continuum

Solid: mean polarization degreeDotted: standard deviation of

polarization degree

The crosstalk is too small to significantly change Stokes profiles

Page 10: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Cross talk VóQ, VóUPlage data

Kuhn et al. 1994Assume small Q <=> U cross talk,

because of the comparison with Hinode!!"# = !$%! − $%$%! − &'$%!%!"# = %$%! − (!!"#'!"# = '$%! − )!!"#

a, b, c, d are not homogeneous spatially => Small cross talksa b

c d

Time (# of data)

Page 11: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Crosstalk• Collados 2003 SPIE– Apply to the umbra of a spot, where the magnetic

field ensures (together with the large magnetic sensitivity of near infrared spectral limes) that the sigma and pi components are completely split apart

– Results of umbra data taken by GREGOR/GRIS(Si�10827 Å can not be split completely to σ and π on the umbrae)

• !!→# = #! $ #"!! % !"

= 0.009

• !!→& = &! $ &"!! % !"

= −0.004• (' − (( = !#→! )' + )( + !&→! +' + +(

– !!→#=0.01– !$→# = −0.09

Reply to a comment

Page 12: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Linear polarization

• Method 20190828_{2,3,4}.pro1. Fit Stokes-I of He 10830 with a gaussian to get

line center wavelength and line width2. Fixing the line center wavelength and the line

width, Stokes-Q/Ic & -U/Ic profiles are fitted with a gaussian=> Stokes-Q&U amplitude

• Noise: standard deviation of Q/Ic & U/Ic at continuum

New slide

Page 13: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Thin: noiseThick: Pl= !! + #!

Noise improve withspatial binning

Signal in-depends on spatial binning

Data with bin 2will be inverted

Oct. 3 -005plage between sunspots, AR 12723

Polarization angleis homogeneous

New slide

Page 14: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Linear polarization in He�1083 nm

1. Fit Stokes-I of He 10830 with a gaussian to get line center wavelength and line width

2. Fixing the line center wavelength and the line width, Stokes-Q/Ic & -U/Ic profiles are fitted with a gaussian=> Stokes-Q&U amplitude Linear polarization aligns

along radial direction

Plage data

Page 15: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Hanle diagram

• ! = 0°• Heliocentric coordinates (693”, −237”)

Solar radius 958”=> % = 49.8°

• +Q = NS direction=> * = 108°

!!0°, 180°10°, 170°20°, 160°30°, 150°40°, 140°50°, 130°60°, 120°70°, 110°80°, 100°90°

ObservationNo 90° ambiguity!! ≈ 90°

& = 0.7∆+ = 11 km/s1 = 200 G

HAZEL

Page 16: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

HAZEL inversion

Page 17: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

HAZELLine of sight angles

• Line of sight angles– ! = 49.8°– ( = 0°– * = 90° (+Q is parallel to limb)

• Rotation of coordinate

XY

θQ

θQ

!! = tan"# &' = −18.9°QQ’

.′0′ = cos 2!! sin 2!!

−sin 2!! cos 2!!.0

Heliocentric coordinates (X,Y)=(693”, −237”)Solar radius 958”

Page 18: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

HAZEL parameters

PhotosphereFunctions of optical depth

• BX in line of sight reference frame(LOS frame)

• BY in LOS frame

• BZ in LOS frame

• Temperature

• Electron pressure

• Microturbulent velocity

• Macroscopic velocity

• Filling factor = 1 (fixed)

Chromosphere• B

• Θ! inclination angle in LOS frame

• Φ! azimuth angle in LOS frame

• Optical depth

• Doppler velocity

• Doppler width

• Damping parameter

• Height of slab = 3.0” (fixed)

• Enhancement factor = 1 (fixed)

• Filling factor = 1 (fixed)

Page 19: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

HAZELThe axis of coordinates for B field

This figure is for chromospheric B fieldX axis for photospheric B is parallel to e1

Test #1 #2γ 90� 0�BX (Si) −69 G −318 GBY (SI) 318 G −69 GΘ! (He) 51� 51�Φ! (He) 349� −10�

Page 20: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

HAZEL, B in chromosphereambiguities

Sequence of HAZEL1. HAZEL inversion => result !!, Θ"! , Φ"! , ⋯2. HAZEL inversion with fixing Φ" = Φ"! + 180°3. HAZEL inversion with fixing Φ" = Φ"! + 90°4. HAZEL inversion with fixing Φ" = Φ"! − 90°

5. Calculate mean square of residual in Q and U within a wavelength range of line center � line width

MS = 10#1$

2$ − 3it$noise$

%

6. Best solution has minimum of MS& +MS'7. Select solution that has

MS& −MS&"()* < 1 & MS' −MS'"()* < 1

Page 21: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Fitting

Stokes I, Si Stokes Q, Si

Stokes U, Si Stokes V, Si

Stokes I, HeStokes Q, He

Stokes U, He Stokes V, He

Page 22: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

line of sight component of B

Photo BLOS Chromo BLOS

Page 23: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

B in photosphere

Θ! Φ!

!"#° ambiguityX for photospheric BX is parallel to limb

Θ! ≈ 45°90�

Correct the result

Page 24: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

B in photosphere, 180° ambiguity 2 solutions

!!

!! "!

"!

LOS

! = 49.8°

Θ! ≈ 45°

Sun center

Correct the result

Page 25: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Distributions of the 2 solutionsGRIS + HAZEL2

New slide

Page 26: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Comparison with HMI SHARP180° ambiguity is resolved using a minimum energy method (Metcalf 1994; Leka et al. 2009)

The vertical solution of GRIS inversion is consistent with HMI

New slide

HMI: a Lambert Cylindrical Equal-Area projection

Page 27: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Comparison with HMI SHARP180° ambiguity is resolved using a minimum energy method (Metcalf 1994; Leka et al. 2009)

Spatial sampling: 0.2”, 0.1” for GRIS (diffraction limit 0.2”), 0.5” for HMI

New slide

Page 28: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Small horizontal B field over plage• Hinode/SOT SP reveals there is horizontal B fields over plage region

– Lifetime 1 – 10 min– Randomly located spatially and temporally– Not associated with existing long lived B field and appear in a region of

insignificant vertical fields– First appears inside the granule, subsequently moves to the inter-

granule lane, and finally disappears in the inter-granular lane

Ishikawa & Tsuneta 2009

New slide

Page 29: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Short summary of the 180�ambiguity

• Two solutions of the 180° ambiguity is classified under vertical and horizontal B fields– Vertical solution is consistent with HMI SHARP CEA– For magnetic elements, which is larger than 1”, the

vertical solution is good– There is a possibility that small horizontal B field exists

• Because it is randomly located spatially and temporally (and it is anisotropic), it is impossible to solve 180�ambiguity for the small B field in a region near the solar limb

• I will “discuss” results in terms of spatial scale

New slide

Page 30: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

B in photosphere

0° radial direction

#! $!B

Page 31: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

B in chromosphereInclination & azimuth angle

Dotted: The solutions of the 4 inversionSolid : The selected solution No !"° ambiguity (see Hanle diagram)

No 90° ambiguity

van Vleck angle

LOS

& = 49.8°

&! ≈ 90°

Sun center

Θ! Φ!

&! .!

No /0"° ambiguity• 3 cos &!" − 1 sinΘ!" =

3 cos &′!" − 1 sinΘ′!"• &! ≈ 90°, Θ! ≈ 45°Þ −0.5 ≈ 3 cos &′!" − 1 sinΘ′!" (1)• Φ! should be ≈ 180°

=> Θ′! ≈ &#! + & (2)• Substitute (2) into (1) => &#! ≈ 80°Stokes V rejects this solution

4 candidates

Page 32: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

B in chromosphere

0° radial direction

#! $!B

Page 33: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

New slideContinuum images taken by HMI, IRIS spectra, and GRIS are spatially aligned carefully

Page 34: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Discussion:Is it true that chromospheric B is horizontal?• Orientation of liner polarization is consistent with Hinode for Sunspot data

– Axis of coordinates can be OK– Cross talk between Stokes Q & U can be OK

• Cross talks among IQUVs are small– 2 methods: Kuhn et al. 1994 & Collados 2003 SPIE

• Horizontal B field is expected from Hanle diagram– Hanle diagram seems to be appropriate

• If it was not true, assumptions of the HAZEL would be wrong– Radiation field?

• The assumption can be OK, because it is granule pattern, which is spatially homogeneous in continuum radiation

– 1 slab atmospheric model?

Page 35: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Discussion: radiation field

• If B is normal to the solar surface in the chromosphere, the Hanle effect is small and linear polarization is produced by anisotropy of radiation field– Radiation irradiating from side should be stronger

than that from vertical– The anisotropy should be homogeneous in the

field of view and oriented to the radial direction on the solar surface

New slide

Page 36: May 5 2020 Tetsu Anananan/shed/... · 5/5/2020  · HMI: a Lambert Cylindrical Equal-Area projection, Not heliocentric coordinates. Cross talk IóQ,U,V 0.01-0.02% Plagedata Polarization

Comment

• Solarsoft: rot_xy enables us to track plage– JSOC��

• GRIS�HMI�"������–������contrast

• Ishikawa et al.#�����$���)%���(���!&(� �(�� �

• Radiative flux�He intensity����'��