ALBERT

Description: Aims. This paper describes the Polarimetric and Helioseismic Imager on the Solar Orbiter mission (SO/PHI), the first magnetograph and helioseismology instrument to observe the Sun from outside the Sun-Earth line. It is the key instrument meant to address the top-level science question: How does the solar dynamo work and drive connections between the Sun and the heliosphere? SO/PHI will also play an important role in answering the other top-level science questions of Solar Orbiter, while hosting the potential of a rich return in further science. Methods. SO/PHI measures the Zeeman effect and the Doppler shift in the Fe I 617.3 nm spectral line. To this end, the instrument carries out narrow-band imaging spectro-polarimetry using a tunable LiNbO3 Fabry-Perot etalon, while the polarisation modulation is done with liquid crystal variable retarders. The line and the nearby continuum are sampled at six wavelength points and the data are recorded by a 2k × 2k CMOS detector. To save valuable telemetry, the raw data are reduced on board, including being inverted under the assumption of a Milne-Eddington atmosphere, although simpler reduction methods are also available on board. SO/PHI is composed of two telescopes; one, the Full Disc Telescope, covers the full solar disc at all phases of the orbit, while the other, the High Resolution Telescope, can resolve structures as small as 200 km on the Sun at closest perihelion. The high heat load generated through proximity to the Sun is greatly reduced by the multilayer-coated entrance windows to the two telescopes that allow less than 4% of the total sunlight to enter the instrument, most of it in a narrow wavelength band around the chosen spectral line. Results. SO/PHI was designed and built by a consortium having partners in Germany, Spain, and France. The flight model was delivered to Airbus Defence and Space, Stevenage, and successfully integrated into the Solar Orbiter spacecraft. A number of innovations were introduced compared with earlier space-based spectropolarimeters, thus allowing SO/PHI to fit into the tight mass, volume, power and telemetry budgets provided by the Solar Orbiter spacecraft and to meet the (e.g. thermal) challenges posed by the mission’s highly elliptical orbit.

Description: Context. Many moons have been detected around planets in our Solar System, but none has been detected unambiguously around any of the confirmed extrasolar planets. Aims. We test the feasibility of a supervised convolutional neural network to classify photometric transit light curves of planet-host stars and identify exomoon transits, while avoiding false positives caused by stellar variability or instrumental noise. Methods. Convolutional neural networks are known to have contributed to improving the accuracy of classification tasks. The network optimization is typically performed without studying the effect of noise on the training process. Here we design and optimize a 1D convolutional neural network to classify photometric transit light curves. We regularize the network by the total variation loss in order to remove unwanted variations in the data features. Results. Using numerical experiments, we demonstrate the benefits of our network, which produces results comparable to or better than the standard network solutions. Most importantly, our network clearly outperforms a classical method used in exoplanet science to identify moon-like signals. Thus the proposed network is a promising approach for analyzing real transit light curves in the future.

Description: Context. The spatial power spectrum of supergranulation does not fully characterize the underlying physics of turbulent convection. For example, it does not describe the non-Gaussianity in the horizontal flow divergence. Aims. Our aim is to statistically characterize the spatial pattern of solar supergranulation beyond the power spectrum. The next-order statistic is the bispectrum. It measures correlations of three Fourier components and is related to the nonlinearities in the underlying physics. It also characterizes how a skewness in the dataset is generated by the coupling of three Fourier components. Methods. We estimated the bispectrum of supergranular horizontal surface divergence maps that were obtained using local correlation tracking (LCT) and time-distance helioseismology (TD) from one year of data from the helioseismic and magnetic imager on-board the solar dynamics observatory starting in May 2010. Results. We find significantly nonzero and consistent estimates for the bispectrum using LCT and TD. The strongest nonlinearity is present when the three coupling wave vectors are at the supergranular scale. These are the same wave vectors that are present in regular hexagons, which have been used in analytical studies of solar convection. At these Fourier components, the bispectrum is positive, consistent with the positive skewness in the data and consistent with supergranules preferentially consisting of outflows surrounded by a network of inflows. We use the bispectral estimates to generate synthetic divergence maps that are very similar to the data. This is done by a model that consists of a Gaussian term and a weaker quadratic nonlinear component. Using this method, we estimate the fraction of the variance in the divergence maps from the nonlinear component to be of the order of 4–6%. Conclusions. We propose that bispectral analysis is useful for understanding the dynamics of solar turbulent convection, for example for comparing observations and numerical models of supergranular flows. This analysis may also be useful to generate synthetic flow fields.

Description: Context. The tilt of solar active regions described by Joy’s law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes what we observe as Joy’s law, acting on east-west flows as they rise towards the surface. Aims. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we can constrain the origins of Joy’s law. Methods. We measured the tilt angle of the primary bipoles in 153 emerging active regions (EARs) in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey. We used line-of-sight magnetic field measurements averaged over 6 h to define the polarities and measure the tilt angle up to four days after emergence. Results. We find that at the time of emergence the polarities are on average aligned east-west, and that neither the separation nor the tilt depends on latitude. We do find, however, that EARs at higher latitudes have a faster north-south separation speed than those closer to the equator at the emergence time. After emergence, the tilt angle increases and Joy’s law is evident about two days later. The scatter in the tilt angle is independent of flux until about one day after emergence, when we find that higher-flux regions have a smaller scatter in tilt angle than lower-flux regions. Conclusions. Our finding that active regions emerge with an east-west alignment is consistent with earlier observations, but is still surprising since thin flux tube models predict that tilt angles of rising flux tubes are generated below the surface. Previously reported tilt angle relaxation of deeply anchored flux tubes can be largely explained by the change in east-west separation. We conclude that Joy’s law is caused by an inherent north-south separation speed present when the flux first reaches the surface, and that the scatter in the tilt angle is consistent with buffeting of the polarities by supergranulation.

Description: Context. Solar Rossby waves have only recently been unambiguously identified in Helioseimsic and Magnetic Imager (HMI) and Michelson Doppler Imager maps of flows near the solar surface. So far this has not been done with the Global Oscillation Network Group (GONG) ground-based observations, which have different noise properties. Aims. We use 17 years of GONG++ data to identify and characterize solar Rossby waves using ring-diagram helioseismology. We compare directly with HMI ring-diagram analysis. Methods. Maps of the radial vorticity were obtained for flows within the top 2 Mm of the surface for 17 years of GONG++ data. The data were corrected for systematic effects including the annual periodicity related to the B0 angle. We then computed the Fourier components of the radial vorticity of the flows in the co-rotating frame. We performed the same analysis on the HMI data that overlap in time. Results. We find that the solar Rossby waves have measurable amplitudes in the GONG++ sectoral power spectra for azimuthal orders between m = 3 and m = 15. The measured mode characteristics (frequencies, lifetimes, and amplitudes) from GONG++ are consistent with the HMI measurements in the overlap period from 2010 to 2018 for m ≤ 9. For higher-m modes the amplitudes and frequencies agree within two sigmas. The signal-to-noise ratio of modes in GONG++ power spectra is comparable to those of HMI for 8 ≤ m ≤ 11, but is lower by a factor of two for other modes. Conclusions. The GONG++ data provide a long and uniform data set that can be used to study solar global-scale Rossby waves from 2001.

Description: Context. Global-scale equatorial Rossby waves have recently been unambiguously identified on the Sun. Like solar acoustic modes, Rossby waves are probes of the solar interior. Aims. We study the latitude and depth dependence of the Rossby wave eigenfunctions. Methods. By applying helioseismic ring-diagram analysis and granulation tracking to observations by HMI aboard SDO, we computed maps of the radial vorticity of flows in the upper solar convection zone (down to depths of more than 16 Mm). The horizontal sampling of the ring-diagram maps is approximately 90 Mm (∼7.5°) and the temporal sampling is roughly 27 hr. We used a Fourier transform in longitude to separate the different azimuthal orders m in the range 3 ≤ m ≤ 15. At each m we obtained the phase and amplitude of the Rossby waves as functions of depth using the helioseismic data. At each m we also measured the latitude dependence of the eigenfunctions by calculating the covariance between the equator and other latitudes. Results. We conducted a study of the horizontal and radial dependences of the radial vorticity eigenfunctions. The horizontal eigenfunctions are complex. As observed previously, the real part peaks at the equator and switches sign near ±30°, thus the eigenfunctions show significant non-sectoral contributions. The imaginary part is smaller than the real part. The phase of the radial eigenfunctions varies by only ±5° over the top 15 Mm. The amplitude of the radial eigenfunctions decreases by about 10% from the surface down to 8 Mm (the region in which ring-diagram analysis is most reliable, as seen by comparing with the rotation rate measured by global-mode seismology). Conclusions. The radial dependence of the radial vorticity eigenfunctions deduced from ring-diagram analysis is consistent with a power law down to 8 Mm and is unreliable at larger depths. However, the observations provide only weak constraints on the power-law exponents. For the real part, the latitude dependence of the eigenfunctions is consistent with previous work (using granulation tracking). The imaginary part is smaller than the real part but significantly nonzero.

Description: Context. There is a wide discrepancy in current estimates of the strength of convection flows in the solar interior obtained using different helioseismic methods applied to observations from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. The cause for these disparities is not known. Aims. As one step in the effort to resolve this discrepancy, we aim to characterize the multi-ridge fitting code for ring-diagram helioseismic analysis that is used to obtain flow estimates from local power spectra of solar oscillations. Methods. We updated the multi-ridge fitting code developed by Greer et al. (2014, Sol. Phys., 289, 2823) to solve several problems we identified through our inspection of the code. In particular, we changed the (1) merit function to account for the smoothing of the power spectra, (2) model for the power spectrum, and (3) noise estimates. We used Monte Carlo simulations to generate synthetic data and to characterize the noise and bias of the updated code by fitting these synthetic data. Results. The bias in the output fit parameters, apart from the parameter describing the amplitude of the p-mode resonances in the power spectrum, is below what can be measured from the Monte-Carlo simulations. The amplitude parameters are underestimated; this is a consequence of choosing to fit the logarithm of the averaged power. We defer fixing this problem as it is well understood and not significant for measuring flows in the solar interior. The scatter in the fit parameters from the Monte-Carlo simulations is well-modeled by the formal error estimates from the code. Conclusions. We document and demonstrate a reliable multi-ridge fitting method for ring-diagram analysis. The differences between the updated fitting results and the original results are less than one order of magnitude and therefore we suspect that the changes will not eliminate the aforementioned orders-of-magnitude discrepancy in the amplitude of convective flows in the solar interior.

Description: Context. While the Solar System contains about 20 times more moons than planets, no moon has been confirmed around any of the thousands of extrasolar planets discovered so far. Considering the large computational load required for the statistical vetting of exomoon candidates in a star–planet–moon framework, tools for an uncomplicated identification of the most promising exomoon candidates could be beneficial to streamline follow-up studies. Aims. Here we study three exomoon indicators that emerge if well-established planet-only models are fitted to a planet–moon transit light curve: transit timing variations (TTVs), transit duration variations (TDVs), and apparent planetary transit radius variations (TRVs). We re-evaluate under realistic conditions the previously proposed exomoon signatures in the TTV and TDV series. Methods. We simulated light curves of a transiting exoplanet with a single moon, taking into account stellar limb darkening, orbital inclinations, planet–moon occultations, and noise from both stellar granulation and instrumental effects. These model light curves were then fitted with a planet-only transit model whilst pretending there were no moon, and we explored the resulting TTV, TDV, and TRV series for evidence of the moon. Results. The previously described ellipse in the TTV-TDV diagram of an exoplanet with a moon emerges only for high-density moons. However, low-density moons distort the sinusoidal shapes of the TTV and the TDV series due to their photometric contribution to the combined planet–moon transit. Sufficiently large moons can nevertheless produce periodic apparent TRVs of their host planets that could be observable. We find that Kepler and PLATO have similar performances in detecting the exomoon-induced TRV effect around simulated bright (mV = 8) stars. Although these stars are rare in the Kepler sample, they will be abundant in the PLATO sample. Moreover, PLATO’s higher cadence yields a stronger TTV signal. We detect substantial TRVs of the Saturn-sized planet Kepler-856 b although an exomoon could only ensure Hill stability in a very narrow orbital range. Conclusions. The periodogram of the sequence of transit radius measurements can indicate the presence of a moon. The TTV and TDV series of exoplanets with moons could be more complex than previously assumed. We propose that TRVs could be a more promising means to identify exomoons in large exoplanet surveys.