Source author record

Alexander G. Kosovichev

Alexander G. Kosovichev appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2026arXiv

Helioseismic Evidence That the Solar Dynamo Originates near the Tachocline

The exact location of the solar dynamo remains uncertain--whether it operates primarily in the near-surface shear layer, throughout the entire convection zone, or near the tachocline, a region of sharp transition in the solar rotation, located at the base of the convection zone, approximately 200,000 km beneath the surface. Various studies have supported each of these possibilities. Notably, the solar magnetic "butterfly" diagram and the pattern of zonal flows ("torsional oscillations") exhibit strikingly similar characteristics, suggesting a link between magnetic field evolution and solar flows. Since magnetic fields cannot be measured directly in the deep solar interior, torsional oscillations and rotation gradients are employed as diagnostic proxies. Our analysis reveals that the gradient of rotation displays "butterfly"-like behavior near the tachocline, which is similar to the magnetic butterfly diagram at the surface. This result supports the idea that the solar dynamo has a deep-seated origin, likely operating either near the tachocline or throughout the convection zone, thereby disfavoring the recent scenario of a shallow, near-surface dynamo. This finding may also have important implications for understanding how stellar dynamos operate in general.

preprint2026arXiv

Simulations of Global Solar Convection with a Fully Compressible CHORUS++ Code

Fluid-dynamics simulations of global solar convection are a critically important tool for assessing the dynamics of the solar interior. However, simulation studies with a fully compressible hydrodynamics code are not yet common. The CHORUS++ code solves robustly and efficiently the fully compressible hydrodynamics equations using a compact local spectral method and semi-unstructured grid system. Using the CHORUS++ code, we simulate the solar interior plasma flows from 0.7 to 0.99 of the solar radius using the values of the solar total luminosity and the sidereal rotation rate. In this paper, we analyze the simulated global flow structures before the statistically stable state to assess the compressibility of the plasma flows obtained with the fully compressible hydrodynamic code. The divergence of mass flux with the compressible model is overall small in the examined state before reaching the equilibrium state of the convection, which implies that the differences between the fully compressible flows and those obtained with the anelastic, incompressible, or linear-equation models are small in the simulated inner part of the convection zone. Although a fully-relaxed stationary convection has not been achieved yet in the examined state, the model qualitatively reproduces the solar-type differential rotation. Simulations for longer periods are needed to achieve the system relaxation state. This work, assessing the early phase of the hydrodynamic evolution of solar convection, is our first step toward a better understanding of the nature of solar convection and the dynamo processes.

preprint2025arXiv

3D Radiative MHD Modeling of Particle Beam Heating of the Solar Atmosphere

While solar flares are primarily associated with enhanced ultraviolet and X-ray emission, a subset of flares exhibit significant continuum brightening in visible light and are classified as white-light flares (WLFs). Despite extensive observational and modeling efforts, the physical mechanisms responsible for the compact, short-lived photospheric brightenings in WLF kernels observed during the impulsive phase of solar flares remain uncertain. Thick-target electron-beam models typically deposit energy in the upper chromosphere, and their ability to reproduce the magnitude and spatial localization of photospheric continuum enhancements observed in white-light flare kernels remains an open question. To investigate the role of realistic atmospheric structuring and multidimensional transport in flare energy deposition, we perform three-dimensional radiative MHD simulations of electron-beam heating using the StellarBox code for beam fluxes of $10^{12}$ erg\,s$^{-1}$\,cm$^{-2}$ and low-energy cutoffs of 10--25\,keV. We then compute Fe\,I 6173\,Å~Stokes profiles using the RH 1.5D radiative transfer code for direct comparison with Helioseismic and Magnetic Imager (HMI) observations. The simulations produce strong upper-chromospheric heating, multiple shock fronts, and continuum enhancements up to a factor of 2.5 relative to pre-flare levels, comparable to continuum enhancements observed during strong X-class white-light flares. Comparison with one-dimensional RADYN simulations highlights the influence of fine-scale structuring on flare dynamics and continuum emission that arises in three-dimensional geometry.

preprint2022arXiv

Constraining Global Solar Models through Helioseismic Analysis

Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this controversy and develop a more robust understanding of global flow regimes in the solar interior, we apply a "forward-modeling" approach to the analysis of helioseismic signatures of meridional circulation profiles obtained from numerical simulations. We employ the global acoustic modeling code GALE to simulate the propagation of acoustic waves through regimes of mean mass flows generated by global hydrodynamic and magnetohydrodynamic models: EULAG, the Pencil Code, and the Rayleigh code. These models are used to create synthetic dopplergram data products, used as inputs for local time-distance helioseismology techniques. Helioseismic travel-time signals from solutions obtained through global numerical simulations are compared directly with inferences from solar observations, in order to set additional constraints on global model parameters in a direct way. We show that even though these models are able to replicate solar-like differential rotation, the resulting rotationally-constrained convection develops a multi-cell global meridional circulation profile that is measurably inconsistent with local time-distance inferences of solar observations. However, we find that the development of rotationally-unconstrained convection close to the model surface is able to maintain solar-like differential rotation, while having a significant impact on the helioseismic travel-time signal, replicating solar observations within one standard deviation of the error due to noise.

preprint2022arXiv

Detection of Travel Time Anisotropy from Subsurface Horizontal Magnetic Fields

A time-distance measurement technique is derived to isolate phase travel time anisotropy caused by subsurface horizontal magnetic fields, and a method which uses the measured anisotropy to estimate the field's orientation is also derived. A simulation of acoustic waves propagating in a uniform, inclined magnetic field with solar background structure is used to verify the derived technique. Then, the procedure is applied to a numerical simulation of a sunspot, for which the subsurface state is known, to provide context for the results obtained from the study of several sunspots observed by the Helioseismic and Magnetic Imager. Significant anisotropies are detected, on the order of one minute, and the subsurface field's azimuth is estimated and compared with the azimuth of the surface magnetic field. In all cases, the subsurface azimuth is found to be well-aligned with that of the surface, and the results from the numerical simulation are used to interpret features in the detected travel time anisotropy.

preprint2022arXiv

Forward Modeling Helioseismic Signatures of One- and Two-cell Meridional Circulation

Using a 3D global solver of the linearized Euler equations, we model acoustic oscillations over background velocity flow fields of single-cell meridional circulation with deep and shallow return flows as well as double-cell meridional circulation with strong and weak reversals. The velocities are generated using a mean-field hydrodynamic and dynamo model -- moving through the regimes with minimal parameter changes; counter-rotation near the base of the tachocline is induced by sign inversion of the non-diffusive action of turbulent Reynolds stresses ($Λ$-effect) due to the radial inhomogeneity of the Coriolis number. By mimicking the stochastic excitation of resonant modes in the convective interior, we simulate realization noise present in solar observations. Using deep-focusing to analyze differences in travel-time signatures between the four regimes, as well as comparing to solar observations, we show that current helioseismology techniques may offer important insights about the location and strength of the return flow, however, that it may not currently be possible to definitively distinguish between profiles of single-cell or double-cell meridional circulation.

preprint2022arXiv

Spatial Spectrum of Solar Convection from Helioseismic Data: Flow Scales and Time Variations

We analyze spectral properties of solar convection in the range of depths from 0 to 19~Mm using subsurface flow maps obtained by the time-distance heiioseismology analysis of solar-oscillation data from the Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO) from May 2010 to September 2020. The results reveal a rapid increase of the horizontal flow scales with the depth, from supergranulation to giant-cell scales, and support the evidence of large-scale convection, previously detected by tracking the motion of supergranular cells on the surface. The total power of convective flows correlates with the solar activity cycle. During the solar maximum, the total power decreases in shallow subsurface layers and increases in the deeper layers.

preprint2021arXiv

Evolution of Subsurface Zonal and Meridional Flows in Solar Cycle 24 from Helioseismological Data

The results of determinations of the azimuthal and meridional velocities by time-distance helioseismology from Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO) from May 2010 to September 2020 at latitudes from -60° to +60° and depths to about 19 Mm below the photosphere are used to analyze spatiotemporal variations of the solar differential rotation and meridional circulation. The pattern of torsional oscillations, or latitudinal belts of alternating `fast' and `slow' zonal flows migrating from high latitudes towards the equator, is found to extend in the time--latitude diagrams over the whole time interval. The oscillation period is comparable with a doubled solar-activity-cycle and can be described as an extended solar cycle. The zonal-velocity variations are related to the solar-activity level, the local-velocity increases corresponding to the sunspot-number increases and being localized at latitudes where the strongest magnetic fields are recorded. The dramatic growth of the zonal velocities in 2018 appears to be a precursor of the beginning of activity Cycle 25. The strong symmetrization of the zonal-velocity field by 2020 can be considered another precursor. The general pattern of poleward meridional flows is modulated by latitudinal variations that are similar to the extended-solar-cycle behavior of the zonal flows. During the activity maximum, these variations are superposed with a higher harmonic corresponding to meridional flows converging to the spot-formation latitudes. Our results indicate that variations of both the zonal and meridional flows exhibit the extended solar-cycle behavior, which is an intrinsic feature of the solar dynamo.

preprint2020arXiv

Estimation of Key Sunquake Parameters through Hydrodynamic Modeling and Cross-Correlation Analysis

Sunquakes are one of the more distinct secondary phenomena related to solar flares, where energy deposition in the lower layers of the Sun's atmosphere excites acoustic waves easily visible in Helioseismic and Magnetic Imager (HMI) dopplergrams. We explore two possible sources of sunquakes in the context of the electron beam hypothesis: an instantaneous transfer of momentum and a gradual applied force due to flare eruption. We model the sunquake excitation and compare with five observed sunquake events using a cross-correlation analysis. We find that at least half the events studied are consistent with the electron beam hypothesis and estimate the energy required to excite the sunquakes to be within the range determined by previous studies.

preprint2020arXiv

Helioseismic Modeling of Background Flows

We present a 3-dimensional (3D) numerical solver of the linearized compressible Euler equations (GALE -- Global Acoustic Linearized Euler), used to model acoustic oscillations throughout the solar interior. The governing equations are solved in conservation form on a fully global spherical mesh ($0 \le ϕ\le 2π$, $0 \le θ\le π$, $0 \le r \le R_{\odot}$) over a background state generated by the standard Solar Model S. We implement an efficient pseudo-spectral computational method to calculate the contribution of the compressible material derivative dyad to internal velocity perturbations, computing oscillations over arbitrary 3D background velocity fields. This model offers a foundation for a "forward-modeling" approach, using helioseismology techniques to explore various regimes of internal mass flows. We demonstrate the efficacy of the numerical method presented in this paper by reproducing observed solar power spectra, showing rotational splitting due to differential rotation, and applying local helioseismology techniques to measure travel times created by a simple model of single-cell meridional circulation.

preprint2016arXiv

A brief history of the solar diameter measurements: a critical quality assessment of the existing data

The size of the diameter of the Sun has been debated for a very long time. First tackled by the Greek astronomers from a geometric point of view, an estimate, although incorrect, has been determined, not truly called into question for several centuries. The French school of astronomy, under the impetus of Mouton and Picard in the XVIIth century can be considered as a pioneer in this issue. It was followed by the German school at the end of the XIXth century whose works led to a canonical value established at 959".63 (second of arc). A number of ground-based observations has been made in the second half of the XIXth century leading to controversial results mainly due to the difficulty to disentangle between the solar and atmospheric effects. Dedicated space measurements yield to a very faint dependence of the solar diameter with time. New studies over the entire radiation spectrum lead to a clear relationship between the solar diameter and the wavelength, reflecting the height at which the lines are formed. Thus the absolute value of the solar diameter, which is a reference for many astrophysical applications, must be stated according to the wavelength. Furthermore, notable features of the Near Sub-Surface Layer (NSSL), called the leptocline, can be established in relation to the solar limb variations, mainly through the shape asphericities coefficients. The exact relationship has not been established yet, but recent studies encourage further in-depth investigations of the solar subsurface dynamics, both observationally and by numerical MHD simulations.

preprint2016arXiv

Local Helioseismology of Emerging Active Regions: A Case Study

Local helioseismology provides a unique opportunity to investigate the subsurface structure and dynamics of active regions and their effect on the large-scale flows and global circulation of the Sun. We use measurements of plasma flows in the upper convection zone, provided by the Time-Distance Helioseismology Pipeline developed for analysis of solar oscillation data obtained by Helioseismic and Magnetic Imager (HMI) on Solar Dynamics Observatory (SDO), to investigate the subsurface dynamics of emerging active region NOAA 11726. The active region emergence was detected in deep layers of the convection zone about 12 hours before the first bipolar magnetic structure appeared on the surface, and 2 days before the emergence of most of the magnetic flux. The speed of emergence determined by tracking the flow divergence with depth is about 1.4 km/s, very close to the emergence speed in the deep layers. As the emerging magnetic flux becomes concentrated in sunspots local converging flows are observed beneath the forming sunspots. These flows are most prominent in the depth range 1-3 Mm, and remain converging after the formation process is completed. On the larger scale converging flows around active region appear as a diversion of the zonal shearing flows towards the active region, accompanied by formation of a large-scale vortex structure. This process occurs when a substantial amount of the magnetic flux emerged on the surface, and the converging flow pattern remains stable during the following evolution of the active region. The Carrington synoptic flow maps show that the large-scale subsurface inflows are typical for active regions. In the deeper layers (10-13 Mm) the flows become diverging, and surprisingly strong beneath some active regions. In addition, the synoptic maps reveal a complex evolving pattern of large-scale flows on the scale much larger than supergranulation.

preprint2016arXiv

Reconstruction of Solar Subsurfaces by Local Helioseismology

Local helioseismology has opened new frontiers in our quest for understanding of the internal dynamics and dynamo on the Sun. Local helioseismology reconstructs subsurface structures and flows by extracting coherent signals of acoustic waves traveling through the interior and carrying information about subsurface perturbations and flows, from stochastic oscillations observed on the surface. The initial analysis of the subsurface flow maps reconstructed from the 5 years of SDO/HMI data by time-distance helioseismology reveals the great potential for studying and understanding of the dynamics of the quiet Sun and active regions, and the evolution with the solar cycle. In particular, our results show that the emergence and evolution of active regions are accompanied by multi-scale flow patterns, and that the meridional flows display the North-South asymmetry closely correlating with the magnetic activity. The latitudinal variations of the meridional circulation speed, which are probably related to the large-scale converging flows, are mostly confined in shallow subsurface layers. Therefore, these variations do not necessarily affect the magnetic flux transport. The North-South asymmetry is also pronounced in the variations of the differential rotation ("torsional oscillations"). The calculations of a proxy of the subsurface kinetic helicity density show that the helicity does not vary during the solar cycle, and that supergranulation is a likely source of the near-surface helicity.

preprint2015arXiv

Detection of Fast-Moving Waves Propagating Outward along Sunspots' Radial Direction in the Photosphere

Helioseismic and magnetohydrodynamic waves are abundant in and above sunspots. Through cross-correlating oscillation signals in the photosphere observed by the SDO/HMI, we reconstruct how waves propagate away from virtual wave sources located inside a sunspot. In addition to the usual helioseismic wave, a fast-moving wave is detected traveling along the sunspot's radial direction from the umbra to about 15 Mm beyond the sunspot boundary. The wave has a frequency range of 2.5 - 4.0 mHz with a phase velocity of 45.3 km/s, substantially faster than the typical speeds of Alfven and magnetoacoustic waves in the photosphere. The observed phenomenon is consistent with a scenario of that a magnetoacoustic wave is excited at approximately 5 Mm beneath the sunspot, and its wavefront travels to and sweeps across the photosphere with a speed higher than the local magnetoacoustic speed. The fast-moving wave, if truly excited beneath the sunspot's surface, will help open a new window to study the internal structure and dynamics of sunspots.

preprint2015arXiv

Dynamics of Turbulent Convection and Convective Overshoot in a Moderate Mass Star

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.

preprint2015arXiv

Dynamo Effects Near The Transition from Solar to Anti-Solar Differential Rotation

Numerical MHD simulations play increasingly important role for understanding mechanisms of stellar magnetism. We present simulations of convection and dynamos in density-stratified rotating spherical fluid shells. We employ a new 3D simulation code for the solution of a physically consistent anelastic model of the process with a minimum number of parameters. The reported dynamo simulations extend into a "buoyancy-dominated" regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients and strong anti-solar differential rotation develops as a result. We find that the self-generated magnetic fields, despite being relatively weak, are able to reverse the direction of differential rotation from anti-solar to solar-like. We also find that convection flows in this regime are significantly stronger in the polar regions than in the equatorial region, leading to non-oscillatory dipole-dominated dynamo solutions, and to concentration of magnetic field in the polar regions. We observe that convection has different morphology in the inner and at the outer part of the convection zone simultaneously such that organized geostrophic convection columns are hidden below a near-surface layer of well-mixed highly-chaotic convection. While we focus the attention on the buoyancy-dominated regime, we also demonstrate that conical differential rotation profiles and persistent regular dynamo oscillations can be obtained in the parameter space of the rotation-dominated regime even within this minimal model.

preprint2015arXiv

Properties of chromospheric evaporation and plasma dynamics of a solar flare from IRIS observations

Dynamics of hot chromospheric plasma of solar flares is a key to understanding of mechanisms of flare energy release and particle acceleration. A moderate M1.0 class flare of 12 June, 2014 (SOL2014-06-12T21:12) was simultaneously observed by NASA's Interface Region Imaging Spectrograph (IRIS), other spacecraft, and also by New Solar Telescope (NST) at the BBSO. This paper presents the first part of our investigation focused on analysis of the IRIS data. Our analysis of the IRIS data in different spectral lines reveals strong redshifted jet-like flow with the speed of ~100 km/s of the chromospheric material before the flare. Strong nonthermal emission of the C II k 1334.5 A line, formed in the chromosphere-corona transition region, is observed at the beginning of the impulsive phase in several small (with a size of ~1 arcsec) points. It is also found that the C II k line is redshifted across the flaring region before, during and after the impulsive phase. A peak of integrated emission of the hot (1.1 MK) plasma in the Fe XXI 1354.1 A line is detected approximately 5 minutes after the integrated emission peak of the lower temperature C II k. A strong blueshift of the Fe XXI line across the flaring region corresponds to evaporation flows of the hot chromospheric plasma with a speed of 50 km/s. Additional analysis of the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) data supports the idea that the upper chromospheric dynamics observed by IRIS has features of "gentle" evaporation driven by heating of the solar chromosphere by accelerated electrons and by a heat flux from the flare energy release site.

preprint2015arXiv

Simulations of Stellar Magnetoconvection using the Radiative MHD Code `StellarBox'

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.

preprint2014arXiv

Helioseismic Constraints and Paradigm Shift in Solar Dynamo

Helioseismology provides important constraints for the solar dynamo problem. However, the basic properties and even the depth of the dynamo process, which operates also in other stars, are unknown. Most of the dynamo models suggest that the toroidal magnetic field that emerges on the surface and forms sunspots is generated near the bottom of the convection zone, in the tachocline. However, there is a number of theoretical and observational problems with justifying the deep-seated dynamo models. This leads to the idea that the subsurface angular velocity shear may play an important role in the solar dynamo. Using helioseismology measurements of the internal rotation and meridional circulation, we investigate a mean-field MHD model of dynamo distributed in the bulk of the convection zone but shaped in a near-surface layer. We show that if the boundary conditions at the top of the dynamo region allow the large-scale toroidal magnetic fields to penetrate into the surface, then the dynamo wave propagates along the isosurface of angular velocity in the subsurface shear layer, forming the butterfly diagram in agreement with the Parker-Yoshimura rule and solar-cycle observations. Unlike the flux-transport dynamo models, this model does not depend on the transport of magnetic field by meridional circulation at the bottom of the convection zone, and works well when the meridional circulation forms two cells in radius, as recently indicated by deep-focus time-distance helioseismology analysis of the SDO/HMI and SOHO/MDI data. We compare the new dynamo model with various characteristics if the solar magnetic cycles, including the cycle asymmetry (Waldmeier's relations) and magnetic `butterfly' diagrams.

preprint2014arXiv

Sunquakes and starquakes

In addition to well-known mechanisms of excitation of solar and stellar oscillations by turbulent convection and instabilities, the oscillations can be excited by an impulsive localized force caused by the energy release in solar and stellar flares. Such oscillations have been observed on the Sun (`sunquakes'), and created a lot of interesting discussions about physical mechanisms of the impulsive excitation and their relationship to the flare physics. The observation and theory have shown that most of a sunquake's energy is released in high-degree, high-frequency p modes. In addition, there have been reports on helioseismic observations of low-degree modes excited by strong solar flares. Much more powerful flares observed on other stars can cause `starquakes' of substantially higher amplitude. Observations of such oscillations can provide new asteroseismic information and also constraints on mechanisms of stellar flares. The basic properties of sunquakes and initial attempts to detect flare-excited oscillations in Kepler short-cadence data are discussed.

preprint2014arXiv

Using MHD simulations to model H-alpha and UV spectral lines for interpretation of IRIS and NST data

We present results of non-LTE modeling of H-alpha 6563 A and Mg II k&h 2796 A and 2803 A lines. This modeling is important for interpretation of coordinated observations from the recently launched NASA's IRIS mission and from the New Solar Telescope at Big Bear Solar Observatory. Among available codes for the non-LTE modeling, the RH code is chosen as the most appropriate for modeling of the line profiles. The most suitable Hydrogen and Magnesium atomic models are selected by performing several tests of the code. The influence of the ionization degree on the line profiles is also studied. Radiative-MHD simulations of the solar atmosphere, obtained with the Bifrost code, are used as input data for calculation of synthetic spectra of the H-alpha and Mg II h&k lines for particular locations evolving with time. The spectral line variations reveal the presence of MHD waves in the simulation results. We construct oscillation power spectra of the line intensity for different wavelength, and compare these with the corresponding height-dependent power spectra of atmospheric parameters from the simulations. We find correlations between the power spectra of intensities of the line profiles at certain wavelengths and the power spectra of the atmospheric parameters at the tau-unity heights for these wavelengths. These results provide a new diagnostic method of chromospheric oscillations; however, larger amounts of data are needed to confirm these correlations.

preprint2012arXiv

Helioseismic Detection of Emerging Magnetic Flux

Investigating the properties of magnetic flux emergence is one of the most important problems of solar physics. In this study we present a newly developed deep-focus time-distance measurement scheme which is able to detect strong emerging flux events in the deep solar interior, before the flux becomes visible on the surface. We discuss in detail the differences between our method and previous methods, and demonstrate step-by-step how the signal-to-noise (S/N) ratio is increased. The method is based on detection of perturbations in acoustic phase travel times determined from cross-covariances of solar oscillations observed on the surface. We detect strong acoustic travel-time reductions of an order of 12 - 16 seconds at a depth of 42 - 75 Mm. These acoustic anomalies are detected 1 - 2 days before high peaks in the photospheric magnetic flux rate implying that the average emerging speed is 0.3 - 0.6 km/s. The results of this work contribute to our understanding of solar magnetism and benefit space weather forecasting.

preprint2012arXiv

Solar wave-field simulation for testing prospects of helioseismic measurements of deep meridional flows

The meridional flow in the Sun is an axisymmetric flow that is generally poleward directed at the surface, and is presumed to be of fundamental importance in the generation and transport of magnetic fields. Its true shape and strength, however, is debated. We present a numerical simulation of helioseismic wave propagation in the whole solar interior in the presence of a prescribed, stationary, single-cell, deep meridional circulation serving as a test-bed for helioseismic measurement techniques. A deep-focusing time-distance helioseismology technique is applied to the artificial data showing that it can in fact be used to measure the effects of the meridional flow very deep in the solar convection zone. It is shown that the ray-approximation which is commonly used for interpretation of helioseismology measurements remains a reasonable approximation even for the very long distances between 12 and 42 degrees corresponding to depths between 52 and 195 Mm considered here. From the measurement noise we extrapolate that on the order of a full solar cycle may be needed to probe the flow all the way to the base of the convection zone.

preprint2012arXiv

Verification of the Travel Time Measurement Technique and the Helioseismic Inversion Procedure for Sound Speed Using Artificial Data

We performed 3D numerical simulations of the solar surface wave field for the quiet Sun and for three models with different localized sound-speed variations in the interior with: (i) deep, (ii) shallow, and (iii) two-layer structures. We used simulated data generated by two different codes which use the same standard solar model as a background model, but utilize two different integration techniques and use different models of stochastic wave excitation. Acoustic travel times were measured from all data sets using the time-distance helioseismology technique and compared with the ray theory predictions, frequently used for helioseismic travel-time inversions. It is found that the measured travel-time shifts agree well with the ray theory in both cases with and without phase-speed filtering for the shallow and deep perturbations. This testing verifies the whole measuring-filtering-inversion procedure for sound-speed anomalies inside the Sun. It is shown, that the phase-speed filtering, frequently used to improve the signal-to-noise ratio does not introduce significant systematic errors. Results of the sound-speed inversion procedure show good agreement with the background sound-speed profiles in all cases. Due to its smoothing nature, the inversion procedure overestimates sound speed variations in areas with sharp gradients of the sound-speed profile.

preprint2011arXiv

Advances in Global and Local Helioseismology: an Introductory Review

Helioseismology studies the structure and dynamics of the Sun's interior by observing oscillations on the surface. These studies provide information about the physical processes that control the evolution and magnetic activity of the Sun. In recent years, helioseismology has made substantial progress towards the understanding of the physics of solar oscillations and the physical processes inside the Sun, thanks to observational, theoretical and modeling efforts. In addition to the global seismology of the Sun based on measurements of global oscillation modes, a new field of local helioseismology, which studies oscillation travel times and local frequency shifts, has been developed. It is capable of providing 3D images of the subsurface structures and flows. The basic principles, recent advances and perspectives of global and local helioseismology are reviewed in this article.

preprint2011arXiv

Helioseismic response to X2.2 solar flare of February 15, 2011

The X2.2-class solar flare of February 15, 2011, produced a powerful sunquake event, representing a helioseismic response to the flare impact in the solar photosphere, which was observed with the HMI instrument on the Solar Dynamics Observatory (SDO). The impulsively excited acoustic waves formed a compact wavepacket traveling through the solar interior and appearing on the surface as expanding wave ripples. The initial flare impacts were observed in the form of compact and rapid variations of the Doppler velocity, line-of-sight magnetic field and continuum intensity. These variations formed a typical two-ribbon flare structure, and are believed to be associated with thermal and hydrodynamic effects of high-energy particles heating the lower atmosphere. The analysis of the SDO/HMI and X-ray data from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) shows that the helioseismic waves were initiated by the photospheric impact in the early impulsive phase, observed prior to the hard X-ray (50-100 keV) impulse, and were probably associated with atmospheric heating by relatively low-energy electrons (~6-50 keV) and heat flux transport. The impact caused a short motion in the sunspot penumbra prior to the appearance of the helioseismic wave. It is found that the helioseismic wave front traveling through a sunspot had a lower amplitude and was significantly delayed relative to the front traveling outside the spot. These observations open new perspectives for studying the flare photospheric impacts and for using the flare-excited waves for sunspot seismology.

preprint2010arXiv

Detection of Supergranulation Alignment in Polar Regions of the Sun by Helioseismology

We report on a new phenomenon of `alignment' of supergranulation cells in the polar regions of the Sun. Recent high-resolution datasets obtained by the Solar Optical Telescope onboard the Hinode satellite enabled us to investigate supergranular structures in high-latitude regions of the Sun. We have carried out a local helioseismology time-distance analysis of the data, and detected acoustic travel-time variations due to the supergranular flows. The supergranulation cells in both the north and south polar regions show systematic alignment patterns in the north-south direction. The south-pole datasets obtained in a month-long Hinode campaign indicate that the supergranulation alignment property may be quite common in the polar regions. We also discuss the latitudinal dependence of the supergranulation cell sizes; the data show that the east-west cell size decreases towards higher latitudes.

preprint2010arXiv

Signatures of Emerging Subsurface Structures in Acoustic Power Maps

We show that under certain conditions, subsurface structures in the solar interior can alter the average acoustic power observed at the photosphere above them. By using numerical simulations of wave propagation, we show that this effect is large enough for it to be potentially used for detecting emerging active regions before they appear on the surface. In our simulations, simplified subsurface structures are modeled as regions with enhanced or reduced acoustic wave speed. We investigate the dependence of the acoustic power above a subsurface region on the sign, depth, and strength of the wave speed perturbation. Observations from the Solar and Heliospheric Observatory/Michelson Doppler Imager (SOHO/MDI) prior and during the emergence of NOAA active region 10488 are used to test the use of acoustic power as a potential precursor of magnetic flux emergence.

preprint2010arXiv

Towards Waveform Heliotomography: Observing Interactions of Helioseismic Waves with a Sunspot

We investigate how helioseismic waves that originate from effective point sources interact with a sunspot. These waves are reconstructed from observed stochastic wavefields on the Sun by cross-correlating photospheric Doppler-velocity signals. We select the wave sources at different locations relative to the sunspot, and investigate the p- and f-mode waves separately. The results reveal a complicated picture of waveform perturbations caused by the wave interaction with the sunspot. In particular, it is found that for waves originating from outside of the sunspot, p-mode waves travel across the sunspot with a small amplitude reduction and slightly higher speed, and wave amplitude and phase get mostly restored to the quiet-Sun values after passing the sunspot. The f-mode wave experiences some amplitude reduction passing through the sunspot, and the reduced amplitude is not recovered after that. The wave-propagation speed does not change before encountering the sunspot and inside the sunspot, but the wavefront becomes faster than the reference wave after passing through the sunspot. For waves originating from inside the sunspot umbra, both f- and p-mode waves show significant amplitude reductions and faster speed during all courses of propagation. A comparison of positive and negative time lags of cross-correlation functions shows an apparent asymmetry in the waveform changes for both the f- and p-mode waves. We suggest that the waveform variations of the helioseismic waves interacting with a sunspot found in this article can be used for developing a method of waveform heliotomography, similar to the waveform tomography of the Earth.

preprint2009arXiv

High Resolution Helioseismic Imaging of Subsurface Structures and Flows of A Solar Active Region Observed by Hinode

We analyze a solar active region observed by the Hinode CaII H line using the time-distance helioseismology technique, and infer wave-speed perturbation structures and flow fields beneath the active region with a high spatial resolution. The general subsurface wave-speed structure is similar to the previous results obtained from SOHO/MDI observations. The general subsurface flow structure is also similar, and the downward flows beneath the sunspot and the mass circulations around the sunspot are clearly resolved. Below the sunspot, some organized divergent flow cells are observed, and these structures may indicate the existence of mesoscale convective motions. Near the light bridge inside the sunspot, hotter plasma is found beneath, and flows divergent from this area are observed. The Hinode data also allow us to investigate potential uncertainties caused by the use of phase-speed filter for short travel distances. Comparing the measurements with and without the phase-speed filtering, we find out that inside the sunspot, mean acoustic travel times are in basic agreement, but the values are underestimated by a factor of 20-40% inside the sunspot umbra for measurements with the filtering. The initial acoustic tomography results from Hinode show a great potential of using high-resolution observations for probing the internal structure and dynamics of sunspots.

preprint2008arXiv

Large-scale horizontal flows in the solar photosphere. IV. On the vertical structure of large-scale horizontal flows

In the recent papers, we introduced a method utilised to measure the flow field. The method is based on the tracking of supergranular structures. We did not precisely know, whether its results represent the flow field in the photosphere or in some sub-photospheric layers. In this paper, in combination with helioseismic data, we are able to estimate the depths in the solar convection envelope, where the detected large-scale flow field is well represented by the surface measurements. We got a clear answer to question what kind of structures we track in full-disc Dopplergrams. It seems that in the quiet Sun regions the supergranular structures are tracked, while in the regions with the magnetic field the structures of the magnetic field are dominant. This observation seems obvious, because the nature of Doppler structures is different in the magnetic regions and in the quiet Sun. We show that the large-scale flow detected by our method represents the motion of plasma in layers down to ~10 Mm. The supergranules may therefore be treated as the objects carried by the underlying large-scale velocity field.

preprint2007arXiv

Initial Helioseismic Observations by Hinode/SOT

Results from initial helioseismic observations by Solar Optical Telescope onboard Hinode are reported. It has been demonstrated that intensity oscillation data from Broadband Filter Imager can be used for various helioseismic analyses. The k-omega power spectra, as well as corresponding time-distance cross-correlation function that promises high-resolution time-distance analysis below 6-Mm travelling distance, were obtained for G-band and CaII-H data. Subsurface supergranular patterns have been observed from our first time-distance analysis. The results show that the solar oscillation spectrum is extended to much higher frequencies and wavenumbers, and the time-distance diagram is extended to much shorter travel distances and times than they were observed before, thus revealing great potential for high-resolution helioseismic observations from Hinode.

preprint2007arXiv

Observations of Sunspot Oscillations in G band and Ca II H line with Solar Optical Telescope on Hinode

Exploiting high-resolution observations made by the Solar Optical Telescope onboard Hinode, we investigate the spatial distribution of power spectral density of oscillatory signal in and around NOAA active region 10935. The G-band data show that in the umbra the oscillatory power is suppressed in all frequency ranges. On the other hand, in Ca II H intensity maps oscillations in the umbra, so-called umbral flashes, are clearly seen with the power peaking around 5.5 mHz. The Ca II H power distribution shows the enhanced elements with the spatial scale of the umbral flashes over most of the umbra but there is a region with suppressed power at the center of the umbra. The origin and property of this node-like feature remain unexplained.

preprint2006arXiv

Validation of Time-Distance Helioseismology by Use of Realistic Simulations of Solar Convection

Recent progress in realistic simulations of solar convection have given us an unprecedented opportunity to evaluate the robustness of solar interior structures and dynamics obtained by methods of local helioseismology. We present results of testing the time-distance method using realistic simulations. By computing acoustic wave propagation time and distance relations for different depths of the simulated data, we confirm that acoustic waves propagate into the interior and then turn back to the photosphere. This demonstrates that in the numerical simulations properties of acoustic waves (p-modes) are similar to the solar conditions, and that these properties can be analyzed by the time-distance technique. For the surface gravity waves (f-mode), we calculate perturbations of their travel times, caused by localized downdrafts, and demonstrate that the spatial pattern of these perturbations (representing so-called sensitivity kernels) is similar to the patterns obtained from the real Sun, displaying characteristic hyperbolic structures. We then test the time-distance measurements and inversions by calculating acoustic travel times from a sequence of vertical velocities at the photosphere of the simulated data, and inferring a mean 3D flow fields by performing inversion based on the ray approximation. The inverted horizontal flow fields agree very well with the simulated data in subsurface areas up to 3 Mm deep, but differ in deeper areas. Due to the cross-talk effects between the horizontal divergence and downward flows, the inverted vertical velocities are significantly different from the mean convection velocities of the simulation dataset.