quantitative analysis of observations of flux emergence
DESCRIPTION
Quantitative Analysis of Observations of Flux Emergence . More on Doppler shifts and flux emergence along PILs ! Still more constraints on photospheric electric fields due to flux emergence! Global energetic implications of flux emergence!. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/1.jpg)
Quantitative Analysis of Observations of Flux Emergence by Brian Welsch1, George Fisher1, Yan Li1, and Xudong Sun2
1Space Sciences Lab, UC-Berkeley; 2Stanford University
1. More on Doppler shifts and flux emergence along PILs!
2. Still more constraints on photospheric electric fields due to flux emergence!
3. Global energetic implications of flux emergence!
![Page 2: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/2.jpg)
a) emergence of new flux (increases total abs. flux)
b) vertical transport of cur-rents in already-emerged flux
NB: This does not increase total unsigned photospheric flux.
NB: New flux only emerges along polarity inversion lines!
Fan
& G
ibso
n 20
07
J. Ch
en a
nd J.
Kra
ll, "A
ccel
erati
on o
f cor
onal
m
ass e
jecti
ons,”
JGR
108,
141
0 (2
003)
A preliminary: When I say flux emergence, I mean “increases total unsigned flux” --- a), not b)!
![Page 3: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/3.jpg)
Magnetogram evolution can be used to estimate electric or velocity fields, E or v, in the magnetogram layer.
• “Component methods” derive v or Eh from the normal component of the ideal induction equation,
Bz/t = -c[ h x Eh ]z= [ x (v x B) ]z
• But the vector induction equation can place additional constraints on E:
B/t = -c( x E)= x (v x B),where I assume the ideal Ohm’s Law,* so v <---> E:
E = -(v x B)/c ==> E·B = 0
*One can instead use E = -(v x B)/c + R, if some model resistivity R is assumed.(I assume R might be a function of B or J or ??, but is not a function of E.)
![Page 4: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/4.jpg)
While tB provides more information about E than tBz alone, it still does not fully determine E.
Note that: tBh also depends upon vertical derivatives in Eh, which single-height magnetograms do not fully constrain.
But most importantly:Faraday’s law only relates tB to the curl of E, not E itself; the gauge electric field is unconstrained by tB.
==> Even multiple-height magnetograms won’t fix this!
(As noted, Ohm’s law & Doppler data are extra constraints.)
![Page 5: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/5.jpg)
In principle, E or v can be used to quantify the flux of magnetic energy into the corona.
• The Poynting flux of magnetic energy into the corona depends upon E =-(v x B)/c:
dU/dt = ∫ dA Sz = c ∫ dA (E x B)z /4π
• Also, coupling of the coronal field Bcor to the photospheric Bph implies E or v can form boundary conditions for data-driven, time-dependent simulations of Bcor.
![Page 6: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/6.jpg)
As George noted, emerging flux might have ~no inductive signature at the emergence site, so Doppler data can help!
Schematic illustration of flux emergence in a bipolar magnetic region, viewed in cross-section normal to the polarity inversion line (PIL).
Note the strong signature of the field change at the edges of the region, while the field change at the PIL is zero.
![Page 7: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/7.jpg)
Problem: We must remove the convective blueshift from measured Doppler shifts!
HMI’s Doppler shifts are not absolutely calibrated! (Helioseismology uses time differences in Doppler shifts.)
There’s a well-known intensity-blueshift correlation, because rising plasma (which is hotter) is brighter (see, e.g., Gray 2009; Hamilton and Lester 1999; or talk to P. Scherrer).
Because magnetic fields suppress convection, lines are redshifted in magnetized regions.
From Gray (2009): Bisectors for 13 spectral lines on the Sun are shown on an absolute velocity scale. The dots indicate the lowest point on the bisectors. (The dashed bisector is for λ6256.)
Lines formed deeper in the atmosphere, where convective upflows are present, are blue-shifted.
![Page 8: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/8.jpg)
How does this bias look in HMI data?
An automated method (Welsch & Li 2008) identified PILs in a subregion of AR 11117. Note predominance of redshifts.
![Page 9: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/9.jpg)
Fortunately, changes in LOS flux can be related to PIL Doppler shifts multiplied by transverse field strengths!
From Faraday’s law,
Since flux can only emerge or submerge at a PIL, €
Φt
= −c da (∇ ×∫ E) ⋅ ˆ n = da (∇ ×∫ (v × B)) ⋅ ˆ n
€
ΦPIL
Δt= Δ
PIL∑ ˆ l ⋅(vLOS× Btransverse)
From LOS m’gram:
Summed Dopplergram and transverse field along PIL pixels.
(Eqn. 2)
€
ΦLOS
Δt= Δ
near PIL∑ x 2[BLOS (t f ) − BLOS (ti)]
In the absence of errors, ΔΦLOS/Δt =ΔΦPIL/Δt.
(Eqn. 1)
![Page 10: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/10.jpg)
To be concrete, consider the change in negative flux within the black, dashed line: ΔΦLOS/Δt should match ΔΦPIL/Δt from Eh (or vz) on the PIL.
ΔΦLOS/ΔtΔΦPIL/Δt
NB: The analysis here applies only near disk center!
![Page 11: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/11.jpg)
We can use this constraint to calibrate the bias in the zero-velocity v0 in observed Doppler shifts!
From Eqn. 2, a bias velocity v0 implies
€
Φapparent
Δt= ΔΦPIL
Δt+ ΔΦv 0
Δt= ΔΦPIL
Δt+ v0 (Δ
PIL∑ ˆ l × Btransverse) ⋅ ˆ l
:= “magnetic length” of PIL
But ΔΦLOS should match ΔΦPIL, so we can solve for v0:
€
v0 =ΔΦapparent − ΔΦLOS
Δt BL= ΔΦbias
Δt BL(Eqn. 3)
NB: v0 should be the SAME for ALL PILs ==> solve statistically!
€
BL
![Page 12: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/12.jpg)
We have solved for v0 using data from > 100 PILs in three successive magnetogram triplets from Fig. 3.
Incorporating error bars on estimates of ΔΦ’s is still a work in progress; ideally, we’d estimate v0 using total least squares.
![Page 13: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/13.jpg)
The inferred offset velocity v0 can be used to correct Doppler shifts along PILs.
![Page 14: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/14.jpg)
Conclusions, #1
• We have demonstrated a method to correct the bias velocity v0 in HMI’s Doppler velocities from convective-blueshifts.
• Discrepancies between ΔΦLOS/Δt & ΔΦPIL/Δt in blueshift-corrected data can arise from departures from ideality --- e.g., Pariat et al. 2004.
• Hence, the method can also be used to identify the effects of magnetic diffusivity.
![Page 15: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/15.jpg)
Future Work, #1
• Incorporate uncertainties in both ΔΦ/Δt and magnetic lengths, to estimate v0 via total least squares (TLS).– Errors are present in both magnetic lengths and
ΔΦ/Δt.
• Investigate cases with large discrepancies between ΔΦLOS/Δt & ΔΦPIL/Δt for evidence of non-ideal dynamics.
![Page 16: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/16.jpg)
Outline, revisited…
1. More on Doppler shifts and flux emergence along PILs!
2. Still more constraints on photospheric electric fields due to flux emergence!
3. Global energetic implications of flux emergence!
![Page 17: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/17.jpg)
The electric field E estimated from magnetogram evolution is not necessarily consistent with ΔΦ/Δt!
• Neither Bz/t nor B/t depend directly upon B, so electric fields inferred to match time differences alone can yield incorrect dynamics.
• We used ΔΦLOS/Δt = ΔΦPIL/Δt to fix v0 – but not E!
• So it is not necessarily true that ΔΦPIL/Δt from -c( x E) is consistent with ΔΦ/Δt -- even with Doppler data!
![Page 18: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/18.jpg)
Consistency can be checked with synthetic magnetograms from our “favorite” MHD simulation...
• PILs can be identified as pixel boundaries between opposite-polarity pixels.
• The change in flux due to emergence can be calculated from Σ Eemrg·dL along all PILs, where cEemrg = vzBh x z ̂
BzBz/t
![Page 19: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/19.jpg)
How do the methods stack up?Method Flux Emergence Rate [Mx s-1]
ANMHD, ΔΦLOS/Δt 1.47 x 1017
ANMHD, ΔΦPIL/Δt 1.41 x 1017
PTD+Ohm’s Law, ΔΦPIL/Δt 1.07 x 1017
PTD+FLCT+Ohm’s, ΔΦPIL/Δt 1.12 x 1017
PTD+Doppler+Ohm’s, ΔΦPIL/Δt 1.31 x 1017
PTD+Doppler+FLCT+Ohm’s, ΔΦPIL/Δt 1.37 x 1017
Clearly, for the ANMHD data, incorporating Doppler data substantially improves the rate of flux emergence in the estimated E field.
Still, incorporating Doppler data does not ensure the estimated E field is consistent with the observed rate of flux emergence.
![Page 20: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/20.jpg)
Away from disk center, the PTD+Doppler approach can be still be used, but does not properly capture emergence.
• Away from disk center, Doppler shifts along PILs of the line-of-sight (LOS) field still constrains E --- but both the normal and horizontal components are constrained.
• Consistency of Eh with the change in radial magnetic flux therefore imposes an independent constraint.
![Page 21: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/21.jpg)
Conclusions, #2
• Flux emergence ΔΦ/Δt places an additional, integral constraint on photospheric electric fields on PILs , so ΔΦPIL/Δt matches.
• Methods of estimating photospheric electric fields that do not incorporate this constraint are inconsistent with observations.
![Page 22: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/22.jpg)
Future Work, #2• Most (all?) methods of estimating horizontal electric
fields at the photosphere are already “inductive:” they match the constraint on Eh from ΔBz/Δt
• Non-inductive electric fields should be added onto Eh to enforce consistency with ΔΦ/Δt. How to do so must be investigated!
• Uniform corrections to either Eh or vz (not equivalent!) could be made.• Better: correct Eh in regions with emergence!
![Page 23: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/23.jpg)
Outline, revisited…
1. More on Doppler shifts and flux emergence along PILs!
2. Still more constraints on photospheric electric fields due to flux emergence!
3. Global energetic implications of flux emergence!
![Page 24: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/24.jpg)
The hypothetical coronal magnetic field with lowest energy is current-free, or “potential.”
• For a given coronal field Bcor, the coronal magnetic energy is:
U dV (Bcor · Bcor)/8.
• The lowest energy coronal field would have current J = 0, and Ampére says 4πJ/c = x B, so x Bmin
= 0.
• Since Bmin is curl-free, Bmin = -; and since B⋅ min = 0 = 2, the
Neumann condition from photospheric Bnormal determines .
Umin dV (Bmin · Bmin
)/8
• The difference Ufree = [U – Umin] is “free” energy stored in the
corona, which can be suddenly released in flares or CMEs.
![Page 25: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/25.jpg)
Flux emergence can induce currents on separatrices between old & new flux, even if the emerged flux is current-free.
Hale’s Law implies that new flux is typically positioned favorably to reconnect with old flux.
ti
tf
Within one hemisphere: Trans-equatorial:
Not a new idea! See, e.g., Hayvaerts et al. 1977
![Page 26: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/26.jpg)
We know the pre- and post- emergence minimum energy states.
• The pre-emergence minimum energy is Upre dV (Bpre
· Bpre )/8, where
x Bpre = 0,
and Bn(ti)⎮∂V = n · Bpre⎮∂V
• The post-emergence minimum energy is Upost dV (Bpost
· Bpost )/8, where
x Bpost = 0,
and Bn(tf)⎮∂V = n · Bpost⎮∂V = n · (Bpre+ Bemrg) ⎮∂V, and Bemrg is the emerged flux at the surface.
![Page 27: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/27.jpg)
From Bemrg, we can also estimate the minimum energy introduced into the corona.
• The min. energy field added by the new flux is:Uadd dV (Badd
· Badd )/8, where
x Badd = 0,
and (Bn(tf) - Bn(ti))⎮∂V = n · Bemrg⎮∂V = n · Badd⎮∂V
• New flux emerged must be separated from previously emerged flux; it’s easiest to do so for entirely new active regions.
![Page 28: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/28.jpg)
Even without knowing Bcor, we can still, in some cases, place bounds on the free energy from flux emergence!
• Linearity in potential fields implies the actual post-emergence energy satisfies:
Uactual ≥ Upost = dV (Bpost· Bpost )/8
= dV (Bpre+ Badd)·(Bpre+ Badd )/8
• But the actual post-emergence energy also satisfies: Uactual ≥ dV (Bpre· Bpre
)/8 + dV (Badd· Badd )/8
• The difference between these bounds is: Uinteraction dV (Bpre· Badd
)/4
• This interaction energy can be positive or negative.
![Page 29: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/29.jpg)
Negative interaction energies imply that reconnection between new and old flux is energetically favorable.
In the Hale’s law examples that I showed earlier, the interaction energy is more than 1031 ergs.
ti
tf
Within one hemisphere: Trans-equatorial:
To repeat: This interaction energy entirely neglects the presence of any electric currents in either the old or emerging flux systems; it arises entirely from the global magnetic configuration.
![Page 30: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/30.jpg)
Conclusions, #3
• Flux emergence can drastically alter the global potential magnetic field (the PFSS field), the large-scale corona’s minimum energy state.
• For negative interaction energies, a bound on free energy introduced by flux emergence is:
Ufree ≥ - dV (Bpre· Badd )/4
![Page 31: Quantitative Analysis of Observations of Flux Emergence](https://reader033.vdocument.in/reader033/viewer/2022051821/56816253550346895dd29d37/html5/thumbnails/31.jpg)
Future Work, #3
• Relationships between free energy implied from the observed evolution in PFSS fields due to flux emergence and energy release (flare activity, CMEs) should be investigated.
• I plan to put an undergrad on the case this fall!