Researcher profile

S. K. Solanki

S. K. Solanki contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

24 published item(s)

preprint2025arXiv

Formation of chromospheric Fe I lines in the near ultraviolet in 1D atmospheres

In the near ultraviolet (NUV) part of the solar spectrum, there are several Fe I lines with very broad profiles, typical of chromospheric lines. These lines are largely unexplored due to the lack of high-resolution data in this region. This changed with the successful Sunrise III flight in 2024, when spectro-polarimetric data were recorded with high spatial, spectral, and temporal resolution covering a large variety of solar targets. The aim of this work is to investigate the formation of the lines and lay the groundwork for further studies. We compute the spectrum of the lines at 358.12 nm, 371.99 nm, 406.36 nm, and 407.17 nm emerging from the standard 1D FAL-atmospheres using the non-local thermodynamic equilibrium (NLTE) radiative transfer code RH. We find that the lines are affected by overionization in the wings, but have line cores strongly affected by scattering. The line cores form well into the chromosphere in the tested atmosphere models except the colder FALX model where the line core forms in the temperature minimum (which lies at traditional chromospheric heights). In the presence of a vertical magnetic field, the Stokes $V$ signal is multi-lobed. The profile can be decomposed into two broad photospheric lobes and two sharper lobes forming in the flanks of the chromospheric line core. We have investigated the properties of four lines of Fe I that form in the lower chromosphere. The results provide a basic understanding of the formation of the lines, which will be useful for later analysis of formation in 3D magnetohydrodynamic simulations and an eventual investigation into their diagnostic potential.

preprint2023arXiv

Focused Space Weather Strategy for Securing Earth, and Human Exploration of the Moon and Mars

This white paper recognizes gaps in observations that will, when addressed, much improve solar radiation hazard and geomagnetic storm forecasting. Radiation forecasting depends on observations of the entire "Solar Radiation Hemisphere" that we will define. Mars exploration needs strategic placement of radiation-relevant observations. We also suggest an orbital solution that will improve geomagnetic storm forecasting through improved in situ and solar/heliospheric remote sensing.

preprint2022arXiv

Chromospheric Extension of the MURaM Code

Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines.

preprint2022arXiv

Faculae cancel out on the surfaces of active Suns

Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres.

preprint2022arXiv

Magnetized supersonic downflows in the chromosphere. A statistical study using the He I 10830 Å lines

The chromosphere above active regions (ARs) on the Sun hosts magnetized supersonic downflows. Studies of these supersonic downflows help to decipher the magnetic fine structure and dynamics of the chromosphere. We perform a statistical analysis of the magnetized supersonic downflows in a number of ARs and survey their characteristics. We analyze spectro-polarimetric scans of parts of 13 ARs obtained in the infrared He I 10830 Å triplet formed in the upper chromosphere recorded with the GREGOR Infrared Spectrograph (GRIS) mounted at the GREGOR solar telescope. We retrieve the line-of-sight velocities and the magnetic field vector using the HeLIx+ inversion code that assumes Milne-Eddington atmospheres. We find magnetized supersonic downflows in all the ARs, with larger area coverage by such flows in ARs observed during their emerging phase. The fact that supersonic downflows were detected in all scans, though they cover only a small fraction, 0.2--6.4%, of the observed field-of-view, suggests that they are a common phenomenon in the upper chromospheres of ARs. The supersonic downflows are found to be associated with many AR features such as pores, sunspot umbrae, sunspot penumbrae, light bridges, plages, He I loops as part of arch filament systems characteristic of emerging fields, and filaments. Although several mechanisms are identified to be causing the supersonic downflows, by far the most common one appears to be the draining of plasma along the legs of rising magnetic loops. The loops mainly drain into forming pores. The line-of-sight velocities of the supersonic downflows reach up to 49 kms-1 and the velocity distribution shows multiple populations. Almost 92% of these supersonic downflows coexist with a subsonic flow component. The weaker, more horizontal fields associated with the supersonic component suggests that it is formed above the subsonic component.

preprint2022arXiv

Predictions of Astrometric Jitter for Sun-like Stars. III. Fast Rotators

A breakthrough in exoplanet detections is foreseen with the unprecedented astrometric measurement capabilities offered by instrumentation aboard Gaia space observatory. Besides, astrometric discoveries of exoplanets are expected from the planned space mission, Small-JASMINE. In this setting, the present series of papers focuses on estimating the effect of magnetic activity of G2V-type host stars on the astrometric signal. This effect interferes with the astrometric detections of Earth-mass planets. While the first two papers considered stars rotating at the solar rotation rate, this paper focuses on stars having solar effective temperature and metallicity but rotating faster than the Sun, and consequently more active. By simulating the distribution of active regions on such stars using the Flux Emergence And Transport model, we show that the contribution of magnetic activity to the astrometric measurements becomes increasingly significant with increasing rotation rates. We further show that the jitter for the most variable periodic Kepler stars is high enough to be detected by Gaia. Furthermore, due to a decrease in the facula-to-spot area ratio for more active stars, the magnetic jitter is found to be spot-dominated for rapid rotators. Our simulations of the astrometric jitter has the potential to aid the interpretation of data from Gaia and upcoming space astrometry missions.

preprint2022arXiv

The magnetic drivers of campfires seen by the Polarimetric and Helioseismic Imager (PHI) on Solar Orbiter

The Extreme Ultraviolet Imager (EUI) on board the Solar Orbiter (SO) spacecraft observed small extreme ultraviolet (EUV) bursts, termed campfires, that have been proposed to be brightenings near the apexes of low-lying loops in the quiet-Sun atmosphere. The underlying magnetic processes driving these campfires are not understood. During the cruise phase of SO and at a distance of 0.523\,AU from the Sun, the Polarimetric and Helioseismic Imager on Solar Orbiter (SO/PHI) observed a quiet-Sun region jointly with SO/EUI, offering the possibility to investigate the surface magnetic field dynamics underlying campfires at a spatial resolution of about 380~km. In 71\% of the 38 isolated events, campfires are confined between bipolar magnetic features, which seem to exhibit signatures of magnetic flux cancellation. The flux cancellation occurs either between the two main footpoints, or between one of the footpoints of the loop housing the campfire and a nearby opposite polarity patch. In one particularly clear-cut case, we detected the emergence of a small-scale magnetic loop in the internetwork followed soon afterwards by a campfire brightening adjacent to the location of the linear polarisation signal in the photosphere, that is to say near where the apex of the emerging loop lays. The rest of the events were observed over small scattered magnetic features, which could not be identified as magnetic footpoints of the campfire hosting loops. The majority of campfires could be driven by magnetic reconnection triggered at the footpoints, similar to the physical processes occurring in the burst-like EUV events discussed in the literature. About a quarter of all analysed campfires, however, are not associated to such magnetic activity in the photosphere, which implies that other heating mechanisms are energising these small-scale EUV brightenings.

preprint2022arXiv

The on-ground data reduction and calibration pipeline for SO/PHI-HRT

The ESA/NASA Solar Orbiter space mission has been successfully launched in February 2020. Onboard is the Polarimetric and Helioseismic Imager (SO/PHI), which has two telescopes, a High Resolution Telescope (HRT) and the Full Disc Telescope (FDT). The instrument is designed to infer the photospheric magnetic field and line-of-sight velocity through differential imaging of the polarised light emitted by the Sun. It calculates the full Stokes vector at 6 wavelength positions at the Fe I 617.3 nm absorption line. Due to telemetry constraints, the instrument nominally processes these Stokes profiles onboard, however when telemetry is available, the raw images are downlinked and reduced on ground. Here the architecture of the on-ground pipeline for HRT is presented, which also offers additional corrections not currently available on board the instrument. The pipeline can reduce raw images to the full Stokes vector with a polarimetric sensitivity of $10^{-3}\cdot I_{c}$ or better.

preprint2021arXiv

A solar coronal loop in a box: Energy generation and heating

Coronal loops are the basic building block of the upper solar atmosphere. Comprehending how these are energized, structured, and evolve is key to understanding stellar coronae. Here we investigate how the energy to heat the loop is generated by photospheric magneto-convection, transported into the upper atmosphere, and how the internal structure of a coronal loop forms. In a 3D magnetohydrodynamics (MHD) model, we study an isolated coronal loop rooted with both footpoints in a shallow layer within the convection zone using the MURaM code. To resolve its internal structure, we limited the computational domain to a rectangular box containing a single coronal loop as a straightened magnetic flux tube. Field-aligned heat conduction, gray radiative transfer in the photosphere and chromosphere, and optically thin radiative losses in the corona were taken into account. The footpoints were allowed to interact self-consistently with the granulation surrounding them. The loop is heated by a Poynting flux that is self-consistently generated through small-scale motions within individual magnetic concentrations in the photosphere. Turbulence develops in the upper layers of the atmosphere as a response to the footpoint motions. We see little sign of heating by large-scale braiding of magnetic flux tubes from different photospheric concentrations at a given footpoint. The synthesized emission, as it would be observed by the Atmospheric Imaging Assembly or the X-ray Telescope, reveals transient bright strands that form in response to the heating events. Overall, our model roughly reproduces the properties and evolution of the plasma as observed within coronal loops. With this model we can build a coherent picture of how the energy flux to heat the upper atmosphere is generated near the solar surface and how this process drives and governs the heating and dynamics of a coronal loop.

preprint2021arXiv

Observing and modelling the young solar analogue EK Draconis: starspot distribution, elemental abundances, and evolutionary status

Observations and modelling of stars with near-solar masses in their early phases of evolution is critical for a better understanding of how dynamos of solar-type stars evolve. We examine the chemical composition and the spot distribution of the pre-main-sequence solar analogue EK Dra. Using spectra from the HERMES Spectrograph (La Palma), we obtain the abundances of 23 elements with respect to the solar ones, which lead to a $[{\rm Fe/H}]=0.03$, with significant overabundance of Li and Ba. The s-process elements Sr, Y, and Ce are marginally overabundant, while Co, Ni, Cu, Zn are marginally deficient compared to solar abundances. The overabundance of Ba is most likely due to the assumption of depth-independent microturbulent velocity. Li abundance is consistent with the age and the other abundances may indicate distinct initial conditions of the pre-stellar nebula. We estimate a mass of 1.04 $M_\odot$ and an age of $27^{+11}_{-8}$\,Myr using various spectroscopic and photometric indicators. We study the surface distribution of dark spots, using 17 spectra collected during 15 nights using the CAFE Spectrograph (Calar Alto). We also conduct flux emergence and transport (FEAT) simulations for EK Dra's parameters and produce 15-day-averaged synoptic maps of the likely starspot distributions. Using Doppler imaging, we reconstruct the surface brightness distributions for the observed spectra and FEAT simulations, which show overall agreement for polar and mid-latitude spots, while in the simulations there is a lack of low-latitude spots compared to the observed image. We find indications that cross-equatorial extensions of mid-latitude spots can be artefacts of the less visible southern-hemisphere activity.

preprint2021arXiv

The Dimmest State of the Sun

How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm-2 to a rise of 6.3 Wm-2. Consequently, the exact contribution by solar forcing to the rise in global temperatures over the past centuries remains inconclusive. Adopting a novel approach based on state-of-the-art solar imagery and numerical simulations, we establish the TSI level of the Sun when it is in its least-active state to be 2.0 +/- 0.7 Wm-2 below the 2019 level. This means TSI could not have risen since the Maunder minimum by more than this amount, thus restricting the possible role of solar forcing in global warming.

preprint2020arXiv

Connecting measurements of solar and stellar brightness variations

Comparing solar and stellar brightness variations is hampered by the difference in spectral passbands used in observations as well as by the possible difference in the inclination of their rotation axes from the line of sight. We calculate the rotational variability of the Sun as it would be measured in passbands used for stellar observations. In particular, we consider the filter systems used by the CoRoT, $Kepler$, TESS, and $Gaia$ space missions. We also quantify the effect of the inclination of the rotation axis on the solar rotational variability. We employ the Spectral And Total Irradiance REconstructions (SATIRE) model to calculate solar brightness variations in different filter systems as observed from the ecliptic plane. We then combine the simulations of the surface distribution of the magnetic features at different inclinations using a surface flux transport model (SFTM) with the SATIRE calculations to compute the dependence of the variability on the inclination. For an ecliptic-bound observer, the amplitude of the solar rotational variability, as observed in the total solar irradiance (TSI) is 0.68 mmag (averaged over solar cycles 21-24). We obtained corresponding amplitudes in the $Kepler$ (0.74 mmag), CoRoT (0.73 mmag), TESS (0.62 mmag), $Gaia~ $ (0.74 mmag), $Gaia~ G_{RP}$ (0.62 mmag), and ), $Gaia~ G_{BP}$ (0.86 mmag) passbands. Decreasing the inclination of the rotation axis decreases the rotational variability. For a sample of randomly inclined stars, the variability is on average 15% lower in all filter systems considered in this work. This almost compensates for the difference in the amplitudes of the variability in TSI and $Kepler$ passbands, making the amplitudes derived from the TSI records an ideal representation of the solar rotational variability for comparison to $Kepler$ stars with unknown inclinations.

preprint2020arXiv

Connecting the Wilson depression to the magnetic field of sunspots

Context: In sunspots, the geometric height of continuum optical depth unity is depressed compared to the quiet Sun. This so-called Wilson depression is caused by the Lorentz force of the strong magnetic field inside the spots. However, it is not understood in detail yet, how the Wilson depression is related to the strength and geometry of the magnetic field or to other properties of the sunspot. Aims: We aim to study the dependence of the Wilson depression on the properties of the magnetic field of the sunspots and how exactly the magnetic field contributes to balancing the Wilson depression with respect to the gas pressure of the surroundings of the spots. Methods: Our study is based on 24 spectropolarimetric scans of 12 individual sunspots performed with Hinode. We derived the Wilson depression for each spot using both, a recently developed method that is based on minimizing the divergence of the magnetic field, and an approach developed earlier that enforces an equilibrium between the gas pressure and the magnetic pressure inside the spot and the gas pressure in the quiet Sun, thus neglecting the influence of the curvature force. We then performed a statistical analysis by comparing the Wilson depression resulting from the two techniques with each other and by relating them to various parameters of the sunspots, such as their size or the strength of the magnetic field. Results: We find that the Wilson depression becomes larger for spots with a stronger magnetic field, but not as much as one would expect from the increased magnetic pressure. This suggests that the curvature integral provides an important contribution to the Wilson depression, particularly for spots with a weak magnetic field. Our results indicate that the geometry of the magnetic field in the penumbra is different between spots with different strengths of the average umbral magnetic field.

preprint2020arXiv

Detection of the strongest magnetic field in a sunspot light bridge

Traditionally, the strongest magnetic fields on the Sun have been measured in sunspot umbrae. More recently, however, much stronger fields have been measured at the ends of penumbral filaments carrying the Evershed and counter-Evershed flows. Superstrong fields have also been reported within a light bridge separating two umbrae of opposite polarities. We aim to accurately determine the strengths of the strongest fields in a light bridge using an advanced inversion technique and to investigate their detailed structure. We analyze observations from the spectropolarimeter on board the Hinode spacecraft of the active region AR 11967. The thermodynamic and magnetic configurations are obtained by inverting the Stokes profiles using an inversion scheme that allows multiple height nodes. Both the traditional 1D inversion technique and the so-called 2D coupled inversions, which take into account the point spread function of the Hinode telescope, are used. We find a compact structure with an area of 32.7 arcsec$^2$ within a bipolar light bridge with field strengths exceeding 5 kG, confirming the strong fields in this light bridge reported in the literature. Two regions associated with downflows of $\sim$5 km s$^{-1}$ harbor field strengths larger than 6.5 kG, covering a total area of 2.97 arcsec$^2$. The maximum field strength found is 8.2 kG, which is the largest ever observed field in a bipolar light bridge up to now.

preprint2020arXiv

Effect of metallicity on the detectability of rotational periods in solar-like stars

Accurate determination of stellar rotation periods is important for estimating stellar ages as well as for understanding stellar activity and evolution. While for about thirty thousand stars in the Kepler field rotation periods can be determined, there are over hundred thousand stars, especially with low photometric variability and irregular pattern of variations, for which rotational periods are unknown. Here, we investigate the effect of metallicity on the detectability of rotation periods. This is done by synthesising light curves of hypothetical stars, which are identical to our Sun, with the exception of the metallicity. These light curves are then used as an input to the period determination algorithms. We find that the success rate for recovering the rotation signal has a minimum close to the solar metallicity value. This can be explained by the compensation effect of facular and spot contributions. In addition, selecting solar-like stars with near-solar effective temperature, near solar photometric variability, and with metallicity between M/H = -0.35 and M/H = 0.35 from the Kepler sample, we analyse the fraction of stars for which rotational periods have been detected as a function of metallicity. In agreement with our theoretical estimate we found a local minimum for the detection fraction close to the solar metallicity. We further report rotation periods of 87~solar-like Kepler stars for the first time.

preprint2020arXiv

No universal connection between the vertical magnetic field and the umbra-penumbra boundary in sunspots

Context. It has been reported that the boundary between the umbra and the penumbra of sunspots occurs at a canonical value of the strength of the vertical magnetic field, independently of the size of the spot. This critical field strength is interpreted as to be the threshold for the onset of magnetoconvection. Aims. Here we investigate the reasons why this criterion, also called the Jurčák criterion in the literature, does not always identify the boundary between umbra and penumbra. Methods. We perform a statistical analysis of 23 sunspots observed with Hinode/SOT. We compare the properties of the continuum intensity and the vertical magnetic field between filaments and spines and how they vary between spots of different sizes. Results. We find that the inner boundary of the penumbra is not related to a universal value of the vertical magnetic field. The properties of spines and filaments vary between spots of different sizes. Both components are darker in larger spots and the spines exhibit stronger vertical magnetic field. These variations of the properties of filaments and spines with spot size are also the reason for the reported invariance of the averaged vertical magnetic field at 50% of the mean continuum intensity. Conclusions. The formation of filaments and the onset of magnetoconvection are not related to a canonical value of the strength of the vertical magnetic field. Such a seemingly unique magnetic field strength is rather an effect of the filling factor of spines and penumbral filaments.

preprint2020arXiv

On the Magnetic Nature of an Exploding Granule as Revealed by Sunrise/IMaX

We study the photospheric evolution of an exploding granule observed in the quiet Sun at high spatial ($0.3^{\prime\prime}$) and temporal (31.5 s) resolution by the imaging magnetograph Sunrise/IMaX in June 2009. These observations show that the exploding granule is cospatial to a magnetic flux emergence event occurring at mesogranular scale (up to 12 Mm$^{2}$ area). Using a modified version of the SIR code for inverting the IMaX spectropolarimetric measurements, we obtain information about the magnetic configuration of this photospheric feature. In particular, we find evidence of highly inclined emerging fields in the structure, carrying a magnetic flux content up to $4 \times 10^{18}$ Mx. The balance between gas and magnetic pressure in the region of flux emergence, compared with a very quiet region of the Sun, indicates that the additional pressure carried by the emerging flux increases by about 5% the total pressure and appears to allow the granulation to be modified, as predicted by numerical simulations. The overall characteristics suggest that a multi-polar structure emerges into the photosphere, resembling an almost horizontal flux sheet. This seems to be associated with exploding granules. Finally, we discuss the origin of such flux emergence events.

preprint2020arXiv

Power spectra of solar brightness variations at different inclinations

Magnetic features on the surfaces of cool stars cause variations of their brightness. Such variations have been extensively studied for the Sun. Recent planet-hunting space telescopes allowed measuring brightness variations in hundred thousands of other stars. The new data posed the question of how typical is the Sun as a variable star. Putting solar variability into the stellar context suffers, however, from the bias of solar observations being made from its near-equatorial plane, whereas stars are observed at all possible inclinations. We model solar brightness variations at timescales from days to years as they would be observed at different inclinations. In particular, we consider the effect of the inclination on the power spectrum of solar brightness variations. The variations are calculated in several passbands routinely used for stellar measurements. We employ the Surface Flux Transport Model (SFTM) to simulate the time-dependent spatial distribution of magnetic features on both near- and far-sides of the Sun. This distribution is then used to calculate solar brightness variations following the SATIRE (Spectral And Total Irradiance REconstruction) approach. We have quantified the effect of the inclination on solar brightness variability at timescales down to a day. Thus, our results allow making solar brightness records directly comparable to those obtained by the planet-hunting space telescopes. Furthermore, we decompose solar brightness variations into the components originating from the solar rotation and from the evolution of magnetic features.

preprint2020arXiv

Power spectrum of turbulent convection in the solar photosphere

The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law.

preprint2020arXiv

Reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars"

This is our reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars" by T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gomez. We hope that both the comment and our reply lead to fruitful discussions which of the two presented scenarios is more likely.

preprint2020arXiv

Rotation periods from the inflection point in the power spectrum of stellar brightness variations: II. The Sun

Young and active stars generally have regular, almost sinusoidal, patterns of variability attributed to their rotation, while the majority of older and less active stars, including the Sun, have more complex and non-regular light-curves which do not have clear rotational-modulation signals. Consequently, the rotation periods have been successfully determined only for a small fraction of the Sun-like stars observed by transit-based planet-hunting missions, such as CoRoT, Kepler, and TESS. This suggests that only a small fraction of such systems have been properly identified as solar-like analogs. We apply a new method for determining rotation periods of low-activity stars, like the Sun. The method is based on calculating the gradient of the power spectrum (GPS) of stellar brightness variations and identifying a tell-tale inflection point in the spectrum. The rotation frequency is then proportional to the frequency of that inflection point. In this paper test this GPS method against available photometric records of the Sun. We apply GPS, autocorrelation functions, Lomb-Scargle periodograms, and wavelet analyses to the total solar irradiance (TSI) time series obtained from the Total Irradiance Monitor (TIM) on the Solar Radiation and Climate Experiment (SORCE) and the Variability of solar IRradiance and Gravity Oscillations (VIRGO) experiment on the SOlar and Heliospheric Observatory (SoHO) missions. We analyse the performance of all methods at various levels of solar activity. We show that the GPS method returns accurate values of solar rotation independently of the level of solar activity. In particular, it performs well during periods of high solar activity, when TSI variability displays an irregular pattern and other methods fail. Our results suggest that the GPS method can successfully determine the rotational periods of stars with both regular and non-regular light-curves.

preprint2020arXiv

Simulations Show that Vortex Flows could Heat the Chromosphere in Solar Plage

The relationship between vortex flows at different spatial scales and their contribution to the energy balance in the chromosphere is not yet fully understood. We perform three-dimensional (3D) radiation-magnetohydrodynamic (MHD) simulations of a unipolar solar plage region at a spatial resolution of 10 km using the MURaM code. We use the swirling-strength criterion that mainly detects the smallest vortices present in the simulation data. We additionally degrade our simulation data to smooth-out the smaller vortices, so that also the vortices at larger spatial scales can be detected. Vortex flows at various spatial scales are found in our simulation data for different effective spatial resolutions. We conclude that the observed large vortices are likely clusters of much smaller ones that are not yet resolved by observations. We show that the vertical Poynting flux decreases rapidly with reduced effective spatial resolutions and is predominantly carried by the horizontal plasma motions rather than vertical flows. Since the small-scale horizontal motions or the smaller vortices carry most of the energy, the energy transported by vortices deduced from low resolution data is grossly underestimated. In full resolution simulation data, the Poynting flux contribution due to vortices is more than adequate to compensate for the radiative losses in plage, indicating their importance for chromospheric heating.

preprint2020arXiv

The influence of NLTE effects in Fe I lines on an inverted atmosphere I. 6301 A and 6302 A lines formed in 1D NLTE

Ultraviolet over-ionisation of iron atoms in the solar atmosphere leads to deviations in their level populations from the Saha-Boltzmann statistics. This causes their line profiles to form in Non-Local Thermodynamic Equilibrium (NLTE) conditions. While inverting such profiles to determine atmospheric parameters, the NLTE effects are often neglected and deviations from LTE are compensated for by tweaking other quantities. We investigate how the routinely employed LTE inversion of iron lines formed in NLTE under- or over-estimates atmospheric quantities such as temperature (T), line-of-sight velocity (v_LOS), magnetic field strength (B) and inclination (gamma) while the previous papers have focused mainly on T. We synthesize the Stokes profiles of Fe I 6301.5 A and 6302.5 A lines in both LTE and NLTE using a snapshot of a 3D MHD simulation. The profiles are then inverted in LTE. By considering the atmosphere inferred from inversion of LTE profiles to be the fiducial model, we compare atmosphere from the inversion of NLTE profiles with it. Any differences observed are attributed to NLTE effects. Neglecting the NLTE effects introduces errors in the inverted atmosphere. While the errors in T can go up to 13%, in v_LOS and B the errors can be as high as 50% or more. We find these errors to be present at all three inversion nodes. Importantly, they survive degradation from spatial averaging of the profiles. We give an overview of how the neglect of NLTE effects influences the values of T, v_LOS, B and gamma determined by inverting Fe I 6300 A line pair, as observed, e.g., by Hinode. Errors are found at the sites of granules, intergranular lanes, magnetic elements, basically in every region susceptible to NLTE effects. For an accurate determination of atmospheric quantities and their stratification, it is therefore important to take account of NLTE effects.

preprint2019arXiv

Inflection point in the power spectrum of stellar brightness variations. I. The model

Context. Considerable effort has been put into using light curves observed by space telescopes such as CoRoT, Kepler and TESS for determining stellar rotation periods. While rotation periods of active stars can be reliably determined, the light curves of many older and less active stars (e.g. stars similar to the Sun) are quite irregular, which hampers determination of their periods. Aims. We examine the factors causing the irregularities in stellar brightness variations and develop a method for determining rotation periods of low activity stars with irregular light curves. Methods. We extend the Spectral And Total Irradiance Reconstruction (SATIRE) approach for modelling solar brightness variations to Sun-like stars. We calculate the power spectra of stellar brightness variations for various combinations of parameters defining the surface configuration and evolution of stellar magnetic features. Results. The short lifetime of spots in comparison to the stellar rotation period as well as the interplay between spot and facular contributions to brightness variations of stars with near solar activity cause irregularities in their light curves. The power spectra of such stars often lack a peak associated with the rotation period. Nevertheless, the rotation period can still be determined by measuring the period where the concavity of the power spectrum plotted in the log-log scale changes sign, i.e. by identifying the position of the inflection point. Conclusions. The inflection point of the (log-log) power spectrum is found to be a new diagnostic for stellar rotation periods that is shown to work even in cases where the power spectrum shows no peak at the rotation rate.