We review recent advances and results in enhancing and developing helioseismic analysis methods and in solar data assimilation. In the first part of this paper we will focus on selected developments in time-distance and global helioseismology. In the second part, we review the application of data assimilation methods on solar data. Relating solar surface observations as well as helioseismic proxies with solar dynamo models by means of the techniques from data assimilation is a promising new approach to explore and to predict the magnetic activity cycle of the Sun.

At intermediate and high degree $l$, solar p- and f modes can be considered as surface waves. Using variational principle, we derive an integral expression for the group velocities of the surface waves in terms of adiabatic eigenfunctions of normal modes, and address the benefits of using group-velocity measurements as a supplementary diagnostic tool in solar seismology. The principal advantage of using group velocities, when compared with direct analysis of the oscillation frequencies, comes from their smaller sensitivity to the uncertainties in the near-photospheric layers. We address some numerical examples where group velocities are used to reveal inconsistencies between the solar models and the seismic data. Further, we implement the group-velocity measurements to the calibration of the specific entropy, helium abundance Y and heavy-element abundance Z in the adiabatically-stratified part of the solar convective envelope, using different recent versions of the equation of state. The results are in close agreement with our earlier measurements based on more sophisticated analysis of the solar oscillation frequencies (Vorontsov et al. 2013, MNRAS 430, 1636). These results bring further support to the downward revision of the solar heavy-element abundances in recent spectroscopic measurements.

We describe a global parametric model for the observed power spectra of solar oscillations of intermediate and low degree. A physically motivated parameterization is used as a substitute for a direct description of mode excitation and damping as these mechanisms remain poorly understood. The model is targeted at the accurate fitting of power spectra coming from Doppler velocity measurements and uses an adaptive response function that accounts for both the vertical and horizontal components of the velocity field on the solar surface and for possible instrumental and observational distortions. The model is continuous in frequency, can easily be adapted to intensity measurements and extends naturally to the analysis of high-frequency pseudo modes (interference peaks at frequencies above the atmospheric acoustic cutoff).

Context. The use of ratios of small to large separations as a diagnostic of stellar interiors. Aims. To demonstrate that model fitting by comparing observed and model separation ratios at the same n values is in error, and to present a correct procedure. Methods. Theoretical analysis using phase shifts and numerical models. Results. We show that the separation ratios of stellar models with the same interior structure, but different outer layers, are not the same when compared at the same n values, but are the same when evaluated at the same frequencies by interpolation. The separation ratios trace the phase shift differences as a function of frequency not of n. We give examples from model fitting where the ratios at the same n values agree within the error estimates, but do not agree when evaluated at the same frequencies and the models do not have the same interior structure. The correct procedure is to compare observed ratios with those of models interpolated to the observed frequencies.

We report our first results on comparing the variations of the solar internal rotation with solar activity, as predicted by non-linear solar dynamo modelling, with helioseismic measurements using the SOHO MDI data.

Analyses of recent helioseismic data indicate that the dynamical regimes at the base of the convection zone can be different from those observed at the top, having either significantly shorter periods or non--periodic behaviour. Recently spatiotemporal fragmentation/bifurcation has been proposed as a dynamical mechanism to account for the multi-mode behaviour that is possibly observed in the solar convection zone, without requiring separate physical mechanisms with different time scales at different depths. Here we study the robustness of this mechanism with respect to changes to the zero order rotation profile, motivated by the uncertainties of and differences between the various reductions of the helioseimic data. We find that spatiotemporal fragmentation is a common feature of the reductions investigated

We investigate the possibility of observing very small amplitude low frequency solar oscillations with the proposed laser interferometer space antenna LISA. For frequencies below $\sim 2\times 10^{-4}$ Hz the dominant contribution is from the near zone time dependent gravitational quadrupole moments associated with the normal modes of oscillation. For frequencies $\nu$ above $\sim 3\times 10^{-4}$ Hz the dominant contribution is from gravitational radiation generated by the quadrupole oscillations which is larger than the Newtonian signal by a factor $\sim (2 \pi r \nu/ c)^4$, where $r$ is the distance to the Sun, and $c$ is the velocity of light. The low order solar quadrupole pressure and gravity oscillation modes have not yet been detected above the solar background by helioseismic velocity and intensity measurements. Our estimates of the amplitudes needed to give a detectable signal on a LISA type space laser interferometer imply surface velocity amplitudes on the sun of the order of 1-10 mm/sec in the frequency range $1 - 5~10^{-4}$Hz. Such surface velocities are below the current sensitivity limits on helioseismic measurements. If modes exist with frequencies and amplitudes in this range they could be detected with a LISA type laser interferometer.