-
SOFIA/HAWC+ Far-Infrared Polarimetric Large-Area CMZ Exploration (FIREPLACE) Survey III: Full Survey Data Set
Authors:
Dylan Paré,
Natalie O. Butterfield,
David T. Chuss,
Jordan A. Guerra,
Jeffrey I. Iuliano,
Kaitlyn Karpovich,
Mark R. Morris,
Edward Wollack
Abstract:
We present the second data release (DR2) of the Far-Infrared Polarimetric Large-Area CMZ Exploration (FIREPLACE) survey. This survey utilized the Stratospheric Observatory for Infrared Astronomy (SOFIA) High-resolution Airborne Wideband Camera plus (HAWC+) instrument at 214 $μ$m (E-band) to observe dust polarization throughout the Central Molecular Zone (CMZ) of the Milky Way. DR2 consists of obse…
▽ More
We present the second data release (DR2) of the Far-Infrared Polarimetric Large-Area CMZ Exploration (FIREPLACE) survey. This survey utilized the Stratospheric Observatory for Infrared Astronomy (SOFIA) High-resolution Airborne Wideband Camera plus (HAWC+) instrument at 214 $μ$m (E-band) to observe dust polarization throughout the Central Molecular Zone (CMZ) of the Milky Way. DR2 consists of observations that were obtained in 2022 covering the region of the CMZ extending roughly from the Brick to the Sgr C molecular clouds (corresponding to a roughly 1$^{\circ}$ $\times$ 0.75$^{\circ}$ region of the sky). We combine DR2 with the first FIREPLACE data release to obtain full coverage of the CMZ (a 1.5$^{\circ}$ $\times$0.75$^{\circ}$ region of the sky). After applying total and polarized intensity significance cuts on the full FIREPLACE data set we obtain $\rm\sim$65,000 Nyquist-sampled polarization pseudovectors. The distribution of polarization pseudovectors confirms a bimodal distribution in the CMZ magnetic field orientations, recovering field components that are oriented predominantly parallel or perpendicular to the Galactic plane. These magnetic field orientations indicate possible connections between the previously observed parallel and perpendicular distributions. We also inspect the magnetic fields toward a set of prominent CMZ molecular clouds (the Brick, Three Little Pigs, 50 km s$\rm^{-1}$, Circum-nuclear Disk, CO 0.02-0.02, 20 km s$\rm^{-1}$, and Sgr C), revealing spatially varying magnetic fields that generally trace the morphologies of the clouds. We find evidence that compression from stellar winds and shear from tidal forces are prominent mechanisms influencing the structure of the magnetic fields observed within the clouds.
△ Less
Submitted 29 April, 2024; v1 submitted 10 January, 2024;
originally announced January 2024.
-
SOFIA/HAWC+ Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) II: Detection of a Magnetized Dust Ring in the Galactic Center
Authors:
Natalie O. Butterfield,
Jordan A. Guerra,
David T. Chuss,
Mark R. Morris,
Dylan Pare,
Edward J. Wollack,
Allison H. Costa,
Matthew J. Hankins,
Johannes Staguhn,
Ellen Zweibel
Abstract:
We present the detection of a magnetized dust ring (M0.8-0.2) in the Central Molecular Zone (CMZ) of the Galactic Center. The results presented in this paper utilize the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey (i.e., FIREPLACE I; Butterfield et al. 2023). The FIREPLACE survey is a 214 $μ$m polarimetic survey of the Galactic Center usi…
▽ More
We present the detection of a magnetized dust ring (M0.8-0.2) in the Central Molecular Zone (CMZ) of the Galactic Center. The results presented in this paper utilize the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey (i.e., FIREPLACE I; Butterfield et al. 2023). The FIREPLACE survey is a 214 $μ$m polarimetic survey of the Galactic Center using the SOFIA/HAWC+ telescope. The M0.8-0.2 ring is a region of gas and dust that has a circular morphology with a central depression. The dust polarization in the M0.8-0.2 ring implies a curved magnetic field that traces the ring-like structure of the cloud. We posit an interpretation in which an expanding shell compresses and concentrates the ambient gas and magnetic field. We argue that this compression results in the strengthening of the magnetic field, as we infer from the observations toward the interior of the ring.
△ Less
Submitted 29 April, 2024; v1 submitted 3 January, 2024;
originally announced January 2024.
-
SOFIA/HAWC+ Far-InfraRed Polarimetric Large Area CMZ Exploration (FIREPLACE) Survey I: General Results from the Pilot Program
Authors:
Natalie O. Butterfield,
David T. Chuss,
Jordan A. Guerra,
Mark R. Morris,
Dylan Pare,
Edward J. Wollack,
C. Darren Dowell,
Matthew J. Hankins,
Kaitlyn Karpovich,
Javad Siah,
Johannes Staguhn,
Ellen Zweibel
Abstract:
We present the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey. The survey was taken using the 214-micron band of the HAWC+ instrument with the SOFIA telescope (19.6$'$ resolution; 0.7 pc). In this first data release we present dust polarization observations covering a ~0.5$°$ region of the Galactic Center's Central Molecular Zone (CMZ), appr…
▽ More
We present the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey. The survey was taken using the 214-micron band of the HAWC+ instrument with the SOFIA telescope (19.6$'$ resolution; 0.7 pc). In this first data release we present dust polarization observations covering a ~0.5$°$ region of the Galactic Center's Central Molecular Zone (CMZ), approximately centered on the Sgr B2 complex. We detect ~25,000 Nyquist-sampled polarization pseudovectors, after applying the standard SOFIA cuts for minimum signal-to-noise in fractional polarization and total intensity of 3 and 200, respectively. Analysis of the magnetic field orientation suggests a bimodal distribution in the field direction. This bimodal distribution shows enhancements in the distribution of field directions for orientations parallel and perpendicular to the Galactic plane, which is suggestive of a CMZ magnetic field configuration with polodial and torodial components. Furthermore, a detailed analysis of individual clouds included in our survey (i.e., Sgr B2, Sgr B2-NW, Sgr B2-Halo, Sgr B1, and Clouds-E/F) shows these clouds have fractional polarization values of 1--10% at 214-micron, with most of the emission having values $<$5%. A few of these clouds (i.e., Sgr B2, Clouds-E/F) show relatively low fractional polarization values toward the cores of the cloud, with higher fractional polarization values toward the less dense periphery. We also observe higher fractional polarization towards compact HII regions which could indicate an enhancement in the grain alignment in the dust surrounding these sources.
△ Less
Submitted 4 December, 2023; v1 submitted 2 June, 2023;
originally announced June 2023.
-
The Strength of the Sheared Magnetic Field in the Galactic's Circum-Nuclear Disk
Authors:
J. A. Guerra,
E. Lopez-Rodriguez,
D. T. Chuss,
N. O. Butterfield,
J. T. Schmelz
Abstract:
Recent high-resolution 53-$μ$m polarimetric observations from SOFIA/HAWC+ have revealed the inferred plane-of-the-sky magnetic field (B-field) orientation in the Galactic center's Circum-Nuclear Disk (CND). The B-field is mostly aligned with the steamers of ionized material falling onto Sgr A* at large, differential velocities (shear). In such conditions, estimating the B-field strength with the `…
▽ More
Recent high-resolution 53-$μ$m polarimetric observations from SOFIA/HAWC+ have revealed the inferred plane-of-the-sky magnetic field (B-field) orientation in the Galactic center's Circum-Nuclear Disk (CND). The B-field is mostly aligned with the steamers of ionized material falling onto Sgr A* at large, differential velocities (shear). In such conditions, estimating the B-field strength with the ``classical" Davis-Chandrasekhar-Fermi (DCF) method does not provide accurate results. We derive a ``modified'' DCF method by solving the ideal MHD equations from first principles considering the effects of a large-scale, shear flow on the propagation of a fast magnetosonic wave. In the context of the DCF approximation, both the value of the shear and its Laplacian affect the inferred B-field strength. Using synthetic polarization data from MHD simulations for a medium dominated by shear flows, we find that the ``classical'' DCF determines B-field strengths only within $>50$\% of the true value where the ``modified" DCF results are improved significantly ($\sim$3-22\%). Applying our ``modified'' DCF method to the CND revealed B-field strengths of 1 - 16 mG in the northern arm, 1 - 13 mG in the eastern arm, and 3 - 27 mG in the western arm at spatial scales $\lesssim1$ pc, with median values of $5.1\pm0.8$, $4.0\pm1.2$, and $8.5\pm2.3$ mG, respectively. The balance between turbulent gas energy (kinetic plus hydrostatic) and turbulent magnetic energy densities suggest that, along the magnetic-field-flow direction, magnetic effects become less dominant as the shear flow increases and weakens the B-field via magnetic convection. Our results indicate that the transition from magnetically to gravitationally dominated accretion of material onto Sgr A* starts at distances $\sim$ 1 pc.
△ Less
Submitted 10 June, 2023; v1 submitted 13 April, 2023;
originally announced April 2023.
-
Tracing PAH Emission in $λ$-Orionis Using COBE/DIRBE Data
Authors:
David T. Chuss,
Brandon S. Hensley,
Alan J. Kogut,
Jordan A. Guerra,
Hayley C. Nofi,
Javad Siah
Abstract:
We use archival COBE/DIRBE data to construct a map of polycyclic aromatic hydrocarbon (PAH) emission in the $λ$-Orionis region. The presence of the 3.3 $μ$m PAH feature within the DIRBE 3.5 $μ$m band and the corresponding lack of significant PAH spectral features in the adjacent DIRBE bands (1.25, 2.2, and 4.9 $μ$m) enable estimation of the PAH contribution to the 3.5 $μ$m data. Having the shortes…
▽ More
We use archival COBE/DIRBE data to construct a map of polycyclic aromatic hydrocarbon (PAH) emission in the $λ$-Orionis region. The presence of the 3.3 $μ$m PAH feature within the DIRBE 3.5 $μ$m band and the corresponding lack of significant PAH spectral features in the adjacent DIRBE bands (1.25, 2.2, and 4.9 $μ$m) enable estimation of the PAH contribution to the 3.5 $μ$m data. Having the shortest wavelength of known PAH features, the 3.3 $μ$m feature probes the smallest PAHs, which are also the leading candidates for carriers of anomalous microwave emission (AME). We use this map to investigate the association between the AME and the emission from PAH molecules. We find that the spatial correlation in $λ$-Orionis is higher between AME and far-infrared dust emission (as represented by the DIRBE 240 $μ$m map) than it is between our PAH map and AME. This finding, in agreement with previous studies using PAH features at longer wavelengths, is in tension with the hypothesis that AME is due to spinning PAHs. However, the expected correlation between mid-infrared and microwave emission could potentially be degraded by different sensitivities of each emission mechanism to local environmental conditions even if PAHs are the carriers of both.
△ Less
Submitted 7 November, 2022; v1 submitted 18 August, 2022;
originally announced August 2022.
-
SOFIA observations of 30 Doradus: II -- Magnetic fields and large scale gas kinematics
Authors:
Le Ngoc Tram,
Lars Bonne,
Yue Hu,
Enrique Lopez-Rodriguez,
Jordan A. Guerra,
Pierre Lesaffre,
Antoine Gusdorf,
Thiem Hoang,
Min-Young Lee,
Alex Lazarian,
B-G Andersson,
Simon Coude,
Archana Soam,
William D. Vacca,
Hyeseung Lee,
Michael Gordon
Abstract:
The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B-fiel…
▽ More
The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B-fields) of 30 Doradus inferred from the far-infrared polarimetric observations by SOFIA/HAWC+ at 89, 154, and 214$\,μ$m. The B-field morphology is complex, showing bending structures around R$\,$136. In addition, we use high spectral and angular resolution [\textsc{CII}] observations from SOFIA/GREAT and CO(2-1) from APEX. The kinematic structure of the region correlates with the B-field morphology and shows evidence of multiple expanding shells. Our B-field strength maps, estimated using the Davis-Chandrasekhar-Fermi method and structure-function, show variations across the cloud within a maximum of 600, 450, and 350$\,μ$G at 89, 154, and 214$\,μ$m, respectively. We estimated that the majority of the 30 Doradus clouds are sub-critical and sub-Alfvénic. The probability distribution function of the gas density shows that the turbulence is mainly compressively driven, while the plasma beta parameter indicates supersonic turbulence. We show that the B-field is sufficient to hold the cloud structure integrity under feedback from R$\,$136. We suggest that supersonic compressive turbulence enables the local gravitational collapse and triggers a new generation of stars to form. The velocity gradient technique (VGT) using [\textsc{CII}] and CO(2-1) is likely to confirm these results.
△ Less
Submitted 12 December, 2022; v1 submitted 24 May, 2022;
originally announced May 2022.
-
The Flare Likelihood and Region Eruption Forecasting (FLARECAST) Project: Flare forecasting in the big data & machine learning era
Authors:
M. K. Georgoulis,
D. S. Bloomfield,
M. Piana,
A. M. Massone,
M. Soldati,
P. T. Gallagher,
E. Pariat,
N. Vilmer,
E. Buchlin,
F. Baudin,
A. Csillaghy,
H. Sathiapal,
D. R. Jackson,
P. Alingery,
F. Benvenuto,
C. Campi,
K. Florios,
C. Gontikakis,
C. Guennou,
J. A. Guerra,
I. Kontogiannis,
V. Latorre,
S. A. Murray,
S. -H. Park,
S. von Stachelski
, et al. (3 additional authors not shown)
Abstract:
The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) differe…
▽ More
The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) different ML techniques, also on equal footing, to optimize the immense Big Data parameter space created by these many predictors; third, the establishment of a robust, three-pronged communication effort oriented toward policy makers, space-weather stakeholders and the wider public. FLARECAST pledged to make all its data, codes and infrastructure openly available worldwide. The combined use of 170+ properties (a total of 209 predictors are now available) in multiple ML algorithms, some of which were designed exclusively for the project, gave rise to changing sets of best-performing predictors for the forecasting of different flaring levels. At the same time, FLARECAST reaffirmed the importance of rigorous training and testing practices to avoid overly optimistic pre-operational prediction performance. In addition, the project has (a) tested new and revisited physically intuitive flare predictors and (b) provided meaningful clues toward the transition from flares to eruptive flares, namely, events associated with coronal mass ejections (CMEs). These leads, along with the FLARECAST data, algorithms and infrastructure, could help facilitate integrated space-weather forecasting efforts that take steps to avoid effort duplication. In spite of being one of the most intensive and systematic flare forecasting efforts to-date, FLARECAST has not managed to convincingly lift the barrier of stochasticity in solar flare occurrence and forecasting: solar flare prediction thus remains inherently probabilistic.
△ Less
Submitted 12 May, 2021;
originally announced May 2021.
-
Ensemble Forecasting of Major Solar Flares: Methods for Combining Models
Authors:
Jordan A. Guerra,
Sophie A. Murray,
D. Shaun Bloomfield,
Peter T. Gallagher
Abstract:
One essential component of operational space weather forecasting is the prediction of solar flares. With a multitude of flare forecasting methods now available online it is still unclear which of these methods performs best, and none are substantially better than climatological forecasts. Space weather researchers are increasingly looking towards methods used by the terrestrial weather community t…
▽ More
One essential component of operational space weather forecasting is the prediction of solar flares. With a multitude of flare forecasting methods now available online it is still unclear which of these methods performs best, and none are substantially better than climatological forecasts. Space weather researchers are increasingly looking towards methods used by the terrestrial weather community to improve current forecasting techniques. Ensemble forecasting has been used in numerical weather prediction for many years as a way to combine different predictions in order to obtain a more accurate result. Here we construct ensemble forecasts for major solar flares by linearly combining the full-disk probabilistic forecasts from a group of operational forecasting methods (ASAP, ASSA, MAG4, MOSWOC, NOAA, and MCSTAT). Forecasts from each method are weighted by a factor that accounts for the method's ability to predict previous events, and several performance metrics (both probabilistic and categorical) are considered. It is found that most ensembles achieve a better skill metric (between 5\% and 15\%) than any of the members alone. Moreover, over 90\% of ensembles perform better (as measured by forecast attributes) than a simple equal-weights average. Finally, ensemble uncertainties are highly dependent on the internal metric being optimized and they are estimated to be less than 20\% for probabilities greater than 0.2. This simple multi-model, linear ensemble technique can provide operational space weather centres with the basis for constructing a versatile ensemble forecasting system -- an improved starting point to their forecasts that can be tailored to different end-user needs.
△ Less
Submitted 1 August, 2020;
originally announced August 2020.
-
Far-Infrared Polarization Spectrum of the OMC-1 Star-Forming Region
Authors:
Joseph M. Michail,
Peter C. Ashton,
Marc G. Berthoud,
David T. Chuss,
C. Darren Dowell,
Jordan A. Guerra,
Doyal A. Harper,
Giles Novak,
Fabio P. Santos,
Javad Siah,
Ezra Sukay,
Aster Taylor,
Le Ngoc Tram,
John E. Vaillancourt,
Edward J. Wollack
Abstract:
We analyze the wavelength dependence of the far-infrared polarization fraction toward the OMC-1 star forming region using observations from HAWC+/SOFIA at 53, 89, 154, and 214 $μ$m. We find that the shape of the far-infrared polarization spectrum is variable across the cloud and that there is evidence of a correlation between the slope of the polarization spectrum and the average line-of-sight tem…
▽ More
We analyze the wavelength dependence of the far-infrared polarization fraction toward the OMC-1 star forming region using observations from HAWC+/SOFIA at 53, 89, 154, and 214 $μ$m. We find that the shape of the far-infrared polarization spectrum is variable across the cloud and that there is evidence of a correlation between the slope of the polarization spectrum and the average line-of-sight temperature. The slope of the polarization spectrum tends to be negative (falling toward longer wavelengths) in cooler regions and positive or flat in warmer regions. This is very similar to what was discovered in $ρ$ Oph A via SOFIA polarimetry at 89 and 154 $μ$m. Like the authors of this earlier work, we argue that the most natural explanation for our falling spectra is line-of-sight superposition of differing grain populations, with polarized emission from the warmer regions and less-polarized emission from the cooler ones. In contrast with the earlier work on $ρ$ Oph A, we do not find a clear correlation of polarization spectrum slope with column density. This suggests that falling spectra are attributable to variations in grain alignment efficiency in a heterogeneous cloud consistent with radiative torques theory. Alternative explanations in which variations in grain alignment efficiency are caused by varying gas density rather than by varying radiation intensity are disfavored.
△ Less
Submitted 24 December, 2020; v1 submitted 1 August, 2020;
originally announced August 2020.
-
Maps of Magnetic Field Strength in the OMC-1 using HAWC+ FIR Polarimetric data
Authors:
Jordan A. Guerra,
David T. Chuss,
C. Darren Dowell,
Martin Houde,
Joseph M. Michail,
Javad Siah,
Edward J. Wollack
Abstract:
Far-infrared (FIR) dust polarimetry enables the study of interstellar magnetic fields via tracing of the polarized emission from dust grains that are partially aligned with the direction of the field. The advent of high quality polarimetric data has permitted the use of statistical methods to extract both the direction and magnitude of the magnetic field. In this work, the Davis-Chandrasekhar-Ferm…
▽ More
Far-infrared (FIR) dust polarimetry enables the study of interstellar magnetic fields via tracing of the polarized emission from dust grains that are partially aligned with the direction of the field. The advent of high quality polarimetric data has permitted the use of statistical methods to extract both the direction and magnitude of the magnetic field. In this work, the Davis-Chandrasekhar-Fermi technique is used to make maps of the plane-of-sky (POS) component of the magnetic field in the Orion Molecular Cloud (OMC-1) by combining polarization maps at 53, 89, 154 and 214 \micron\ from HAWC+/SOFIA with maps of density and velocity dispersion. In addition, maps of the local dispersion of polarization angles are used in conjuction with Zeeman measurements to estimate a map of the strength of the line-of-sight (LOS) component of the field. Combining these maps, information about the three-dimensional magnetic field configuration (integrated along the line-of-sight) is inferred over the OMC-1 region. POS magnetic field strengths of up to 2 mG are observed near the BN/KL object, while the OMC-1 bar shows strengths of up to a few hundred $μ$G. These estimates of the magnetic field components are used to produce maps of the mass-to-magnetic flux ratio ($M/Φ$) -- a metric for probing the conditions for star formation in molecular clouds -- and determine regions of sub- and super-criticality in OMC-1. Such maps can provide invaluable input and comparison to MHD simulations of star formation processes in filamentary structures of molecular clouds.
△ Less
Submitted 28 December, 2020; v1 submitted 9 July, 2020;
originally announced July 2020.
-
Which Photospheric Characteristics are Most Relevant to Active-Region Coronal Mass Ejections?
Authors:
Ioannis Kontogiannis,
Manolis K. Georgoulis,
Jordan A. Guerra,
Sung-Hong Park,
D. Shaun Bloomfield
Abstract:
We investigate the relation between characteristics of coronal mass ejections and parameterizations of the eruptive capability of solar active regions widely used in solar flare prediction schemes. These parameters, some of which are explored for the first time, are properties related to topological features, namely, magnetic polarity inversion lines (MPILs) that indicate large amounts of stored n…
▽ More
We investigate the relation between characteristics of coronal mass ejections and parameterizations of the eruptive capability of solar active regions widely used in solar flare prediction schemes. These parameters, some of which are explored for the first time, are properties related to topological features, namely, magnetic polarity inversion lines (MPILs) that indicate large amounts of stored non-potential (i.e. free) magnetic energy. We utilize the Space Weather Database of Notifications, Knowledge, Information (DONKI) and the Large Angle and Spectrometric Coronograph (LASCO) databases to find flare-associated coronal mass ejections and their kinematic characteristics while properties of MPILs are extracted from Helioseismic and Magnetic Imager (HMI) vector magnetic-field observations of active regions to extract the properties of source-region MPILs. The correlation between all properties and the characteristics of CMEs ranges from moderate to very strong. More significant correlations hold particularly for fast CMEs, which are most important in terms of adverse space-weather manifestations. Non-neutralized currents and the length of the main MPIL exhibit significantly stronger correlations than the rest of the properties. This finding supports a causal relationship between coronal mass ejections and non-neutralized electric currents in highly sheared, conspicuous MPILs. In addition, non-neutralized currents and MPIL length carry distinct, independent information as to the eruptive potential of active regions. The combined total amount of non-neutralized electric currents and the length of the main polarity inversion line, therefore, reflect more efficiently than other parameters the eruptive capacity of solar active regions and the CME kinematic characteristics stemming from these regions.
△ Less
Submitted 13 September, 2019;
originally announced September 2019.
-
SOFIA/HAWC+ traces the magnetic fields in NGC 1068
Authors:
E. Lopez-Rodriguez,
C. D. Dowell,
T. J. Jones,
D. A. Harper,
M. Berthoud,
D. Chuss,
D. A. Dale,
J. A. Guerra,
R. T. Hamilton,
L. W. Looney,
J. M. Michail,
R. Nikutta,
G. Novak,
F. P. Santos,
K. Sheth,
J. Siah,
J. Staguhn,
I. W. Stephens,
K. Tassis,
C. Q. Trinh,
D. Ward-Thompson,
M. Werner,
E. J. Wollack,
E. Zweibel
Abstract:
We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $μ$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and…
▽ More
We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $μ$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $μ$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.
△ Less
Submitted 15 November, 2019; v1 submitted 15 July, 2019;
originally announced July 2019.
-
The far-infrared polarization spectrum of Rho Ophiuchi A from HAWC+/SOFIA observations
Authors:
Fabio P. Santos,
David T. Chuss,
C. Darren Dowell,
Martin Houde,
Leslie W. Looney,
Enrique Lopez Rodriguez,
Giles Novak,
Derek Ward-Thompson,
Marc Berthoud,
Daniel A. Dale,
Jordan A. Guerra,
Ryan T. Hamilton,
Shaul Hanany,
Doyal A. Harper,
Thomas K. Henning,
Terry Jay Jones,
Alex Lazarian,
Joseph M. Michail,
Mark R. Morris,
Johannes Staguhn,
Ian W. Stephens,
Konstantinos Tassis,
Christopher Q. Trinh,
Eric Van Camp,
C. G. Volpert
, et al. (1 additional authors not shown)
Abstract:
We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $μ$m) and D (154 $μ$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_…
▽ More
We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $μ$m) and D (154 $μ$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.
△ Less
Submitted 19 July, 2019; v1 submitted 2 May, 2019;
originally announced May 2019.
-
HAWC+/SOFIA Multiwavelength Polarimetric Observations of OMC-1
Authors:
David T. Chuss,
B-G Andersson,
John Bally,
Jessie L. Dotson,
C. Darren Dowell,
Jordan A. Guerra,
Doyal A. Harper,
Martin Houde,
Terry Jay Jones,
A. Lazarian,
Enrique Lopez Rodriguez,
Joseph M. Michail,
Mark R. Morris,
Giles Novak,
Javad Siah,
Johannes Staguhn,
John E. Vaillancourt,
C. G. Volpert,
Michael Werner,
Edward J. Wollack,
Dominic J. Benford,
Marc Berthoud,
Erin G. Cox,
Richard Crutcher,
Daniel A. Dale
, et al. (13 additional authors not shown)
Abstract:
We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometr…
▽ More
We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.
△ Less
Submitted 8 January, 2019; v1 submitted 18 October, 2018;
originally announced October 2018.
-
Photospheric Shear Flows in Solar Active Regions and Their Relation to Flare Occurrence
Authors:
Sung-Hong Park,
Jordan A. Guerra,
Peter T. Gallagher,
Manolis K. Georgoulis,
D. Shaun Bloomfield
Abstract:
Solar active regions (ARs) that produce major flares typically exhibit strong plasma shear flows around photospheric magnetic polarity inversion lines (MPILs). It is therefore important to quantitatively measure such photospheric shear flows in ARs for a better understanding of their relation to flare occurrence. Photospheric flow fields were determined by applying the Differential Affine Velocity…
▽ More
Solar active regions (ARs) that produce major flares typically exhibit strong plasma shear flows around photospheric magnetic polarity inversion lines (MPILs). It is therefore important to quantitatively measure such photospheric shear flows in ARs for a better understanding of their relation to flare occurrence. Photospheric flow fields were determined by applying the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) method to a large data set of 2,548 co-aligned pairs of AR vector magnetograms with 12-min separation over the period 2012-2016. From each AR flow-field map, three shear-flow parameters were derived corresponding to the mean (<S>), maximum (S_max) and integral (S_sum) shear-flow speeds along strong-gradient, strong-field MPIL segments. We calculated flaring rates within 24 hr as a function of each shear-flow parameter, and also investigated the relation between the parameters and the waiting time (τ) until the next major flare (class M1.0 or above) after the parameter observation. In general, it is found that the larger S_sum an AR has, the more likely it is for the AR to produce flares within 24 hr. It is also found that among ARs which produce major flares, if one has a larger value of S_sum then τ generally gets shorter. These results suggest that large ARs with widespread and/or strong shear flows along MPILs tend to not only be more flare productive, but also produce major flares within 24 hr or less.
△ Less
Submitted 20 July, 2018;
originally announced July 2018.
-
Testing and Improving a Set of Morphological Predictors of Flaring Activity
Authors:
Ioannis Kontogiannis,
Manolis K. Georgoulis,
Sung-Hong Park,
Jordan A Guerra
Abstract:
Efficient prediction of solar flares relies on parameters that quantify the eruptive capability of solar active regions. Several such quantitative predictors have been proposed in the literature, inferred mostly from photospheric magnetograms and/or white-light observations. Two of them are the Ising energy and the sum of the total horizontal magnetic field gradient. The former has been developed…
▽ More
Efficient prediction of solar flares relies on parameters that quantify the eruptive capability of solar active regions. Several such quantitative predictors have been proposed in the literature, inferred mostly from photospheric magnetograms and/or white-light observations. Two of them are the Ising energy and the sum of the total horizontal magnetic field gradient. The former has been developed from line-of-sight magnetograms, while the latter uses sunspot detections and characteristics, based on continuum images. Aiming to include these parameters in an automated prediction scheme, we test their applicability on regular photospheric magnetic field observations provided by the Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). We test their efficiency as predictors of flaring activity on a representative sample of active regions and investigate possible modifications of these quantities. The Ising energy appears to be an efficient predictor, and the efficiency is even improved if it is modified to describe interacting magnetic partitions or sunspot umbrae. The sum of the horizontal magnetic field gradient appears to be slightly more promising than the three variations of the Ising energy we implement in this article. The new predictors are also compared with two very promising predictors: the effective connected magnetic field strength and the total unsigned non-neutralized current. Our analysis shows that the efficiency of morphological predictors depends on projection effects in a nontrivial way. All four new predictors are found useful for inclusion in an automated flare forecasting facility, such as the Flare Likelihood and Region Eruption Forecasting (FLARECAST), but their utility, among others, will ultimately be determined by the validation effort underway in the framework of the FLARECAST project.
△ Less
Submitted 17 July, 2018;
originally announced July 2018.
-
Connecting Coronal Mass Ejections to their Solar Active Region Sources: Combining Results from the HELCATS and FLARECAST Projects
Authors:
Sophie A. Murray,
Jordan A. Guerra,
Pietro Zucca,
Sung-Hong Park,
Eoin P. Carley,
Peter T. Gallagher,
Nicole Vilmer,
Volker Bothmer
Abstract:
Coronal mass ejections (CMEs) and other solar eruptive phenomena can be physically linked by combining data from a multitude of ground-based and space-based instruments alongside models, however this can be challenging for automated operational systems. The EU Framework Package 7 HELCATS project provides catalogues of CME observations and properties from the Helio- spheric Imagers onboard the two…
▽ More
Coronal mass ejections (CMEs) and other solar eruptive phenomena can be physically linked by combining data from a multitude of ground-based and space-based instruments alongside models, however this can be challenging for automated operational systems. The EU Framework Package 7 HELCATS project provides catalogues of CME observations and properties from the Helio- spheric Imagers onboard the two NASA/STEREO spacecraft in order to track the evolution of CMEs in the inner heliosphere. From the main HICAT catalogue of over 2,000 CME detections, an automated algorithm has been developed to connect the CMEs observed by STEREO to any corresponding solar flares and active region (AR) sources on the solar surface. CME kinematic properties, such as speed and angular width, are compared with AR magnetic field properties, such as magnetic flux, area, and neutral line characteristics. The resulting LOWCAT catalogue is also compared to the extensive AR property database created by the EU Horizon 2020 FLARECAST project, which provides more complex magnetic field parameters derived from vector magnetograms. Initial statistical analysis has been undertaken on the new data to provide insight into the link between flare and CME events, and characteristics of eruptive ARs. Warning thresholds determined from analysis of the evolution of these parameters is shown to be a useful output for operational space weather purposes. Parameters of particular interest for further analysis include total unsigned flux, vertical current, and current helicity. The automated method developed to create the LOWCAT catalogue may also be useful for future efforts to develop operational CME forecasting.
△ Less
Submitted 17 March, 2018;
originally announced March 2018.
-
Forecasting Solar Flares Using Magnetogram-based Predictors and Machine Learning
Authors:
Kostas Florios,
Ioannis Kontogiannis,
Sung-Hong Park,
Jordan A. Guerra,
Federico Benvenuto,
D. Shaun Bloomfield,
Manolis K. Georgoulis
Abstract:
We propose a forecasting approach for solar flares based on data from Solar Cycle 24, taken by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) mission. In particular, we use the Space-weather HMI Active Region Patches (SHARP) product that facilitates cut-out magnetograms of solar active regions (AR) in the Sun in near-real-time (NRT), taken over a five-year…
▽ More
We propose a forecasting approach for solar flares based on data from Solar Cycle 24, taken by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) mission. In particular, we use the Space-weather HMI Active Region Patches (SHARP) product that facilitates cut-out magnetograms of solar active regions (AR) in the Sun in near-real-time (NRT), taken over a five-year interval (2012 - 2016). Our approach utilizes a set of thirteen predictors, which are not included in the SHARP metadata, extracted from line-of-sight and vector photospheric magnetograms. We exploit several Machine Learning (ML) and Conventional Statistics techniques to predict flares of peak magnitude >M1 and >C1, within a 24 h forecast window. The ML methods used are multi-layer perceptrons (MLP), support vector machines (SVM) and random forests (RF). We conclude that random forests could be the prediction technique of choice for our sample, with the second best method being multi-layer perceptrons, subject to an entropy objective function. A Monte Carlo simulation showed that the best performing method gives accuracy ACC=0.93(0.00), true skill statistic TSS=0.74(0.02) and Heidke skill score HSS=0.49(0.01) for >M1 flare prediction with probability threshold 15% and ACC=0.84(0.00), TSS=0.60(0.01) and HSS=0.59(0.01) for >C1 flare prediction with probability threshold 35%.
△ Less
Submitted 17 January, 2018;
originally announced January 2018.
-
Active Region Photospheric Magnetic Properties Derived from Line-of-sight and Radial Fields
Authors:
Jordan A. Guerra,
Sung-H. Park,
Ioannis Kontogiannis,
Peter T. Gallagher,
Manolis Georgoulis,
D. Shaun Bloomfield
Abstract:
The effect of using two representations of the normal-to-surface magnetic field to calculate photospheric measures that are related to active region (AR) potential for flaring is presented. Several AR properties were computed using line-of-sight ($B_{\rm los}$) and spherical-radial ($B_{r}$) magnetograms from the Space-weather HMI Active Region Patch (SHARP) products of the \emph{Solar Dynamics Ob…
▽ More
The effect of using two representations of the normal-to-surface magnetic field to calculate photospheric measures that are related to active region (AR) potential for flaring is presented. Several AR properties were computed using line-of-sight ($B_{\rm los}$) and spherical-radial ($B_{r}$) magnetograms from the Space-weather HMI Active Region Patch (SHARP) products of the \emph{Solar Dynamics Observatory}, characterizing the presence and features of magnetic polarity inversion lines, fractality, and magnetic connectivity of the AR photospheric field. The data analyzed corresponds to $\approx$ 4,000 AR observations, achieved by randomly selecting 25 % of days between September 2012 and May 2016 for analysis at 6-hr cadence. Results from this statistical study include: i) the $B_{r}$ component results in a slight upwards shift of property values in a manner consistent with a field-strength underestimation by the $B_{\rm los}$ component; ii) using the $B_{r}$ component results in significantly lower inter-property correlation in one-third of the cases, implying more independent information about the state of the AR photospheric magnetic field; iii) flaring rates for each property vary between the field components in a manner consistent with the differences in property-value ranges resulting from the components; iv) flaring rates generally increase for higher values of properties, except Fourier spectral power index that has flare rates peaking around a value of $5/3$. These findings indicate that there may be advantages in using $B_{r}$ rather than $B_{\rm los}$ in calculating flare-related AR magnetic properties, especially for regions located far from central meridian.
△ Less
Submitted 19 December, 2017;
originally announced December 2017.
-
Non-neutralized electric currents in solar active regions and flare productivity
Authors:
Ioannis Kontogiannis,
Manolis K. Georgoulis,
Sung-Hong Park,
Jordan A. Guerra
Abstract:
We explore the association of non-neutralized currents with solar flare occurrence in a sizable sample of observations, aiming to show the potential of such currents in solar flare prediction. We use the regularly produced high quality vector magnetograms by the Helioseismic Magnetic Imager and more specifically the Space weather HMI Active Region Patches (SHARP). Through a newly established metho…
▽ More
We explore the association of non-neutralized currents with solar flare occurrence in a sizable sample of observations, aiming to show the potential of such currents in solar flare prediction. We use the regularly produced high quality vector magnetograms by the Helioseismic Magnetic Imager and more specifically the Space weather HMI Active Region Patches (SHARP). Through a newly established method, that incorporates detailed error analysis, we calculate the non-neutralized currents contained in active regions (AR). Two predictors are produced, namely the total and the maximum unsigned non-neutralized current. Both are tested in AR time-series and a representative sample of point-in-time observations during the interval 2012-2016. The average values of non-neutralized currents in flaring active regions are by more than an order of magnitude higher than in non-flaring regions and correlate very well with the corresponding flare index. The temporal evolution of these parameters appears to be connected to physical processes, such as flux emergence and/or magnetic polarity inversion line formation that are associated with increased solar flare activity. Using Bayesian inference of flaring probabilities, it is shown that the total unsigned non-neutralized current outperforms significantly the total unsigned magnetic flux and other well established current-related predictors. Thus, it shows good prospects for inclusion in an operational flare forecasting service. We plan to use the new predictor in the framework of the FLARECAST project along with other highly performing predictors.
△ Less
Submitted 23 August, 2017;
originally announced August 2017.
-
Ensemble Forecasting of Major Solar Flares -- First Results
Authors:
J. A. Guerra,
A. Pulkkinen,
V. M. Uritsky
Abstract:
We present the results from the first ensemble prediction model for major solar flares (M and X classes). The primary aim of this investigation is to explore the construction of an ensemble for an initial prototyping of this new concept. Using the probabilistic forecasts from three models hosted at the Community Coordinated Modeling Center (NASA-GSFC) and the NOAA forecasts, we developed an ensemb…
▽ More
We present the results from the first ensemble prediction model for major solar flares (M and X classes). The primary aim of this investigation is to explore the construction of an ensemble for an initial prototyping of this new concept. Using the probabilistic forecasts from three models hosted at the Community Coordinated Modeling Center (NASA-GSFC) and the NOAA forecasts, we developed an ensemble forecast by linearly combining the flaring probabilities from all four methods. Performance-based combination weights were calculated using a Monte-Carlo-type algorithm that applies a decision threshold $P_{th}$ to the combined probabilities and maximizing the Heidke Skill Score (HSS). Using the data for 13 recent solar active regions between years 2012 - 2014, we found that linear combination methods can improve the overall probabilistic prediction and improve the categorical prediction for certain values of decision thresholds. Combination weights vary with the applied threshold and none of the tested individual forecasting models seem to provide more accurate predictions than the others for all values of $P_{th}$. According to the maximum values of HSS, a performance-based weights calculated by averaging over the sample, performed similarly to a equally weighted model. The values $P_{th}$ for which the ensemble forecast performs the best are 25 \% for M-class flares and 15 \% for X-class flares. When the human-adjusted probabilities from NOAA are excluded from the ensemble, the ensemble performance in terms of the Heidke score, is reduced.
△ Less
Submitted 4 September, 2015; v1 submitted 5 April, 2015;
originally announced April 2015.
-
Spatio-temporal scaling of turbulent photospheric line-of-sight magnetic field in active region NOAA 11158
Authors:
Jordan A. Guerra,
Antti A. Pulkkinen,
Vadim M. Uritsky,
Seiji Yashiro
Abstract:
We study structure and dynamics of the turbulent photospheric magnetic field in the active region NOAA 11158 by characterizing the spatial and temporal scaling properties of the line-of-sight (LOS) component. Using high-resolution high-cadence LOS magnetograms from SDO/HMI, we measured power-law exponents $α$ and $β$ describing wavenumber- ($k$) and frequency-domain ($f$) Fourier power spectra, re…
▽ More
We study structure and dynamics of the turbulent photospheric magnetic field in the active region NOAA 11158 by characterizing the spatial and temporal scaling properties of the line-of-sight (LOS) component. Using high-resolution high-cadence LOS magnetograms from SDO/HMI, we measured power-law exponents $α$ and $β$ describing wavenumber- ($k$) and frequency-domain ($f$) Fourier power spectra, respectively, and investigated their evolution during the passage of the active region through the field of view of HMI. Flaring active region NOAA 11158 produces an average one-dimensional spatial power spectral density that follows approximately a $k^{-2}$ power law -- a spectrum that suggests parallel MHD fluctuations in an anisotropic turbulent medium. In addition, we found that values of $α$ capture systematic changes in the configuration of LOS photospheric magnetic field during flaring activity in the corona. Position-dependent values of the temporal scaling exponent $β$ showed that, on average, the core of the active region scales with $β>$ 3 surrounded by a diffusive region with an approximately $f^{-2}$-type spectrum. Our results indicate that only about 1 - 3 \% of the studied LOS photospheric magnetic flux displays $β\approxα$, implying that Taylor's hypothesis of frozen-in-flow turbulence is typically invalid for this scalar in presence of turbulent photospheric flows. In consequence, both spatial and temporal variations of the plasma and magnetic field must be included in a complete description of the turbulent evolution of active regions.
△ Less
Submitted 19 November, 2014; v1 submitted 24 February, 2014;
originally announced February 2014.
