We investigate the parameters of global solar p-mode oscillations, namely damping width $\Gamma$, amplitude $A$, mean squared velocity $\langle v^2\rangle$, energy $E$, and energy supply rate $\mathrm{d}E/\mathrm{d} t$, derived from two solar cycles' worth (1996 - 2018) of Global Oscillation Network Group (GONG) time series for harmonic degrees $l=0 - 150$. We correct for the effect of fill factor, apparent solar radius, and spurious jumps in the mode amplitudes. We find that the amplitude of the activity related changes of $\Gamma$ and $A$ depends on both frequency and harmonic degree of the modes, with the largest variations of $\Gamma$ for modes with $2400\,\mu\text{Hz}\le \nu \le 3300\,\mu\text{Hz}$ and $31\le l \le 60$ with a min-to-max variation of $26.6\pm0.3\%$ and of $A$ for modes with $2400\,\mu\text{Hz}\le\nu\le 3300\,\mu\text{Hz}$ and $61\le l \le 100$ with a min-to-max variation of $27.4\pm0.4\%$. The level of correlation between the solar radio flux $F_{10.7}$ and mode parameters also depends on mode frequency and harmonic degree. As a function of mode frequency, the mode amplitudes are found to follow an asymmetric Voigt profile with $\nu_{\text{max}}=3073.59\pm0.18\,\mu\text{Hz}$. From the mode parameters, we calculate physical mode quantities and average them over specific mode frequency ranges. This way, we find that the mean squared velocities $\langle v^2\rangle$ and energies $E$ of p modes are anti-correlated with the level of activity, varying by $14.7\pm0.3\%$ and $18.4\pm0.3\%$, respectively, and that the mode energy supply rates show no significant correlation with activity. With this study we expand previously published results on the temporal variation of solar p-mode parameters. Our results will be helpful to future studies of the excitation and damping of p modes, i.e., the interplay between convection, magnetic field, and resonant acoustic oscillations.

The pattern of migrating zonal flow bands associated with the solar cycle, known as the torsional oscillation, has been monitored with continuous global helioseismic observations by the Global Oscillations Network Group, together with those made by the Michelson Doppler Imager onboard the Solar and Heliosepheric Observatory and its successor the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory, since 1995, giving us nearly two full solar cycles of observations. We report that the flows now show traces of the mid-latitude acceleration that is expected to become the main equatorward-moving branch of the zonal flow pattern for Cycle 25. Based on the current position of this branch, we speculate that the onset of widespread activity for Cycle 25 is unlikely to be earlier than the middle of 2019.

Using continuous observations for 22 years from ground-based network GONG and space-borne instruments MDI onboard {\it SoHO} and HMI onboard {\it SDO}, we report both global and local properties of the convection zone and their variations with time.

We examine the frequency shifts in low-degree helioseismic modes from the Birmingham Solar-Oscillations Network (BiSON) covering the period from 1985 - 2016, and compare them with a number of global activity proxies well as a latitudinally-resolved magnetic index. As well as looking at frequency shifts in different frequency bands, we look at a parametrization of the shift as a cubic function of frequency. While the shifts in the medium- and highfrequency bands are very well correlated with all of the activity indices (with the best correlation being with the 10.7 cm radio flux), we confirm earlier findings that there appears to have been a change in the frequency response to activity during solar cycle 23, and the low frequency shifts are less correlated with activity in the last two cycles than they were in Cycle 22. At the same time, the more recent cycles show a slight increase in their sensitivity to activity levels at medium and higher frequencies, perhaps because a greater proportion of activity is composed of weaker or more ephemeral regions. This lends weight to the speculation that a fundamental change in the nature of the solar dynamo may be in progress.

The solar-cyle variation of acoustic mode frequencies has a frequency dependence related to the inverse mode inertia. The discrepancy between model predictions and measured oscillation frequencies for solar and solar-type stellar acoustic modes includes a significant frequency-dependent term known as the surface term that is also related to the inverse mode inertia. We parametrize both the surface term and the frequency variations for low-degree solar data from Birmingham Solar-Oscillations Network (BiSON) and medium-degree data from the Global Oscillations Network Group (GONG) using the mode inertia together with cubic and inverse frequency terms. We find that for the central frequency of rotationally split multiplets the cubic term dominates both the average surface term and the temporal variation, but for the medium-degree case the inverse term improves the fit to the temporal variation. We also examine the variation of the even-order splitting coefficients for the medium-degree data and find that, as for the central frequency, the latitude-dependent frequency variation, which reflects the changing latitudinal distribution of magnetic activity over the solar cycle, can be described by the combination of a cubic and an inverse function of frequency scaled by inverse mode inertia. The results suggest that this simple parametrization could be used to assess the activity-related frequency variation in solar-like asteroseismic targets.

Continuous high-cadence and high-spatial resolution Dopplergrams allow us to study sub-surface dynamics that may be further extended to explore precursors of visible solar activity on the surface. Since the p-mode power is absorbed in the regions of high magnetic field, the inferences in these regions are often presumed to have large uncertainties. In this paper, using the Dopplergrams from space-borne Helioseismic Magnetic Imager (HMI), we compare horizontal flows in a shear layer below the surface and the photospheric layer in and around active regions. The photospheric flows are calculated using local correlation tracking (LCT) method while the ring-diagram (RD) technique of helioseismology is used to infer flows in the sub-photospheric shear layer. We find a strong positive correlation between flows from both methods near the surface. This implies that despite the absorption of acoustic power in the regions of strong magnetic field, the flows inferred from the helioseismology are comparable to those from the surface measurements. However, the magnitudes are significantly different; the flows from the LCT method are smaller by a factor of 2 than the helioseismic measurements. Also, the median difference between direction of corresponding vectors is 49 degree.

Near-surface flows measured by the ring-diagram technique of local helioseismology show structures that persist over multiple rotations. We examine these phenomena using data from the {\em Global Oscillation Network Group} (GONG) and the {\em Helioseismic and Magnetic Imager} (HMI) and show that a correlation analysis of the structures can be used to estimate the rotation rate as a function of latitude, giving a result consistent with the near-surface rate from global helioseismology and slightly slower than that obtained from a similar analysis of the surface magnetic field strength. At latitudes of 60$^{\circ}$ and above the HMI flow data reveal a strong signature of a two-sided zonal flow structure. This signature may be related to recent reports of "giant cells" in solar convection.

We compare the average photospheric current helicity $H_c$, photospheric twist parameter $\alpha$ (a well-known proxy for the full relative magnetic helicity), and subsurface kinetic helicity $H_k$ for 194 active regions observed between 2006-2013. We use 2440 Hinode photospheric vector magnetograms, and the corresponding subsurface fluid velocity data derived from GONG (2006-2012) and HMI (2010-2013) Dopplergrams. We find a significant hemispheric bias in all three parameters. The subsurface kinetic helicity is preferentially positive in the southern hemisphere and negative in the northern hemisphere. The photospheric current helicity and the $\alpha$ parameter have the same bias for strong fields ($|B|>1000$ G) and no significant bias for weak fields (100 G $<|B|<500$ G). We find no significant region-by-region correlation between the subsurface kinetic helicity and either the strong-field current helicity or $\alpha$. Subsurface fluid motions of a given handedness correspond to photospheric helicities of both signs in approximately equal numbers. However, common variations appear in annual averages of these quantities over all regions. Furthermore, in a subset of 77 regions we find significant correlations between the temporal profiles of the subsurface and photospheric helicities. In these cases, the sign of the linear correlation coefficient matches the sign relationship between the helicities, indicating that the photospheric magnetic field twist is sensitive to the twisting motions below the surface.

The temporal variation of the horizontal velocity in subsurface layers beneath three different types of active regions is studied using the technique of ring diagrams. In this study, we select active regions (ARs) 10923, 10930, 10935 from three consecutive Carrington rotations: AR 10930 contains a fast-rotating sunspot in a strong emerging active region while other two have non-rotating sunspots with emerging flux in AR 10923 and decaying flux in AR 10935. The depth range covered is from the surface to about 12 Mm. In order to minimize the influence of systematic effects, the selection of active and quiet regions is made so that these were observed at the same heliographic locations on the solar disk. We find a significant variation in both components of the horizontal velocity in active regions as compared to quiet regions. The magnitude is higher in emerging-flux regions than in the decaying-flux region, in agreement with earlier findings. Further, we clearly see a significant temporal variation in depth profiles of both zonal and meridional flow components in AR 10930, with the variation in the zonal component being more pronounced. We also notice a significant influence of the plasma motion in areas closest to the rotating sunspot in AR 10930 while areas surrounding the non-rotating sunspots in all three cases are least affected by the presence of the active region in their neighborhood.

In this paper, we compare components of the horizontal flow below the solar surface in and around regions consisting of rotating and non-rotating sunspots. Our analysis suggests that there is a significant variation in both components of the horizontal flow at the beginning of sunspot rotation as compared to the non-rotating sunspot. In most cases, the flows in surrounding areas are relatively small. However, there is a significant influence of the motion on flows in an area closest to the sunspot rotation.

Local helioseismology provides unique information about the subsurface structure and dynamics of sunspots and active regions. However, because of complexity of sunspot regions local helioseismology diagnostics require careful analysis of systematic uncertainties and physical interpretation of the inversion results. We present new results of comparison of the ring-diagram analysis and time-distance helioseismology for active region NOAA 9787, for which a previous comparison showed significant differences in the subsurface sound-speed structure, and discuss systematic uncertainties of the measurements and inversions. Our results show that both the ring-diagram and time-distance techniques give qualitatively similar results, revealing a characteristic two-layer seismic sound-speed structure consistent with the results for other active regions. However, a quantitative comparison of the inversion results is not straightforward. It must take into account differences in the sensitivity, spatial resolution and the averaging kernels. In particular, because of the acoustic power suppression, the contribution of the sunspot seismic structure to the ring-diagram signal can be substantially reduced. We show that taking into account this effect reduces the difference in the depth of transition between the negative and positive sound-speed variations inferred by these methods. Further detailed analysis of the sensitivity, resolution and averaging properties of the local helioseismology methods is necessary for consolidation of the inversion results. It seems to be important that both methods indicate that the seismic structure of sunspots is rather deep and extends to at least 20 Mm below the surface, putting constraints on theoretical models of sunspots.

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

We show here a component of the meridional circulation developing at medium-high latitudes (40-50 degrees) before the new solar cycle starts. Like the torsional oscillation of the zonal flows, this extra circulation seems to precede the onset of magnetic activity at the solar surface and move slowly towards lower latitudes. However, the behavior of this component differs from that of the torsional oscillation regarding location and convergence towards the equator at the end of the cycle. The observation of this component before the magnetic regions appear at the solar surface has only been possible due to the prolonged solar minimum. The results could settle the discussion as to whether the extra component of the meridional circulation around the activity belts, which has been known for some time, is or is not an effect of material motions around the active regions.

The availability of continuous helioseismic data for two consecutive solar minima has provided a unique opportunity to study the changes in the solar interior that might have led to this unusual minimum. We present preliminary analysis of inter mediate-degree mode frequencies in the 3 mHz band during the current period of minimal solar activity and show that the mode frequencies are significantly lower than those during the previous activity minimum. Our analysis do not show any signature of the beginning of cycle 24 till the end of 2008. In addition, the zonal and meridional flow patterns inferred from inverting frequencies also hint for a delayed onset of a new cycle. The estimates of travel time are higher than the previous minimum confirming a relatively weak solar activity during the current minimum.

We examine the evolution of the zonal flow pattern in the upper solar convection zone during the current extended solar minimum, and compare it with that during the previous minimum. The results suggest that a configuration matching that at the previous minimum was reached during 2008, but that the flow band corresponding to the new cycle has been moving more slowly towards the equator than was observed in the previous cycle, resulting in a gradual increase in the apparent length of the cycle during the 2007 -- 2008 period. The current position of the lower-latitude fast-rotating belt corresponds to that seen around the onset of activity in the previous cycle.

Temporal variations of the subsurface meridional flow with the solar cycle have been reported by several authors. The measurements are typically averaged over periods of time during which surface magnetic activity existed in the regions were the velocities are calculated. The present work examines the possible contamination of these measurements due to the extra velocity fields associated with active regions plus the uncertainties in the data obtained where strong magnetic fields are present. We perform a systematic analysis of more than five years of GONG data and compare meridional flows obtained by ring-diagram analysis before and after removing the areas of strong magnetic field. The overall trend of increased amplitude of the meridional flow towards solar minimum remains after removal of large areas associated with surface activity. We also find residual circulation toward the active belts that persist even after the removal of the surface magnetic activity, suggesting the existence of a global pattern or longitudinally-located organized flows.

We study the variation of the frequency splitting coefficients describing the solar asphericity in both GONG and MDI data, and use these data to investigate temporal sound-speed variations as a function of both depth and latitude during the period from 1995-2000 and a little beyond. The temporal variations in even splitting coefficients are found to be correlated to the corresponding component of magnetic flux at the solar surface. We confirm that the sound-speed variations associated with the surface magnetic field are superficial. Temporally averaged results show a significant excess in sound speed around 0.92 solar radii and latitude of 60 degrees.

Helioseismic observations have detected small temporal variations of the rotation rate below the solar surface corresponding to the so-called `torsional oscillations' known from Doppler measurements of the surface. These appear as bands of slower and faster than average rotation moving equatorward. Here we establish, using complementary helioseismic observations over four years from the GONG network and from the MDI instrument on board SOHO, that the banded flows are not merely a near-surface phenomenon: rather they extend downward at least 60 Mm (some 8% of the total solar radius) and thus are evident over a significant fraction of the nearly 200 Mm depth of the solar convection zone.