High cadence observations of solar activity (active regions, flares, filaments) in the H$\alpha$ line were performed at Meudon and Haute Provence Observatories from 1956 to 2004. More than 7 million images were recorded, mainly on 35 mm films. After a review of the scientific interest of solar surveys at high temporal resolution and the historical background, we describe the new instrument which will operate automatically in 2020 at the Calern station of the C\^{o}te d'Azur observatory (1270 m). It will replace the former heliographs with improved cadence, seeing and time coverage. We summarize the capabilities of the optical design and present new scientific perspectives in terms of flare onset and Moreton wave detection.

The recent claimed detection of g modes was obtained from the analysis of the power spectrum of the time series of round-trip travel time of p modes. The goal of this paper is to reproduce these results on which the claims are based for confirming or invalidating the detection of g modes with the method used to make the claims. We computed the time series of round-trip travel time using the procedure given in Fossat et al. (2017), and used different variations of the times series for comparison. We used the recently calibrated GOLF data (published in Paper I) with different sampling, different photomultipliers, different length of data for reproducing the analysis. We also correlated the power spectrum with an asymptotic model of g-mode frequencies in a similar manner to Fossat and Schmider (2018). We devised a scheme for optimising the correlation both for pure noise and for the GOLF data. We confirm the analysis performed in Fossat et al. (2017) but draw different conclusions. Their claims of detection of g modes cannot be confirmed when changing parameters such as sampling interval, length of time series, or photomultipliers. Other instrument such as GONG and BiSON do not confirm their detection. We also confirm the analysis performed in Fossat and Schmider (2018), but again draw different conclusions. For GOLF, the correlation of the power spectrum with the asymptotic model of g-mode frequencies for $l=1$ and $l=2$ show a high correlation at lag=0 and at lag corresponding to the rotational splitting $\nu_l$, but the same occurs for pure noise due to the large number of peaks present in the model. In addition, other very different parameters defining the asymptotic model also provide a high correlation at these lags. We conclude that the detection performed in Fossat and Schmider (2018) is an artefact of the methodology.

In this work we study in details the influence of pure astronomical refraction on solar metrologic measurements made from ground-based full disk imagery and provide the tools for correcting the measurements and estimating the associated uncertainties.For a given standard atmospheric model, we first use both analytical and numerical methods in order to test the validity of the commonly or historically used approximations of the differential effect of refraction as a function of zenith distance. For a given refraction model, we provide the exact formulae for correcting solar radius measurements at any heliographic angle and for any zenith distance. Then, using solar images recorded in the near infrared between 2011 and 2016, we show that these corrections can be applied up to $70\deg$ using the usual approximate formulae and can be extended up to $80\deg$ of zenith distance provided that a standard atmospheric model and a full numerical integration of the refraction integral are used. We also provide estimates of the absolute uncertainties associated with the differential refraction corrections and show that approximate formulae can be used up to $80\deg$ of zenith distance for computing these uncertainties. For a given instrumental setup and the knowledge of the uncertainties associated with local weather records, this can be used to fix the maximum zenith distance one can observe depending on the required astrometric accuracy.

The recent claims of g-mode detection have restarted the search for these potentially extremely important modes. These claims can be reassessed in view of the different data sets available from the SoHO instruments and ground-based instruments. We produce a new calibration of the GOLF data with a more consistent p-mode amplitude and a more consistent time shift correction compared to the time series used in the past. The calibration of 22 years of GOLF data is done with a simpler approach that uses only the predictive radial velocity of the SoHO spacecraft as a reference. Using p modes, we measure and correct the time shift between ground- and space-based instruments and the GOLF instrument. The p-mode velocity calibration is now consistent to within a few percent with other instruments. The remaining time shifts are within $\pm$ 5 s for 99.8% of the time series. Accepted by A&A on June 29, 2018

PICARD is a scientific space mission dedicated to the study of the solar variability origin. A French micro-satellite will carry an imaging telescope for measuring the solar diameter, limb shape and solar oscillations, and two radiometers for measuring the total solar irradiance and the irradiance in five spectral domains, from ultraviolet to infrared. The mission is planed to be launched in 2009 for a 3-year duration. This article presents the PICARD Payload Data Centre, which role is to collect, process and distribute the PICARD data. The Payload Data Centre is a joint project between laboratories, space agency and industries. The Belgian scientific policy office funds the industrial development and future operations under the European Space Agency program. The development is achieved by the SPACEBEL Company. The Belgian operation centre is in charge of operating the PICARD Payload Data Centre. The French space agency leads the development in partnership with the French scientific research centre, which is responsible for providing all the scientific algorithms. The architecture of the PICARD Payload Data Centre (software and hardware) is presented. The software system is based on a Service Oriented Architecture. The host structure is made up of the basic functions such as data management, task scheduling and system supervision including a graphical interface used by the operator to interact with the system. The other functions are mission-specific: data exchange (acquisition, distribution), data processing (scientific and non-scientific processing) and managing the payload (programming, monitoring). The PICARD Payload Data Centre is planned to be operated for 5 years. All the data will be stored into a specific data centre after this period.

We present the identification of very low frequency g modes in the asymptotic regime and two important parameters that have long been waited for: the core rotation rate, and the asymptotic equidistant period spacing of these g modes. The GOLF instrument on board the SOHO space observatory has provided two decades of full-disk helioseismic data. In the present study, we search for possible collective frequency modulations that are produced by periodic changes in the deep solar structure. Such modulations provide access to only very low frequency g modes, thus allowing statistical methods to take advantage of their asymptotic properties. For oscillatory periods in the range between 9 and nearly 48 hours, almost 100 g modes of spherical harmonic degree 1 and more than 100 g modes of degree 2 are predicted. They are not observed individually, but when combined, they unambiguouslyprovide their asymptotic period equidistance and rotational splittings, in excellent agreement with the requirements of the asymptotic approximations. Previously, p-mode helioseismology allowed the g-mode period equidistance parameter $P_0$ to be bracketed inside a narrow range, between approximately 34 and 35 minutes. Here, $P_0$ is measured to be 34 min 01 s, with a 1 s uncertainty. The previously unknown g-mode splittings have now been measured from a non-synodic reference with very high accuracy, and they imply a mean weighted rotation of 1277 $\pm$ 10 nHz (9-day period) of their kernels, resulting in a rapid rotation frequency of 1644 $\pm$ 23 nHz (period of one week) of the solar core itself, which is a factor 3.8 $\pm$ 0.1 faster than the rotation of the radiative envelope. The g modes are known to be the keys to a better understanding of the structure and dynamics of the solar core. Their detection with these precise parameters will certainly stimulate a new era of research in this field.

While magnetic field has only second order effects on global sound waves, it can affect the mode frequencies and flow fields inferred using local helioseismology. A strong localized field can bias our results when one try to infer global scale flows such as the meridional circulation. On the other hand, one important object of research in this area is to determine whether one can detect changes in frequencies or sub-surface flow fields that could be a precursor of surface magnetic events. In this paper we review recent development in this area focusing on results obtained with the ring-diagram analysis technique.

Developed at the Observatoire de la C\^ote d'Azur (OCA) within the framework of the PICARD space mission (Thuillier et al., 2006) and with support from the french spatial agency (CNES), MISOLFA (Moniteur d'Images Solaires Franco-Alg\'erien) is a new generation of daytime turbulence monitor. Its objective is to measure both the spatial and temporal turbulence parameters in order to quantify their effects on the solar diameter measurements that will be made from ground using the qualification model of the SODISM (SOlar Diameter Imager and Surface Mapper) instrument onboard PICARD. The comparison of simultaneous images from ground and space should allow us, with the help of the solar monitor, to find the best procedure possible to measure solar diameter variations from ground on the long term. MISOLFA is now installed at the Calern facility of OCA and PICARD is scheduled to be launched in 2010. We present here the principles of the instrument and the first results obtained on the characteristics of the turbulence observed at Calern observatory using this monitor while waiting for the launch of the space mission.

Continuum intensity images from PICARD/SODISM and SDO/HMI covering a 209 days period in 2011 are analyzed in order to extract mode parameters for spherical harmonics up to $l=100$. SODISM helioseismology signal is affected by the low orbit of PICARD and by important gaps and CCD persistence effects. SODISM intensity signal has a lower signal to noise ratio and duty cycle than HMI and less modes were successfully fitted over the same period. A comparison of the rotation profiles obtained from both sets of continuum images shows however that the results remain compatible within one standard deviation of HMI formal errors.

The Solar Diameter Imager and Surface Mapper (SODISM) on board the PICARD space mission provides wide-field images of the photosphere and chromosphere of the Sun in five narrow pass bands (centered at 215.0, 393.37, 535.7, 607.1, and 782.2 nm). PICARD is a space mission, which was successfully launched on 15 June 2010 into a Sun synchronous dawn-dusk orbit. It represents a European asset aiming at collecting solar observations that can serve to estimate some of the inputs to Earth climate models. The scientific payload consists of the SODISM imager and of two radiometers, SOVAP (SOlar VAriability PICARD) and PREMOS (PREcision MOnitor Sensor), which carry out measurements that allow estimating the Total Solar Irradiance (TSI) and the Solar Spectral Irradiance (SSI) from the middle ultraviolet to the red. The SODISM telescope monitors solar activity continuously. It thus produces images that can also feed SSI reconstruction models. Further, the objectives of SODISM encompass the probing of the interior of the Sun via helioseismic analysis of observations in intensity (on the solar disc and at the limb), and via astrometric investigations at the limb. The latter addresses especially the spectral dependence of the radial limb shape, and the temporal evolution of the solar diameter and asphericity. After a brief review of its original science objectives, this paper presents the detailed design of the SODISM instrument, its expected performance, and the scheme of its flight operations. Some observations with SODISM are presented and discussed.

PICARD is a CNES micro-satellite launched in June 2010 (Thuillier at al. 2006). Its main goal is to measure the solar shape, total and spectral irradiance during the ascending phase of the activity cycle. The SODISM telescope onboard PICARD also allows us to conduct a program for helioseismology in intensity at 535.7 nm (Corbard et al. 2008). One-minute cadence low-resolution full images are available for a so-called medium-$l$ program, and high-resolution images of the limb recorded every 2 minutes are used to study mode amplification near the limb in the perspective of g-mode search. First analyses and results from these two programs are presented here.

Context. We compare solar differential rotation of subphotospheric layers derived from local helioseismology analysis of GONG++ dopplergrams and the one derived from tracing small bright coronal structures (SBCS) using EIT/SOHO images for the period August 2001 - December 2006, which correspond to the declining phase of solar cycle 23. Aims. The study aims to find a relationship between the rotation of the SBCS and the subphotospheric angular velocity. The northsouth asymmetries of both rotation velocity measurements are also investigated. Methods. Subphotospheric differential rotation was derived using ring-diagram analysis of GONG++ full-disk dopplergrams of 1 min cadence. The coronal rotation was derived by using an automatic method to identify and track the small bright coronal structures in EIT full-disk images of 6 hours cadence. Results. We find that the SBCS rotate faster than the considered upper subphotospheric layer (3Mm) by about 0.5 deg/day at the equator. This result joins the results of several other magnetic features (sunspots, plages, faculae, etc.) with a higher rotation than the solar plasma. The rotation rate latitudinal gradients of the SBCS and the subphotospheric layers are very similar. The SBCS motion shows an acceleration of about 0.005 deg/day/month during the declining phase of solar cycle 23, whereas the angular velocity of subsurface layers does not display any evident variation with time, except for the well known torsional oscillation pattern. Finally, both subphotospheric and coronal rotations of the southern hemisphere are predominantly larger than those of the northern hemisphere. At latitudes where the north-south asymmetry of the angular velocity increases (decreases) with activity for the SBCS, it decreases (increases) for subphotospheric layers.

The POLAR Investigation of the Sun (POLARIS) mission uses a combination of a gravity assist and solar sail propulsion to place a spacecraft in a 0.48 AU circular orbit around the Sun with an inclination of 75 degrees with respect to solar equator. This challenging orbit is made possible by the challenging development of solar sail propulsion. This first extended view of the high-latitude regions of the Sun will enable crucial observations not possible from the ecliptic viewpoint or from Solar Orbiter. While Solar Orbiter would give the first glimpse of the high latitude magnetic field and flows to probe the solar dynamo, it does not have sufficient viewing of the polar regions to achieve POLARIS' primary objective : determining the relation between the magnetism and dynamics of the Sun's polar regions and the solar cycle.

Context. The most recent determination of the solar chemical composition, using a time-dependent, 3D hydrodynamical model of the solar atmosphere, exhibits a significant decrease of C, N, O abundances compared to their previous values. Solar models that use these new abundances are not consistent with helioseismological determinations of the sound speed profile, the surface helium abundance and the convection zone depth. Aims. We investigate the effect of changes of solar abundances on low degree p-mode and g-mode characteristics which are strong constraints of the solar core. We consider particularly the increase of neon abundance in the new solar mixture in order to reduce the discrepancy between models using new abundances and helioseismology. Methods. The observational determinations of solar frequencies from the GOLF instrument are used to test solar models computed with different chemical compositions. We consider in particular the normalized small frequency spacings in the low degree p-mode frequency range. Results. Low-degree small frequency spacings are very sensitive to changes in the heavy-element abundances, notably neon. We show that by considering all the seismic constraints, including the small frequency spacings, a rather large increase of neon abundance by about (0.5 +/- 0.05)dex can be a good solution to the discrepancy between solar models that use new abundances and low degree helioseismology, subject to adjusting slightly the solar age and the highest abundances. We also show that the change in solar abundances, notably neon, considerably affects g-mode frequencies, with relative frequency differences between the old and the new models higher than 1.5%

The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the Sun-Earth relationship. Thus obtaining a complete MHD solar picture is a clear objective for the next decades and it requires complementary observations of the dynamics of the radiative zone. For that ambitious goal, space prototypes are being developed to improve gravity mode detection. The Sun is unique to progress on the topology of deep internal magnetic fields and to understand the complex mechanisms which provoke photospheric and coronal magnetic changes and possible longer cycles important for human life. We propose the following roadmap in Europe to contribute to this "impressive" revolution in Astronomy and in our Sun-Earth relationship: SoHO (1995-2007), PICARD (2008-2010), DynaMICS (2009-2017) in parallel to SDO (2008-2017) then a world-class mission located at the L1 orbit or above the solar poles.

Of all the fundamental parameters of the Sun (diameter, mass, temperature...), the gravitational multipole moments (of degree l and order m) that determine the solar moments of inertia, are still poorly known. However, at the first order (l=2), the quadrupole moment is relevant to many astrophysical applications. It indeed contributes to the relativistic perihelion advance of planets, together with the post-Newtonian (PN) parameters; or to the precession of the orbital plane about the Sun polar axis, the latter being unaffected by the purely relativistic PN contribution. Hence, a precise knowledge of the quadrupole moment is necessary for accurate orbit determination, and alternatively, to obtain constraints on the PN parameters. Moreover, the successive gravitational multipole moments have a physical meaning: they describe deviations from a purely spherical mass distribution. Thus, their precise determination gives indications on the solar internal structure. Here, we explain why it is difficult to compute these parameters, how to derive the best values, and how they will be determined in a near future by means of space experiments.

From the analysis of low-order GOLF+MDI sectoral modes and LOWL data (l > 3), we derive the solar radial rotation profile assuming no latitudinal dependance in the solar core. These low-order acoustic modes contain the most statistically significant information about rotation of the deepest solar layers and should be least influenced by internal variability associated with the solar dynamo. After correction of the sectoral splittings for their contamination by the rotation of the higher latitudes, we obtain a flat rotation profile down to 0.2 solar radius.

We report quantitative analysis of the radial gradient of solar angular velocity at depths down to about 15 Mm below the solar surface for latitudes up to 75 degree using the Michelson Doppler Imager (MDI) observations of surface gravity waves (f modes) from the Solar and Heliospheric Observatory (SoHO). A negative outward gradient of around -400 nHz/R, equivalent to logarithmic gradient of the rotation frequency with respect to radius which is very close to -1, is found to be remarkably constant between the equator and 30 degree of latitude. Above 30 degree it decreases in absolute magnitude to a very small value at around 50 degree. At higher latitudes the gradient may reverse its sign: if so this reversal takes place in a thin layer extending only 5 Mm beneath the visible surface, as evidenced by the most superficial modes (with degrees l>250). The signature of the torsional oscillations is seen in this layer, but no other significant temporal variations of the gradient and value of the rotation rate there are found.

High quality observations of the low-degree acoustic modes (p-modes) exist for almost two complete solar cycles using the solar spectrophotometer Mark-I, located at the Observatorio del Teide (Tenerife, Spain) and operating now as part of the Birmingham Solar Oscillations Network (BiSON). We have performed a Fourier analysis of 30 calibrated time-series of one year duration covering a total period of 15 years between 1984 and 1999. Applying different techniques to the resulting power spectra, we study the signature of the solar activity changes on the low-degree p-modes. We show that the variation of the central frequencies and the total velocity power (TVP) changes. A new method of simultaneous fit is developed and a special effort has been made to study the frequency-dependence of the frequency shift. The results confirm a variation of the central frequencies of acoustic modes of about 450 nHz, peak-to-peak, on average for low degree modes between 2.5 and 3.7 mHz. The TVP is anti-correlated with the common activity indices with a decrease of about 20% between the minimum and the maximum of solar cycle 22. The results are compared with those obtained for intermediate degrees, using the LOWL data. The frequency shift is found to increase with the degree with a weak l-dependence similar to that of the inverse mode mass. This verifies earlier suggestions that near surface effects are predominant.

High quality observations of the low degree p-modes exist for almost two complete solar cycles using the solar spectrophotometer Mark-I, located and operating at the Observatorio del Teide (Tenerife, Spain). In this work, the observations available have been re-analyzed over a much wider time interval than before. We analyze the time variation of the yearly frequency shift and its frequency dependence. This information will be used in order to average annual power spectra by removing the effect of the solar cycle. Using this average power spectrum, a new estimate of the rotational splittings is attempted.

We present an analysis of the frequency shift and the even terms of the frequency splitting coefficients carried out using six years of LOWL data, starting in 1994. The temporal variations, and their dependences with the frequency and degree are addressed. The results are consistent with previous analysis.

The solar tachocline, located at the interface between the latitude-dependent rotation of the convection zone and the rigid radiative interior, presents high gradients of angular velocity which are of particular interest for the models of the solar dynamo and angular momentum transport. Furthermore, latitudinal and temporal variations of the tachocline parameters, if any, are also of particular interest in order to constrain models. We present a review of some of the theories of the tachocline and their predictions that may be tested by helioseismology. We describethe methods for inferring the tachocline parameters from observations and the associated difficulties. A review of results previously obtained is given and an analysis of the new 6 years database of LOWL observations is presented which yields no compelling evidence of variations or general trend of the tachocline parameters during the ascending phase of the current solar cycle (1994-2000).

Recent studies have established that peaks in solar oscillation power spectra are not Lorentzian in shape, but have a distinct asymmetry. Fitting a symmetric Lorentzian profile to the peaks therefore produces a shift in frequency of the modes. Accurate determination of low-frequency modes is essential to infer the structure of the solar core by inversion of the mode frequencies. In this paper we investigate how the changes in frequencies of low-degree modes obtained by fitting symmetric and asymmetric peak profiles change the inferred properties of the solar core. We use data obtained by the Global Oscillations at Low Frequencies (GOLF) project on board the SoHO spacecraft. Two different solar models and inversion procedures are used to invert the data to determine the sound speed in the solar core. We find that for a given set of modes no significant difference in the inferred sound-speed results from taking asymmetry into account when fitting the low-degree modes.

Inversions of rotational splittings have shown that there exists at the base of the solar convection zone a region called the tachocline in which high radial gradients of the rotation rate occur. The usual linear regularization methods tend to smooth out any high gradients in the solution, and may not be appropriate for the study of this zone. In this paper we use, in the helioseismic context of rotation inversions, regularization methods that have been developed for edge-preserving regularization in computed imaging. It is shown from Monte-Carlo simulations that this approach can lead directly to results similar to those reached by linear inversions which however required some assumptions on the shape of the transition in order to be deconvolved. The application of this method to LOWL data leads to a very thin tachocline. From the discussions on the parameters entering the inversion and the Monte-Carlo simulations, our conclusion is that the tachocline width is very likely below 0.05R_sun which lowers our previous estimate of 0.05+/- 0.03R_sun obtained from the same dataset (Corbard et al. 1998).

Inversions of rotational splittings have shown that the surface layers and the so-called solar tachocline at the base of the convection zone are regions in which high radial gradients of the rotation rate occur. The usual regularization methods tend to smooth out every high gradients in the solution and may not be appropriate for the study of a zone like the tachocline. In this paper we use nonlinear regularization methods that are developed for edge-preserving regularization in computed imaging (e.g. Blanc-F\'eraud et al. 1995) and we apply them in the helioseismic context of rotational inversions.

The low degree splittings obtained from one year of GOLF data analysis are combined with the MDI medium-l 144-day splittings in order to infer the solar internal rotation as a function of the radius down to $0.2R_\odot$. Several inverse methods are applied to the same data and the uncertainties on the solution as well as the resolution reachable are discussed. The results are compared with the one obtained from the low degree splittings estimated from GONG network.

Helioseismic inversions, carried out for several years on various ground-based and spatial observations, have shown that the solar rotation rate presents two principal regimes: a quasi-rigid rotation in the radiative interior and a latitude-dependent rotation in the whole convection zone. The thin layer, named solar tachocline, between these two regimes is difficult to infer through inverse techniques because of the ill-posed nature of the problem that requires regularization techniques which, in their global form, tend to smooth out any high gradient in the solution. Thus, most of the previous attempts to study the rotation profile of the solar tachocline have been carried out through forward modeling. In this work we show that some appropriate inverse techniques can also be used and we compare the ability of three 1D inverse techniques combined with two automatic strategies for the choice of the regularization parameter, to infer the solar tachocline profile in the equatorial plane. Our work, applied on LOWL (LOWL is an abbreviation for low degree denoted by L) two years dataset, argue in favor of a very sharp (0.05+/-0.03R_sun) transition zone located at 0.695+/-0.005R_sun which is in good agreement with the previous forward analysis carried out on Global Oscillations Network Group (GONG), Big Bear Solar Observatory (BBSO) and LOWL datasets.