Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory/Helioseismic and Magnetic Imager Time?–?Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a?time?–?distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time?–?distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a?minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference travel-time perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.     

Observational Evidence for a Double-Helix Structure in CMEs and Magnetic Clouds

We compare recent observations of a solar eruptive prominence as seen in extreme-UV light on 30 March 2010 by the Solar Dynamics Observatory (SDO) with the multi-tube model for interplanetary magnetic clouds (Osherovich, Fainberg, Stone, Geophys. Res. Lett. 26, 2597, 1999). Our model is based on an exact analytical solution of the plasma equilibrium with magnetic force balanced by a gradient of scalar gas pressure. Topologically, this solution describes two magnetic helices with opposite magnetic polarity embedded in a cylindrical magnetic flux tube that creates magnetic flux inequality between the two helices by enhancing one helix and suppressing the other. The magnetic field in this model is continuous everywhere and has a finite magnetic energy per unit length of the tube. These configurations have been introduced as MHD bounded states (Osherovich, Soln. Dannye 5, 70, 1975). Apparently, the SDO observations depict two non-equal magnetically interacting helices described by this analytical model. We consider magnetic and thermodynamic signatures of multiple magnetic flux ropes inside the same magnetic cloud, using in situ observations. The ratio of magnetic energy density to bulk speed solar wind energy density has been defined as a solar wind quasi-invariant (QI). We analyze the structure of the QI profile to probe the topology of the internal structure of magnetic clouds. From the superposition of 12 magnetically isolated clouds observed by Ulysses, we have found that the corresponding QI is consistent with our double helix model.     

The Evolution of Barbs of a Polar Crown Filament Observed by SDO

From 16 to 21 August 2010, a northern (???N60) polar crown filament was observed by Solar Dynamics Observatory (SDO). Employing the six-day SDO/AIA data, we identify 69 barbs, and select 58 of them, which appeared away from the western solar limb (???W60), as our sample. We systematically investigate the evolution of filament barbs. Three different types of apparent formation of barbs are detected, including i)?the convergence of surrounding moving plasma condensations, comprised 55.2?% of our sample, ii)?the flows of plasma condensations from the filament, comprised 37.9?%, and iii)?the plasma injections from the neighboring brightening regions, comprised 6.9?%. We also find three different ways that barb disappear, involving: i)?bi-lateral movements (44.8?%), and ii)?outflowing of barb plasma (27.6?%) results in the disappearance of a barb, as well as iii)?disappearance of a barb is associated with a neighboring brightening (27.6?%). The evolution of the magnetic fields, e.g. emergence and cancellation of magnetic flux, may cause the formation or disappearance of the barb magnetic structures. Barbs exchange plasma condensations with the surrounding atmosphere, filament, and nearby brightenings, leading to the increase or drainage of barb material. Furthermore, we find that all the barbs undergo oscillations. The average oscillation period, amplitude, and velocity are 30?min, 2.4?Mm, and 5.7?km?s^?1, respectively. Besides the oscillations, 21 (36?%) barbs manifested sideward motions having an average speed of 0.45?km?s^?1. Small-scale wave-like propagating disturbances caused by small-scale brightenings are detected, and the barb oscillations associated with these disturbances are also found. We propose that the kinematics of barbs are influenced or even caused by the evolution of the neighboring photospheric magnetic fields.     

Very Fast and Accurate Azimuth Disambiguation of Vector Magnetograms

We present a method for fast and accurate azimuth disambiguation of vector magnetogram data regardless of the location of the analyzed region on the solar disk. The direction of the transverse field is determined with the principle of minimum deviation of the field from the reference (potential) field. The new disambiguation (NDA) code is examined on the well-known models of Metcalf et al. (Solar Phys. 237, 267, 2006) and Leka et al. (Solar Phys. 260, 83, 2009), and on an artificial model based on the observed magnetic field of AR 10930 (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We compare Hinode/SOT-SP vector magnetograms of AR 10930 disambiguated with three codes: the NDA code, the nonpotential magnetic-field calculation (NPFC: Georgoulis, Astrophys. J. Lett. 629, L69, 2005), and the spherical minimum-energy method (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We then illustrate the performance of NDA on SDO/HMI full-disk magnetic-field observations. We show that our new algorithm is more than four times faster than the fastest algorithm that provides the disambiguation with a satisfactory accuracy (NPFC). At the same time, its accuracy is similar to that of the minimum-energy method (a very slow algorithm). In contrast to other codes, the NDA code maintains high accuracy when the region to be analyzed is very close to the limb.     

On the Origin of Pulsations of Sub-THz Emission from Solar Flares

We propose a model to explain fast pulsations in sub-THz emission from solar flares. The model is based on the approach of a flaring loop as an equivalent electric circuit and explains the pulse-repetition rate, the high-quality factor, Q≥10³, low modulation depth, pulse synchronism at different frequencies, and the dependence of the pulse-repetition rate on the emission flux, observed by Kaufmann et al. (Astrophys. J. 697, 420, 2009). We solved the nonlinear equation for electric current oscillations using a Van der Pol method and found the steady-state value for the amplitude of the current oscillations. Using the pulse rate variation during the flare on 4 November 2003, we found a decrease of the electric current from 1.7×10¹² A in the flare maximum to 4×10¹⁰ A just after the burst. Our model is consistent with the plasma mechanism of sub-THz emission suggested recently by Zaitsev, Stepanov, and Melnikov (Astron. Lett. 39, 650, 2013).     

Transient Artifacts in a Flare Observed by the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory

The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) provides a new tool for the systematic observation of white-light flares, including Doppler and magnetic information as well as continuum. In our initial analysis of the highly impulsive $\mathrm{\gamma}$ -ray flare SOL2010-06-12T00:57 (Martínez Oliveros et al., Solar Phys. 269, 269, 2011), we reported the signature of a strong blueshift in the two footpoint sources. Concerned that this might be an artifact due to aliasing peculiar to the HMI instrument, we undertook a comparative analysis of Global Oscillation Network Group (GONG++) observations of the same flare, using the PArametric Smearing Correction ALgorithm (PASCAL) algorithm to correct for artifacts caused by variations in atmospheric smearing. This analysis confirms the artifactual nature of the apparent blueshift in the HMI observations, finding weak redshifts at the footpoints instead. We describe the use of PASCAL with GONG++ observations as a complement to the SDO observations and discuss constraints imposed by the use of HMI far from its design conditions. With proper precautions, these data provide rich information on flares and transients.     

Preprocessing the Photospheric Vector Magnetograms for an NLFFF Extrapolation Using a Potential-Field Model and an Optimization Method

Numerical reconstruction/extrapolation of the coronal nonlinear force-free magnetic field (NLFFF) usually takes the photospheric vector magnetogram as input at the bottom boundary. The magnetic field observed at the photosphere, however, contains a force that is in conflict with the fundamental assumption of the force-free model. It also contains measurement noise, which hinders the practical computation. Wiegelmann, Inhester, and Sakurai (Solar Phys. 233, 215, 2006) have proposed to preprocess the raw magnetogram to remove the force and noise to provide better input for NLFFF modeling. In this paper we develop a new code of magnetogram preprocessing that is consistent with our extrapolation method CESE–MHD–NLFFF (Jiang, Feng, and Xiang in Astrophys. J. 755, 62, 2012; Jiang and Feng in Astrophys. J. 749, 135, 2012a). Based on the magnetic-splitting rule that a magnetic field can be split into a potential-field part and a non-potential part, we split the magnetogram and dealt with the two parts separately. The preprocessing of the magnetogram’s potential part is based on a numerical potential-field model, and the non-potential part is preprocessed using the similar optimization method of Wiegelmann, Inhester, and Sakurai (2006). The code was applied to the SDO/HMI data, and results show that the method can remove the force and noise efficiently and improve the extrapolation quality.     

Torsional Oscillations in a Global Solar Dynamo

We characterize and analyze rotational torsional oscillations developing in a large-eddy magnetohydrodynamical simulation of solar convection (Ghizaru, Charbonneau, and Smolarkiewicz, Astrophys. J. Lett. 715, L133, 2010; Racine et al., Astrophys. J. 735, 46, 2011) producing an axisymmetric, large-scale, magnetic field undergoing periodic polarity reversals. Motivated by the many solar-like features exhibited by these oscillations, we carry out an analysis of the large-scale zonal dynamics. We demonstrate that simulated torsional oscillations are not driven primarily by the periodically varying large-scale magnetic torque, as one might have expected, but rather via the magnetic modulation of angular-momentum transport by the large-scale meridional flow. This result is confirmed by a straightforward energy analysis. We also detect a fairly sharp transition in rotational dynamics taking place as one moves from the base of the convecting layers to the base of the thin tachocline-like shear layer formed in the stably stratified fluid layers immediately below. We conclude by discussing the implications of our analyses with regard to the mechanism of amplitude saturation in the global dynamo operating in the simulation, and speculate on the possible precursor value of torsional oscillations for the forecast of solar-cycle characteristics.     

Absolute magnitude calibration for red clump stars

We combined the (K _s , J?K _s ) data in Laney et al. (Mon. Not. R. Astron. Soc. 419:1637, 2012) with the V apparent magnitudes and trigonometric parallaxes taken from the Hipparcos catalogue and used them to fit the $M_{K_{s}}$ absolute magnitude to a linear polynomial in terms of V?K _s colour. The mean and standard deviation of the absolute magnitude residuals, ?0.001 and 0.195 mag, respectively, estimated for 224 red clump stars in Laney et al. (2012) are (absolutely) smaller than the corresponding ones estimated by the procedure which adopts a mean $M_{K_{s}}=-1.613~\mbox{mag}$ absolute magnitude for all red clump stars, ?0.053 and 0.218 mag, respectively. The statistics estimated by applying the linear equation to the data of 282 red clump stars in Alves (Astrophys. J. 539:732, 2000) are larger, $\Delta M_{K_{s}}=0.209$ and σ=0.524 mag, which can be explained by a different absolute magnitude trend, i.e. condensation along a horizontal distribution.     

Inclusion of sdBs in evolutionary population synthesis for binary stellar populations and the application: the determination of photo-z and galaxy morphology

Subdwarf B stars (sdBs) can significantly change the ultraviolet spectra of populations at age t～1 Gyr, and have been even included in the evolutionary population synthesis (EPS) models by Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007). In this study we present the spectral energy distributions (SEDs) of binary stellar populations (BSPs) by combining the EPS models of Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007) and those of the Yunnan group (Zhang et al. in Astron. Astrophys. 415:117, 2004; Mon. Not. R. Astron. Soc. 357:1088, 2005), which have included various binary interactions (except sdBs) in EPS models. This set of SEDs is available upon request from the authors. Using this set of SEDs of BSPs we build the spectra of Burst, E, S0–Sd and Irr types of galaxies by using the package of Bruzual and Charlot (Mon. Not. R. Astron. Soc. 344:1000, 2003). Combined with the photometric data (filters and magnitudes), we obtain the photometric redshifts and morphologies of 1502 galaxies by using the Hyperz code of Bolzonella et al. (Astron. Astrophys. 363:476, 2000). This sample of galaxies is obtained by removing those objects, mismatched with the SDSS/DR7 and GALEX/DR4, from the catalogue of Fukugita et al. (Astron. J. 134:579, 2007). By comparison the results with the SDSS spectroscopic redshifts and the morphological index of Fukugita et al. (Astron. J. 134:579, 2007), we find that the photo-z fluctuate with the SDSS spectroscopic redshifts, while the Sa–Sc galaxies in the catalogue of Fukugita et al. (Astron. J. 134:579, 2007) are classified earlier as Burst-E galaxies.     

Influence of the coorbital resonance on the rotation of the Trojan satellites of Saturn

The Cassini spacecraft collects high resolution images of the Saturnian satellites and reveals the surface of these new worlds. Tiscareno et?al. succeeded to determine the Epimetheus rotation from the Cassini Imaging Science Subsystem data, initiating studies on the rotation of Epimetheus and its companion Janus (Tiscareno et?al., Icarus 204:254?C261, 2009; Noyelles, Icarus 207:887?C902, 2010; Robutel et?al., Icarus 211:758?C769, 2011). Especially, Epimetheus is characterized by its horseshoe shape orbit and the presence of the swap has to be introduced explicitly into rotational models. During its journey in the Saturnian system, Cassini spacecraft accumulates the observational data of the other satellites and it will be possible to determine the rotational parameters of several of them. To prepare these future observations, we built rotational models of the coorbital (also called Trojan) satellites Telesto, Calypso, Helene, and Polydeuces, in addition to Janus and Epimetheus. Indeed, Telesto and Calypso orbit around the L ₄ and L ₅ Lagrange points of Saturn-Tethys while Helene and Polydeuces are coorbital of Dione. The goal of this study is to understand how the departure from the Keplerian motion induced by the perturbations of the coorbital body, influences the rotation of these satellites. To this aim, we introduce explicitly the perturbation in the rotational equations by using the formalism developed by érdi (Celest Mech 15:367?C383, 1977) to represent the coorbital motions, and so we describe the rotational motion of the coorbitals, Janus and Epimetheus included, in compact form.     

On the Structure and Evolution of a Polar Crown Prominence/Filament System

Polar crown prominences, that partially circle the Sun’s poles between 60° and 70° latitude, are made of chromospheric plasma. We aim to diagnose the 3D dynamics of a polar crown prominence using high-cadence EUV images from the Solar Dynamics Observatory (SDO)/AIA at 304, 171, and 193 Å and the Ahead spacecraft of the Solar Terrestrial Relations Observatory (STEREO-A)/EUVI at 195 Å. Using time series across specific structures, we compare flows across the disk in 195 Å with the prominence dynamics seen on the limb. The densest prominence material forms vertical columns that are separated by many tens of Mm and connected by dynamic bridges of plasma that are clearly visible in 304/171 Å two-colour images. We also observe intermittent but repetitious flows with velocity 15 km?s^?1 in the prominence that appear to be associated with EUV bright points on the solar disk. The boundary between the prominence and the overlying cavity appears as a sharp edge. We discuss the structure of the coronal cavity seen both above and around the prominence. SDO/HMI and GONG magnetograms are used to infer the underlying magnetic topology. The evolution and structure of the prominence with respect to the magnetic field seems to agree with the filament-linkage model.     

Magneto-Acoustic Energetics Study of the Seismically Active Flare of 15 February 2011

Multi-wavelength studies of energetic solar flares with seismic emissions have revealed interesting common features between them. We studied the first GOES X-class flare of Solar Cycle 24, as detected by the Solar Dynamics Observatory (SDO). For context, seismic activity from this flare (SOL2011-02-15T01:55-X2.2, in NOAA AR 11158) has been reported by Kosovichev (Astrophys. J. Lett., 734, L15, 2011) and Zharkov et?al. (Astrophys. J. Lett., 741, L35, 2011). Based on Dopplergram data from the Helioseismic and Magnetic Imager (HMI), we applied standard methods of local helioseismology in order to identify the seismic sources in this event. RHESSI hard X-ray data are used to check the correlation between the location of the seismic sources and the particle-precipitation sites in during the flare. Using HMI magnetogram data, the temporal profile of fluctuations in the photospheric line-of-sight magnetic field is used to estimate the magnetic-field change in the region where the seismic signal was observed. This leads to an estimate of the work done by the Lorentz-force transient on the photosphere of the source region. In this instance, this is found to be a significant fraction of the acoustic energy in the attendant seismic emission, suggesting that Lorentz forces can contribute significantly to the generation of sunquakes. However, there are regions in which the signature of the Lorentz force is much stronger, but from which no significant acoustic emission emanates.     

Piecewise mass flows within a solar prominence observed by the New Vacuum Solar Telescope

The material of solar prominences is often observed in a state of flowing. These mass flows (MF) are important and useful for us to understand the internal structure and dynamics of prominences. In this paper, we present a high resolution H\(\alpha \) observation of MFs within a quiescent solar prominence. From the observation, we find that the plasma primarily has a circular motion and a downward motion separately in the middle section and legs of the prominence, which creates a piecewise mass flow along the observed prominence. Moreover, the observation also shows a clear displacement of MF’s velocity peaks in the middle section of the prominence. All of these provide us with a detailed record of MFs within a solar prominence and show a new approach to detecting the physical properties of prominence.     

Nonsingular charged analogues of Schwarzschild’s interior solution

The problem of finding nonsingular charged analogue of Schwarzschild’s interior solutions has been reduced to that of finding a monotonically decreasing function f. The models are discussed in generality by imposing reality condition on f. It is shown that the physical solutions are possible only for surface density to central density ratio greater than or equal to 2/3 i.e. $\frac{\rho_{a}}{\rho_{0}}\ge2/3$ . The unphysical nature of solutions with linear equation state has been proved. A generalization procedure has been utilized to generalize solutions by Guilfoyle (1999). Recently found solutions by Gupta and Kumar (2005a, 2005b, 2005c) are generalized by taking particular form of f and seen to have higher mass and more stable. The maximum mass is found to be 1.59482 M _Θ . The models have been found to be stable once the physical requirements are established due to mass to radius less than 4/9, total charge to total mass ratio less than 1 and redshift quite low.     

Capture probability in the 3:1 mean motion resonance with Jupiter: an application to the Vesta family

We study the capture and crossing probabilities in the 3:1 mean motion resonance with Jupiter for a small asteroid that migrates from the inner to the middle Main Belt under the action of the Yarkovsky effect. We use an algebraic mapping of the averaged planar restricted three-body problem based on the symplectic mapping of Hadjidemetriou (Celest Mech Dyn Astron 56:563–599, 1993), adding the secular variations of the orbit of Jupiter and non-symplectic terms to simulate the migration. We found that, for fast migration rates, the captures occur at discrete windows of initial eccentricities whose specific locations depend on the initial resonant angles, indicating that the capture phenomenon is not probabilistic. For slow migration rates, these windows become narrower and start to accumulate at low eccentricities, generating a region of mutual overlap where the capture probability tends to 100 %, in agreement with the theoretical predictions for the adiabatic regime. Our simulations allow us to predict the capture probabilities in both the adiabatic and non-adiabatic cases, in good agreement with results of Gomes (Celest Mech Dyn Astron 61:97–113, 1995) and Quillen (Mon Not RAS 365:1367–1382, 2006). We apply our model to the case of the Vesta asteroid family in the same context as Roig et al. (Icarus 194:125–136, 2008), and found results indicating that the high capture probability of Vesta family members into the 3:1 mean motion resonance is basically governed by the eccentricity of Jupiter and its secular variations.     

Stability of triangular equilibrium points in the elliptic restricted problem of three bodies with radiating and triaxial primaries

This paper studies the stability of infinitesimal motions about the triangular equilibrium points in the elliptic restricted three body problem assuming bigger primary as a source of radiation and the smaller one a triaxial rigid body. The perturbation technique developed by Bennet (Icarus 4:177, 1965b) has been used for determination of characteristic exponents. This technique is based on Floquet’s Theory for determination of characteristic exponents in the system with periodic coefficients. The results of the study are analytical and numerical expressions are simulated for the transition curves bounding the region of stability in the μ–e plane, accurate to O(e ²). The unstable region is found to be divided into three parts. The effect of radiation parameter is significant. For small values of e, the results are in favor with the numerical analysis of Danby (Astron. J. 69:166, 1964), Bennet (Icarus 4:177, 1965b), Alfriend and Rand (AIAA J. 6:1024, 1969). The effect of radiation pressure is significant than the oblateness and triaxiality of the primaries.     

Prominence Cavity Regions Observed Using SWAP 174 Å Filtergrams and Simultaneous Eclipse Flash Spectra

SWAP images from PROBA2 taken at 174 Å in the Fe ix/x lines are compared with simultaneous slitless flash spectra obtained during the solar total eclipse of 11 July 2010. Myriad faint low-excitation emission lines together with the He i and He ii Paschen α chromospheric lines are recorded on eclipse spectra where regions of limb prominences are obtained with space-borne imagers. We analyzed a deep flash spectrum obtained by summing 80 individual spectra to evaluate the intensity modulations of the continuum. Intensity deficits are observed and measured at the prominences boundaries in both eclipse and SWAP images. The prominence cavities interpreted as a relative depression of plasma density, produced inside the corona surrounding the prominences, and some intense heating occurring in these regions, are discussed. Photometric measurements are shown at different scales and different, spectrally narrow, intervals for both the prominences and the coronal background.     

Seismology of a Sunspot Atmosphere

Two competing theories of sunspot oscillations are discussed. It is pointed out that the normal mode (eigenoscillations) theory is in contradiction with a number of observations. The reasons for this are discussed. The revised filter theory of the three-minute sunspot oscillations is outlined. It is shown that the reason for the occurrence of the multipassband filter for the slow waves is the interference that appears from the multilayer structure of the sunspot atmosphere. In contrast with Zhugzhda and Locans (Sov. Astron. Lett. 7, 25?–?27, 1981) it is shown that along with the Fabry?–?Perot chromospheric passband the cutoff frequency passband and a number of the high-frequency passbands occur. The effect of the nonlinearity of the sunspot oscillations in the upper chromosphere and the transition region is taken into account. The spectra of the distinct empirical models of the sunspot atmosphere are explored. An example of the interpretation of the sunspot oscillations based on the revised filter theory is presented. Only the filter theory can explain the complicated behavior of the oscillations across the sunspot. The observations provide evidence of the nonuniformity of the sunspot atmosphere.     

Long-Term Variations of Solar Differential Rotation and Sunspot Activity: Revisited

Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential-rotation law obtained by Makarov, Tlatov, and Callebaut (Solar Phys. 170, 373, 1997), Javaraiah, Bertello, and Ulrich (Astrophys. J. 626, 579, 2005a; Solar Phys. 232, 25, 2005b), and Javaraiah et al. (Solar Phys. 257, 61, 2009). Our results indicate that the solar-surface-rotation rate at the Equator (indicated by the A-parameter of the standard solar-rotation law) shows a secular decrease since Cycle 12 onwards, given by about 1?–?1.5×10^?3 (deg?day^?1?year^?1). The B-parameter of the standard differential-rotation law seems to also show a secular decrease since Cycle 12 onwards, but of weak statistical significance. The rotation rate averaged over latitudes 0^°?–?40^° does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since Cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation on long-term scales.