Determination of the Coronal and Interplanetary Magnetic Field Strength and Radial Profiles from Large-Scale Photospheric Magnetic Fields

We propose a new model for the magnetic field at different distances from the Sun during different phases of the solar cycle. The model depends on the observed large-scale non-polar (\({\pm}\, 55^{\circ }\)) photospheric magnetic field and on the magnetic field measured at polar regions from \(55^{\circ }\) N to \(90^{\circ }\) N and from \(55^{\circ }\) S to \(90^{\circ }\) S, which are the visible manifestations of cyclic changes in the toroidal and poloidal components of the global magnetic field of the Sun. The modeled magnetic field is determined as the superposition of the non-polar and polar photospheric magnetic field and considers cycle variations. The agreement between the model predictions and magnetic fields derived from direct in situ measurements at different distances from the Sun, obtained with different methods and at different solar activity phases, is quite satisfactory. From a comparison of the magnetic fields as observed and calculated from the model at 1 AU, we conclude that the model magnetic field variations adequately explain the main features of the interplanetary magnetic field (IMF) radial, \(B_{\mathrm{x}}\), component cycle evolution at Earth’s orbit. The modeled magnetic field averaged over a Carrington rotation (CR) correlates with the IMF \(B_{\mathrm{x}}\) component also averaged over a CR at Earth’s orbit with a coefficient of 0.691, while for seven CR-averaged data, the correlation reaches 0.81. The radial profiles of the modeled magnetic field are compared with those of already existing models. In contrast to existing models, ours provides realistic magnetic-field radial distributions over a wide range of heliospheric distances at different cycle phases, taking into account the cycle variations of the solar toroidal and poloidal magnetic fields. The model is a good approximation of the cycle behavior of the magnetic field in the heliosphere. In addition, the decrease in the non-polar and polar photospheric magnetic fields is shown. Furthermore, the magnetic field during solar cycle maxima and minima decreased from Cycle 21 to Cycle 24. This implies that both the toroidal and poloidal components, and therefore the solar global magnetic field, decreased from Cycle 21 to Cycle 24.     

Recurring coronal holes and their rotation rates during the solar cycles 22–24

Coronal holes (CHs) play a significant role in making the Earth geo-magnetically active during the declining and minimum phases of the solar cycle. In this study, we analysed the evolutionary characteristics of the Recurring CHs from the year 1992 to 2016. The extended minimum of Solar Cycle 23 shows unusual characteristics in the number of persistent coronal holes in the mid- and low-latitude regions of the Sun. Carrington rotation maps of He 10830 Å and EUV 195 Å observations are used to identify the Coronal holes. The latitude distribution of the RCHs shows that most of them are appeared between \(\pm 20^{\circ }\) latitudes. In this period, more number of recurring coronal holes appeared in and around \(100^{\circ }\) and \(200^{\circ }\) Carrington longitudes. The large sized coronal holes lived for shorter period and they appeared close to the equator. From the area distribution over the latitude considered, it shows that more number of recurring coronal holes with area \(<10^{21}~\mbox{cm}^{2}\) appeared in the southern latitude close to the equator. The rotation rates calculated from the RCHs appeared between \(\pm 60^{\circ }\) latitude shows rigid body characteristics. The derived rotational profiles of the coronal holes show that they have anchored to a depth well below the tachocline of the interior, and compares well with the helioseismology results.     

Analysis of the exoplanet containing system Kepler-13

We have applied the close binary system analysis program WinFitter, with its physically detailed fitting function, to an intensive study of the complex multiple system Kepler-13 using photometry data from all 13 short cadence quarters downloaded from the NASA Exoplanet Archive (NEA) (http://exoplanetarchive.ipac.caltech.edu). The data-point error of our normalized, phase-sequenced and binned (380 points per bin: 0.00025 phase interval) flux values, at 14 ppm, allows the model’s specification for the mean reference flux level of the system to a precision better than 1 ppm. Our photometrically derived values for the mass and radius of KOI13.01 are \(6.8\pm0.6~\mbox{M}_{\mathrm{J}}\) and \(1.44\pm0.04~\mbox{R}_{\mathrm{J}}\). The star has a radius of \(1.67\pm0.05~\mbox{R}_{\odot}\). Our modelling sets the mean of the orbital inclination \(i\) at \(94.35\pm0.14^{\circ}\), with the star’s mean precession angle \(\phi_{p}\)—\(49.1\pm5.0^{\circ}\) and obliquity \(\theta_{o}\)\(67.9 \pm 3.0^{\circ}\), though there are known ambiguities about the sense in which such angles are measured.Our findings did not confirm secular variation in the transit modelling parameters greater than their full correlated errors, as argued by previous authors, when each quarter’s data was best-fitted with a determinable parameter set without prejudice. However, if we accept that most of the parameters remain the same for each transit, then we could confirm a small but steady diminution in the cosine of the orbital inclination over the 17 quarter timespan. This is accompanied by a slight increase of the star’s precession angle (less negative), but with no significant change in the obliquity of its spin axis. There are suggestions of a history of strong dynamical interaction with a highly distorted planet rotating in a 3:2 resonance with its revolution, together with a tidal lag of \(\sim30~\mbox{deg}\). The mean precessional period is derived to be about 1000 y, but at the present time the motion of the star’s rotation axis appears to be supporting the gravitational torque, rather than providing the balance against it that would be expected over long periods of time.The planet has a small but detectable backwarming effect on the star, which helps to explain the difference in brightness just after transit and just before occultation eclipses. In assessing these findings it is recognized that sources of uncertainty remain, notably with possible inherent micropulsational effects, variations from other components of the multiple star, stellar activity, differential rotation and the neglect of higher order terms (than \(r_{1}^{5}\)) in the fitting function, where \(r_{1}\) is the ratio of the radius of the star to the mean orbital separation of planet and host star.     

Preferred design and error analysis for the future dedicated deep-space Mars-SST satellite gravity mission

Because the precise measurement of the Martian gravitational field plays a significant role in the future Mars exploration program, the future dedicated Mars satellite-to-satellite tracking (Mars-SST) gravity mission in China is investigated in detail for producing the next generation of the Mars gravity field model with high accuracy. Firstly, a new semi-numerical synthetical error model of the cumulative Martian geoid height influenced by the major error sources of the space-borne instruments is precisely established and efficiently verified. Secondly, the deep space network in combination with the satellite-to-satellite tracking in the low-low (DSN-SST-LL) mode is a preferred design owing to the high precision determination of the gravity maps, the low technical complexity of the satellite system and the successful experiences with the Earth’s Gravity Recovery and Climate Experiment (GRACE) projects and the lunar Gravity Recovery and Interior Laboratory (GRAIL) program. Finally, the future twin Mars-SST satellites plan to adopt the optimal matching accuracy indices of the satellite-equipped sensors (e.g., \(10^{-7}\) m/s in the inter-satellite range-rate from the interferometric laser ranging system (ILRS), 35 m in the orbital position tracked by the DSN and \(3\times 10^{-11}\) m/s² in the non-conservative force from the drag-free control system (DFCS)) and the preferred orbital parameters (e.g., the orbital altitude of \(100\pm 50\) km and the inter-satellite range of \(50\pm 10\) km).     

Eruptive Solar Prominence at 37 GHz

On 27 June 2012, an eruptive solar prominence was observed in the extreme ultraviolet (EUV) and radio wavebands. At the Aalto University Metsähovi Radio Observatory (MRO) it was observed at 37 GHz. It was the first time that the MRO followed a radio prominence with dense sampling in the millimetre wavelengths. This prompted us to study the connection of the 37 GHz event with other wavelength domains. At 37 GHz, the prominence was tracked to a height of around \(1.6~\mathrm{R}_{\odot}\), at which the loop structure collapsed. The average velocity of the radio prominence was \(55 \pm 6~\mbox{km}\,\mbox{s}^{-1}\). The brightness temperature of the prominence varied between \(800 \pm 100\) K and \(3200 \pm 100\) K. We compared our data with the Solar Dynamic Observatory (SDO)/Atmospheric Imaging Assembly (AIA) instrument’s 304 Å EUV data, and found that the prominence behaves very similarly in both wavelengths. The EUV data also reveal flaring activity nearby the prominence. We present a scenario in which this flare works as a trigger that causes the prominence to move from a stable stage to an acceleration stage.     

Influence of the Atmospheric Model on Hanle Diagnostics

We clarify the uncertainty in the inferred magnetic field vector via the Hanle diagnostics of the hydrogen Lyman-\(\upalpha\) line when the stratification of the underlying atmosphere is unknown. We calculate the anisotropy of the radiation field with plane-parallel semi-empirical models under the nonlocal thermal equilibrium condition and derive linear polarization signals for all possible parameters of magnetic field vectors based on an analytical solution of the atomic polarization and Hanle effect. We find that the semi-empirical models of the inter-network region (FAL-A) and network region (FAL-F) show similar degrees of anisotropy in the radiation field, and this similarity results in an acceptable inversion error (e.g., \({\sim}\, 40~\mbox{G}\) instead of 50 G in field strength and \({\sim}\,100^{\circ}\) instead of \(90^{\circ}\) in inclination) when FAL-A and FAL-F are swapped. However, the semi-empirical models of FAL-C (averaged quiet-Sun model including both inter-network and network regions) and FAL-P (plage regions) yield an atomic polarization that deviates from all other models, which makes it difficult to precisely determine the magnetic field vector if the correct atmospheric model is not known (e.g., the inversion error is much larger than 40% of the field strength; \({>}\,70~\mbox{G}\) instead of \(50~\mbox{G}\)). These results clearly demonstrate that the choice of model atmosphere is important for Hanle diagnostics. As is well known, one way to constrain the average atmospheric stratification is to measure the center-to-limb variation of the linear polarization signals. The dependence of the center-to-limb variations on the atmospheric model is also presented in this paper.     

Tidal synchronization of close-in satellites and exoplanets. III. Tidal dissipation revisited and application to Enceladus

This paper deals with a new formulation of the creep tide theory (Ferraz-Mello in Celest Mech Dyn Astron 116:109, 2013—Paper I) and with the tidal dissipation predicted by the theory in the case of stiff bodies whose rotation is not synchronous but is oscillating around the synchronous state with a period equal to the orbital period. We show that the tidally forced libration influences the amount of energy dissipated in the body and the average perturbation of the orbital elements. This influence depends on the libration amplitude and is generally neglected in the study of planetary satellites. However, they may be responsible for a 27% increase in the dissipation of Enceladus. The relaxation factor necessary to explain the observed dissipation of Enceladus (\(\gamma =1.2{-}3.8\times 10^{-7}\ \mathrm{s}^{-1}\)) has the expected order of magnitude for planetary satellites and corresponds to the viscosity \(0.6{-}1.9 \times 10^{14}\) Pa s, which is in reasonable agreement with the value recently estimated by Efroimsky (Icarus 300:223, 2018) (\(0.24 \times 10^{14}\) Pa s) and with the value adopted by Roberts and Nimmo (Icarus 194:675, 2008) for the viscosity of the ice shell (\(10^{13}{-}10^{14}\) Pa s). For comparison purposes, the results are extended also to the case of Mimas and are consistent with the negligible dissipation and the absence of observed tectonic activity. The corrections of some mistakes and typos of paper II (Ferraz-Mello in Celest Mech Dyn Astron 122:359, 2015) are included at the end of the paper.     

A study of the geomagnetic indices asymmetry based on the interplanetary magnetic field polarities

Data of geomagnetic indices (aa, Kp, Ap, and Dst) recorded near 1 AU over the period 1967–2016, have been studied based on the asymmetry between the interplanetary magnetic field (IMF) directions above and below of the heliospheric current sheet (HCS). Our results led to the following conclusions: (i) Throughout the considered period, 31 random years (62%) showed apparent asymmetries between Toward (\(\mathbf{T}\)) and Away (\(\mathbf{A}\)) polarity days and 19 years (38%) exhibited nearly a symmetrical behavior. The days of \(\mathbf{A}\) polarity predominated over the \(\mathbf{T}\) polarity days by 4.3% during the positive magnetic polarity epoch (1991–1999). While the days of \(\mathbf{T}\) polarity exceeded the days of \(\mathbf{A}\) polarity by 5.8% during the negative magnetic polarity epoch (2001–2012). (ii) Considerable yearly North–South (N–S) asymmetries of geomagnetic indices observed throughout the considered period. (iii) The largest toward dominant peaks for \(aa\) and \(Ap\) indices occurred in 1995 near to minimum of solar activity. Moreover, the most substantial away dominant peaks for \(aa\) and \(Ap\) indices occurred in 2003 (during the descending phase of the solar cycle 23) and in 1991 (near the maximum of solar activity cycle) respectively. (iv) The N–S asymmetry of \(Kp\) index indicated a most significant away dominant peak occurred in 2003. (v) Four of the away dominant peaks of Dst index occurred at the maxima of solar activity in the years 1980, 1990, 2000, and 2013. The largest toward dominant peak occurred in 1991 (at the reversal of IMF polarity). (vi) The geomagnetic indices (aa, Ap, and \(Kp\)) all have northern dominance during positive magnetic polarity epoch (1971–1979), while the asymmetries shifts to the southern solar hemisphere during negative magnetic polarity epoch (2001–2012).     

Dynamics of Large-Scale Solar-Wind Streams Obtained by the Double Superposed Epoch Analysis: 3. Deflection of the Velocity Vector

This work is a continuation of our previous articles (Yermolaev et al. in J. Geophys. Res.120, 7094, 2015 and Yermolaev et al. in Solar Phys.292, 193, 2017), which describe the average temporal profiles of interplanetary plasma and field parameters in large-scale solar-wind (SW) streams: corotating interaction regions (CIRs), interplanetary coronal mass ejections (ICMEs, including both magnetic clouds (MCs) and ejecta), and sheaths as well as interplanetary shocks (ISs). Changes in the longitude angle, \(\varphi\), in CIRs from ?2 to \(2^{\circ}\) agree with earlier results (e.g. Gosling and Pizzo, 1999). We have also analyzed the average temporal profiles of the bulk velocity angles in sheaths and ICMEs. We have found that the angle \(\varphi\) in ICMEs changes from 2 to \(-2^{\circ}\), while in sheaths it changes from ?2 to \(2^{\circ}\) (similar to the change in CIRs), i.e. the angle in CIRs and sheaths deflects in the opposite sense to ICMEs. When averaging the latitude angle \(\vartheta\) on all the intervals of the chosen SW types, the angle \(\vartheta\) is almost constant at \({\sim}\,1^{\circ}\). We made for the first time a selection of SW events with increasing and decreasing \(\vartheta\) and found that the average \(\vartheta\) temporal profiles in the selected events have the same “integral-like” shape as for \(\varphi\). The difference in \(\varphi\) and \(\vartheta\) average profiles is explained by the fact that most events have increasing profiles for the angle in the ecliptic plane as a result of solar rotation, while for the angle in the meridional plane, the numbers of events with increasing and decreasing profiles are equal.     

Constraints on dissipation in the deep interiors of Ganymede and Europa from tidal phase-lags

Jupiter’s satellites are subject to strong tidal forces which result in variations of the gravitational potential and deformations of the satellites’ surfaces on the diurnal tidal cycle. Such variations are described by the Love numbers \(k_2\) and \(h_2\) for the tide-induced potential variation due to internal mass redistribution and the radial surface displacement, respectively. The phase-lags \( \phi _{k_2}\) and \( \phi _{h_2}\) of these complex numbers contain information about the rheological and dissipative states of the satellites. Starting from interior structure models and assuming a Maxwell rheology to compute the tidal deformation, we calculate the phase-lags in application to Ganymede and Europa. For both satellites we assume a decoupling of the outer ice-shell from the deep interior by a liquid subsurface water ocean. We show that, in this case, the phase-lag difference \(\varDelta \phi = \phi _{k_2}- \phi _{h_2}\) can provide information on the rheological and thermal state of the deep interiors if the viscosities of the deeper layers are small. In case of Ganymede, phase-lag differences can reach values of a few degrees for high-pressure ice viscosities \({<}10^{14}\) Pa s and would indicate a highly dissipative state of the deep interior. In this case \(\varDelta \phi \) is dominated by dissipation in the high-pressure ice layer rather than dissipation within the ice-I shell. These phase lags would be detectable from spacecraft in orbit around the satellite. For Europa \(\varDelta \phi \) could reach values exceeding \(20^\circ \) and phase-lag measurements could help distinguish between (1) a hot dissipative silicate mantle which would in thermal equilibrium correspond to a very thin outer ice-I shell and (2) a cold deep interior implying that dissipation would mainly occur in a thick (several tens of km) outer ice-I shell. These measurements are highly relevant for ESA’s Jupiter Icy Moons Explorer (JUICE) and NASA’s Europa Multiple Flyby Mission, both targeted for the Jupiter system.     

Subsurface Zonal and Meridional Flow During Cycles 23 and 24

We study the solar-cycle variation of subsurface flows from the surface to a depth of 16 Mm. We have used ring-diagram analysis to analyze Dopplergrams obtained with the Michelson Doppler Imager (MDI) Dynamics Program, the Global Oscillation Network Group (GONG), and the Helioseismic and Magnetic Imager (HMI) instrument. We combined the zonal and meridional flows from the three data sources and scaled the flows derived from MDI and GONG to match those from HMI observations. In this way, we derived their temporal variation in a consistent manner for Solar Cycles 23 and 24. We have corrected the measured flows for systematic effects that vary with disk positions. Using time-depth slices of the corrected subsurface flows, we derived the amplitudes and times of the extrema of the fast and slow zonal and meridional flows during Cycles 23 and 24 at every depth and latitude. We find an average difference between maximum and minimum amplitudes of \(8.6 \pm0.4~\mbox{m}\,\mbox{s}^{-1}\) for the zonal flows and \(7.9 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the meridional flows associated with Cycle 24 averaged over a depth range from 2 to 12 Mm. The corresponding values derived from GONG data alone are \(10.5 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the zonal and \(10.8 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the meridional flow. For Cycle 24, the flow patterns are precursors of the magnetic activity. The timing difference between the occurrence of the flow pattern and the magnetic one increases almost linearly with increasing latitude. For example, the fast zonal and meridional flow appear \(2.1 \pm 0.6\) years and \(2.5\pm 0.6\) years, respectively, before the magnetic pattern at \(30^{\circ}\) latitude in the northern hemisphere, while in the southern hemisphere, the differences are \(3.2 \pm 1.2\) years and \(2.6 \pm 0.6\) years. The flow patterns of Cycle 25 are present and have reached \(30^{\circ}\) latitude. The amplitude differences of Cycle 25 are about 22% smaller than those of Cycle 24, but are comparable to those of Cycle 23. Moreover, polynomial fits of meridional flows suggest that equatorward meridional flows (counter-cells) might exist at about \(80^{\circ}\) latitude except during the declining phase of the solar cycle.     

The Pleiades apex and its kinematical structure

A study of cluster characteristics and internal kinematical structure of the middle-aged Pleiades open star cluster is presented. The individual star apexes and various cluster kinematical parameters including the velocity ellipsoid parameters are determined using both Hipparcos and Gaia data. Modern astrometric parameters were taken from the Gaia Data Release 1 (DR1) in combination with the Radial Velocity Experiment Fifth Data Release (DR5). The necessary set of parameters including parallaxes, proper motions and radial velocities are used for \(n=17\) stars from Gaia DR1+RAVE DR5 and for \(n=19\) stars from the Hipparcos catalog using SIMBAD data base. Single stars are used to improve accuracy by eliminating orbital movements. RAVE DR5 measurements were taken only for the stars with the radial velocity errors not exceeding \(2~\mbox{km}/\mbox{s}\). For the Pleiades stars taken from Gaia, we found mean heliocentric distance as \(136.8 \pm 6.4\) pc, and the apex position is calculated as: \(A_{CP}=92^{\circ }.52\pm 1^{\circ }.72\), \(D_{CP}=-42^{\circ }.28\pm 2^{\circ }.56\) by the convergent point method and \(A_{0}=95^{\circ }.59\pm 2^{\circ }.30\) and \(D_{0}=-50^{\circ }.90\pm 2^{\circ }.04\) using AD-diagram method (\(n=17\) in both cases). The results are compared with those obtained historically before the Gaia mission era.     

Accretion flow dynamics during 1999 outburst of XTE J1859+226—modeling of broadband spectra and constraining the source mass

We examine the dynamical behavior of accretion flow around XTE J1859+226 during the 1999 outburst by analyzing the entire outburst data (～166 days) from RXTE Satellite. Towards this, we study the hysteresis behavior in the hardness intensity diagram (HID) based on the broadband (3–150 keV) spectral modeling, spectral signature of jet ejection and the evolution of Quasi-periodic Oscillation (QPO) frequencies using the two-component advective flow model around a black hole. We compute the flow parameters, namely Keplerian accretion rate (\({\dot{m}}_{d}\)), sub-Keplerian accretion rate (\({\dot{m}}_{h}\)), shock location (\(r_{s}\)) and black hole mass (\(M_{\mathit{bh}}\)) from the spectral modeling and study their evolution along the q-diagram. Subsequently, the kinetic jet power is computed as \(L^{\mathrm{obs}}_{\mathrm{jet}} \sim3\mbox{--}6 \times10^{37}~\mbox{erg}\,\mbox{s}^{-1}\) during one of the observed radio flares which indicates that jet power corresponds to 8–16% mass outflow rate from the disc. This estimate of mass outflow rate is in close agreement with the change in total accretion rate (～14%) required for spectral modeling before and during the flare. Finally, we provide a mass estimate of the source XTE J1859+226 based on the spectral modeling that lies in the range of 5.2–7.9 \(M_{\odot}\) with 90% confidence.     

Study and application of the resonant secular dynamics beyond Neptune

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (Celest Mech Dyn Astron, doi: 10.1007/s10569-016-9700-5, 2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or \(2012\text {VP}_{113}\), which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at \(a\in [100;300]\) AU, \(q\in [50;70]\) AU and \(I\in [30;50]^{\circ }\). In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.     

Investigation of the Geoeffectiveness of Disk-Centre Full-Halo Coronal Mass Ejections

We studied the occurrence and characteristics of geomagnetic storms associated with disk-centre full-halo coronal mass ejections (DC-FH-CMEs). Such coronal mass ejections (CMEs) can be considered as the most plausible cause of geomagnetic storms. We selected front-side full-halo coronal mass ejections detected by the Large Angle and Spectrometric Coronagraph onboard the Solar and Heliospheric Observatory (SOHO/LASCO) from the beginning of 1996 till the end of 2015 with source locations between solar longitudes E10 and W10 and latitudes N20 and S20. The number of selected CMEs was 66 of which 33 (50%) were deduced to be the cause of 30 geomagnetic storms with \(\mathrm{Dst} \leq- 50~\mbox{nT}\). Of the 30 geomagnetic storms, 26 were associated with single disk-centre full-halo CMEs, while four storms were associated, in addition to at least one disk-centre full-halo CME, also with other halo or wide CMEs from the same active region. Thirteen of the 66 CMEs (20%) were associated with 13 storms with \(-100~\mbox{nT} < \mbox{Dst} \leq- 50~\mbox{nT}\), and 20 (30%) were associated with 17 storms with \(\mbox{Dst}\leq- 100~\mbox{nT}\). We investigated the distributions and average values of parameters describing the DC-FH-CMEs and their interplanetary counterparts encountering Earth. These parameters included the CME sky-plane speed and direction parameter, associated solar soft X-ray flux, interplanetary magnetic field strength, \(B_{t}\), southward component of the interplanetary magnetic field, \(B_{s}\), solar wind speed, \(V_{sw}\), and the \(y\)-component of the solar wind electric field, \(E_{y}\). We found only a weak correlation between the Dst of the geomagnetic storms associated with DC-FH-CMEs and the CME sky-plane speed and the CME direction parameter, while the correlation was strong between the Dst and all the solar wind parameters (\(B_{t}\), \(B_{s}\), \(V_{sw}\), \(E_{y}\)) measured at 1 AU. We investigated the dependences of the properties of DC-FH-CMEs and the associated geomagnetic storms on different phases of solar cycles and the differences between Solar Cycles 23 and 24. In the rise phase of Solar Cycle 23 (SC23), five out of eight DC-FH-CMEs were geoeffective (\(\mbox{Dst} \leq- 50~\mbox{nT}\)). In the corresponding phase of SC24, only four DC-FH-CMEs were observed, three of which were nongeoeffective (\(\mbox{Dst} > - 50~\mbox{nT}\)). The largest number of DC-FH-CMEs occurred at the maximum phases of the cycles (21 and 17, respectively). Most of the storms with \(\mbox{Dst}\leq- 100~\mbox{nT}\) occurred at or close to the maximum phases of the cycles. When comparing the storms during epochs of corresponding lengths in Solar Cycles 23 and 24, we found that during the first 85 months of Cycle 23 the geoeffectiveness rate of the disk-centre full-halo CMEs was 58% with an average minimum value of the Dst index of \(- 146~\mbox{nT}\). During the corresponding epoch of Cycle 24, only 35% of the disk-centre full-halo CMEs were geoeffective with an average value of Dst of \(- 97~\mbox{nT}\).     

Local stellar kinematics from RAVE data—VIII. Effects of the Galactic disc perturbations on stellar orbits of red clump stars

We aim to probe the dynamic structure of the extended Solar neighborhood by calculating the radial metallicity gradients from orbit properties, which are obtained for axisymmetric and non-axisymmetric potential models, of red clump (RC) stars selected from the RAdial Velocity Experiment’s Fourth Data Release. Distances are obtained by assuming a single absolute magnitude value in near-infrared, i.e. \(M_{Ks}=-1.54\pm0.04\) mag, for each RC star. Stellar orbit parameters are calculated by using the potential functions: (i) for the MWPotential2014 potential, (ii) for the same potential with perturbation functions of the Galactic bar and transient spiral arms. The stellar age is calculated with a method based on Bayesian statistics. The radial metallicity gradients are evaluated based on the maximum vertical distance (\(z_{max}\)) from the Galactic plane and the planar eccentricity (\(e_{p}\)) of RC stars for both of the potential models. The largest radial metallicity gradient in the \(0< z_{max} \leq0.5\) kpc distance interval is \(-0.065\pm0.005~\mbox{dex}\,\mbox{kpc}^{-1}\) for a subsample with \(e_{p}\leq0.1\), while the lowest value is \(-0.014\pm0.006~\mbox{dex}\,\mbox{kpc}^{-1}\) for the subsample with \(e_{p}\leq0.5\). We find that at \(z_{max}>1\) kpc, the radial metallicity gradients have zero or positive values and they do not depend on \(e_{p}\) subsamples. There is a large radial metallicity gradient for thin disc, but no radial gradient found for thick disc. Moreover, the largest radial metallicity gradients are obtained where the outer Lindblad resonance region is effective. We claim that this apparent change in radial metallicity gradients in the thin disc is a result of orbital perturbation originating from the existing resonance regions.     

LOFAR Observations of Fine Spectral Structure Dynamics in Type IIIb Radio Bursts

Solar radio emission features a large number of fine structures demonstrating great variability in frequency and time. We present spatially resolved spectral radio observations of type IIIb bursts in the 30?–?80 MHz range made by the Low Frequency Array (LOFAR). The bursts show well-defined fine frequency structuring called “stria” bursts. The spatial characteristics of the stria sources are determined by the propagation effects of radio waves; their movement and expansion speeds are in the range of \((0.1\,\mbox{--}\,0.6)c\). Analysis of the dynamic spectra reveals that both the spectral bandwidth and the frequency drift rate of the striae increase with an increase of their central frequency. The striae bandwidths are in the range of \({\approx}\,(20\,\mbox{--}\,100)\) kHz and the striae drift rates vary from zero to \({\approx}\,0.3~\mbox{MHz}\,\mbox{s}^{-1}\). The observed spectral characteristics of the stria bursts are consistent with the model involving modulation of the type III burst emission mechanism by small-amplitude fluctuations of the plasma density along the electron beam path. We estimate that the relative amplitude of the density fluctuations is of \(\Delta n/n\sim10^{-3}\), their characteristic length scale is less than 1000 km, and the characteristic propagation speed is in the range of \(400\,\mbox{--}\,800~\mbox{km}\,\mbox{s}^{-1}\). These parameters indicate that the observed fine spectral structures could be produced by propagating magnetohydrodynamic waves.     

Addendum to: Strength of the Solar Coronal Magnetic Field – A Comparison of Independent Estimates Using Contemporaneous Radio and White-Light Observations

This addendum uses an alternate fit for the electron density distribution \(N(r)\) (see Figure 1) and estimates the coronal magnetic field using the new model. We find that the estimates of the magnetic field are in close agreement using both the models.

We have fit the \(N(r)\) distribution obtained from STEREO-A/COR1 and SOHO/LASCO-C2 using a fifth-order polynomial (see Figure 1). The expression can be written as

$$\begin{aligned} N_{\text{cor}}(r) &= 1.43 \times 10^{9} r^{-5} - 1.91 \times 10^{9} r^{-4} + 1.07 \times 10^{9} r^{-3} - 2.87 \times 10^{8} r^{-2} \\ &\quad {} + 3.76 \times 10^{7} r^{-1} - 1.91 \times 10^{6} , \end{aligned}$$

(1)

where \(N_{\text{cor}}(r)\) is in units of cm^?3 and \(r\) is in units of \(\mathrm{R}_{\odot}\). The background coronal electron density is enhanced by a factor of 5.5 at 2.63 \(\mathrm{R}_{\odot}\) during the coronal mass ejection (CME). The estimated coronal magnetic field strength (\(B\)) using radio data indicates that \(B(r) \approx(0.51\text{\,--\,}0.48) \pm 0.02\ \mathrm{G}\) in the range \(r \approx2.65\text{\, --\,}2.82\ \mathrm{R}_{\odot}\). The field strengths for STEREO-A/COR1 and SOHO/LASCO-C2 are ≈?0.32 G at \(r \approx 3.11\ \mathrm{R}_{\odot}\) and ≈?0.12 G at \(r \approx 4.40\ \mathrm{R}_{\odot}\), respectively.

     

Interstellar communication: The colors of optical SETI

It has recently been argued from a laser engineering point of view that there are only a few magic colors for optical SETI. These are primarily the Nd:YAG line at \(1{,}064\,\hbox {nm}\) and its second harmonic (532.1 nm). Next best choices would be the sum frequency and/or second harmonic generation of Nd:YAG and Nd:YLF laser lines, 393.8 nm (near Fraunhofer CaK), 656.5 nm (\(\hbox {H}\alpha \)) and 589.1 nm (NaD2). In this paper, we examine the interstellar extinction, atmospheric transparency and scintillation, as well as noise conditions for these laser lines. For strong signals, we find that optical wavelengths are optimal for distances \(d\lesssim \,\hbox {kpc}\). Nd:YAG at \(\lambda =1{,}064\,\hbox {nm}\) is a similarly good choice, within a factor of two, under most conditions and out to \(d\lesssim 3\,\hbox {kpc}\). For weaker transmitters, where the signal-to-noise ratio with respect to the blended host star is relevant, the optimal wavelength depends on the background source, such as the stellar type. Fraunhofer spectral lines, while providing lower stellar background noise, are irrelevant in most use cases, as they are overpowered by other factors. Laser-pushed spaceflight concepts, such as “Breakthrough Starshot”, would produce brighter and tighter beams than ever assumed for OSETI. Such beamers would appear as naked eye stars out to kpc distances. If laser physics has already matured and converged on the most efficient technology, the laser line of choice for a given scenario (e.g., Nd:YAG for strong signals) can be observed with a narrow filter to dramatically reduce background noise, allowing for large field-of-view observations in fast surveys.     

Turbulence and Heating in the Flank and Wake Regions of a Coronal Mass Ejection

As a coronal mass ejection (CME) passes, the flank and wake regions are typically strongly disturbed. Various instruments, including the Large Angle and Spectroscopic Coronagraph (LASCO), the Atmospheric Imaging Assembly (AIA), and the Coronal Multi-channel Polarimeter (CoMP), observed a CME close to the east limb on 26 October 2013. A hot (\({\approx}\,10~\mbox{MK}\)) rising blob was detected on the east limb, with an initial ejection flow speed of \({\approx}\, 330~\mbox{km}\,\mbox{s}^{-1}\). The magnetic structures on both sides and in the wake of the CME were strongly distorted, showing initiation of turbulent motions with Doppler-shift oscillations enhanced from \({\approx}\, \pm 3~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\, \pm 15~\mbox{km}\,\mbox{s}^{-1}\) and effective thermal velocities from \({\approx}\,30~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\,60~\mbox{km}\,\mbox{s}^{-1}\), according to the CoMP observations at the Fe?xiii line. The CoMP Doppler-shift maps suggest that the turbulence behaved differently at various heights; it showed clear wave-like torsional oscillations at lower altitudes, which are interpreted as the antiphase oscillation of an alternating red/blue Doppler shift across the strands at the flank. The turbulence seems to appear differently in the channels of different temperatures. Its turnover time was \({\approx}\,1000\) seconds for the Fe 171 Å channel, while it was \({\approx}\,500\) seconds for the Fe 193 Å channel. Mainly horizontal swaying rotations were observed in the Fe 171 Å channel, while more vertical vortices were seen in the Fe 193 Å channel. The differential-emission-measure profiles in the flank and wake regions have two components that evolve differently: the cool component decreased over time, evidently indicating a drop-out of cool materials due to ejection, while the hot component increased dramatically, probably because of the heating process, which is suspected to be a result of magnetic reconnection and turbulence dissipation. These results suggest a new turbulence-heating scenario of the solar corona and solar wind.