Active region NOAA 12673 was extremely volatile in 2017 September, producing many solar flares, including the largest of solar cycle 24, an X9.3 flare of 06 September 2017. It has been reported that this flare produced a number of sunquakes along the flare ribbon (Sharykin & Kosovichev 2018; Zhao & Chen 2018). We have used co-temporal and co-spatial Helioseismic and Magnetic Imager (HMI) line-of-sight (LOS) and Swedish 1-m Solar Telescope observations to show evidence of the chromospheric response to these sunquakes. Analysis of the Ca II 8542 \AA\space line profiles of the wavefronts revealed that the crests produced a strong blue asymmetry, whereas the troughs produced at most a very slight red asymmetry. We used the combined HMI, SST datasets to create time-distance diagrams and derive the apparent transverse velocity and acceleration of the response. These velocities ranged from 4.5 km s$^{-1}$ to 29.5 km s$^{-1}$ with a constant acceleration of 8.6 x 10$^{-3}$ km s$^{-2}$. We employed NICOLE inversions, in addition to the Center-of-Gravity (COG) method to derive LOS velocities ranging 2.4 km s$^{-1}$ to 3.2 km s$^{-1}$. Both techniques show that the crests are created by upflows. We believe that this is the first chromospheric signature of a flare induced sunquake.

Sunquakes are created by the hydrodynamic response of the lower atmosphere to a sudden deposition of energy and momentum. In this study we investigate a sunquake that occurred in NOAA active region 11675 on 17 February 2013. Observations of the corona, chromosphere and photosphere are brought together for the first time with a non-linear force-free model of the active region's magnetic field in order to probe the magnetic environment in which the sunquake was initiated. We find that the sunquake was associated with the destabilization of a flux rope and an associated M-class GOES flare. Active region 11675 was in its emergence phase at the time of the sunquake and photospheric motions caused by the emergence heavily modified the flux rope and its associated quasi-separatrix layers, eventually triggering the flux rope's instability. The flux rope was surrounded by an extended envelope of field lines rooted in a small area at the approximate position of the sunquake. We argue that the configuration of the envelope, by interacting with the expanding flux rope, created a "magnetic lens" that may have focussed energy in one particular location the photosphere, creating the necessary conditions for the initiation of the sunquake.

Using a summary curve of two eigen vectors of solar magnetic field oscillations derived with Principal Components Analysis (PCA) from synoptic maps for solar cycles 21-24 as a proxy of solar activity, we extrapolate this curve backwards three millennia revealing 9 grand cycles lasting 350-400 years each. The summary curve shows a remarkable resemblance to the past sunspot and terrestrial activity: grand minima - Maunder Minimum (1645-1715 AD), Wolf minimum (1280-1350 AD), Oort minimum (1010-1050 AD) and Homer minimum (800-900 BC); grand maxima - modern warm period (1990-2015), medieval warm period (900-1200 AD), Roman warm period (400-10 BC) and others. We verify the extrapolated activity curve by the pre-telescope observations of large sunspots with naked eye, by comparing the observed and simulated butterfly diagrams for Maunder Minimum (MM), by a maximum of the terrestrial temperature and extremely intense terrestrial auroras seen in the past grand cycle occurred in 14-16 centuries. We confirm the occurrence of upcoming Modern grand minimum in 2020-2053, which will have a shorter duration (3 cycles) and, thus, higher solar activity compared to MM. We argue that Sporer minimum (1450-1550) derived from the increased abundances of isotopes 14C and 10Be is likely produced by a strong increase of the terrestrial background radiation caused by the galactic cosmic rays of powerful supernovae.

The presence of flare related acoustic emission (sunquakes) in some flares represents a severe challenge to our current understanding of flare energy transport processes. We present a comparison of new spectral observations from Hinode's EUV imaging Spectrometer (EIS) and the Interface Region Imaging Spectrograph (IRIS) of the atmosphere above a sunquake, and compare them to the spectra observed in a part of the flaring region with no acoustic signature. Evidence for the sunquake is determined using both time-distance and acoustic holography methods, and we find that, unlike many previous sunquake detections, the signal is rather dispersed, but that the time-distance and 6 and 7 mHz sources converge at the same spatial location. We also see some evidence for different evolution at different frequencies, with an earlier peak at 7 mHz than at 6 mHz. Using spectroscopic measurements we find that in this location at the time of the 7 mHz peak the spectral emission is significantly more intense, shows larger velocity shifts and substantially broader profiles than in the location with no sunquake, and that there is a good correlation between blue-shifted, hot coronal, hard X-ray (HXR) and red-shifted chromospheric emission, consistent with the idea of a strong downward motion driven by rapid heating by non-thermal electrons and the formation of chromospheric shocks. Exploiting the diagnostic potential of the Mg II triplet lines, we also find evidence for a single, large temperature increase deep in the atmosphere, consistent with this scenario. The time of the 6 mHz and time-distance peak signal coincides with a secondary peak in the energy release process, but in this case we find no evidence of HXR emission in the quake location, but very broad spectral lines, strongly shifted to the red, indicating the possible presence of a significant flux of downward propagating Alfven waves.

Several mechanisms have been proposed to explain the transient seis- mic emission, i.e., sunquakes, from some solar flares. Some theories associate high-energy electrons and/or white-light emission with sunquakes. High-energy charged particles and their subsequent heating of the photosphere and/or chro- mosphere could induce acoustic waves in the solar interior. We carried out a correlative study of solar flares with emission in hard-X rays (HXRs), enhanced continuum emission at 6173{\AA}, and transient seismic emission. We selected those flares observed by RHESSI (Reuven Ramaty High Energy Solar Spectroscopic Imager) with a considerable flux above 50 keV between January 1, 2010 and June 26, 2014. We then used data from the Helioseismic and Magnetic Imager onboard the Solar Dynamic Observatory (SDO/HMI) to search for excess visible continuum emission and new sunquakes not previously reported. We found a total of 18 sunquakes out of 75 investigated. All of the sunquakes were associated with a enhancement of the visible continuum during the flare time. Finally, we calculated a coefficient of correlation for a set of dichotomic variables related to these observations. We found a strong correlation between two of the standard helioseismic detection techniques, and between sunquakes and visible continuum enhancements. We discuss the phenomenological connectivity between these physical quantities and the observational difficulties of detecting seismic signals and excess continuum radiation.

We present observational evidence for the presence of MHD waves in the solar photosphere deduced from SOHO/MDI Dopplergram velocity observations. The magneto-acoustic perturbations are observed as acoustic power enhancement in the sunspot umbra at high-frequency bands in the velocity component perpendicular to the magnetic field. We use numerical modelling of wave propagation through localised non-uniform magnetic field concentration along with the same filtering procedure as applied to the observations to identify the observed waves. Guided by the results of the numerical simulations we classify the observed oscillations as magneto-acoustic waves excited by the trapped sub-photospheric acoustic waves. We consider the potential application of the presented method as a diagnostic tool for magnetohelioseismology.

We present results of our preliminary analysis of acoustically active X-class flare of September 7, 2011. We report two acoustic sources detected via acoustic holography and verified by finding a ridge in time-distance diagrams. We compare the directional information extracted from time-distance and acoustic holography, showing a good agreement in this case. We report that the direction where amplitude of the wave-front is the largest lies through the strong magnetic field and sunspot, suggesting that absorption of the acoustic wave power by magnetic field can be ruled out as a wave anisotropy mechanism in this case.

Acoustic waves generated by a point source in stratified plasma are considered in this paper. Analytical parametric solution for monochromatic source is derived for plane-parallel polytrope model of the solar interior. The solution is used to gain insight into the properties of the generated wavefront as a function of excitation frequency and depth. A slowly varying pressure perturbation moving in upper layers of solar photosphere with supersonic speed is also considered. It is shown to excite acoustic waves putting certain restrictions upon their geometry of the generated wavefront. The results are discussed in relation to flare generated sunquakes.

The first near-side X-class flare of the Solar Cycle 24 occurred in February 2011 and produced a very strong seismic response in the photosphere. One sunquake was reported by Kosovichev (2011) followed by the discovery of a second sunquake by Zharkov et al (2011). The flare had a two-ribbon structure and was associated with a flux rope eruption and a halo coronal mass ejection (CME) as reported in the CACTus catalogue. Following the discovery of the second sunquake and the spatial association of both sources with the locations of the feet of the erupting flux rope (Zharkov et al 2011) we present here a more detailed analysis of the observed photospheric changes in and around the seismic sources. These sunquakes are quite unusual, taking place early in the impulsive stage of the flare, with the seismic sources showing little hard X-ray (HXR) emission, and strongest X-ray emission sources located in the flare ribbons. We present a directional time--distance diagram computed for the second source, which clearly shows a ridge corresponding to the travelling acoustic wave packet and find that the quake at the second source happened about 45 seconds to one minute earlier than the first source. Using acoustic holography we report different frequency responses of the two sources. We find strong downflows at both seismic locations and a supersonic horizontal motion at the second site of acoustic wave excitation.

We present an analysis of the 15 February 2011 X-class solar flare, previously reported to produce the first sunquake in solar cycle 24 (Kosovichev 2011). Using acoustic holography, we confirm the first, and report a second, weaker, seismic source associated with this flare. We find that the two sources are located at either end of a sigmoid which indicates the presence of a flux rope. Contrary to the majority of previously reported sunquakes, the acoustic emission precedes the peak of major hard X-ray (HXR) sources by several minutes. Furthermore, the strongest hard X-ray footpoints derived from RHESSI data are found to be located away from the seismic sources in the flare ribbons. We account for these discrepancies within the context of a phenomenological model of a flux rope eruption and accompanying two-ribbon flare. We propose that the sunquakes are triggered at the footpoints of the erupting flux rope at the start of the flare impulsive phase and eruption onset, while the main hard X-ray sources appear later at the footpoints of the flare loops formed under the rising flux rope. Possible implications of this scenario for the theoretical interpretation of the forces driving sunquakes are discussed.

The first observations of seismic responses to solar flares were carried out using time-distance (TD) and holography techniques applied to SOHO/MDI Dopplergrams obtained from space and un-affected by terrestrial atmospheric disturbances. However, the ground-based network GONG is potentially a very valuable source of sunquake observations, especially in cases where space observations are unavailable. In this paper we present updated technique for pre-processing of GONG observations for application of subjacent vantage holography. Using this method and TD diagrams we investigate several sunquakes observed in association with M and X-class solar flares and compare the outcomes with those reported earlier using MDI data. In both GONG and MDI datasets, for the first time, we also detect the TD ridge associated with the September 9, 2001 flare. Our results show reassuringly positive identification of sunquakes from GONG data that can provide further information about the physics of seismic processes associated with solar flares.

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.

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.

Aims. In this paper we analyse numerically the propagation and dispersion of acoustic waves in the solar-like sub-photosphere with localised non-uniform magnetic field concentrations, mimicking sunspots with various representative magnetic field configurations. Methods. Numerical simulations of wave propagation through the solar sub-photosphere with a localised magnetic field concentration are carried out using SAC, which solves the MHD equations for gravitationally stratified plasma. The initial equilibrium density and pressure stratifications are derived from a standard solar model. Acoustic waves are generated by a source located at the height approximately corresponding to the visible surface of the Sun. We analyse the response of vertical velocity to changes in the interior due to magnetic field at the level corresponding to the visible solar surface, by the means of local time-distance helioseismology. Results. The results of numerical simulations of acoustic wave propagation and dispersion in the solar sub-photosphere with localised magnetic field concentrations of various types are presented. Time-distance diagrams of the vertical velocity perturbation at the level corresponding to the visible solar surface show that the magnetic field perturbs and scatters acoustic waves and absorbs the acoustic power of the wave packet. For the weakly magnetised case the effect of magnetic field is mainly thermodynamic, since the magnetic field changes the temperature stratification. However, we observe the signature of slow magnetoacoustic mode, propagating downwards, for the strong magnetic field cases.

We investigate the emergence of Active Region NOAA 10790 by means of time--distance helioseismology. Shallow regions of increased sound speed at the location of increased magnetic activity are observed, with regions becoming deeper at the locations of sunspot pores. We also see a long-lasting region of decreased sound speed located underneath the region of the flux emergence, possibly relating to a temperature perturbation due to magnetic quenching of eddy diffusivity, or to a dense flux tube. We detect and track an object in the subsurface layers of the Sun characterised by increased sound speed which could be related to emerging magnetic flux and thus obtain a provisional estimate of the speed of emergence of around $1 {\rm km s^{-1}}$.

The 3 seismic sources S1, S2 and S3 detected from MDI dopplergrams using the time-distance diagram technique are presented with the locations, areas and vertical and horizontal velocities of the visible wave displacements. Within the datacube of 120 Mm the horizontal velocities and the wave propagation times slightly vary from source to source. The momenta and start times measured from the TD diagrams in the sources S1-S3 are compared with those delivered to the photosphere by different kinds of high energy particles with the parameters deduced from hard X-ray and $\gamma$-ray emission as well as by the hydrodynamic shocks caused by these particles. The energetic protons (power laws combined with quasi-thermal ones, or jets) are shown to deliver momentum high enough and to form the hydrodynamic shocks deeply in a flaring atmosphere that allows them to be delivered to the photosphere through much shorter distances and times. Then the seismic waves observed in the sources S2 and S3 can be explained by the momenta produced by hydrodynamic shocks which are caused by mixed proton beams and jets occurring nearly simultaneously with the third burst of hard X-ray (HXR) and $\gamma$-ray emission in the loops with footpoints in the locations of these sources. The seismic wave in the source S1, delayed by 4 and 2 minutes from the first and second HXR bursts, respectively, is likely to be associated with a hydrodynamic shock occurring in this loop from precipitation of a very powerful and hard electron beam with higher energy cutoff mixed with quasi-thermal protons generated by either of these 2 bursts.

We present a comparative seismic study of conditions around and beneath isolated sunspots. Using the European Grid of Solar Observations' Solar Feature Catalogue of sunspots derived from SOHO/MDI continuum and magnetogram data, 1996-2005, we identify a set of isolated sunspots by checking that within a Carrington Rotation there were no other spots detected in the vicinity. We then use level-2 tracked MDI Dopplergrams available from SOHO website to investigate wave-speed perturbations of such sunspots using time-distance helioseismology.