-
Forecasting SEP Events During Solar Cycles 23 and 24 Using Interpretable Machine Learning
Authors:
Spiridon Kasapis,
Irina N. Kitiashvili,
Paul Kosovich,
Alexander G. Kosovichev,
Viacheslav M. Sadykov,
Patrick O'Keefe,
Vincent Wang
Abstract:
Prediction of the Solar Energetic Particle (SEP) events garner increasing interest as space missions extend beyond Earth's protective magnetosphere. These events, which are, in most cases, products of magnetic reconnection-driven processes during solar flares or fast coronal-mass-ejection-driven shock waves, pose significant radiation hazards to aviation, space-based electronics, and particularly,…
▽ More
Prediction of the Solar Energetic Particle (SEP) events garner increasing interest as space missions extend beyond Earth's protective magnetosphere. These events, which are, in most cases, products of magnetic reconnection-driven processes during solar flares or fast coronal-mass-ejection-driven shock waves, pose significant radiation hazards to aviation, space-based electronics, and particularly, space exploration. In this work, we utilize the recently developed dataset that combines the Solar Dynamics Observatory/Helioseismic and Magnetic Imager's (SDO/HMI) Space weather HMI Active Region Patches (SHARP) and the Solar and Heliospheric Observatory/Michelson Doppler Imager's (SoHO/MDI) Space Weather MDI Active Region Patches (SMARP). We employ a suite of machine learning strategies, including Support Vector Machines (SVM) and regression models, to evaluate the predictive potential of this new data product for a forecast of post-solar flare SEP events. Our study indicates that despite the augmented volume of data, the prediction accuracy reaches 0.7 +- 0.1, which aligns with but does not exceed these published benchmarks. A linear SVM model with training and testing configurations that mimic an operational setting (positive-negative imbalance) reveals a slight increase (+ 0.04 +- 0.05) in the accuracy of a 14-hour SEP forecast compared to previous studies. This outcome emphasizes the imperative for more sophisticated, physics-informed models to better understand the underlying processes leading to SEP events.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
Predicting the Emergence of Solar Active Regions Using Machine Learning
Authors:
Spiridon Kasapis,
Irina N. Kitiashvili,
Alexander G. Kosovichev,
John T. Stefan,
Bhairavi Apte
Abstract:
To create early warning capabilities for upcoming Space Weather disturbances, we have selected a dataset of 61 emerging active regions, which allows us to identify characteristic features in the evolution of acoustic power density to predict continuum intensity emergence. For our study, we have utilized Doppler shift and continuum intensity observations from the Helioseismic and Magnetic Imager (H…
▽ More
To create early warning capabilities for upcoming Space Weather disturbances, we have selected a dataset of 61 emerging active regions, which allows us to identify characteristic features in the evolution of acoustic power density to predict continuum intensity emergence. For our study, we have utilized Doppler shift and continuum intensity observations from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). The local tracking of 30.66 x 30.66-degree patches in the vicinity of active regions allowed us to trace the evolution of active regions starting from the pre-emergence state. We have developed a machine learning model to capture the acoustic power flux density variations associated with upcoming magnetic flux emergence. The trained Long Short-Term Memory (LSTM) model is able to predict 5 hours ahead whether, in a given area of the solar surface, continuum intensity values will decrease. The performed study allows us to investigate the potential of the machine learning approach to predict the emergence of active regions using acoustic power maps as input.
△ Less
Submitted 13 February, 2024;
originally announced February 2024.
-
Time Series of Magnetic Field Parameters of Merged MDI and HMI Space-Weather Active Region Patches as Potential Tool for Solar Flare Forecasting
Authors:
Paul A. Kosovich,
Viacheslav M. Sadykov,
Alexander G. Kosovichev,
Spiridon Kasapis,
Irina N. Kitiashvili,
Patrick M. O'Keefe,
Aatiya Ali,
Vincent Oria,
Samuel Granovsky,
Chun Jie Chong,
Gelu M. Nita
Abstract:
Solar flare prediction studies have been recently conducted with the use of Space-Weather MDI (Michelson Doppler Imager onboard Solar and Heliospheric Observatory) Active Region Patches (SMARP) and Space-Weather HMI (Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory) Active Region Patches (SHARP), which are two currently available data products containing magnetic field character…
▽ More
Solar flare prediction studies have been recently conducted with the use of Space-Weather MDI (Michelson Doppler Imager onboard Solar and Heliospheric Observatory) Active Region Patches (SMARP) and Space-Weather HMI (Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory) Active Region Patches (SHARP), which are two currently available data products containing magnetic field characteristics of solar active regions. The present work is an effort to combine them into one data product, and perform some initial statistical analyses in order to further expand their application in space weather forecasting. The combined data are derived by filtering, rescaling, and merging the SMARP with SHARP parameters, which can then be spatially reduced to create uniform multivariate time series. The resulting combined MDI-HMI dataset currently spans the period between April 4, 1996, and December 13, 2022, and may be extended to a more recent date. This provides an opportunity to correlate and compare it with other space weather time series, such as the daily solar flare index or the statistical properties of the soft X-ray flux measured by the Geostationary Operational Environmental Satellites (GOES). Time-lagged cross-correlation indicates that a relationship may exist, where some magnetic field properties of active regions lead the flare index in time. Applying the rolling window technique makes it possible to see how this leader-follower dynamic varies with time. Preliminary results indicate that areas of high correlation generally correspond to increased flare activity during the peak solar cycle.
△ Less
Submitted 22 February, 2024; v1 submitted 10 January, 2024;
originally announced January 2024.
-
The Puzzling Structure of Solar Convection: Window into the Dynamo
Authors:
Nicholas A. Featherstone,
Kyle C. Augustson,
Jonathan M. Aurnou,
Catherine Blume,
Benjamin P. Brown,
Nicholas Brummell,
Keaton J. Burns,
Michael A. Calkins,
Maria Camisassa,
Mausumi Dikpati,
Yuhong Fan,
J. R. Fuentes,
Gustavo Guerrero,
Bradley W. Hindman,
Keith Julien,
Irina N. Kitiashvili,
Lydia Korre,
Daniel Lecoanet,
Bhishek Manek,
Loren Matilsky,
Mark Miesch,
Nicholas J. Nelson,
Jeffrey S. Oishi,
Whitney T. Powers,
Matthias Rempel
, et al. (3 additional authors not shown)
Abstract:
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective…
▽ More
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements.
△ Less
Submitted 15 May, 2023;
originally announced May 2023.
-
Non-zero phase-shifts of acoustic waves in the lower solar atmosphere measured from realistic simulations and their role in local helioseismology
Authors:
M. Waidele,
Junwei Zhao,
I. N. Kitiashvili
Abstract:
Previous studies analyzing the evanescent nature of acoustic waves in the lower solar atmosphere, up to 300\,km above the photosphere, have shown an unexpected phase shift of an order of 1\,s between different heights. Those studies investigated the spectral line \ion{Fe}{1} 6173.3\,Å, commonly used for helioseismic measurements. Such phase-shifts can contribute to a misinterpretation of the measu…
▽ More
Previous studies analyzing the evanescent nature of acoustic waves in the lower solar atmosphere, up to 300\,km above the photosphere, have shown an unexpected phase shift of an order of 1\,s between different heights. Those studies investigated the spectral line \ion{Fe}{1} 6173.3\,Å, commonly used for helioseismic measurements. Such phase-shifts can contribute to a misinterpretation of the measured travel times in local helioseismology, complicating inferences of, e.g., the deep meridional flow. In this study, we carry out phase-shift computations using a simulated, fully radiative, and convective atmosphere from which the \ion{Fe}{1} 6173.3\,Å line is synthesized. The resulting phase-shifts as functions of frequency across multiple heights show non-zero values in evanescent waves, similar to what was found in observational data. Comparing the Doppler-velocities estimated from the synthesized absorption line with the true velocities directly obtained from the simulated plasma motions, we find substantial differences in phase-shifts between the two. This leads us to hypothesize that the non-adiabaticity of the solar atmosphere yields extra phase-shift contributions to Doppler velocities. Finally, computing phase-differences for different viewing angles reveals a systematic center-to-limb variation, similar to what is present in observations. Overall, this study helps to improve our understanding of the physical cause of the helioseismic center-to-limb effect.
△ Less
Submitted 4 April, 2023;
originally announced April 2023.
-
The Random Hivemind: An Ensemble Deep Learner Application to Solar Energetic Particle Prediction Problem
Authors:
Patrick M. O'Keefe,
Viacheslav Sadykov,
Alexander Kosovichev,
Irina N. Kitiashvili,
Vincent Oria,
Gelu M. Nita,
Fraila Francis,
Chun-Jie Chong,
Paul Kosovich,
Aatiya Ali,
Russell D. Marroquin
Abstract:
The application of machine learning and deep learning, including the wide use of non-ensemble, conventional neural networks (CoNN), for predicting various phenomena has become very popular in recent years thanks to the efficiencies and the abilities of these techniques to find relationships in data without human intervention. However, certain CoNN setups may not work on some datasets, especially i…
▽ More
The application of machine learning and deep learning, including the wide use of non-ensemble, conventional neural networks (CoNN), for predicting various phenomena has become very popular in recent years thanks to the efficiencies and the abilities of these techniques to find relationships in data without human intervention. However, certain CoNN setups may not work on some datasets, especially if the parameters passed to it, including model parameters and hyperparameters, are arguably arbitrary in nature and need to continuously be updated with the need to retrain the model. This concern can be partially alleviated by employing committees of neural networks that are identical in terms of input features and architectures, initialized randomly, and "vote" on the decisions made by the committees as a whole. Yet, it is possible for the committee members to "agree" on identical sets of weights and biases for all nodes and edges. Members of these committees also cannot be expanded to accommodate new features and entire committees must therefore be retrained in order to do so. We propose the Random Hivemind (RH) approach, which helps to alleviate this concern by having multiple neural network estimators make decisions based on random permutations of features and prescribing a method to determine the weight of the decision of each individual estimator. The effectiveness of RH is demonstrated through experimentation in the predictions of hazardous Solar Energetic Particle (SEP) events by comparing it to that of using both CoNNs and the aforementioned setup of committees. Our results demonstrate that RH, while having a comparable or better performance than the CoNN and a Committee-based approach, demonstrates a lesser score spread for the individual experiments, and shows promising results with respect to capturing almost every single flare instance leading to SEPs.
△ Less
Submitted 7 May, 2024; v1 submitted 14 March, 2023;
originally announced March 2023.
-
Statistical Study of the Correlation between Solar Energetic Particles and Properties of Active Regions
Authors:
Russell D. Marroquin,
Viacheslav Sadykov,
Alexander Kosovichev,
Irina N. Kitiashvili,
Vincent Oria,
Gelu M. Nita,
Egor Illarionov,
Patrick M. O'Keefe,
Fraila Francis,
Chun-Jie Chong,
Paul Kosovich,
Aatiya Ali
Abstract:
The flux of energetic particles originating from the Sun fluctuates during the solar cycles. It depends on the number and properties of Active Regions (ARs) present in a single day and associated solar activities, such as solar flares and coronal mass ejections (CMEs). Observational records of the Space Weather Prediction Center (SWPC NOAA) enable the creation of time-indexed databases containing…
▽ More
The flux of energetic particles originating from the Sun fluctuates during the solar cycles. It depends on the number and properties of Active Regions (ARs) present in a single day and associated solar activities, such as solar flares and coronal mass ejections (CMEs). Observational records of the Space Weather Prediction Center (SWPC NOAA) enable the creation of time-indexed databases containing information about ARs and particle flux enhancements, most widely known as Solar Energetic Particle events (SEPs). In this work, we utilize the data available for Solar Cycles 21-24, and the initial phase of Cycle 25 to perform a statistical analysis of the correlation between SEPs and properties of ARs inferred from the McIntosh and Hale classifications. We find that the complexity of the magnetic field, longitudinal location, area, and penumbra type of the largest sunspot of ARs are most correlated with the production of SEPs. It is found that most SEPs ($\approx$60\%, or 108 out of 181 considered events) were generated from an AR classified with the 'k' McIntosh subclass as the second component, and these ARs are more likely to produce SEPs if they fall in a Hale class containing $δ$ component. The resulting database containing information about SEP events and ARs is publicly available and can be used for the development of Machine Learning (ML) models to predict the occurrence of SEPs.
△ Less
Submitted 22 March, 2023; v1 submitted 10 March, 2023;
originally announced March 2023.
-
Predicting Solar Proton Events of Solar Cycles 22-24 using GOES Proton & soft X-Ray flux features
Authors:
Aatiya Ali,
Viacheslav Sadykov,
Alexander Kosovichev,
Irina N. Kitiashvili,
Vincent Oria,
Gelu M. Nita,
Egor Illarionov,
Patrick M. O'Keefe,
Fraila Francis,
Chun-Jie Chong,
Paul Kosovich,
Russell D. Marroquin
Abstract:
Solar Energetic Particle (SEP) events and their major subclass, Solar Proton Events (SPEs), can have unfavorable consequences on numerous aspects of life and technology, making them one of the most harmful effects of solar activity. Garnering knowledge preceding such events by studying operational data flows is essential for their forecasting. Considering only Solar Cycle (SC) 24 in our previous s…
▽ More
Solar Energetic Particle (SEP) events and their major subclass, Solar Proton Events (SPEs), can have unfavorable consequences on numerous aspects of life and technology, making them one of the most harmful effects of solar activity. Garnering knowledge preceding such events by studying operational data flows is essential for their forecasting. Considering only Solar Cycle (SC) 24 in our previous study, Sadykov et al. 2021, we found that it may be sufficient to utilize only proton and soft X-ray (SXR) parameters for SPE forecasts. Here, we report a catalog recording $\geq$ 10 MeV $\geq$ 10 particle flux unit SPEs with their properties, spanning SCs 22-24, using NOAA's Geostationary Operational Environmental Satellite flux data. We report an additional catalog of daily proton and SXR flux statistics for this period, employing it to test the application of machine learning (ML) on the prediction of SPEs using a Support Vector Machine (SVM) and eXtreme Gradient Boosting (XGBoost). We explore the effects of training models with data from one and two SCs, evaluating how transferable a model can be across different time periods. XGBoost proved to be more accurate than SVMs for almost every test considered, while outperforming operational SWPC NOAA predictions and a persistence forecast. Interestingly, training done with SC 24 produces weaker TSS and HSS2, even when paired with SC 22 or SC 23, indicating transferability issues. This work contributes towards validating forecasts using long-spanning data -- an understudied area in SEP research that should be considered to verify the cross-cycle robustness of ML-driven forecasts.
△ Less
Submitted 7 November, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
Leptocline as a Shallow Substructure of Near-Surface Shear Layer in 3D Radiative Hydrodynamic Simulations
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray,
Viacheslav M. Sadykov,
Gustavo Guerrero
Abstract:
Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees lati…
▽ More
Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees latitude. The simulation results reveal the development of a shallow 10-Mm deep substructure of the Near-Surface Shear Layer (NSSL), characterized by a strong radial rotational gradient and self-organized meridional flows. This shallow layer ("leptocline") is located in the hydrogen ionization zone associated with enhanced anisotropic overshooting-type flows into a less unstable layer between the H and HeII ionization zones. We discuss current observational evidence of the presence of the leptocline and show that the radial variations of the differential rotation and meridional flow profiles obtained from the simulations in this layer qualitatively agree with helioseismic observations.
△ Less
Submitted 27 August, 2022; v1 submitted 2 March, 2022;
originally announced March 2022.
-
Modeling Stellar Jitter for the Detection of Earth-Mass Exoplanets via Precision Radial Velocity Measurements
Authors:
Samuel Granovsky,
Irina N. Kitiashvili,
Alan Wray
Abstract:
The detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity measurements. Increasing radial velocity precision requires developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking ste…
▽ More
The detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity measurements. Increasing radial velocity precision requires developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking stellar variability and instrumental uncertainties into account. In this work, we investigate and quantify stellar disturbances of the planet-hosting solar-type star HD121504 from 3D radiative modeling obtained with the StellarBox code. The model has been used for determining statistical properties of the turbulent plasma and obtaining synthetic spectroscopic observations for several Fe I lines at different locations on the stellar disk to mimic high-resolution spectroscopic observations.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
3D Modeling of the Structure and Dynamics of a Main-Sequence F-type Star
Authors:
Irina N. Kitiashvili,
Alan A. Wray
Abstract:
Current state-of-the-art computational modeling makes it possible to build realistic models of stellar convection zones and atmospheres that take into account chemical composition, radiative effects, ionization, and turbulence. The standard 1D mixing-length-based evolutionary models are not able to capture many physical processes of the stellar interior dynamics. Mixing-length models provide an in…
▽ More
Current state-of-the-art computational modeling makes it possible to build realistic models of stellar convection zones and atmospheres that take into account chemical composition, radiative effects, ionization, and turbulence. The standard 1D mixing-length-based evolutionary models are not able to capture many physical processes of the stellar interior dynamics. Mixing-length models provide an initial approximation of stellar structure that can be used to initialize 3D radiative hydrodynamics simulations which include realistic modeling of turbulence, radiation, and other phenomena.
In this paper, we present 3D radiative hydrodynamic simulations of an F-type main-sequence star with 1.47 solar mass. The computational domain includes the upper layers of the radiation zone, the entire convection zone, and the photosphere. The effects of stellar rotation is modeled in the f-plane approximation. These simulations provide new insight into the properties of the convective overshoot region, the dynamics of the near-surface, highly turbulent layer, and the structure and dynamics of granulation. They reveal anti-solar-type differential rotation and latitudinal dependence of the tachocline location.
△ Less
Submitted 26 July, 2021;
originally announced July 2021.
-
Radiation Data Portal: Integration of Radiation Measurements at the Aviation Altitudes and Solar-Terrestrial Environment Observations
Authors:
Viacheslav M Sadykov,
Irina N Kitiashvili,
W Kent Tobiska,
Madhulika Guhathakurta
Abstract:
The impact of radiation dramatically increases at high altitudes in the Earth's atmosphere and in space. Therefore, monitoring and access to radiation environment measurements are critical for estimating the radiation exposure risks of aircraft and spacecraft crews and the impact of space weather disturbances on electronics. Addressing these needs requires convenient access to multi-source radiati…
▽ More
The impact of radiation dramatically increases at high altitudes in the Earth's atmosphere and in space. Therefore, monitoring and access to radiation environment measurements are critical for estimating the radiation exposure risks of aircraft and spacecraft crews and the impact of space weather disturbances on electronics. Addressing these needs requires convenient access to multi-source radiation environment data and enhancement of visualization and search capabilities. The Radiation Data Portal represents an interactive web-based application for search and visualization of in-flight radiation measurements. The Portal enhances the exploration capabilities of various properties of the radiation environment and provides measurements of dose rates along with information on space weather-related conditions. The Radiation Data Portal back-end is a MySQL relational database that contains the radiation measurements obtained from the Automated Radiation Measurements for Aerospace Safety (ARMAS) device and the soft X-ray and proton flux measurements from the Geostationary Operational Environmental Satellite (GOES). The implemented Application Programming Interface (API) and Python routines allow a user to retrieve the database records without interaction with the web interface. As a use case of the Radiation Data Portal, we examine ARMAS measurements during an enhancement of the Solar Proton (SP) fluxes, known as Solar Proton Events (SPEs), and compare them to measurements during SP-quiet periods.
△ Less
Submitted 12 March, 2021;
originally announced March 2021.
-
Compression of Solar Spectroscopic Observations: a Case Study of Mg II k Spectral Line Profiles Observed by NASA's IRIS Satellite
Authors:
Viacheslav M Sadykov,
Irina N Kitiashvili,
Alberto Sainz Dalda,
Vincent Oria,
Alexander G Kosovichev,
Egor Illarionov
Abstract:
In this study we extract the deep features and investigate the compression of the Mg II k spectral line profiles observed in quiet Sun regions by NASA's IRIS satellite. The data set of line profiles used for the analysis was obtained on April 20th, 2020, at the center of the solar disc, and contains almost 300,000 individual Mg II k line profiles after data cleaning. The data are separated into tr…
▽ More
In this study we extract the deep features and investigate the compression of the Mg II k spectral line profiles observed in quiet Sun regions by NASA's IRIS satellite. The data set of line profiles used for the analysis was obtained on April 20th, 2020, at the center of the solar disc, and contains almost 300,000 individual Mg II k line profiles after data cleaning. The data are separated into train and test subsets. The train subset was used to train the autoencoder of the varying embedding layer size. The early stopping criterion was implemented on the test subset to prevent the model from overfitting. Our results indicate that it is possible to compress the spectral line profiles more than 27 times (which corresponds to the reduction of the data dimensionality from 110 to 4) while having a 4 DN average reconstruction error, which is comparable to the variations in the line continuum. The mean squared error and the reconstruction error of even statistical moments sharply decrease when the dimensionality of the embedding layer increases from 1 to 4 and almost stop decreasing for higher numbers. The observed occasional improvements in training for values higher than 4 indicate that a better compact embedding may potentially be obtained if other training strategies and longer training times are used. The features learned for the critical four-dimensional case can be interpreted. In particular, three of these four features mainly control the line width, line asymmetry, and line dip formation respectively. The presented results are the first attempt to obtain a compact embedding for spectroscopic line profiles and confirm the value of this approach, in particular for feature extraction, data compression, and denoising.
△ Less
Submitted 16 June, 2021; v1 submitted 12 March, 2021;
originally announced March 2021.
-
Connecting Atmospheric Properties and Synthetic Emission of Shock Waves Using 3D RMHD Simulations of Quiet Sun
Authors:
Viacheslav M. Sadykov,
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray
Abstract:
We analyze the evolution of shock waves in high-resolution 3D radiative MHD simulations of the quiet Sun and their synthetic emission characteristics. The simulations model the dynamics of a 12.8x12.8x15.2 Mm quiet-Sun region (including a 5.2 Mm layer of the upper convection zone and a 10 Mm atmosphere from the photosphere to corona) with an initially uniform vertical magnetic field of 10 G, natur…
▽ More
We analyze the evolution of shock waves in high-resolution 3D radiative MHD simulations of the quiet Sun and their synthetic emission characteristics. The simulations model the dynamics of a 12.8x12.8x15.2 Mm quiet-Sun region (including a 5.2 Mm layer of the upper convection zone and a 10 Mm atmosphere from the photosphere to corona) with an initially uniform vertical magnetic field of 10 G, naturally driven by convective flows. We synthesize the Mg II and C II spectral lines observed by the IRIS satellite and EUV emission observed by the SDO/AIA telescope. Synthetic observations are obtained using the RH1.5D radiative transfer code and temperature response functions at both the numerical and instrumental resolutions. We found that the Doppler velocity jumps of the C II 1334.5 A IRIS line and a relative enhancement of the emission in the 335 A SDO/AIA channel are the best proxies for the enthalpy deposited by shock waves into the corona (with Kendall's $τ$ correlation coefficients of 0.59 and 0.38, respectively). The synthetic emission of the lines and extreme ultraviolet passbands are correlated with each other during the shock wave propagation. All studied shocks are mostly hydrodynamic (i.e., the magnetic energy carried by horizontal fields is $\leq{}$2.6% of the enthalpy for all events) and have Mach numbers > 1.0-1.2 in the low corona. The study reveals the possibility of diagnosing energy transport by shock waves into the solar corona, as well as their other properties, by using IRIS and SDO/AIA sensing observations.
△ Less
Submitted 12 March, 2021; v1 submitted 13 August, 2020;
originally announced August 2020.
-
Global Evolution of Solar Magnetic Fields and Prediction of Activity Cycles
Authors:
Irina N. Kitiashvili
Abstract:
Prediction of solar activity cycles is challenging because physical processes inside the Sun involve a broad range of multiscale dynamics that no model can reproduce and because the available observations are highly limited and cover mostly surface layers. Helioseismology makes it possible to probe solar dynamics in the convective zone, but variations in differential rotation and meridional circul…
▽ More
Prediction of solar activity cycles is challenging because physical processes inside the Sun involve a broad range of multiscale dynamics that no model can reproduce and because the available observations are highly limited and cover mostly surface layers. Helioseismology makes it possible to probe solar dynamics in the convective zone, but variations in differential rotation and meridional circulation are currently available for only two solar activity cycles. It has been demonstrated that sunspot observations, which cover over 400 years, can be used to calibrate the Parker-Kleeorin-Ruzmaikin dynamo model, and that the Ensemble Kalman Filter (EnKF) method can be used to link the modeled magnetic fields to sunspot observations and make reliable predictions of a following activity cycle. However, for more accurate predictions, it is necessary to use actual observations of the solar magnetic fields, which are available only for the last four solar cycles. In this paper I briefly discuss the influence of the limited number of available observations on the accuracy of EnKF estimates of solar cycle parameters, the criteria to evaluate the predictions, and application of synoptic magnetograms to the prediction of solar activity.
△ Less
Submitted 10 March, 2020;
originally announced March 2020.
-
Resolving Power of Asteroseismic Inversion of the Kepler Legacy Sample
Authors:
Alexander G. Kosovichev,
Irina N. Kitiashvili
Abstract:
The Kepler Asteroseismic Legacy Project provided frequencies, separation ratios, error estimates, and covariance matrices for 66 Kepler main sequence targets. Most of the previous analysis of these data was focused on fitting standard stellar models. We present results of direct asteroseismic inversions using the method of optimally localized averages (OLA), which effectively eliminates the surfac…
▽ More
The Kepler Asteroseismic Legacy Project provided frequencies, separation ratios, error estimates, and covariance matrices for 66 Kepler main sequence targets. Most of the previous analysis of these data was focused on fitting standard stellar models. We present results of direct asteroseismic inversions using the method of optimally localized averages (OLA), which effectively eliminates the surface effects and attempts to resolve the stellar core structure. The inversions are presented for various structure properties, including the density stratification and sound speed. The results show that the mixed modes observed in post-main sequence F-type stars allow us to resolve the stellar core structure and reveal significant deviations from the evolutionary models obtained by the grid-fitting procedure to match the observed oscillation frequencies.
△ Less
Submitted 23 February, 2020;
originally announced February 2020.
-
Effects of Observational Data Shortage on Accuracy of Global Solar Activity Forecast
Authors:
Irina N. Kitiashvili
Abstract:
Building a reliable forecast of solar activity is a long-standing problem that requires to accurately describe past and current global dynamics. However, synoptic observations of magnetic fields and subsurface flows became available relatively recently. In this paper, we present an investigation of effects of short observational data series on accuracy of solar cycle prediction. This analysis is p…
▽ More
Building a reliable forecast of solar activity is a long-standing problem that requires to accurately describe past and current global dynamics. However, synoptic observations of magnetic fields and subsurface flows became available relatively recently. In this paper, we present an investigation of effects of short observational data series on accuracy of solar cycle prediction. This analysis is performed using the annual sunspot number time-series applied to the Parker-Kleeorin-Ruzmaikin dynamo model and employing the Ensemble Kalman Filter (EnKF) data assimilation method. The testing of the cycle prediction accuracy is performed for the last six cycles (from Solar Cycle 19 to 24) by sequentially shortening the observational data series that are used for prediction of a target cycle, and evaluating the prediction accuracy according to specified criteria. According to the analysis, reliable activity predictions can be made using relatively short time-series of the sunspot number. It demonstrated that even two cycles of available observations allow us to obtain reasonable forecasts.
△ Less
Submitted 21 September, 2020; v1 submitted 25 January, 2020;
originally announced January 2020.
-
Application of Synoptic Magnetograms to Global Solar Activity Forecast
Authors:
Irina N. Kitiashvili
Abstract:
Synoptic magnetograms provide us with knowledge about the evolution of magnetic fields on the solar surface and present important information for forecasting future solar activity. In this work, poloidal and toroidal magnetic field components derived from synoptic magnetograms are assimilated, using the Ensemble Kalman Filter method, into a mean-field dynamo model based on Parker's migratory dynam…
▽ More
Synoptic magnetograms provide us with knowledge about the evolution of magnetic fields on the solar surface and present important information for forecasting future solar activity. In this work, poloidal and toroidal magnetic field components derived from synoptic magnetograms are assimilated, using the Ensemble Kalman Filter method, into a mean-field dynamo model based on Parker's migratory dynamo theory complemented by magnetic helicity conservation. It was found that the predicted toroidal field is in good agreement with observations for almost the entire following solar cycle. However, poloidal field predictions agree with observations only for the first 2 - 3 years of the predicted cycle. The results indicate that the upcoming Solar Maximum of Cycle 25 (SC25) is expected to be weaker than the current Cycle 24. The model results show that a deep extended solar activity minimum is expected during 2019 - 2021, and that the next solar maximum will occur in 2024 - 2025. The sunspot number at the maximum will be about 50 with an error estimate of 15 - 30 %. The maximum will likely have a double peak or show extended periods (for 2 - 2.5 years) of high activity. According to the hemispheric prediction results, SC25 will start in 2020 in the Southern hemisphere, and will have a maximum in 2024 with a sunspot number of about 28. In the Northern hemisphere the cycle will be delayed for about 1 year (with an error of $\pm 0.5$ year), and reach a maximum in 2025 with a sunspot number of about 23.
△ Less
Submitted 2 October, 2019;
originally announced October 2019.
-
Response of SDO/HMI observables to heating of the solar atmosphere by precipitating high-energy electrons
Authors:
Viacheslav M Sadykov,
Alexander G Kosovichev,
Irina N Kitiashvili,
Graham S Kerr
Abstract:
We perform an analysis of the line-of-sight (LOS) observables of the SDO/HMI for models of the solar atmosphere heated by precipitating high-energy electrons during solar flares. The radiative hydrodynamic (RADYN) flare models are obtained from the F-CHROMA database. The Stokes profiles for the Fe 6173 A line observed by SDO/HMI are calculated using the radiative transfer code RH1.5D assuming stat…
▽ More
We perform an analysis of the line-of-sight (LOS) observables of the SDO/HMI for models of the solar atmosphere heated by precipitating high-energy electrons during solar flares. The radiative hydrodynamic (RADYN) flare models are obtained from the F-CHROMA database. The Stokes profiles for the Fe 6173 A line observed by SDO/HMI are calculated using the radiative transfer code RH1.5D assuming statistical equilibrium for atomic level populations and imposing uniform background vertical magnetic field of various strength. The SDO/HMI observing sequence and LOS data processing pipeline algorithm are applied to derive the observables (continuum intensity, line depth, Doppler velocity, LOS magnetic field). Our results reveal that the strongest deviations of the observables from the actual spectroscopic line parameters are found for the model with a total energy deposited of $E_{total}=1.0\times{}10^{12}$ erg cm$^{-2}$, injected with power-law spectral index of $δ=3$ above a low-energy cutoff of $E_{c}=25$ keV. The magnitudes of the velocity and magnetic field deviations depend on the imposed magnetic field, and can reach 0.35 km/s for LOS velocities, 90 G for LOS magnetic field, and 3\% for continuum enhancement for the 1000 G imposed LOS magnetic field setup. For $E_{total}\geq{}3.0\times{}10^{11}$ erg cm$^{-2}$ models, the velocity and magnetic field deviations are most strongly correlated with the energy flux carried by $\sim$50 keV electrons, and the continuum enhancement is correlated with the synthesized $\sim$55-60 keV hard X-ray photon flux. The relatively low magnitudes of perturbations of the observables and absence of magnetic field sign reversals suggest that the considered radiative hydrodynamic beam heating models augmented with the uniform vertical magnetic field setups cannot explain strong transient changes found in the SDO/HMI observations.
△ Less
Submitted 1 March, 2020; v1 submitted 25 June, 2019;
originally announced June 2019.
-
Reconstructing Extreme Space Weather from Planet Hosting Stars
Authors:
V. S. Airapetian,
V. Adibekyan,
M. Ansdell,
D. Alexander,
T. Bastian,
S. Boro Saikia,
A. S. Brun,
O. Cohen,
M. Cuntz,
W. Danchi,
J. Davenport,
J. DeNolfo,
R. DeVore,
C. F. Dong,
J. J. Drake,
K. France,
F. Fraschetti,
K. Herbst,
K. Garcia-Sage,
M. Gillon,
A. Glocer,
J. L. Grenfell,
G. Gronoff,
N. Gopalswamy,
M. Guedel
, et al. (58 additional authors not shown)
Abstract:
The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditi…
▽ More
The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
Origin of Deep Acoustic Sources Associated with Solar Magnetic Structures
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Nagi N. Mansour,
Alan A. Wray,
Timothy A. Sandstrom
Abstract:
It is generally accepted that solar acoustic (p) modes are excited by near-surface turbulent motions, in particular, by downdrafts and interacting vortices in intergranular lanes. Recent analysis of Solar Dynamics Observatory data by (Zhao et al., 2015) revealed fast-moving waves around sunspots, which are consistent with magnetoacoustic waves excited approximately 5 Mm beneath the sunspot. We ana…
▽ More
It is generally accepted that solar acoustic (p) modes are excited by near-surface turbulent motions, in particular, by downdrafts and interacting vortices in intergranular lanes. Recent analysis of Solar Dynamics Observatory data by (Zhao et al., 2015) revealed fast-moving waves around sunspots, which are consistent with magnetoacoustic waves excited approximately 5 Mm beneath the sunspot. We analyzed 3D radiative MHD simulations of solar magnetoconvection with a self-organized pore-like magnetic structure, and identified more than 600 individual acoustic events both inside and outside of this structure. By performing a case-by-case study, we found that surrounding the magnetic structure, acoustic sources are associated with downdrafts. Their depth correlates with downdraft speed and magnetic fields. The sources often can be transported into deeper layers by downdrafts. The wave front shape, in the case of a strong or inclined downdraft, can be stretched along the downdraft. Inside the magnetic structure, excitation of acoustic waves is driven by converging flows. Frequently, strong converging plasma streams hit the structure boundaries, causing compressions in its interior that excite acoustic waves. Analysis of the depth distribution of acoustic events shows the strongest concentration at 0.2 - 1 Mm beneath the surface for the outside sources and 2.5 - 3 Mm for the excitation event inside the structure.
△ Less
Submitted 14 October, 2018;
originally announced October 2018.
-
Statistical Properties of Soft X-ray emission of solar flares
Authors:
Viacheslav M Sadykov,
Alexander G Kosovichev,
Irina N Kitiashvili,
Alexander Frolov
Abstract:
We present a statistical analysis of properties of Soft X-Ray (SXR) emission, plasma temperature (T), and emission measure (EM), derived from GOES observations of flares in 2002-2017. The temperature and emission measures are obtained using the TEBBS algorithm (Ryan et al. 2012), which delivers reliable results together with uncertainties even for weak B-class flare events. More than 96% of flares…
▽ More
We present a statistical analysis of properties of Soft X-Ray (SXR) emission, plasma temperature (T), and emission measure (EM), derived from GOES observations of flares in 2002-2017. The temperature and emission measures are obtained using the TEBBS algorithm (Ryan et al. 2012), which delivers reliable results together with uncertainties even for weak B-class flare events. More than 96% of flares demonstrate a sequential appearance of T, SXR, and EM maxima, in agreement with the expected behavior of the chromospheric evaporation process. The relative number of such flares increases with increasing the SXR flux maximum. The SXR maximum is closer in time to the T maximum for B-class flares than for >=C-class flares, while it is very close to the EM maximum for M- and X-class flares. We define flares as "T-controlled" if the time interval between the SXR and T maxima is at least two times shorter than the interval between the EM and SXR maxima, and as "EM-controlled" if the time interval between the EM and SXR maxima is at least two times shorter than the interval between the SXR and T maxima. For any considered flare class range, the T-controlled events compared to EM-controlled events have: a) higher EM but lower T; b) longer durations and shorter relative growth times; c) longer FWHM and characteristic decay times. Interpretation of these statistical results based on analysis of a single loop dynamics suggests that for flares of the same class range, the T-controlled events can be developed in longer loops than the EM-controlled events.
△ Less
Submitted 12 February, 2019; v1 submitted 12 October, 2018;
originally announced October 2018.
-
2-D and 3-D Models of Convective Turbulence and Oscillations in Intermediate-Mass Main-Sequence Stars
Authors:
Joyce A. Guzik,
T. H. Morgan,
N. J. Nelson,
C. Lovekin,
K. Kosak,
I. N. Kitiashvili,
N. N. Mansour,
A. Kosovichev
Abstract:
We present multidimensional modeling of convection and oscillations in main-sequence stars somewhat more massive than the Sun, using three separate approaches: 1) Using the 3-D planar StellarBox radiation hydrodynamics code to model the envelope convection zone and part of the radiative zone. Our goals are to examine the interaction of stellar pulsations with turbulent convection in the envelope,…
▽ More
We present multidimensional modeling of convection and oscillations in main-sequence stars somewhat more massive than the Sun, using three separate approaches: 1) Using the 3-D planar StellarBox radiation hydrodynamics code to model the envelope convection zone and part of the radiative zone. Our goals are to examine the interaction of stellar pulsations with turbulent convection in the envelope, excitation of acoustic modes, and the role of convective overshooting; 2) Applying the spherical 3-D MHD ASH (Anelastic Spherical Harmonics) code to simulate the core convection and radiative zone. Our goal is to determine whether core convection can excite low-frequency gravity modes, and thereby explain the presence of low frequencies for some hybrid gamma Doradus/delta Scuti variables for which the envelope convection zone is too shallow for the convective blocking mechanism to drive gravity modes; 3) Applying the ROTORC 2-D stellar evolution and dynamics code to calculate evolution with a variety of initial rotation rates and extents of core convective overshooting. The nonradial adiabatic pulsation frequencies of these nonspherical models are calculated using the 2-D pulsation code NRO. We present new insights into pulsations for stars of one to two solar masses gained by multidimensional modeling.
△ Less
Submitted 14 May, 2016;
originally announced May 2016.
-
Dynamics of Turbulent Convection and Convective Overshoot in a Moderate Mass Star
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Nagi N. Mansour,
Alan A. Wray
Abstract:
Continued progress in observational stellar astrophysics requires a deep understanding of the underlying convection dynamics. We present results of realistic 3D radiative hydrodynamic simulations of the outer layers of a moderate mass star (1.47 Msun), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface gr…
▽ More
Continued progress in observational stellar astrophysics requires a deep understanding of the underlying convection dynamics. We present results of realistic 3D radiative hydrodynamic simulations of the outer layers of a moderate mass star (1.47 Msun), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface granulation has a broad range of scales, from 2 to 12 Mm, and that large granules are organized in well-defined clusters, consisting of several granules. Comparison of the mean structure profiles from 3D simulations with the corresponding 1D standard stellar model shows an increase of the stellar radius by ~800 km, as well as significant changes in the thermodynamic structure and turbulent properties of the ionization zones. Convective downdrafts in the intergranular lanes between granulation clusters reach speeds of more than 20 km/s, penetrate through the whole convection zone, hit the radiative zone, and form a 8 Mm thick overshoot layer. Contrary to semi-empirical overshooting models, our results show that the 3D dynamic overshoot region consists of two layers: a nearly adiabatic extension of the convection zone and a deeper layer of enhanced subadiabatic stratification. This layer is formed because of heating caused by the braking of the overshooting convective plumes. This effect has to be taken into account in stellar modeling and the interpretation of asteroseismology data. In particular, we demonstrate that the 3D model can qualitatively explain deviations from the standard solar model found by helioseismology.
△ Less
Submitted 22 December, 2015;
originally announced December 2015.
-
Simulations of Stellar Magnetoconvection using the Radiative MHD Code `StellarBox'
Authors:
Alan A. Wray,
Khalil Bensassi,
Irina N. Kitiashvili,
Nagi N. Mansour,
Alexander G. Kosovichev
Abstract:
Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of their magnetic activity and variability. The goal of this paper is to present a detailed description and test results of a compressible radiative MHD code, `Stell…
▽ More
Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of their magnetic activity and variability. The goal of this paper is to present a detailed description and test results of a compressible radiative MHD code, `StellarBox', specifically developed for simulating the convection zones, surface, and atmospheres of the Sun and moderate-mass stars. The code solves the three-dimensional, fully coupled compressible MHD equations using a fourth-order Padé spatial differentiation scheme and a fourth-order Runge-Kutta scheme for time integration. The radiative transfer equation is solved using the Feautrier method for bi-directional ray tracing and an opacity-binning technique. A specific feature of the code is the implementation of subgrid-scale MHD turbulence models. The data structures are automatically configured, depending on the computational grid and the number of available processors, to achieve good load balancing. We present test results and illustrate the code's capabilities for simulating the granular convection on the Sun and a set of main-sequence stars. The results reveal substantial changes in the near-surface turbulent convection in these stars, which in turn affect properties of the surface magnetic fields. For example, in the solar case initially uniform vertical magnetic fields tend to self-organize into compact (pore-like) magnetic structures, while in more massive stars such structures are not formed and the magnetic field is distributed more-or-less uniformly in the intergranular lanes.
△ Less
Submitted 28 July, 2015;
originally announced July 2015.
-
Realistic modeling of local dynamo processes on the Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
Magnetic fields are usually observed in the quiet Sun as small-scale elements that cover the entire solar surface (the `salt and pepper' patterns in line-of-sight magnetograms). By using 3D radiative MHD numerical simulations we find that these fields result from a local dynamo action in the top layers of the convection zone, where extremely weak 'seed' magnetic fields (e.g., from a $10^{-6}$ G) c…
▽ More
Magnetic fields are usually observed in the quiet Sun as small-scale elements that cover the entire solar surface (the `salt and pepper' patterns in line-of-sight magnetograms). By using 3D radiative MHD numerical simulations we find that these fields result from a local dynamo action in the top layers of the convection zone, where extremely weak 'seed' magnetic fields (e.g., from a $10^{-6}$ G) can locally grow above the mean equipartition field, to a stronger than 2000~G field localized in magnetic structures. Our results reveal that the magnetic flux is predominantly generated in regions of small-scale helical downflows. We find that the local dynamo action takes place mostly in a shallow, about 500~km deep, subsurface layer, from which the generated field is transported into the deeper layers by convective downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and horizontal fields above the photosphere can be explained by small-scale magnetic loops produced by the dynamo. Such small-scale loops play an important role in the structure and dynamics of the solar atmosphere and that their detection in observations is critical for understanding the local dynamo action on the Sun.
△ Less
Submitted 29 June, 2015;
originally announced June 2015.
-
Using Realistic MHD Simulations for Modeling and Interpretation of Quiet-Sun Observations with the Solar Dynamics Observatory Helioseismic and Magnetic Imager
Authors:
Irina N. Kitiashvili,
Sebastien Couvidat,
Andreas Lagg
Abstract:
The solar atmosphere is extremely dynamic, and many important phenomena develop on small scales that are unresolved in observations with the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO). For correct calibration and interpretation of the observations, it is very important to investigate the effects of small-scale structures and dynamics on the HMI observ…
▽ More
The solar atmosphere is extremely dynamic, and many important phenomena develop on small scales that are unresolved in observations with the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO). For correct calibration and interpretation of the observations, it is very important to investigate the effects of small-scale structures and dynamics on the HMI observables, such as Doppler shift, continuum intensity, spectral line depth, and width. We use 3D radiative hydrodynamics simulations of the upper turbulent convective layer and the atmosphere of the Sun, and a spectro-polarimetric radiative transfer code to study observational characteristics of the Fe I 6173A line observed by HMI in quiet-Sun regions. We use the modeling results to investigate the sensitivity of the line Doppler shift to plasma velocity, and also sensitivities of the line parameters to plasma temperature and density, and determine effective line formation heights for observations of solar regions located at different distances from the disc center. These estimates are important for the interpretation of helioseismology measurements. In addition, we consider various center-to-limb effects, such as convective blue-shift, variations of helioseismic travel-times, and the 'concave' Sun effect, and show that the simulations can qualitatively reproduce the observed phenomena, indicating that these effects are related to a complex interaction of the solar dynamics and radiative transfer.
△ Less
Submitted 4 June, 2015; v1 submitted 9 July, 2014;
originally announced July 2014.
-
Spectro-Polarimetric Properties of Small-Scale Plasma Eruptions Driven by Magnetic Vortex Tubes
Authors:
Irina N. Kitiashvili
Abstract:
Highly turbulent nature of convection on the Sun causes strong multi-scale interaction of subsurface layers with the photosphere and chromosphere. According to realistic 3D radiative MHD numerical simulations ubiquitous small-scale vortex tubes are generated by turbulent flows below the visible surface and concentrated in the intergranular lanes. The vortex tubes can capture and amplify magnetic f…
▽ More
Highly turbulent nature of convection on the Sun causes strong multi-scale interaction of subsurface layers with the photosphere and chromosphere. According to realistic 3D radiative MHD numerical simulations ubiquitous small-scale vortex tubes are generated by turbulent flows below the visible surface and concentrated in the intergranular lanes. The vortex tubes can capture and amplify magnetic field, penetrate into chromospheric layers and initiate quasi-periodic flow eruptions that generates Alfvénic waves, transport mass and energy into the solar atmosphere. The simulations revealed high-speed flow patterns, and complicated thermodynamic and magnetic structures in the erupting vortex tubes. The spontaneous eruptions are initiated and driven by strong pressure gradients in the near-surface layers, and accelerated by the Lorentz force in the low chromosphere. In this paper, the simulation data are used to further investigate the dynamics of the eruptions, their spectro-polarimetric characteristics for the Fe I 6301.5 and 6302.5 A spectral lines, and demonstrate expected signatures of the eruptions in the Hinode SP data. We found that the complex dynamical structure of vortex tubes (downflows in the vortex core and upflows on periphery) can be captured by the Stokes I profiles. During an eruption, the ratio of down and upflows can suddenly change, and this effect can be observed in the Stokes V profile. Also, during the eruption the linear polarization signal increases, and this also can be detected with Hinode SP.
△ Less
Submitted 8 July, 2014;
originally announced July 2014.
-
Mechanism of local dynamo action on the Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
In the quiet Sun, magnetic fields are usually observed as small-scale magnetic elements, `salt and pepper', covering the entire solar surface. By using 3D radiative MHD numerical simulations we demonstrate that these fields are a result of local dynamo action in the top layers of the convection zone, where extremely weak `seed' magnetic fields can locally grow above the mean equipartition field (e…
▽ More
In the quiet Sun, magnetic fields are usually observed as small-scale magnetic elements, `salt and pepper', covering the entire solar surface. By using 3D radiative MHD numerical simulations we demonstrate that these fields are a result of local dynamo action in the top layers of the convection zone, where extremely weak `seed' magnetic fields can locally grow above the mean equipartition field (e.g., from a $10^{-6}$ G `seed' field to more than 1000 G magnetic structures). We find that the local dynamo action takes place only in a shallow, about 500 km deep, subsurface layer, from which the generated field is transported into deeper layers by convection downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and the horizontal fields above the photosphere can be explained by multi-scale magnetic loops produced by the dynamo.
△ Less
Submitted 3 December, 2013;
originally announced December 2013.
-
Ubiquitous Solar Eruptions Driven by Magnetized Vortex Tubes
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
S. K. Lele,
N. N. Mansour,
A. A. Wray
Abstract:
The solar surface is covered by high-speed jets transporting mass and energy into the solar corona and feeding the solar wind. The most prominent of these jets have been known as spicules. However, the mechanism initiating these eruptions events is still unknown. Using realistic numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by…
▽ More
The solar surface is covered by high-speed jets transporting mass and energy into the solar corona and feeding the solar wind. The most prominent of these jets have been known as spicules. However, the mechanism initiating these eruptions events is still unknown. Using realistic numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by the Sun's turbulent convection in subsurface layers. The swirling vortex tubes (resembling tornadoes) penetrate into the solar atmosphere, capture and stretch background magnetic field, and push surrounding material up, generating quasiperiodic shocks. Our simulations reveal a complicated high-speed flow patterns, and thermodynamic and magnetic structure in the erupting vortex tubes. We found that the eruptions are initiated in the subsurface layers and are driven by the high-pressure gradients in the subphotosphere and photosphere, and by the Lorentz force in the higher atmosphere layers.
△ Less
Submitted 31 December, 2012;
originally announced January 2013.
-
Turbulent Kinetic Energy Spectra of Solar Convection from NST Observations and Realistic MHD Simulations
Authors:
I. N. Kitiashvili,
V. I. Abramenko,
P. R. Goode,
A. G. Kosovichev,
S. K. Lele,
N. N. Mansour,
A. A. Wray,
V. B. Yurchyshyn
Abstract:
Turbulent properties of the quiet Sun represent the basic state of surface conditions, and a background for various processes of solar activity. Therefore understanding of properties and dynamics of this `basic' state is important for investigation of more complex phenomena, formation and development of observed phenomena in the photosphere and atmosphere. For characterization of the turbulent pro…
▽ More
Turbulent properties of the quiet Sun represent the basic state of surface conditions, and a background for various processes of solar activity. Therefore understanding of properties and dynamics of this `basic' state is important for investigation of more complex phenomena, formation and development of observed phenomena in the photosphere and atmosphere. For characterization of the turbulent properties we compare kinetic energy spectra on granular and sub-granular scales obtained from infrared TiO observations with the New Solar Telescope (Big Bear Solar Observatory) and from 3D radiative MHD numerical simulations ('SolarBox' code). We find that the numerical simulations require a high spatial resolution with 10 - 25 km grid-step in order to reproduce the inertial (Kolmogorov) turbulence range. The observational data require an averaging procedure to remove noise and potential instrumental artifacts. The resulting kinetic energy spectra show a good agreement between the simulations and observations, opening new perspectives for detailed joint analysis of more complex turbulent phenomena on the Sun, and possibly on other stars. In addition, using the simulations and observations we investigate effects of background magnetic field, which is concentrated in self-organized complicated structures in intergranular lanes, and find an increase of the small-scale turbulence energy and its decrease at larger scales due to magnetic field effects.
△ Less
Submitted 22 June, 2012;
originally announced June 2012.
-
Dynamics of Magnetized Vortex Tubes in the Solar Chromosphere
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
We use 3D radiative MHD simulations to investigate the formation and dynamics of small-scale (less than 0.5 Mm in diameter) vortex tubes spontaneously generated by turbulent convection in quiet-Sun regions with initially weak mean magnetic fields. The results show that the vortex tubes penetrate into the chromosphere and substantially affect the structure and dynamics of the solar atmosphere. The…
▽ More
We use 3D radiative MHD simulations to investigate the formation and dynamics of small-scale (less than 0.5 Mm in diameter) vortex tubes spontaneously generated by turbulent convection in quiet-Sun regions with initially weak mean magnetic fields. The results show that the vortex tubes penetrate into the chromosphere and substantially affect the structure and dynamics of the solar atmosphere. The vortex tubes are mostly concentrated in intergranular lanes and are characterized by strong (near sonic) downflows and swirling motions that capture and twist magnetic field lines, forming magnetic flux tubes that expand with height and which attain magnetic field strengths ranging from 200 G in the chromosphere to more than 1 kG in the photosphere. We investigate in detail the physical properties of these vortex tubes, including thermodynamic properties, flow dynamics, and kinetic and current helicities, and conclude that magnetized vortex tubes provide an important path for energy and momentum transfer from the convection zone into the chromosphere.
△ Less
Submitted 26 January, 2012;
originally announced January 2012.
-
Vortex Tubes of Turbulent Solar Convection
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
S. K. Lele,
A. A. Wray
Abstract:
Investigation of the turbulent properties of solar convection is extremely important for understanding the multi-scale dynamics observed on the solar surface. In particular, recent high-resolution observations have revealed ubiquitous vortical structures, and numerical simulations have demonstrated links between vortex tube dynamics and magnetic field organization and have shown the importance of…
▽ More
Investigation of the turbulent properties of solar convection is extremely important for understanding the multi-scale dynamics observed on the solar surface. In particular, recent high-resolution observations have revealed ubiquitous vortical structures, and numerical simulations have demonstrated links between vortex tube dynamics and magnetic field organization and have shown the importance of vortex tube interactions in the mechanisms of acoustic wave excitation on the Sun. In this paper we investigate the mechanisms of the formation of vortex tubes in highly-turbulent convective flows near the solar surface by using realistic radiative hydrodynamic LES simulations. Analysis of data from the simulations indicates two basic processes of vortex tube formation: 1) development of small-scale convective instability inside convective granules, and 2) a Kelvin-Helmholtz type instability of shearing flows in intergranular lanes. Our analysis shows that vortex stretching during these processes is a primary source of generation of small-scale vorticity on the Sun.
△ Less
Submitted 11 January, 2012; v1 submitted 27 December, 2011;
originally announced December 2011.
-
LES of turbulent convection in solar-type stars and formation of large-scale magnetic structures
Authors:
I. Rogachevskii,
N. Kleeorin,
I. N. Kitiashvili,
A. G. Kosovichev,
A. A. Wray,
N. N. Mansour
Abstract:
In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed f…
▽ More
In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed for simulating the upper solar convection zone and lower atmosphere. To this end we derived the energy budget equations, which include the effects of the subgrid-scale turbulence on the Lorentz-force, and implemented the new subgrid-scale turbulence model (TELF-Model) in a three-dimensional nonlinear MHD LES code. Using imposed initial vertical and horizontal uniform magnetic fields in LES with the TELF-Model, we have shown that the magnetic flux tubes formation is started when the initial mean magnetic field is larger than a threshold value (about 100 G). This is in agreement with the theoretical studies by Rogachevskii and Kleeorin (2007). We have determined the vertical profiles of the velocity and magnetic fluctuations, total MHD energy and anisotropy of turbulent magneto-convection, kinetic and current and cross helicities.
△ Less
Submitted 6 February, 2011;
originally announced February 2011.
-
Numerical MHD Simulations of Solar Magnetoconvection and Oscillations in Inclined Magnetic Field Regions
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed ou…
▽ More
The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed outflow in the penumbra. In this paper, we use radiative MHD simulations to investigate the influence of the magnetic field inclination on the power spectrum of vertical velocity oscillations. The results reveal a strong shift of the resonance mode peaks to higher frequencies in the case of a highly inclined magnetic field. The frequency shift for the inclined field is significantly greater than that in vertical field regions of similar strength. This is consistent with the behavior of fast MHD waves.
△ Less
Submitted 30 November, 2010; v1 submitted 24 November, 2010;
originally announced November 2010.
-
Excitation of acoustic waves by vortices in the quiet Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
Five-minutes oscillations is one of the basic properties of solar convection. Observations show mixture of a large number of acoustic wave fronts propagating from their sources. We investigate the process of acoustic waves excitation from the point of view of individual events, by using realistic 3D radiative hydrodynamic simulation of the quiet Sun. The results show that the excitation events are…
▽ More
Five-minutes oscillations is one of the basic properties of solar convection. Observations show mixture of a large number of acoustic wave fronts propagating from their sources. We investigate the process of acoustic waves excitation from the point of view of individual events, by using realistic 3D radiative hydrodynamic simulation of the quiet Sun. The results show that the excitation events are related to dynamics vortex tubes (or swirls) in the intergranular lanes. These whirlpool-like flows are characterized by very strong horizontal velocities (7 - 11 km/s) and downflows (~ 7 km/s), and are accompanied by strong decreases of the temperature, density and pressure at the surface and in a ~ 0.5-1 Mm deep layer below the surface. High-speed whirlpool flows can attract and capture other vortices. According to our simulation results, the processes of the vortex interaction, such as vortex annihilation, can cause the excitation of acoustic waves.
△ Less
Submitted 16 November, 2010;
originally announced November 2010.
-
A precise asteroseismic age and radius for the evolved Sun-like star KIC 11026764
Authors:
T. S. Metcalfe,
M. J. P. F. G. Monteiro,
M. J. Thompson,
J. Molenda-Zakowicz,
T. Appourchaux,
W. J. Chaplin,
G. Dogan,
P. Eggenberger,
T. R. Bedding,
H. Bruntt,
O. L. Creevey,
P. -O. Quirion,
D. Stello,
A. Bonanno,
V. Silva Aguirre,
S. Basu,
L. Esch,
N. Gai,
M. P. Di Mauro,
A. G. Kosovichev,
I. N. Kitiashvili,
J. C. Suarez,
A. Moya,
L. Piau,
R. A. Garcia
, et al. (33 additional authors not shown)
Abstract:
The primary science goal of the Kepler Mission is to provide a census of exoplanets in the solar neighborhood, including the identification and characterization of habitable Earth-like planets. The asteroseismic capabilities of the mission are being used to determine precise radii and ages for the target stars from their solar-like oscillations. Chaplin et al. (2010) published observations of thre…
▽ More
The primary science goal of the Kepler Mission is to provide a census of exoplanets in the solar neighborhood, including the identification and characterization of habitable Earth-like planets. The asteroseismic capabilities of the mission are being used to determine precise radii and ages for the target stars from their solar-like oscillations. Chaplin et al. (2010) published observations of three bright G-type stars, which were monitored during the first 33.5 days of science operations. One of these stars, the subgiant KIC 11026764, exhibits a characteristic pattern of oscillation frequencies suggesting that it has evolved significantly. We have derived asteroseismic estimates of the properties of KIC 11026764 from Kepler photometry combined with ground-based spectroscopic data. We present the results of detailed modeling for this star, employing a variety of independent codes and analyses that attempt to match the asteroseismic and spectroscopic constraints simultaneously. We determine both the radius and the age of KIC 11026764 with a precision near 1%, and an accuracy near 2% for the radius and 15% for the age. Continued observations of this star promise to reveal additional oscillation frequencies that will further improve the determination of its fundamental properties.
△ Less
Submitted 20 October, 2010;
originally announced October 2010.
-
Mechanism of spontaneous formation of stable magnetic structures on the Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
A. A. Wray,
N. N. Mansour
Abstract:
One of the puzzling features of solar magnetism is formation of long-living compact magnetic structures; such as sunspots and pores, in the highly turbulent upper layer of the solar convective zone. We use realistic radiative 3D MHD simulations to investigate the interaction between magnetic field and turbulent convection. In the simulations, a weak vertical uniform magnetic field is imposed in a…
▽ More
One of the puzzling features of solar magnetism is formation of long-living compact magnetic structures; such as sunspots and pores, in the highly turbulent upper layer of the solar convective zone. We use realistic radiative 3D MHD simulations to investigate the interaction between magnetic field and turbulent convection. In the simulations, a weak vertical uniform magnetic field is imposed in a region of fully developed granular convection; and the total magnetic flux through the top and bottom boundaries is kept constant. The simulation results reveal a process of spontaneous formation of stable magnetic structures, which may be a key to understanding of the magnetic self-organization on the Sun and formation of pores and sunspots. This process consists of two basic steps: 1) formation of small-scale filamentary magnetic structures associated with concentrations of vorticity and whirlpool-type motions, and 2) merging of these structures due to the vortex attraction, caused by converging downdrafts around magnetic concentration below the surface. In the resulting large-scale structure maintained by the converging plasma motions, the magnetic field strength reaches ~1.5 kG at the surface and ~6 kG in the interior; and the surface structure resembles solar pores. The magnetic structure remains stable for the whole simulation run of several hours with no sign of decay.
△ Less
Submitted 13 April, 2010;
originally announced April 2010.
-
Explanation of the sea-serpent magnetic structure of sunspot penumbrae
Authors:
I. N. Kitiashvili,
L. R. Bellot Rubio,
A. G. Kosovichev,
N. N. Mansour,
A. Sainz Dalda,
A. A. Wray
Abstract:
Recent spectro-polarimetric observations of a sunspot showed the formation of bipolar magnetic patches in the mid penumbra and their propagation toward the outer penumbral boundary. The observations were interpreted as being caused by sea-serpent magnetic fields near the solar surface (Sainz Dalda & Bellot Rubio 2008). In this Letter, we develop a 3D radiative MHD numerical model to explain the…
▽ More
Recent spectro-polarimetric observations of a sunspot showed the formation of bipolar magnetic patches in the mid penumbra and their propagation toward the outer penumbral boundary. The observations were interpreted as being caused by sea-serpent magnetic fields near the solar surface (Sainz Dalda & Bellot Rubio 2008). In this Letter, we develop a 3D radiative MHD numerical model to explain the sea-serpent structure and the wave-like behavior of the penumbral magnetic field lines. The simulations reproduce the observed behavior, suggesting that the sea-serpent phenomenon is a consequence of magnetoconvection in a strongly inclined magnetic field. It involves several physical processes: filamentary structurization, high-speed overturning convective motions in strong, almost horizontal magnetic fields with partially frozen field lines, and traveling convective waves. The results demonstrate a correlation of the bipolar magnetic patches with high-speed Evershed downflows in the penumbra. This is the first time that a 3D numerical model of the penumbra results in downward directed magnetic fields, an essential ingredient of sunspot penumbrae that has eluded explanation until now.
△ Less
Submitted 26 February, 2010;
originally announced March 2010.
-
The asteroseismic potential of Kepler: first results for solar-type stars
Authors:
W. J. Chaplin,
T. Appourchaux,
Y. Elsworth,
R. A. Garcia,
G. Houdek,
C. Karoff,
T. S. Metcalfe,
J. Molenda-Zakowicz,
M. J. P. F. G. Monteiro,
M. J. Thompson,
T. M. Brown,
J. Christensen-Dalsgaard,
R. L. Gilliland,
H. Kjeldsen,
W. J. Borucki,
D. Koch,
J. M. Jenkins,
J. Ballot,
S. Basu,
M. Bazot,
T. R. Bedding,
O. Benomar,
A. Bonanno,
I. M. Brandao,
H. Bruntt
, et al. (83 additional authors not shown)
Abstract:
We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies a…
▽ More
We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe.
△ Less
Submitted 18 January, 2010; v1 submitted 4 January, 2010;
originally announced January 2010.
-
Traveling Waves of Magnetoconvection and the Origin of the Evershed Effect in Sunspots
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray,
Nagi N. Mansour
Abstract:
Discovered in 1909 the Evershed effect represents strong mass outflows in sunspot penumbra, where the magnetic field of sunspots is filamentary and almost horizontal. These flows play important role in sunspots and have been studied in detail using large ground-based and space telescopes, but the basic understanding of its mechanism is still missing. We present results of realistic numerical sim…
▽ More
Discovered in 1909 the Evershed effect represents strong mass outflows in sunspot penumbra, where the magnetic field of sunspots is filamentary and almost horizontal. These flows play important role in sunspots and have been studied in detail using large ground-based and space telescopes, but the basic understanding of its mechanism is still missing. We present results of realistic numerical simulations of the Sun's subsurface dynamics, and argue that the key mechanism of this effect is in non-linear magnetoconvection that has properties of traveling waves in the presence of strong, highly inclined magnetic field. The simulations reproduce many observed features of the Evershed effect, including the high-speed "Evershed clouds", the filamentary structure of the flows, and the non-stationary quasi-periodic behavior. The results provide a synergy of previous theoretical models and lead to an interesting prediction of a large-scale organization of the outflows.
△ Less
Submitted 23 April, 2009;
originally announced April 2009.
-
Prediction of Sunspot Cycles by Data Assimilation Method
Authors:
I. N. Kitiashvili,
A. G. Kosovichev
Abstract:
Despite the known general properties of the solar cycles, a reliable forecast of the 11-year sunspot number variations is still a problem. The difficulties are caused by the apparent chaotic behavior of the sunspot numbers from cycle to cycle and by the influence of various turbulent dynamo processes, which are far from understanding. For predicting the solar cycle properties we make an initial…
▽ More
Despite the known general properties of the solar cycles, a reliable forecast of the 11-year sunspot number variations is still a problem. The difficulties are caused by the apparent chaotic behavior of the sunspot numbers from cycle to cycle and by the influence of various turbulent dynamo processes, which are far from understanding. For predicting the solar cycle properties we make an initial attempt to use the Ensemble Kalman Filter (EnKF), a data assimilation method, which takes into account uncertainties of a dynamo model and measurements, and allows to estimate future observational data. We present the results of forecasting of the solar cycles obtained by the EnKF method in application to a low-mode nonlinear dynamical system modeling the solar $αΩ$-dynamo process with variable magnetic helicity. Calculations of the predictions for the previous sunspot cycles show a reasonable agreement with the actual data. This forecast model predicts that the next sunspot cycle will be significantly weaker (by $\sim 30%$) than the previous cycle, continuing the trend of low solar activity.
△ Less
Submitted 11 February, 2009; v1 submitted 21 July, 2008;
originally announced July 2008.
-
Nonlinear Solar Dynamo Model with Magnetic Helicity
Authors:
I. N. Kitiashvili,
A. G. Kosovichev
Abstract:
Variations of the sunspot number are important indicators of the solar activity cycles. The sunspot formation is a result of a dynamo process inside the Sun, which is far from being understood. We use simple dynamical models of the dynamo process to simulate the magnetic field evolution and investigate general properties of the sunspot number variations during the solar cycles. We have found tha…
▽ More
Variations of the sunspot number are important indicators of the solar activity cycles. The sunspot formation is a result of a dynamo process inside the Sun, which is far from being understood. We use simple dynamical models of the dynamo process to simulate the magnetic field evolution and investigate general properties of the sunspot number variations during the solar cycles. We have found that the classical Parker's model with a standard kinetic helicity quenching cannot represent the typical profiles of the solar-cycle variations of the sunspot number, and also does not give chaotic solutions. For modeling of the solar cycle properties we use a nonlinear dynamo model of Kleeorin and Ruzmaikin (1982), which takes into account dynamics of the turbulent magnetic helicity. We have obtained a series of periodic and chaotic solutions for different layers of the convective zone. The solutions qualitatively reproduce some basic observational features of the solar cycle properties, in particular, the relationship between the growth time and the cycle amplitude. Also, on the longer time scale the dynamo model with the magnetic helicity has intermittent solutions, which may be important for modeling long-term variations of the solar cycles.
△ Less
Submitted 11 February, 2009; v1 submitted 20 July, 2008;
originally announced July 2008.