Modes of magnetic field generation in the low-mode αΩ-dynamo model with dynamic regulation of the α-effect by the field energy

В работе используется маломодовая модель αΩ-динамо для моделирования режимов генерации магнитного поля при незначительных изменениях поля скорости вязкой жидкости. В рамках этой модели интенсивность α-эффекта регулируется процессом с памятью, который вводится в магнитогидродинамическую систему (МГД-система) как аддитивная поправка в виде функционала Z(t) от энергии поля. В качестве ядра J(t) функционала Z(t) выбрана функция, определяющая затухающие колебания с варьируемым коэффициентом затухания и постоянной частотой затухания, принятой равной единице. Исследование поведения магнитного поля проводится на больших временных масштабах, поэтому для численных расчётов используется перемасштабированная и обезразмеренная МГД-система, где в качестве единицы времени принято время диссипации магнитного поля (104 лет). Управляющими параметрами системы выступают число Рейнольдса и амплитуда α-эффекта, в которых заложена информация о крупномасштабном и турбулентном генераторах. Результаты численного моделирования режимов генерации магнитного поля при различных значениях коэффициента затухания и постоянной частоте затухания отражены на фазовой плоскости управляющих параметров. В работе исследуется вопрос о динамике изменения картины на фазовой плоскости в зависимости от значения коэффициента затухания. Проводится сравнение с результатами, полученными ранее при постоянной интенсивности α-эффекта и при изменении интенсивности α — эффекта, которое определялось функционалом Z(t) с показательным ядром и аналогичными значениями коэффициента затухания. In this paper, we use a low-mode αΩ-dynamo model to simulate the modes of magnetic field generation with insignificant changes in the velocity field of a viscous fluid. Within the framework of this model, an additive correction is introduced into the magnetohydrodynamic system to control the intensity of the α-effect in the form of a function Z(t) from the field energy. As the kernel J(t) of the function Z(t) is chosen the function that determines damped oscillations with the different values of the damping coefficient and a constant damping frequency taken equal to one. The study of the magnetic field behavior is carried out on a large time scales, therefore, for numerical calculations, a rescaled and dimensionless MHD-system is used, where the time of the magnetic field dissipation (104 years) is accepted as the unit of time. The main parameters of the system are the Reynolds number and the amplitude of the α-effect, which contains information about the large-scale and turbulent generators, respectively. According to the results of numerical simulation, an increase in the values of the damping coefficient is characterized an increase in the inhibition effect of the process Z(t) on the α-effect and decrease of the magnetic field divergence region on the plane of the main parameters.

Subcritical dynamo and hysteresis in a Babcock-Leighton type kinematic dynamo model

In the Sun and Sun-like stars, it is believed that cycles of the large-scale magnetic field are produced due to the existence of differential rotation and helicity in the plasma flows in their convection zones (CZs). Hence, it is expected that for each star, there is a critical dynamo number for the operation of a large-scale dynamo. As a star slows down, it is expected that the large-scale dynamo ceases to operate above a critical rotation period. In our study, we explore the possibility of the operation of the dynamo in the subcritical region using the Babcock–Leighton type kinematic dynamo model. In some parameter regimes, we find that the dynamo shows hysteresis behavior, i.e., two dynamo solutions are possible depending on the initial parameters—decaying solution if starting with weak field and strong oscillatory solution (subcritical dynamo) when starting with a strong field. However, under large fluctuations in the dynamo parameter, the subcritical dynamo mode is unstable in some parameter regimes. Therefore, our study supports the possible existence of subcritical dynamo in some stars which was previously demonstrated in a mean-field dynamo model with distributed α and MHD turbulent dynamo simulations.

Modeling of solar magnetic field using the kinematic-gravitational ion dynamo model

In present article the kinematic-gravitational ion dynamo model accounting for influence of tidal forces on electric currents in ionized substances is applied to modeling of the magnetic field of the Sun. Estimates of currents and field values obtained using a seven-layer model indicate that tidal forces influence is not insignificant. A correlation method for assessment of the Sun’s polarity was created and applied for a detailed analysis of the polarity of magnetic field of the Sun in the 21 and 22 cycles.

Solar evolution and extrema: current state of understanding of long-term solar variability and its planetary impacts

The activity of stars such as the Sun varies over timescales ranging from the very short to the very long—stellar and planetary evolutionary timescales. Experience from our solar system indicates that short-term, transient events such as stellar flares and coronal mass ejections create hazardous space environmental conditions that impact Earth-orbiting satellites and planetary atmospheres. Extreme events such as stellar superflares may play a role in atmospheric mass loss and create conditions unsuitable for life. Slower, long-term evolutions of the activity of Sun-like stars over millennia to billions of years result in variations in stellar wind properties, radiation flux, cosmic ray flux, and frequency of magnetic storms. This coupled evolution of star-planet systems eventually determines planetary and exoplanetary habitability. The Solar Evolution and Extrema (SEE) initiative of the Variability of the Sun and Its Terrestrial Impact (VarSITI) program of the Scientific Committee on Solar-Terrestrial Physics (SCOSTEP) aimed to facilitate and build capacity in this interdisciplinary subject of broad interest in astronomy and astrophysics. In this review, we highlight progress in the major themes that were the focus of this interdisciplinary program, namely, reconstructing and understanding past solar activity including grand minima and maxima, facilitating physical dynamo-model-based predictions of future solar activity, understanding the evolution of solar activity over Earth’s history including the faint young Sun paradox, and exploring solar-stellar connections with the goal of illuminating the extreme range of activity that our parent star—the Sun—may have displayed in the past, or may be capable of unleashing in the future.

Hereditary low-mode dynamo model

В данной статье рассматривается модель динамо в виде двумерной динамической системы в интегро-дифференциальной форме. В модели реализован стабилизирующий генерацию поля механизм обратной связи в виде подавления α-эффекта функционалом сверточного типа от актуальных и предыдущих значений спиральности и энергии. Наличие этого механизма подавления вводит в модель эредитарность (память). Для модели была построена численная схема ввиде совмещение разностных схем для дифференциальной и интегральной части, двухступенчатый неявный методы Рунге-Кутты и метод трапеций соотвественно. Так же были рассмотрены и графически представлены динамические режимы нашей модели. This article discusses a dynamo model in the form of a two-dimensional dynamical system in integro-differential form. The model implements a stabilizing polarization generator in the form of suppression of the a effect of convolutional type functional from current and previous helicity and energy values. The presence of this suppression mechanism introduces hereditarity (memory) into the model. For modeling, a digital scheme was constructed in the form of a combination of difference schemes for the differential and integral parts, a twostep implicit Runge-Kutta method and a trapezium method, respectively. We also reviewed and graphically presented the dynamic modes of our model.

A secular variation candidate model for IGRF-13 based on Swarm data and ensemble inverse geodynamo modelling

This paper describes the design of a candidate secular variation model for the 13th generation of the International Geomagnetic Reference Field. This candidate is based upon the integration of an ensemble of 100 numerical models of the geodynamo between epochs 2019.0 and 2025.0. The only difference between each ensemble member lies in the initial condition that is used for the numerical integration, all other control parameters being fixed. An initial condition is defined as follows: an estimate of the magnetic field and its rate-of-change at the core surface for 2019.0 is obtained from a year (2018.5–2019.5) of vector Swarm data. This estimate (common to all ensemble members) is subject to prior constraints: the statistical properties of the numerical dynamo model for the main geomagnetic field and its secular variation, and prescribed covariances for the other sources. One next considers 100 three-dimensional core states (in terms of flow, buoyancy and magnetic fields) extracted at different discrete times from a dynamo simulation that is not constrained by observations, with the time distance between each state exceeding the dynamo decorrelation time. Each state is adjusted (in three dimensions) in order to take the estimate of the geomagnetic field and its rate-of-change for 2019.0 into account. This methodology provides 100 different initial conditions for subsequent numerical integration of the dynamo model up to epoch 2025.0. Focussing on the 2020.0–2025.0 time window, we use the median average rate-of-change of each Gauss coefficient of the ensemble and its statistics to define the geomagnetic secular variation over that time frame and its uncertainties.

The magnetic helicity density patterns from non-axisymmetric solar dynamo

We study the helicity density patterns which can result from the emerging bipolar regions. Using the relevant dynamo model and the magnetic helicity conservation law we find that the helicity density patterns around the bipolar regions depend on the configuration of the ambient large-scale magnetic field, and in general they show a quadrupole distribution. The position of this pattern relative to the equator can depend on the tilt of the bipolar region. We compute the time–latitude diagrams of the helicity density evolution. The longitudinally averaged effect of the bipolar regions shows two bands of sign for the density distributions in each hemisphere. Similar helicity density patterns are provided by the helicity density flux from the emerging bipolar regions subjected to surface differential rotation.

MEDIUM-TERM OSCILLATIONS OF THE SOLAR ACTIVITY

In addition to the well-known 11-year cycle, longer and shorter characteristic periods can be isolated in variations of the parameters of helio-geophysical activity. Periods of about 36 and 60 years were revealed in variations of the geomagnetic activity and an approximately 60-year periodicity, in the evolution of correlation between the pressure in the lower atmosphere and the solar activity. Similar periods are observed in the cyclonic activity. Such periods in the parameters of the solar activity are difficult to identify because of a limited database available; however, they are clearly visible in variations of the asymmetry of the sunspot activity in the northern and southern solar hemispheres. In geomagnetic variations, one can also isolate oscillations with the characteristic periods of 5-6 years (QSO) and 2-3 years (QBO). We have considered 5-6-year periodicities (about half the main cycle) observed in variations of the sunspot numbers and the intensity of the dipole component of the solar magnetic field. A comparison with different magnetic dynamo models allowed us to determine the possible origin of these oscillations. A similar result can be reproduced in a dynamo model with nonlinear parameter variations. In this case, the activity cycle turns out to be anharmonic and contains other periodicities in addition to the main one. As a result of the study, we conclude that the 5-6-year activity variations are related to the processes of nonlinear saturation of the dynamo in the solar interior. Quasi-biennial oscillations are actually separate pulses related little to each other. Therefore, the methods of the spectral analysis do not reveal them over large time intervals. They are a direct product of local fields, are generated in the near-surface layers, and are reliably recorded only in the epochs of high solar activity.