-
Constraining binary mergers in AGN disks using the non-observation of lensed gravitational waves
Authors:
Samson H. W. Leong,
Justin Janquart,
Aditya Kumar Sharma,
Paul Martens,
Parameswaran Ajith,
Otto A. Hannuksela
Abstract:
The dense and dynamic environments within active galactic nuclei (AGN) accretion disks may serve as prolific birthplaces for binary black holes (BBHs) and one possible origin for some of the BBHs detected by gravitational-wave (GW) observatories. We show that a considerable fraction of the BBH in AGN disks will be strongly lensed by the central supermassive black hole (SMBH). Thus, the non-observa…
▽ More
The dense and dynamic environments within active galactic nuclei (AGN) accretion disks may serve as prolific birthplaces for binary black holes (BBHs) and one possible origin for some of the BBHs detected by gravitational-wave (GW) observatories. We show that a considerable fraction of the BBH in AGN disks will be strongly lensed by the central supermassive black hole (SMBH). Thus, the non-observation of lensed GW signals can be used to constrain the fraction of BBH binaries residing in AGN disks. The non-detection of lensing with current ${\cal O}(100)$ detections will be sufficient to start placing constraints on the fraction of BBHs living within accretion disks near the SMBH. In the next-generation detectors era, with ${\cal O}(10^5)$ BBH observations and no lensed events, we will be able to rule out most migration traps as dominant birthplaces of BBH mergers; moreover, we will be able to constrain the minimum size of the accretion disk. On the other hand, should AGNs constitute a major formation channel, lensed events from AGNs will become prominent in the future.
△ Less
Submitted 23 August, 2024;
originally announced August 2024.
-
Precise and Accurate Short-term Forecasting of Solar Energetic Particle Events with Multivariate Time Series Classifiers
Authors:
Sumanth A. Rotti,
Berkay Aydin,
Petrus C. Martens
Abstract:
Solar energetic particle (SEP) events are one of the most crucial aspects of space weather that require continuous monitoring and forecasting using robust methods. We demonstrate a proof of concept of using a data-driven supervised classification framework on a multivariate time series data set covering solar cycles 22, 23, and 24. We implement ensemble modeling that merges the results from three…
▽ More
Solar energetic particle (SEP) events are one of the most crucial aspects of space weather that require continuous monitoring and forecasting using robust methods. We demonstrate a proof of concept of using a data-driven supervised classification framework on a multivariate time series data set covering solar cycles 22, 23, and 24. We implement ensemble modeling that merges the results from three proton channels (E$\geq$10 MeV, 50 MeV, and 100 MeV) and the long band X-ray flux (1-8Å) channel from the Geostationary Operational Environmental Satellite missions. Our task is binary classification, such that the aim of the model is to distinguish strong SEP events from nonevents. Here, strong SEP events are those crossing the Space Weather Prediction Center's "S1" threshold of solar radiation storm and proton fluxes below that are weak SEP events. In addition, we consider periods of non-occurrence of SEPs following a flare with magnitudes $\geq$C6.0 to maintain a natural imbalance of sample distribution. In our data set, there are 244 strong SEP events comprising the positive class. There are 189 weak events and 2,460 "SEP-quiet" periods for the negative class. We experiment with summary statistic classifier, one-nearest neighbor and supervised time series forest (STSF), and compare their performances to validate our methods for prediction windows from 5 min up to 60 min. We find STSF to perform better under all circumstances. For an optimal classification threshold of $\approx$0.3 and a 60 min prediction window, we obtain: TSS = 0.850, HSS = 0.878, GSS = 0.783.
△ Less
Submitted 10 August, 2024;
originally announced August 2024.
-
Bispectrum from inflation/bouncing Universe in VCDM
Authors:
Alexander Ganz,
Paul Martens,
Shinji Mukohyama,
Ryo Namba
Abstract:
We discuss the non-linear interactions within the VCDM model, a type II minimally modified gravity model with the same number of degrees of freedom as in General Relativity but not connected to the latter by field redefinitions. During an inflationary phase in the early universe, if the VCDM potential does not modify the slow-roll behavior of the inflaton field, we recover, up to the leading order…
▽ More
We discuss the non-linear interactions within the VCDM model, a type II minimally modified gravity model with the same number of degrees of freedom as in General Relativity but not connected to the latter by field redefinitions. During an inflationary phase in the early universe, if the VCDM potential does not modify the slow-roll behavior of the inflaton field, we recover, up to the leading order, the standard results for the bispectrum in slow-roll inflation. On the other hand, if the VCDM potential becomes dominant, the interactions can strongly deviate and even violate the Maldacena's consistency relation for the local non-Gaussianities. Furthermore, we apply the formalism to the recently introduced bouncing model in VCDM, and show that the bispectrum still respects current observational constraints. Future measurements on non-Gaussianities of the local type should provide the test ground for the model's validity.
△ Less
Submitted 25 July, 2024; v1 submitted 3 July, 2024;
originally announced July 2024.
-
Short-term Classification of Strong Solar Energetic Particle Events using Multivariate Time Series Classifiers
Authors:
Sumanth A. Rotti,
Berkay Aydin,
Petrus C. Martens
Abstract:
Solar energetic particle (SEP) events are one of the most crucial aspects of space weather that require continuous monitoring and forecasting. Their prediction depends on various factors including source eruptions. In the present work, we use the Geostationary Solar Energetic Particle (GSEP) data set covering solar cycles 22, 23, and 24. We develop a framework using time series-based machine learn…
▽ More
Solar energetic particle (SEP) events are one of the most crucial aspects of space weather that require continuous monitoring and forecasting. Their prediction depends on various factors including source eruptions. In the present work, we use the Geostationary Solar Energetic Particle (GSEP) data set covering solar cycles 22, 23, and 24. We develop a framework using time series-based machine learning (ML) models with the aim of developing robust short-term forecasts by classifying SEP events. For this purpose, we introduce an ensemble learning approach that merges the results from univariate time series of three proton channels (10 MeV, 50 MeV, and 100 MeV) and the long band X-ray flux channel from the Geostationary Operational Environmental Satellite (GOES) missions and analyze their performance. We consider three models, namely, time series forest (TSF), supervised time series forest (STSF) and bag of SFA symbols (BOSS). Our study also focuses on understanding and developing confidence in the predictive capabilities of our models. Therefore, we utilize multiple evaluation techniques and metrics. Based on that, we find STSF to perform well in all scenarios. The summary of metrics for the STSF model is as follows: AUC = 0.981; F1-score = 0.960; TSS = 0.919; HSS = 0.920; GSS = 0.852; and MCC = 0.920. The Brier score loss of the STSF model is 0.077. This work lays the foundation for building near-real-time (NRT) short-term SEP event predictions using robust ML methods.
△ Less
Submitted 26 March, 2024;
originally announced March 2024.
-
Magnetic Helicity Signs and Flaring Propensity: Comparing Force-free Parameter with the Helicity signs of Hα Filaments and X-ray Sigmoids
Authors:
V. Aparna,
Manolis K. Georgoulis,
Petrus C. Martens
Abstract:
Sigmoids produce strong eruptive events. Earlier studies have shown that the ICME axial magnetic field Bz can be predicted with some credibility by observing the corresponding filament or the polarity inversion line in the region of eruption and deriving the magnetic field direction from that. Sigmoids are coronal structures often associated with filaments in the sigmoidal region. In this study, f…
▽ More
Sigmoids produce strong eruptive events. Earlier studies have shown that the ICME axial magnetic field Bz can be predicted with some credibility by observing the corresponding filament or the polarity inversion line in the region of eruption and deriving the magnetic field direction from that. Sigmoids are coronal structures often associated with filaments in the sigmoidal region. In this study, firstly we compare filament chirality with sigmoid handedness to observe their correlation. Secondly, we perform non-linear force-free approximations of the coronal magnetic connectivity using photospheric vector magnetograms underneath sigmoids to obtain a weighted-average value of the force-free parameter and to correlate it with filament chirality and the observed coronal sigmoid handedness. Importantly, we find that the sigmoids and their filament counterparts do not always have the same helicity signs. Production of eruptive events by regions that do not have the same signs of helicities is $\sim$3.5 times higher than when they do. A case study of magnetic energy/ helicity evolution in NOAA AR 12473 is also presented.
△ Less
Submitted 25 March, 2024;
originally announced March 2024.
-
Investigating Performance Trends of Simulated Real-time Solar Flare Predictions: The Impacts of Training Windows, Data Volumes, and the Solar Cycle
Authors:
Griffin T. Goodwin,
Viacheslav M. Sadykov,
Petrus C. Martens
Abstract:
This study explores the behavior of machine learning-based flare forecasting models deployed in a simulated operational environment. Using Georgia State University's Space Weather Analytics for Solar Flares benchmark dataset (Angryk et al. 2020a,b), we examine the impacts of training methodology and the solar cycle on decision tree, support vector machine, and multilayer perceptron performance. We…
▽ More
This study explores the behavior of machine learning-based flare forecasting models deployed in a simulated operational environment. Using Georgia State University's Space Weather Analytics for Solar Flares benchmark dataset (Angryk et al. 2020a,b), we examine the impacts of training methodology and the solar cycle on decision tree, support vector machine, and multilayer perceptron performance. We implement our classifiers using three temporal training windows: stationary, rolling, and expanding. The stationary window trains models using a single set of data available before the first forecasting instance, which remains constant throughout the solar cycle. The rolling window trains models using data from a constant time interval before the forecasting instance, which moves with the solar cycle. Finally, the expanding window trains models using all available data before the forecasting instance. For each window, a number of input features (1, 5, 10, 25, 50, 120) and temporal sizes (5, 8, 11, 14, 17, 20 months) were tested. To our surprise, we found that for a 20-month window, skill scores were comparable regardless of the window type, feature count, and classifier selected. Furthermore, reducing the size of this window only marginally decreased stationary and rolling window performance. This implies that, given enough data, a stationary window can be chosen over other window types, eliminating the need for model retraining. Lastly, a moderately strong positive correlation was found to exist between a model's false positive rate and the solar X-ray background flux. This suggests that the solar cycle phase has a considerable influence on forecasting.
△ Less
Submitted 7 February, 2024;
originally announced February 2024.
-
Taming the Beast: Fully Automated Unit Testing with Coyote C++
Authors:
Sanghoon Rho,
Philipp Martens,
Seungcheol Shin,
Yeoneo Kim
Abstract:
In this paper, we present Coyote C++, a fully automated white-box unit testing tool for C and C++. Whereas existing tools have struggled to realize unit test generation for C++, Coyote C++ is able to produce high coverage results from unit test generation at a testing speed of over 10,000 statements per hour. This impressive feat is made possible by the combination of a powerful concolic execution…
▽ More
In this paper, we present Coyote C++, a fully automated white-box unit testing tool for C and C++. Whereas existing tools have struggled to realize unit test generation for C++, Coyote C++ is able to produce high coverage results from unit test generation at a testing speed of over 10,000 statements per hour. This impressive feat is made possible by the combination of a powerful concolic execution engine with sophisticated automated test harness generation. Additionally, the GUI of Coyote C++ displays detailed code coverage visualizations and provides various configuration features for users seeking to manually optimize their coverage results. Combining potent one-click automated testing with rich support for manual tweaking, Coyote C++ is the first automated testing tool that is practical enough to make automated testing of C++ code truly viable in industrial applications.
△ Less
Submitted 4 January, 2024; v1 submitted 2 January, 2024;
originally announced January 2024.
-
Coyote C++: An Industrial-Strength Fully Automated Unit Testing Tool
Authors:
Sanghoon Rho,
Philipp Martens,
Seungcheol Shin,
Yeoneo Kim,
Hoon Heo,
SeungHyun Oh
Abstract:
Coyote C++ is an automated testing tool that uses a sophisticated concolic-execution-based approach to realize fully automated unit testing for C and C++. While concolic testing has proven effective for languages such as C and Java, tools have struggled to achieve a practical level of automation for C++ due to its many syntactical intricacies and overall complexity. Coyote C++ is the first automat…
▽ More
Coyote C++ is an automated testing tool that uses a sophisticated concolic-execution-based approach to realize fully automated unit testing for C and C++. While concolic testing has proven effective for languages such as C and Java, tools have struggled to achieve a practical level of automation for C++ due to its many syntactical intricacies and overall complexity. Coyote C++ is the first automated testing tool to breach the barrier and bring automated unit testing for C++ to a practical level suitable for industrial adoption, consistently reaching around 90% code coverage. Notably, this testing process requires no user involvement and performs test harness generation, test case generation and test execution with "one-click" automation. In this paper, we introduce Coyote C++ by outlining its high-level structure and discussing the core design decisions that shaped the implementation of its concolic execution engine. Finally, we demonstrate that Coyote C++ is capable of achieving high coverage results within a reasonable timespan by presenting the results from experiments on both open-source and industrial software.
△ Less
Submitted 22 October, 2023;
originally announced October 2023.
-
Analysis of SEP events and their possible precursors based on the GSEP Catalog
Authors:
Sumanth A. Rotti,
Petrus C. Martens
Abstract:
Solar energetic particle (SEP) events are one of the most crucial aspects of space weather. Their prediction depends on various factors including the source solar eruptions such as flares and coronal mass ejections (CMEs). The Geostationary Solar Energetic Particle (GSEP) Events catalog was developed as an extensive data set towards this effort for solar cycles 22, 23 and 24. In the present work,…
▽ More
Solar energetic particle (SEP) events are one of the most crucial aspects of space weather. Their prediction depends on various factors including the source solar eruptions such as flares and coronal mass ejections (CMEs). The Geostationary Solar Energetic Particle (GSEP) Events catalog was developed as an extensive data set towards this effort for solar cycles 22, 23 and 24. In the present work, we review and extend the GSEP data set by; (1) adding "weak" SEP events that have proton enhancements from 0.5 to 10 pfu in the E>10 MeV channel, and (2) improving the associated solar source eruptions information. We analyze and discuss spatio-temporal properties such as flare magnitudes, locations, rise times, and speed and width of CMEs. We check for the correlation of these parameters with peak proton fluxes and event fluences. Our study also focuses on understanding feature importance towards the optimal performance of machine learning (ML) models for SEP event forecasting. We implement random forest (RF), extreme gradient boosting (XGBoost), logistic regression (LR) and support vector machines (SVM) classifiers in a binary classification schema. Based on the evaluation of our best models, we find both the flare and CME parameters are requisites to predict the occurrence of an SEP event. This work is a foundation for our further efforts on SEP event forecasting using robust ML methods.
△ Less
Submitted 20 June, 2023;
originally announced June 2023.
-
Spherical scalar collapse in a type-II minimally modified gravity
Authors:
Atabak Fathe Jalali,
Paul Martens,
Shinji Mukohyama
Abstract:
We investigate the spherically-symmetric gravitational collapse of a massless scalar field in the framework of a type-II minimally modified gravity theory called VCDM. This theory propagates only two local physical degrees of freedom supplemented by the so-called instantaneous (or shadowy) mode. Imposing asymptotically flat spacetime in the standard Minkowski time slicing, one can integrate out th…
▽ More
We investigate the spherically-symmetric gravitational collapse of a massless scalar field in the framework of a type-II minimally modified gravity theory called VCDM. This theory propagates only two local physical degrees of freedom supplemented by the so-called instantaneous (or shadowy) mode. Imposing asymptotically flat spacetime in the standard Minkowski time slicing, one can integrate out the instantaneous mode. Consequently, the equations of motion reduce to those in general relativity (GR) with the maximal slicing. Unlike GR, however, VCDM lacks 4D diffeomorphism invariance, and thus one cannot change the time slicing that is preferred by the theory. We then numerically evolve the system to see if and how a black hole forms. For small amplitudes of the initial scalar profile, we find that its collapse does not generate any black hole, singularity or breakdown of the time slicing. For sufficiently large amplitudes, however, the collapse does indeed result in the formation of an apparent horizon in a finite time. After that, the solution outside the horizon is described by a static configuration, i.e. the Schwarzschild geometry with a finite and time-independent lapse function. Inside the horizon, on the other hand, the numerical results indicate that the lapse function keeps decreasing towards zero so that the central singularity is never reached. This implies the necessity for a UV completion of the theory to describe physics inside the horizon. Still, we can conclude that VCDM is able to fully describe the entire time evolution of the Universe outside the black hole horizon without knowledge about such a UV completion.
△ Less
Submitted 5 January, 2024; v1 submitted 18 June, 2023;
originally announced June 2023.
-
Bouncing Cosmology in VCDM
Authors:
Alexander Ganz,
Paul Martens,
Shinji Mukohyama,
Ryo Namba
Abstract:
We construct an asymmetric bouncing scenario within the VCDM model - also known as type-II minimally modified gravity -, a modified gravity theory with two local physical degrees of freedom. The scenario is exempt of any ghost or gradient instability, ad-hoc matching conditions or anisotropic stress issue (BKL instability). It moreover succeeds in generating the cosmological perturbations compatib…
▽ More
We construct an asymmetric bouncing scenario within the VCDM model - also known as type-II minimally modified gravity -, a modified gravity theory with two local physical degrees of freedom. The scenario is exempt of any ghost or gradient instability, ad-hoc matching conditions or anisotropic stress issue (BKL instability). It moreover succeeds in generating the cosmological perturbations compatible with the observations. The scalar spectral index can be adapted by the choice of the equation of state of the matter sector and the form of the VCDM potential leading to an almost scale-invariant power spectrum. Satisfying the CMB bounds on the tensor-to-scalar ratio leads to a blue tensor spectrum.
△ Less
Submitted 27 April, 2023; v1 submitted 27 December, 2022;
originally announced December 2022.
-
Reheating after relaxation of large cosmological constant
Authors:
Paul Martens,
Shinji Mukohyama,
Ryo Namba
Abstract:
We present a cosmological model of an early-time scenario that incorporates a relaxation process of the would-be large vacuum energy, followed by a reheating era connecting to the standard hot big bang universe. Avoiding fine-tuning the cosmological constant is achieved by the dynamics of a scalar field whose kinetic term is modulated by an inverse power of spacetime curvature. While it is at work…
▽ More
We present a cosmological model of an early-time scenario that incorporates a relaxation process of the would-be large vacuum energy, followed by a reheating era connecting to the standard hot big bang universe. Avoiding fine-tuning the cosmological constant is achieved by the dynamics of a scalar field whose kinetic term is modulated by an inverse power of spacetime curvature. While it is at work against radiative corrections to the dark energy, this mechanism alone would wipe out not only the vacuum energy but also all other matter contents. Our present work aims to complete the scenario by exploiting a null-energy-condition violating sector whose energy is eventually transferred to a reheating sector. We provide an explicit example of this process and thus a concrete scenario of the cosmic onset that realizes the thermal history of the Universe with a negligible cosmological constant.
△ Less
Submitted 23 May, 2022;
originally announced May 2022.
-
DeWitt wave function in Hořava-Lifshitz cosmology with tensor perturbation
Authors:
Paul Martens,
Hiroki Matsui,
Shinji Mukohyama
Abstract:
We present a well-tempered DeWitt wave function, which vanishes at the classical big-bang singularity, in Hořava-Lifshitz (HL) cosmology with tensor perturbation, both analytically and numerically. In general relativity, the DeWitt wave function is ill-behaved once the tensor perturbation is taken into account. This is essentially because the amplitude of the perturbation diverges at the singulari…
▽ More
We present a well-tempered DeWitt wave function, which vanishes at the classical big-bang singularity, in Hořava-Lifshitz (HL) cosmology with tensor perturbation, both analytically and numerically. In general relativity, the DeWitt wave function is ill-behaved once the tensor perturbation is taken into account. This is essentially because the amplitude of the perturbation diverges at the singularity and the perturbative expansion completely breaks down. On the other hand, in HL gravity it is known that the higher dimensional operators required by the perturbative renormalizability render the tensor perturbation scale-invariant and regular all the way up to the singularity. In this paper we analytically show that in $d+1$ dimensional HL gravity the DeWitt wave function for tensor perturbation is indeed well-defined around the classical big-bang singularity. Also, we numerically demonstrate the well-behaved DeWitt wave function for tensor perturbation from the big-bang to a finite size of the Universe.
△ Less
Submitted 23 May, 2022;
originally announced May 2022.
-
Stability of the fundamental quasinormal mode in time-domain observations against small perturbations
Authors:
Emanuele Berti,
Vitor Cardoso,
Mark Ho-Yeuk Cheung,
Francesco Di Filippo,
Francisco Duque,
Paul Martens,
Shinji Mukohyama
Abstract:
Black hole spectroscopy with gravitational waves is an important tool to measure the mass and spin of astrophysical black holes and to test their Kerr nature. Next-generation ground- and space-based detectors will observe binary black hole mergers with large signal-to-noise ratios and perform spectroscopy routinely. It was recently shown that small perturbations due, e.g., to environmental effects…
▽ More
Black hole spectroscopy with gravitational waves is an important tool to measure the mass and spin of astrophysical black holes and to test their Kerr nature. Next-generation ground- and space-based detectors will observe binary black hole mergers with large signal-to-noise ratios and perform spectroscopy routinely. It was recently shown that small perturbations due, e.g., to environmental effects (the "flea") to the effective potential governing gravitational-wave generation and propagation in black hole exteriors (the "elephant") can lead to arbitrarily large changes in the black hole's quasinormal spectrum, including the fundamental mode, which is expected to dominate the observed signal. This raises an important question: is the black hole spectroscopy program robust against perturbations? We clarify the physical behavior of time-domain signals under small perturbations in the potential, and we show that changes in the amplitude of the fundamental mode in the prompt ringdown signal are parametrically small. This implies that the fundamental quasinormal mode extracted from the observable time-domain signal is stable against small perturbations. The stability of overtones deserves further investigation.
△ Less
Submitted 24 October, 2022; v1 submitted 17 May, 2022;
originally announced May 2022.
-
Integrated Geostationary Solar Energetic Particle Events Catalog: GSEP
Authors:
Sumanth Rotti,
Berkay Aydin,
Manolis K. Georgoulis,
Petrus C. Martens
Abstract:
We present a catalog of solar energetic particle (SEP) events covering solar cycles 22, 23 and 24. We correlate and integrate three existing catalogs based on Geostationary Operational Environmental Satellite (GOES) integral proton flux data. We visually verified and labeled each event in the catalog to provide a homogenized data set. We have identified a total of 341 SEP events of which 245 cross…
▽ More
We present a catalog of solar energetic particle (SEP) events covering solar cycles 22, 23 and 24. We correlate and integrate three existing catalogs based on Geostationary Operational Environmental Satellite (GOES) integral proton flux data. We visually verified and labeled each event in the catalog to provide a homogenized data set. We have identified a total of 341 SEP events of which 245 cross the space weather prediction center (SWPC) threshold of a significant proton event. The metadata consists of physical parameters and observables concerning the possible source solar eruptions, namely flares and coronal mass ejections for each event. The sliced time series data of each event, along with intensity profiles of proton fluxes in several energy bands, have been made publicly available. This data set enables researchers in machine learning (ML) and statistical analysis to understand the SEPs and the source eruption characteristics useful for space weather prediction.
△ Less
Submitted 8 August, 2022; v1 submitted 25 April, 2022;
originally announced April 2022.
-
The best of both worlds: Using automatic detection and limited human supervision to create a homogenous magnetic catalog spanning four solar cycles
Authors:
A. Munoz-Jaramillo,
Z. A. Werginz,
J. P. Vargas-Acosta,
M. D. DeLuca,
J. C. Windmueller,
J. Zhang,
D. W. Longcope,
D. A. Lamb,
C. E. DeForest,
S. Vargas-Dominguez,
J. W. Harvey,
P. C. H. Martens
Abstract:
Bipolar magnetic regions (BMRs) are the cornerstone of solar variability. They are tracers of the large-scale magnetic processes that give rise to the solar cycle, shapers of the solar corona, building blocks of the large-scale solar magnetic field, and significant contributors to the free-energetic budget that gives rise to flares and coronal mass ejections. Surprisingly, no homogeneous catalog o…
▽ More
Bipolar magnetic regions (BMRs) are the cornerstone of solar variability. They are tracers of the large-scale magnetic processes that give rise to the solar cycle, shapers of the solar corona, building blocks of the large-scale solar magnetic field, and significant contributors to the free-energetic budget that gives rise to flares and coronal mass ejections. Surprisingly, no homogeneous catalog of BMRs exists today, in spite of the existence of systematic measurements of the magnetic field since the early 1970's. The purpose of this work is to address this deficiency by creating a homogenous catalog of BMRs from the 1970's until the present. For this purpose, in this paper we discuss the strengths and weaknesses of the automatic and manual detection of BMRs and how both methods can be combined to form the basis of our Bipolar Active Region Detection (BARD) code and its supporting human supervision module. At present, the BARD catalog contains more than 10,000 unique BMRs tracked and characterized during every day of their observation. Here we also discuss our future plans for the creation of an extended multi-scale magnetic catalog combining the SWAMIS and BARD catalogs.
△ Less
Submitted 22 March, 2022;
originally announced March 2022.
-
Revisiting the Solar Research Cyberinfrastructure Needs: A White Paper of Findings and Recommendations
Authors:
Gelu Nita,
Azim Ahmadzadeh,
Serena Criscuoli,
Alisdair Davey,
Dale Gary,
Manolis Georgoulis,
Neal Hurlburt,
Irina Kitiashvili,
Dustin Kempton,
Alexander Kosovichev,
Piet Martens,
Ryan McGranaghan,
Vincent Oria,
Kevin Reardon,
Viacheslav Sadykov,
Ryan Timmons,
Haimin Wang,
Jason T. L. Wang
Abstract:
Solar and Heliosphere physics are areas of remarkable data-driven discoveries. Recent advances in high-cadence, high-resolution multiwavelength observations, growing amounts of data from realistic modeling, and operational needs for uninterrupted science-quality data coverage generate the demand for a solar metadata standardization and overall healthy data infrastructure. This white paper is prepa…
▽ More
Solar and Heliosphere physics are areas of remarkable data-driven discoveries. Recent advances in high-cadence, high-resolution multiwavelength observations, growing amounts of data from realistic modeling, and operational needs for uninterrupted science-quality data coverage generate the demand for a solar metadata standardization and overall healthy data infrastructure. This white paper is prepared as an effort of the working group "Uniform Semantics and Syntax of Solar Observations and Events" created within the "Towards Integration of Heliophysics Data, Modeling, and Analysis Tools" EarthCube Research Coordination Network (@HDMIEC RCN), with primary objectives to discuss current advances and identify future needs for the solar research cyberinfrastructure. The white paper summarizes presentations and discussions held during the special working group session at the EarthCube Annual Meeting on June 19th, 2020, as well as community contribution gathered during a series of preceding workshops and subsequent RCN working group sessions. The authors provide examples of the current standing of the solar research cyberinfrastructure, and describe the problems related to current data handling approaches. The list of the top-level recommendations agreed by the authors of the current white paper is presented at the beginning of the paper.
△ Less
Submitted 17 March, 2022;
originally announced March 2022.
-
Improved Measurements of the Sun's Meridional Flow and Torsional Oscillation from Correlation tracking on MDI \& HMI magnetograms
Authors:
Sushant S. Mahajan,
David H. Hathaway,
Andrés Muñoz-Jaramillo,
Petrus C. Martens
Abstract:
The Sun's axisymmetric flows, differential rotation and meridional flow, govern the dynamics of the solar magnetic cycle and variety of methods are used to measure these flows, each with its own strengths and weaknesses. Flow measurements based on cross-correlating images of the surface magnetic field have been made since the 1970s which require advanced numerical techniques that are capable of de…
▽ More
The Sun's axisymmetric flows, differential rotation and meridional flow, govern the dynamics of the solar magnetic cycle and variety of methods are used to measure these flows, each with its own strengths and weaknesses. Flow measurements based on cross-correlating images of the surface magnetic field have been made since the 1970s which require advanced numerical techniques that are capable of detecting movements of less than the pixel size in images of the Sun. We have identified several systematic errors in addition to the center-to-limb effect that influence previous measurements of these flows and propose numerical techniques that can minimize these errors by utilizing measurements of displacements at several time-lags. Our analysis of line-of-sight magnetograms from the {\em Michelson Doppler Imager} (MDI) on the ESA/NASA {\em Solar and Heliospheric Observatory} (SOHO) and {\em Helioseismic and Magnetic Imager} (HMI) on the NASA {\em Solar Dynamics Observatory} (SDO) shows long-term variations in the meridional flow and differential rotation over two sunspot cycles from 1996 to 2020. These improved measurements can serve as vital inputs for solar dynamo and surface flux transport simulations.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
Analysis of the Coronal Mass Ejections through Axial Field Direction of Solar Filaments and IMF Bz
Authors:
Kashvi Mundra,
V. Aparna,
Petrus C. H. Martens
Abstract:
In the past, there have been many studies claiming that the effects of geomagnetic storms strongly depends on the orientation of the magnetic-cloud part of the Coronal Mass Ejections (CMEs). Aparna & Martens (2020), using Halo-CME data from 2007-2017, have shown that the magnetic field orientation of filaments at the location where CMEs originate can be effectively used for predicting the onset of…
▽ More
In the past, there have been many studies claiming that the effects of geomagnetic storms strongly depends on the orientation of the magnetic-cloud part of the Coronal Mass Ejections (CMEs). Aparna & Martens (2020), using Halo-CME data from 2007-2017, have shown that the magnetic field orientation of filaments at the location where CMEs originate can be effectively used for predicting the onset of geo-magnetic storms. The purpose of this study is to extend their survey by analyzing the halo-CME data for 1996-2006. The correlation of filament axial direction and their corresponding Bz signatures are used to form a more extensive reasoning for the claims presented by Aparna & Martens before. This study utilizes SOHO EIT 195 Å, MDI magnetogram images, KSO and BBSO H$α$ images for the time period, along with ACE data for inter-planetary magnetic field signatures. Correlating all these, we have found that the trend in Aparna & Martens' study of a high likelihood of the correlation between the axial field direction and Bz orientation, persists for the data between 1996-2006 as well.
△ Less
Submitted 3 November, 2020;
originally announced November 2020.
-
A Catalog of Solar Flare Events Observed by the SOHO/EIT
Authors:
Sumanth A. Rotti,
Petrus C. H. Martens,
Berkay Aydin
Abstract:
We have compiled a catalog of solar flares as observed by the Extreme ultraviolet Imaging Telescope (EIT) aboard the Solar and Heliospheric Observatory (SOHO) spacecraft and the GOES spacecraft over a span from 1997 to 2010. During mid-1998, the cadence of EIT images was revised from two images per day to 12 minutes. However, the low temporal resolution causes significant data gaps in capturing mu…
▽ More
We have compiled a catalog of solar flares as observed by the Extreme ultraviolet Imaging Telescope (EIT) aboard the Solar and Heliospheric Observatory (SOHO) spacecraft and the GOES spacecraft over a span from 1997 to 2010. During mid-1998, the cadence of EIT images was revised from two images per day to 12 minutes. However, the low temporal resolution causes significant data gaps in capturing much of the flaring phenomenon. Therefore, we monitor possible errors in flare detection by flare parameters such as temporal overlap, observational wavelength, and considering full field of view (FOV) images. We consider the GOES flare catalog as the primary source. We describe the technique used to enhance the GOES detected flares using the Extreme Ultraviolet (EUV) image captured by the EIT instrument. In order to detect brightenings, we subtract the images with a maximum cadence of 25 minutes. We have downloaded and analyzed the EIT data via the Virtual Solar Observatory (VSO). This flare dataset from the SOHO/EIT period proves indispensable to the process of the solar flare predictions as the instrument has covered most of Solar Cycle 23.
△ Less
Submitted 10 July, 2020;
originally announced July 2020.
-
Machine Learning in Heliophysics and Space Weather Forecasting: A White Paper of Findings and Recommendations
Authors:
Gelu Nita,
Manolis Georgoulis,
Irina Kitiashvili,
Viacheslav Sadykov,
Enrico Camporeale,
Alexander Kosovichev,
Haimin Wang,
Vincent Oria,
Jason Wang,
Rafal Angryk,
Berkay Aydin,
Azim Ahmadzadeh,
Xiaoli Bai,
Timothy Bastian,
Soukaina Filali Boubrahimi,
Bin Chen,
Alisdair Davey,
Sheldon Fereira,
Gregory Fleishman,
Dale Gary,
Andrew Gerrard,
Gregory Hellbourg,
Katherine Herbert,
Jack Ireland,
Egor Illarionov
, et al. (16 additional authors not shown)
Abstract:
The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, model…
▽ More
The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.
△ Less
Submitted 22 June, 2020;
originally announced June 2020.
-
Roadmap for Reliable Ensemble Forecasting of the Sun-Earth System
Authors:
Gelu Nita,
Rafal Angryk,
Berkay Aydin,
Juan Banda,
Tim Bastian,
Tom Berger,
Veronica Bindi,
Laura Boucheron,
Wenda Cao,
Eric Christian,
Georgia de Nolfo,
Edward DeLuca,
Marc DeRosa,
Cooper Downs,
Gregory Fleishman,
Olac Fuentes,
Dale Gary,
Frank Hill,
Todd Hoeksema,
Qiang Hu,
Raluca Ilie,
Jack Ireland,
Farzad Kamalabadi,
Kelly Korreck,
Alexander Kosovichev
, et al. (22 additional authors not shown)
Abstract:
The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the und…
▽ More
The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.
△ Less
Submitted 29 October, 2018; v1 submitted 19 October, 2018;
originally announced October 2018.
-
The Extended Solar Cycle: Muddying the Waters of Solar/Stellar Dynamo Modeling Or Providing Crucial Observational Constraints?
Authors:
A. K. Srivastava,
S. W. McIntosh,
N. Arge,
D. Banerjee,
E. Cliver,
M. Dikpati,
B. N. Dwivedi,
M. Guhathakurta,
B. B. Karak,
R. J. Leamon,
P. Martens,
S. K. Matthew,
A. Munoz-Jaramillo,
D. Nandi,
A. Norton,
L. Upton,
S. Chatterjee,
R. Mazumder,
Yamini K. Rao,
R. Yadav
Abstract:
In 1844 Schwabe discovered that the number of sunspots increased and decreased over a period of about 11 years, that variation became known as the sunspot cycle. Almost eighty years later, Hale described the nature of the Sun's magnetic field, identifying that it takes about 22 years for the Sun's magnetic polarity to cycle. It was also identified that the latitudinal distribution of sunspots rese…
▽ More
In 1844 Schwabe discovered that the number of sunspots increased and decreased over a period of about 11 years, that variation became known as the sunspot cycle. Almost eighty years later, Hale described the nature of the Sun's magnetic field, identifying that it takes about 22 years for the Sun's magnetic polarity to cycle. It was also identified that the latitudinal distribution of sunspots resembles the wings of a butterfly showing migration of sunspots in each hemisphere that abruptly start at mid-latitudes towards the Sun's equator over the next 11 years. These sunspot patterns were shown to be asymmetric across the equator. In intervening years, it was deduced that the Sun (and sun-like stars) possess magnetic activity cycles that are assumed to be the physical manifestation of a dynamo process that results from complex circulatory transport processes in the star's interior. Understanding the Sun's magnetism, its origin and its variation, has become a fundamental scientific objective \-- the distribution of magnetism, and its interaction with convective processes, drives various plasma processes in the outer atmosphere. In the past few decades, a range of diagnostic techniques have been employed to systematically study finer scale magnetized objects, and associated phenomena. The patterns discerned became known as the ``Extended Solar Cycle'' (ESC). The patterns of the ESC appeared to extend the wings of the activity butterfly back in time, nearly a decade before the formation of the sunspot pattern, and to much higher solar latitudes. In this short review, we describe their observational patterns of the ESC and discuss possible connections to the solar dynamo as we depart on a multi-national collaboration to investigate the origins of solar magnetism through a blend of archived and contemporary data analysis with the goal of improving solar dynamo understanding and modeling.
△ Less
Submitted 19 July, 2018;
originally announced July 2018.
-
On the Prediction of >100 MeV Solar Energetic Particle Events Using GOES Satellite Data
Authors:
Soukaina Filali Boubrahimi,
Berkay Aydin,
Petrus Martens,
Rafal Angryk
Abstract:
Solar energetic particles are a result of intense solar events such as solar flares and Coronal Mass Ejections (CMEs). These latter events all together can cause major disruptions to spacecraft that are in Earth's orbit and outside of the magnetosphere. In this work we are interested in establishing the necessary conditions for a major geo-effective solar particle storm immediately after a major f…
▽ More
Solar energetic particles are a result of intense solar events such as solar flares and Coronal Mass Ejections (CMEs). These latter events all together can cause major disruptions to spacecraft that are in Earth's orbit and outside of the magnetosphere. In this work we are interested in establishing the necessary conditions for a major geo-effective solar particle storm immediately after a major flare, namely the existence of a direct magnetic connection. To our knowledge, this is the first work that explores not only the correlations of GOES X-ray and proton channels, but also the correlations that happen across all the proton channels. We found that proton channels auto-correlations and cross-correlations may also be precursors to the occurrence of an SEP event. In this paper, we tackle the problem of predicting >100 MeV SEP events from a multivariate time series perspective using easily interpretable decision tree models.
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
Data Handling and Assimilation for Solar Event Prediction
Authors:
Petrus C. Martens,
Rafal A. Angryk
Abstract:
The prediction of solar flares, eruptions, and high energy particle storms is of great societal importance. The data mining approach to forecasting has been shown to be very promising. Benchmark datasets are a key element in the further development of data-driven forecasting. With one or more benchmark data sets established, judicious use of both the data themselves and the selection of prediction…
▽ More
The prediction of solar flares, eruptions, and high energy particle storms is of great societal importance. The data mining approach to forecasting has been shown to be very promising. Benchmark datasets are a key element in the further development of data-driven forecasting. With one or more benchmark data sets established, judicious use of both the data themselves and the selection of prediction algorithms is key to developing a high quality and robust method for the prediction of geo-effective solar activity. We review here briefly the process of generating benchmark datasets and developing prediction algorithms.
△ Less
Submitted 4 December, 2017;
originally announced December 2017.
-
Hemispheric Preference and Cyclic Variation of Solar Filament Chirality from 2000 to 2016
Authors:
Soumitra Hazra,
Sushant S. Mahajan,
William Keith Douglas Jr.,
Petrus C. H. Martens
Abstract:
It is well known that solar filaments are features in the solar atmosphere which show a hemispheric preference in their chirality. The hemispheric preference is such that the dextral chirality dominates in the northern hemisphere while the sinistral chirality dominates in the southern. Determining the strength and cyclic variation of the degree of this hemispheric preference, however, is challengi…
▽ More
It is well known that solar filaments are features in the solar atmosphere which show a hemispheric preference in their chirality. The hemispheric preference is such that the dextral chirality dominates in the northern hemisphere while the sinistral chirality dominates in the southern. Determining the strength and cyclic variation of the degree of this hemispheric preference, however, is challenging and tedious and thus needs to be automated. In this paper, we follow Dr. Pietro Bernasconi's algorithm (Bernasconi et al. 2005) to detect filament chirality with two parallel channels of application. The algorithm is applied to H-alpha images by "Advanced Automated Filament Detection and Characterization Code" (AAFDCC) (Bernasconi et al. 2005) or the algorithm is explained to a human and the human determines the chirality of the solar filament. We have conducted this exercise on data during the month of August from years 2000 to 2016 and we found that 83 % of our visually determined filaments follow the hemispheric chirality preference, while 58 % of automatically determined filaments follow it. Our visually compiled results have over 90 % agreement with those of Pevtsov et al. (2003), yet the visually determined chiralities of filaments disagree with automated determinations significantly. We find that the hemispheric preference remained the same between solar cycles 23 and 24 but the preference is very difficult to determine during the solar minimum of 2008-2010 primarily due to the absence of filaments.
△ Less
Submitted 5 October, 2018; v1 submitted 15 November, 2017;
originally announced November 2017.
-
The Faint Young Sun and Faint Young Stars Paradox
Authors:
Petrus C. Martens
Abstract:
The purpose of this paper is to explore a resolution for the Faint Young Sun Paradox that has been mostly rejected by the community, namely the possibility of a somewhat more massive young Sun with a large mass loss rate sustained for two to three billion years. This would make the young Sun bright enough to keep both the terrestrial and Martian oceans from freezing, and thus resolve the paradox.…
▽ More
The purpose of this paper is to explore a resolution for the Faint Young Sun Paradox that has been mostly rejected by the community, namely the possibility of a somewhat more massive young Sun with a large mass loss rate sustained for two to three billion years. This would make the young Sun bright enough to keep both the terrestrial and Martian oceans from freezing, and thus resolve the paradox. It is found that a large and sustained mass loss is consistent with the well observed spin-down rate of Sun-like stars, and indeed may be required for it. It is concluded that a more massive young Sun must be considered a plausible hypothesis.
△ Less
Submitted 3 June, 2017;
originally announced June 2017.
-
Optimal undulating swimming for a single fish-like body and for a pair of interacting swimmers
Authors:
Audrey P. Maertens,
Amy Gao,
Michael S. Triantafyllou
Abstract:
We establish through numerical simulation conditions for optimal undulatory propulsion for a single fish, and for a pair of hydrodynamically interacting fish, accounting for linear and angular recoil. We first employ systematic 2D simulations to identify conditions for minimal propulsive power of a self-propelled fish, and continue with targeted 3D simulations for a danio-like fish. We find that t…
▽ More
We establish through numerical simulation conditions for optimal undulatory propulsion for a single fish, and for a pair of hydrodynamically interacting fish, accounting for linear and angular recoil. We first employ systematic 2D simulations to identify conditions for minimal propulsive power of a self-propelled fish, and continue with targeted 3D simulations for a danio-like fish. We find that the Strouhal number, phase angle between heave and pitch at the trailing edge, and angle of attack are principal parameters. Angular recoil has significant impact on efficiency, while optimized body bending requires maximum bending amplitude upstream of the trailing edge. For 2D simulations, imposing a deformation based on measured displacement for carangiform swimming provides efficiency of 40%, which increases for an optimized profile to 57%; for a 3D fish, the corresponding increase is from 22% to 35%; all at Reynolds number 5000.
Next, we turn to 2D simulation of two hydrodynamically interacting fish. We find that the upstream fish benefits energetically only for small distances. In contrast, the downstream fish can benefit at any position that allows interaction with the upstream wake, provided its body motion is timed appropriately with respect to the oncoming vortices. For an in-line configuration, one body length apart, the optimal efficiency of the downstream fish can increase to 66%; for an offset arrangement it can reach 81%. This proves that in groups of fish, energy savings can be achieved for downstream fish through interaction with oncoming vortices, even when the downstream fish lies directly inside the jet-like flow of an upstream fish.
△ Less
Submitted 4 April, 2016;
originally announced April 2016.
-
The Minimum of Solar Cycle 23: As Deep as It Could Be?
Authors:
Andrés Muñoz-Jaramillo,
Ryan R. Senkpeil,
Dana W. Longcope,
Andrey G. Tlatov,
Alexei A. Pevtsov,
Laura A. Balmaceda,
Edward E. DeLuca,
Petrus C. H. Martens
Abstract:
In this work we introduce a new way of binning sunspot group data with the purpose of better understanding the impact of the solar cycle on sunspot properties and how this defined the characteristics of the extended minimum of cycle 23. Our approach assumes that the statistical properties of sunspots are completely determined by the strength of the underlying large-scale field and have no addition…
▽ More
In this work we introduce a new way of binning sunspot group data with the purpose of better understanding the impact of the solar cycle on sunspot properties and how this defined the characteristics of the extended minimum of cycle 23. Our approach assumes that the statistical properties of sunspots are completely determined by the strength of the underlying large-scale field and have no additional time dependencies. We use the amplitude of the cycle at any given moment (something we refer to as activity level) as a proxy for the strength of this deep-seated magnetic field.
We find that the sunspot size distribution is composed of two populations: one population of groups and active regions and a second population of pores and ephemeral regions. When fits are performed at periods of different activity level, only the statistical properties of the former population, the active regions, is found to vary.
Finally, we study the relative contribution of each component (small-scale versus large-scale) to solar magnetism. We find that when hemispheres are treated separately, almost every one of the past 12 solar minima reaches a point where the main contribution to magnetism comes from the small-scale component. However, due to asymmetries in cycle phase, this state is very rarely reached by both hemispheres at the same time. From this we infer that even though each hemisphere did reach the magnetic baseline, from a heliospheric point of view the minimum of cycle 23 was not as deep as it could have been.
△ Less
Submitted 5 August, 2015;
originally announced August 2015.
-
Small-Scale and Global Dynamos and the Area and Flux Distributions of Active Regions, Sunspot Groups, and Sunspots: A Multi-Database Study
Authors:
Andrés Muñoz-Jaramillo,
Ryan R. Senkpeil,
John C. Windmueller,
Ernest C. Amouzou,
Dana W. Longcope,
Andrey G. Tlatov,
Yury A. Nagovitsyn,
Alexei A. Pevtsov,
Gary A. Chapman,
Angela M. Cookson,
Anthony R. Yeates,
Fraser T. Watson,
Laura A. Balmaceda,
Edward E. DeLuca,
Petrus C. H. Martens
Abstract:
In this work we take advantage of eleven different sunspot group, sunspot, and active region databases to characterize the area and flux distributions of photospheric magnetic structures. We find that, when taken separately, different databases are better fitted by different distributions (as has been reported previously in the literature). However, we find that all our databases can be reconciled…
▽ More
In this work we take advantage of eleven different sunspot group, sunspot, and active region databases to characterize the area and flux distributions of photospheric magnetic structures. We find that, when taken separately, different databases are better fitted by different distributions (as has been reported previously in the literature). However, we find that all our databases can be reconciled by the simple application of a proportionality constant, and that, in reality, different databases are sampling different parts of a composite distribution. This composite distribution is made up by linear combination of Weibull and log-normal distributions -- where a pure Weibull (log-normal) characterizes the distribution of structures with fluxes below (above) $10^{21}$Mx ($10^{22}$Mx). We propose that this is evidence of two separate mechanisms giving rise to visible structures on the photosphere: one directly connected to the global component of the dynamo (and the generation of bipolar active regions), and the other with the small-scale component of the dynamo (and the fragmentation of magnetic structures due to their interaction with turbulent convection). Additionally, we demonstrate that the Weibull distribution shows the expected linear behavior of a power-law distribution (when extended into smaller fluxes), making our results compatible with the results of Parnell et al. (2009).
△ Less
Submitted 10 February, 2015; v1 submitted 23 October, 2014;
originally announced October 2014.
-
Efficiency of Fish Propulsion
Authors:
A. P. Maertens,
M. S. Triantafyllou,
D. K. P. Yue
Abstract:
It is shown that the system efficiency of a self-propelled flexible body is ill-defined unless one considers the concept of quasi-propulsive efficiency, defined as the ratio of the power needed to tow a body in rigid-straight condition over the power it needs for self-propulsion, both measured for the same speed. Through examples we show that the quasi-propulsive efficiency is the only rational no…
▽ More
It is shown that the system efficiency of a self-propelled flexible body is ill-defined unless one considers the concept of quasi-propulsive efficiency, defined as the ratio of the power needed to tow a body in rigid-straight condition over the power it needs for self-propulsion, both measured for the same speed. Through examples we show that the quasi-propulsive efficiency is the only rational non-dimensional metric of the propulsive fitness of fish and fish-like mechanisms. Using two-dimensional viscous simulations and the concept of quasi-propulsive efficiency, we discuss the efficiency two-dimensional undulating foils. We show that low efficiencies, due to adverse body-propulsor hydrodynamic interactions, cannot be accounted for by the increase in friction drag.
△ Less
Submitted 25 September, 2014;
originally announced September 2014.
-
The unusual minimum of sunspot cycle 23 a consequence of Sun's meridional plasma flow variations
Authors:
Dibyendu Nandy,
Andrés Muñoz-Jaramillo,
Petrus C. H. Martens
Abstract:
Direct observations over the past four centuries show that the number of sunspots observed on the Sun's surface vary periodically, going through successive maxima and minima. Following sunspot cycle 23, the Sun went into a prolonged minimum characterized by a very weak polar magnetic field and an unusually large number of days without sunspots. Sunspots are strongly magnetized regions and are gene…
▽ More
Direct observations over the past four centuries show that the number of sunspots observed on the Sun's surface vary periodically, going through successive maxima and minima. Following sunspot cycle 23, the Sun went into a prolonged minimum characterized by a very weak polar magnetic field and an unusually large number of days without sunspots. Sunspots are strongly magnetized regions and are generated by a dynamo mechanism which recreates the solar polar field mediated via plasma flows. Here we report results from kinematic dynamo simulations which demonstrate that a fast meridional flow in the early half of a cycle, followed by a slower flow in the latter half, reproduces both the characteristics of the minimum of sunspot cycle 23 - a large number of spotless days and a relatively weak polar field. Our model predicts that, in general, very deep minima are associated with weak polar fields. Sunspots govern the solar radiative energy and radio flux, and in conjunction with the polar field, modulate the solar wind, heliospheric open flux and consequently cosmic ray flux at Earth; our results therefore provide a causal link between solar internal dynamics and the atypical values of these heliospheric forcing parameters during the recently concluded solar minimum.
△ Less
Submitted 1 March, 2013;
originally announced March 2013.
-
Fast Differential Emission Measure Inversion of Solar Coronal Data
Authors:
Joseph Plowman,
Charles Kankelborg,
Petrus Martens
Abstract:
We present a fast method for reconstructing Differential Emission Measures (DEMs) using solar coronal data. On average, the method computes over 1000 DEMs per second for a sample active region observed by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO), and achieves reduced chi-squared of order unity with no negative emission in all but a few test cases. The high per…
▽ More
We present a fast method for reconstructing Differential Emission Measures (DEMs) using solar coronal data. On average, the method computes over 1000 DEMs per second for a sample active region observed by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO), and achieves reduced chi-squared of order unity with no negative emission in all but a few test cases. The high performance of this method is especially relevant in the context of AIA, which images of order one million solar pixels per second. This paper describes the method, analyzes its fidelity, compares its performance and results with other DEM methods, and applies it to an active region and loop observed by AIA and by the Extreme-ultraviolet Imaging Spectrometer (EIS) on Hinode.
△ Less
Submitted 22 October, 2012; v1 submitted 27 April, 2012;
originally announced April 2012.
-
EUV Analysis of a Quasi-Static Coronal Loop Structure
Authors:
Jason Scott,
P. C. H. Martens,
David McKenzie
Abstract:
Decaying active region 10942 is investigated from 4:00-16:00 UT on February 24, 2007 using a suite of EUV observing instruments. Results from Hinode/EIS, STEREO and TRACE show that although the active region has decayed and no sunspot is present, the physical mechanisms that produce distinguishable loop structures, spectral line broadening, and plasma flows still occur. A coronal loop that appears…
▽ More
Decaying active region 10942 is investigated from 4:00-16:00 UT on February 24, 2007 using a suite of EUV observing instruments. Results from Hinode/EIS, STEREO and TRACE show that although the active region has decayed and no sunspot is present, the physical mechanisms that produce distinguishable loop structures, spectral line broadening, and plasma flows still occur. A coronal loop that appears as a blue-shifted structure in Doppler maps is apparent in intensity images of log(T) = 6.0-6.3 ions. The loop structure is found to be anti-correlated with spectral line broadening generally attributed to nonthermal velocities. This coronal loop structure is investigated physically (temperature, density, geometry) and temporally. Lightcurves created from imaging instruments show brightening and dimming of the loop structure on two different time scales; short pulses of 10-20 min and long duration dimming of 2-4 hours until its disappearance. The coronal loop structure, formed from relatively blue-shifted material that is anti-correlated with spectral line broadening, shows a density of 10^10 to 10^9.3 cm-3 and is visible for longer than characteristic cooling times. The maximum nonthermal spectral line broadenings are found to be adjacent to the footpoint of the coronal loop structure.
△ Less
Submitted 25 October, 2011;
originally announced October 2011.
-
Magnetic Quenching of Turbulent Diffusivity: Reconciling Mixing-length Theory Estimates with Kinematic Dynamo Models of the Solar Cycle
Authors:
Andrés Muñoz-Jaramillo,
Dibyendu Nandy,
Petrus C. H. Martens
Abstract:
The turbulent magnetic diffusivity in the solar convection zone is one of the most poorly constrained ingredients of mean-field dynamo models. This lack of constraint has previously led to controversy regarding the most appropriate set of parameters, as different assumptions on the value of turbulent diffusivity lead to radically different solar cycle predictions. Typically, the dynamo community u…
▽ More
The turbulent magnetic diffusivity in the solar convection zone is one of the most poorly constrained ingredients of mean-field dynamo models. This lack of constraint has previously led to controversy regarding the most appropriate set of parameters, as different assumptions on the value of turbulent diffusivity lead to radically different solar cycle predictions. Typically, the dynamo community uses double step diffusivity profiles characterized by low values of diffusivity in the bulk of the convection zone. However, these low diffusivity values are not consistent with theoretical estimates based on mixing-length theory -- which suggest much higher values for turbulent diffusivity. To make matters worse, kinematic dynamo simulations cannot yield sustainable magnetic cycles using these theoretical estimates. In this work we show that magnetic cycles become viable if we combine the theoretically estimated diffusivity profile with magnetic quenching of the diffusivity. Furthermore, we find that the main features of this solution can be reproduced by a dynamo simulation using a prescribed (kinematic) diffusivity profile that is based on the spatiotemporal geometric-average of the dynamically quenched diffusivity. Here, we provide an analytic fit to the dynamically quenched diffusivity profile, which can be used in kinematic dynamo simulations. Having successfully reconciled the mixing-length theory estimated diffusivity profile with kinematic dynamo models, we argue that they remain a viable tool for understanding the solar magnetic cycle.
△ Less
Submitted 7 July, 2010;
originally announced July 2010.
-
Empirical Modeling of Radiative versus Magnetic Flux for the Sun-as-a-Star
Authors:
Dora Preminger,
Dibyendu Nandy,
Gary Chapman,
Petrus C. H. Martens
Abstract:
We study the relationship between full-disk solar radiative flux at different wavelengths and average solar photospheric magnetic-flux density, using daily measurements from the Kitt Peak magnetograph and other instruments extending over one or more solar cycles. We use two different statistical methods to determine the underlying nature of these flux-flux relationships. First, we use statistical…
▽ More
We study the relationship between full-disk solar radiative flux at different wavelengths and average solar photospheric magnetic-flux density, using daily measurements from the Kitt Peak magnetograph and other instruments extending over one or more solar cycles. We use two different statistical methods to determine the underlying nature of these flux-flux relationships. First, we use statistical correlation and regression analysis and show that the relationships are not monotonic for total solar irradiance and for continuum radiation from the photosphere, but are approximately linear for chromospheric and coronal radiation. Second, we use signal theory to examine the flux-flux relationships for a temporal component. We find that a well-defined temporal component exists and accounts for some of the variance in the data. This temporal component arises because active regions with high magnetic field strength evolve, breaking up into small-scale magnetic elements with low field strength, and radiative and magnetic fluxes are sensitive to different active-region components. We generate empirical models that relate radiative flux to magnetic flux, allowing us to predict spectral-irradiance variations from observations of disk-averaged magnetic-flux density. In most cases, the model reconstructions can account for 85-90% of the variability of the radiative flux from the chromosphere and corona. Our results are important for understanding the relationship between magnetic and radiative measures of solar and stellar variability.
△ Less
Submitted 22 June, 2010;
originally announced June 2010.
-
A Double-Ring Algorithm for Modeling Solar Active Regions: Unifying Kinematic Dynamo Models and Surface Flux-Transport Simulations
Authors:
Andrés Muñoz-Jaramillo,
Dibyendu Nandy,
Petrus C. H. Martens,
Anthony R. Yeates
Abstract:
The emergence of tilted bipolar active regions and the dispersal of their flux, mediated via processes such as diffusion, differential rotation and meridional circulation is believed to be responsible for the reversal of the Sun's polar field. This process (commonly known as the Babcock-Leighton mechanism) is usually modeled as a near-surface, spatially distributed $α$-effect in kinematic mean-fie…
▽ More
The emergence of tilted bipolar active regions and the dispersal of their flux, mediated via processes such as diffusion, differential rotation and meridional circulation is believed to be responsible for the reversal of the Sun's polar field. This process (commonly known as the Babcock-Leighton mechanism) is usually modeled as a near-surface, spatially distributed $α$-effect in kinematic mean-field dynamo models. However, this formulation leads to a relationship between polar field strength and meridional flow speed which is opposite to that suggested by physical insight and predicted by surface flux-transport simulations. With this in mind, we present an improved double-ring algorithm for modeling the Babcock-Leighton mechanism based on active region eruption, within the framework of an axisymmetric dynamo model. Using surface flux-transport simulations we first show that an axisymmetric formulation -- which is usually invoked in kinematic dynamo models -- can reasonably approximate the surface flux dynamics. Finally, we demonstrate that our treatment of the Babcock-Leighton mechanism through double-ring eruption leads to an inverse relationship between polar field strength and meridional flow speed as expected, reconciling the discrepancy between surface flux-transport simulations and kinematic dynamo models.
△ Less
Submitted 10 August, 2010; v1 submitted 22 June, 2010;
originally announced June 2010.
-
Solar Cycle Variation of Magnetic Flux Ropes in a Quasi-Static Coronal Evolution Model
Authors:
A. R. Yeates,
J. A. Constable,
P. C. H. Martens
Abstract:
The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free "potential field" extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helic…
▽ More
The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free "potential field" extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helicity due to flux emergence and large-scale photospheric motions. This helicity concentrates into twisted magnetic flux ropes, which may lose equilibrium and be ejected. Here, we consider how the magnetic structure predicted by this model-in particular the flux ropes-varies over the solar activity cycle, based on photospheric input data from six periods of cycle 23. The number of flux ropes doubles from minimum to maximum, following the total length of photospheric polarity inversion lines. However, the number of flux rope ejections increases by a factor of eight, following the emergence rate of active regions. This is broadly consistent with the observed cycle modulation of coronal mass ejections, although the actual rate of ejections in the simulation is about a fifth of the rate of observed events. The model predicts that, even at minimum, differential rotation will produce sheared, non-potential, magnetic structure at all latitudes.
△ Less
Submitted 24 March, 2010;
originally announced March 2010.
-
Comparison of a Global Magnetic Evolution Model with Observations of Coronal Mass Ejections
Authors:
A. R. Yeates,
G. D. R. Attrill,
Dibyendu Nandy,
D. H. Mackay,
P. C. H. Martens,
A. A. van Ballegooijen
Abstract:
The relative importance of different initiation mechanisms for coronal mass ejections (CMEs) on the Sun is uncertain. One possible mechanism is the loss of equilibrium of coronal magnetic flux ropes formed gradually by large-scale surface motions. In this paper, the locations of flux rope ejections in a recently-developed quasi-static global evolution model are compared with observed CME source…
▽ More
The relative importance of different initiation mechanisms for coronal mass ejections (CMEs) on the Sun is uncertain. One possible mechanism is the loss of equilibrium of coronal magnetic flux ropes formed gradually by large-scale surface motions. In this paper, the locations of flux rope ejections in a recently-developed quasi-static global evolution model are compared with observed CME source locations over a 4.5-month period in 1999. Using EUV data, the low-coronal source locations are determined unambiguously for 98 out of 330 CMEs. Despite the incomplete observations, positive correlation (with coefficient up to 0.49) is found between the distributions of observed and simulated ejections, but only when binned into periods of one month or longer. This binning timescale corresponds to the time interval at which magnetogram data are assimilated into the coronal simulations, and the correlation arises primarily from the large-scale surface magnetic field distribution; only a weak dependence is found on the magnetic helicity imparted to the emerging active regions. The simulations are limited in two main ways: they produce fewer ejections, and they do not reproduce the strong clustering of observed CME sources into active regions. Due to this clustering, the horizontal gradient of radial photospheric magnetic field is better correlated with the observed CME source distribution (coefficient 0.67). Our results suggest that, while the gradual formation of magnetic flux ropes over weeks can account for many observed CMEs, especially at higher latitudes, there exists a second class of CMEs (at least half) for which dynamic active region flux emergence on shorter timescales must be the dominant factor.
△ Less
Submitted 17 December, 2009;
originally announced December 2009.
-
Helioseismic data inclusion in solar dynamo models
Authors:
Andrés Muñoz-Jaramillo,
Dibyendu Nandy,
Petrus C. H. Martens
Abstract:
An essential ingredient in kinematic dynamo models is the velocity field within the solar convection zone. In particular, the differential rotation is now well constrained by helioseismic observations. Helioseismology also gives us information about the depth-dependence of the meridional circulation in the near-surface layers. The typical velocity inputs used in solar dynamo models, however, con…
▽ More
An essential ingredient in kinematic dynamo models is the velocity field within the solar convection zone. In particular, the differential rotation is now well constrained by helioseismic observations. Helioseismology also gives us information about the depth-dependence of the meridional circulation in the near-surface layers. The typical velocity inputs used in solar dynamo models, however, continue to be an analytic fit to the observed differential rotation and a theoretically constructed meridional flow profile that matches only the peak flow speed at the surface. Here we take the first steps towards realistic helioseismic data assimilation, by presenting methodologies for constructing differential rotation and meridional circulation profiles that more closely conform to the observational constraints currently available. We also present simulations driven by the assimilated rotation and four plausible profiles for the internal meridional circulation -- all of which match the helioseismically inferred near-surface depth-dependence, but whose magnitudes are made to vary. We discuss how the results compare with those that are driven by purely analytic fits. Our results indicate that the latitudinal shear of the rotation in the bulk of the solar convection zone plays a more important role, than either the tachocline or surface radial shear, in the induction of toroidal field. We also find that it is the speed of the equatorward counter-flow in the meridional flow at the base of the convection zone, and not how far into the radiative interior it penetrates, that primarily determines the dynamo cycle period. Given that improved helioseismic constraints are expected to be available in the future, our analysis lays the basis for assimilating these data within dynamo models.
△ Less
Submitted 29 June, 2009; v1 submitted 20 November, 2008;
originally announced November 2008.
-
Solar Coronal Structures and Stray Light in TRACE
Authors:
C. E. DeForest,
P. C. H. Martens,
M. J. Wills-Davey
Abstract:
Using the 2004 Venus transit of the Sun to constrain a semi-empirical point-spread function for the TRACE EUV solar telescope, we have measured the effect of stray light in that telescope. We find that 43% of 171A EUV light that enters TRACE is scattered, either through diffraction off the entrance filter grid or through other nonspecular effects. We carry this result forward, via known-PSF deco…
▽ More
Using the 2004 Venus transit of the Sun to constrain a semi-empirical point-spread function for the TRACE EUV solar telescope, we have measured the effect of stray light in that telescope. We find that 43% of 171A EUV light that enters TRACE is scattered, either through diffraction off the entrance filter grid or through other nonspecular effects. We carry this result forward, via known-PSF deconvolution of TRACE images, to identify its effect on analysis of TRACE data. Known-PSF deconvolution by this derived PSF greatly reduces the effect of visible haze in the TRACE 171A images, enhances bright features, and reveals that the smooth background component of the corona is considerably less bright (and hence much more rarefied) than commonly supposed. Deconvolution reveals that some prior conlclusions about the Sun appear to have been based on stray light in the images. In particular, the diffuse background "quiet corona" becomes consistent with hydrostatic support of the coronal plasma; feature contrast is greatly increased, possibly affecting derived parameters such as the form of the coronal heating function; and essentially all existing differential emission measure studies of small features appear to be affected by contamination from nearby features. We speculate on further implications of stray light for interpretation of EUV images from TRACE and similar instruments, and advocate deconvolution as a standard tool for image analysis with future instruments such as SDO/AIA.
△ Less
Submitted 29 August, 2008;
originally announced August 2008.
-
Scaling Laws and Temperature Profiles for Solar and Stellar Coronal Loops with Non-uniform Heating
Authors:
P. C. H. Martens
Abstract:
The bulk of solar coronal radiative loss consists of soft X-ray emission from quasi-static loops at the cores of Active Regions. In order to develop diagnostics for determining the heating mechanism of these loops from observations by coronal imaging instruments, I have developed analytical solutions for the temperature structure and scaling laws of loop strands for a wide range of heating funct…
▽ More
The bulk of solar coronal radiative loss consists of soft X-ray emission from quasi-static loops at the cores of Active Regions. In order to develop diagnostics for determining the heating mechanism of these loops from observations by coronal imaging instruments, I have developed analytical solutions for the temperature structure and scaling laws of loop strands for a wide range of heating functions, including footpoint heating, uniform heating, and heating concentrated at the loop apex. Key results are that the temperature profile depends only weakly on the heating distribution -- not sufficiently to be of significant diagnostic value -- and that the scaling laws survive for this wide range of heating distributions, but with the constant of proportionality in the RTV scaling law ($P_{0}L \thicksim T_{max}^3$) depending on the specific heating function. Furthermore, quasi-static analytical solutions do not exist for an excessive concentration of heating near the loop footpoints, a result in agreement with recent numerical simulations. It is demonstrated that a generalization of the solutions to the case of a strand with a variable diameter leads to only relatively small correction factors in the scaling laws and temperature profiles for constant diameter loop strands. A quintet of leading theoretical coronal heating mechanisms is shown to be captured by the formalism of this paper, and the differences in thermal structure between them may be verified through observations. Preliminary results from full numerical simulations demonstrate that, despite the simplifying assumptions, the analytical solutions from this paper are stable and accurate.
△ Less
Submitted 16 April, 2008; v1 submitted 14 April, 2008;
originally announced April 2008.
-
Multithermal Analysis of a CDS Coronal Loop
Authors:
J. T. Schmelz,
P. C. H. Martens
Abstract:
The observations from 1998 April 20 taken with the Coronal Diagnostics Spectrometer CDS on SOHO of a coronal loop on the limb have shown that the plasma was multi-thermal along each line of sight investigated, both before and after background subtraction. The latter result relied on Emission Measure Loci plots, but in this Letter, we used a forward folding technique to produce Differential Emiss…
▽ More
The observations from 1998 April 20 taken with the Coronal Diagnostics Spectrometer CDS on SOHO of a coronal loop on the limb have shown that the plasma was multi-thermal along each line of sight investigated, both before and after background subtraction. The latter result relied on Emission Measure Loci plots, but in this Letter, we used a forward folding technique to produce Differential Emission Measure curves. We also calculate DEM-weighted temperatures for the chosen pixels and find a gradient in temperature along the loop as a function of height that is not compatible with the flat profiles reported by numerous authors for loops observed with EIT on SOHO and TRACE. We also find discrepancies in excess of the mathematical expectation between some of the observed and predicted CDS line intensities. We demonstrate that these differences result from well-known limitations in our knowledge of the atomic data and are to be expected. We further show that the precision of the DEM is limited by the intrinsic width of the ion emissivity functions that are used to calculate the DEM. Hence we conclude that peaks and valleys in the DEM, while in principle not impossible, cannot be confirmed from the data.
△ Less
Submitted 16 November, 2005;
originally announced November 2005.
-
Cosmological stabilization of moduli with steep potentials
Authors:
R. Brustein,
S. P. de Alwis,
P. Martens
Abstract:
A scenario which overcomes the well-known cosmological overshoot problem associated with stabilizing moduli with steep potentials in string theory is proposed. Our proposal relies on the fact that moduli potentials are very steep and that generically their kinetic energy quickly becomes dominant. However, moduli kinetic energy red-shifts faster than other sources when the universe expands. So, i…
▽ More
A scenario which overcomes the well-known cosmological overshoot problem associated with stabilizing moduli with steep potentials in string theory is proposed. Our proposal relies on the fact that moduli potentials are very steep and that generically their kinetic energy quickly becomes dominant. However, moduli kinetic energy red-shifts faster than other sources when the universe expands. So, if any additional sources are present, even in very small amounts, they will inevitably become dominant. We show that in this case cosmic friction allows the dissipation of the large amount of moduli kinetic energy that is required for the field to be able to find an extremely shallow minimum. We present the idea using analytic methods and verify with some numerical examples.
△ Less
Submitted 20 August, 2004;
originally announced August 2004.
-
Can Streamer Blobs prevent the Buildup of the Interplanetetary Magnetic Field?
Authors:
M. K. van Aalst,
P. C. H. Martens,
A. J. C. Beliën
Abstract:
Coronal Mass Ejections continuously drag closed magnetic field lines away from the Sun, adding new flux to the interplanetary magnetic field (IMF). We propose that the outward-moving blobs that have been observed in helmet streamers are evidence of ongoing, small-scale reconnection in streamer current sheets, which may play an important role in the prevention of an indefinite buildup of the IMF.…
▽ More
Coronal Mass Ejections continuously drag closed magnetic field lines away from the Sun, adding new flux to the interplanetary magnetic field (IMF). We propose that the outward-moving blobs that have been observed in helmet streamers are evidence of ongoing, small-scale reconnection in streamer current sheets, which may play an important role in the prevention of an indefinite buildup of the IMF. Reconnection between two open field lines from both sides of a streamer current sheet creates a new closed field line, which becomes part of the helmet, and a disconnected field line, which moves outward. The blobs are formed by plasma from the streamer that is swept up in the trough of the outward moving field line. We show that this mechanism is supported by observations from SOHO/LASCO. Additionally, we propose a thorough statistical study to quantify the contribution of blob formation to the reduction of the IMF, and indicate how this mechanism may be verified by observations with SOHO/UVCS and the proposed NASA STEREO and ESA Polar Orbiter missions.
△ Less
Submitted 4 December, 1998;
originally announced December 1998.
-
An approximate self-consistent theory of the magnetic field of fluted penumbrae
Authors:
T. Neukirch,
P. C. H. Martens
Abstract:
A self-consistent mathematical description of the magnetic field of fluted sunspot penumbrae is presented. This description is based on an expansion of the nonlinear force-free magnetohydrostatic equations written in cylindrical coordinates. The lowest order solutions are mathematically equivalent to laminated force-free equilibria in Cartesian geometry. The lowest order solutions have no toroid…
▽ More
A self-consistent mathematical description of the magnetic field of fluted sunspot penumbrae is presented. This description is based on an expansion of the nonlinear force-free magnetohydrostatic equations written in cylindrical coordinates. The lowest order solutions are mathematically equivalent to laminated force-free equilibria in Cartesian geometry. The lowest order solutions have no toroidal component of the magnetic field and the magnetic pressure does not vary with azimuth but the solutions allow arbitrary variations of the magnetic field components with azimuth. Explicit solutions are presented which have a realistic radial profile of the magnetic field strength and reproduce the basic features of the observations.
△ Less
Submitted 8 January, 1998; v1 submitted 28 November, 1997;
originally announced November 1997.