Global-scale Rossby waves have recently been unambiguously identified on the Sun. Here we study the latitude and depth dependence of the Rossby wave eigenfunctions. By applying helioseismic ring-diagram analysis and granulation tracking to SDO/HMI observations, we compute maps of the radial vorticity of flows in the upper solar convection zone (down to depths of more than $16$ Mm). We use a Fourier transform in longitude to separate the different azimuthal orders m in the range $3 \le m \le 15$. At each $m$ we obtain the phase and amplitude of the Rossby waves as a function of depth using the helioseismic data. At each $m$ we also measure the latitude dependence of the eigenfunctions by calculating the covariance between the equator and other latitudes. We then study 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 $\pm 30^\circ$, 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 roughly $\pm 5^\circ$ over the top $15$ Mm. The amplitude of the radial eigenfunctions decreases by about $10\%$ from the surface down to $8$ Mm (the region where ring-diagram analysis is most reliable, as seen by comparing with the rotation rate measured by global-mode seismology). 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.

Over the past decades and even centuries, the astronomical community has accumulated a signif-icant heritage of recorded observations of a great many astronomical objects. Those records con-tain irreplaceable information about long-term evolutionary and non-evolutionary changes in our Universe, and their preservation and digitization is vital. Unfortunately, most of those data risk becoming degraded and thence totally lost. We hereby call upon the astronomical community and US funding agencies to recognize the gravity of the situation, and to commit to an interna-tional preservation and digitization efforts through comprehensive long-term planning supported by adequate resources, prioritizing where the expected scientific gains, vulnerability of the origi-nals and availability of relevant infrastructure so dictates. The importance and urgency of this issue has been recognized recently by General Assembly XXX of the International Astronomical Union (IAU) in its Resolution B3: "on preservation, digitization and scientific exploration of his-torical astronomical data". We outline the rationale of this promotion, provide examples of new science through successful recovery efforts, and review the potential losses to science if nothing it done.

The Sun's complex dynamics is controlled by buoyancy and rotation in the convection zone and by magnetic forces in the atmosphere and corona. While small-scale solar convection is well understood, the dynamics of large-scale flows in the solar convection zone is not explained by theory or simulations. Waves of vorticity due to the Coriolis force, known as Rossby waves, are expected to remove energy out of convection at the largest scales. Here we unambiguously detect and characterize retrograde-propagating vorticity waves in the shallow subsurface layers of the Sun at angular wavenumbers below fifteen, with the dispersion relation of textbook sectoral Rossby waves. The waves have lifetimes of several months, well-defined mode frequencies below 200 nHz in a co-rotating frame, and eigenfunctions of vorticity that peak at the equator. Rossby waves have nearly as much vorticity as the convection at the same scales, thus they are an essential component of solar dynamics. We find a transition from turbulence-like to wave-like dynamics around the Rhines scale of angular wavenumber of twenty; this might provide an explanation for the puzzling deficit of kinetic energy at the largest spatial scales.

We analyze the effect of a sunspot in its quiet surroundings applying a helioseismic technique on almost three years of Helioseismic and Magnetic Imager (HMI) observations obtained during solar cycle 24 to further study the sunspot structure below the solar surface. The attenuation of acoustic waves with frequencies lower than 4.2 mHz depends more strongly on the wave direction at a distance of $6^o-7^o$ from the sunspot center. The amplification of higher frequency waves is highest 6$^o$ away from the active region and it is largely independent of the wave's direction. We observe a mean clockwise flow around active regions, which angular speed decreases exponentially with distance and has a coefficient close to $-0.7$ degree$^{-1}$. The observed horizontal flow in the direction of the nearby active region agrees with a large-scale circulation around the sunspot in the shape of cylindrical shell. The center of the shell seems to be centered around 7$^o$ from the sunspot center, where we observe an inflow close to the surface down to $\sim$2 Mm, followed by an outflow at deeper layers until at least 7 Mm.

In order to quantify the influence of magnetic fields on acoustic mode parameters and flows in and around active regions, we analyse the differences in the parameters in magnetically quiet regions nearby an active region (which we call `nearby regions'), compared with those of quiet regions at the same disc locations for which there are no neighboring active regions. We also compare the mode parameters in active regions with those in comparably located quiet regions. Our analysis is based on ring diagram analysis of all active regions observed by HMI during almost five years. We find that the frequency at which the mode amplitude changes from attenuation to amplification in the quiet nearby regions is around 4.2 mHz, in contrast to the active regions, for which it is about 5.1 mHz. This amplitude enhancement (the `acoustic halo effect') is as large as that observed in the active regions, and has a very weak dependence on the wave propagation direction. The mode energy difference in nearby regions also changes from a deficit to an excess at around 4.2 mHz, but averages to zero over all modes. The frequency difference in nearby regions increases with increasing frequency until a point at which the frequency shifts turn over sharply as in active regions. However, this turnover occurs around 4.9 mHz, significantly below the acoustic cutoff frequency. Inverting the horizontal flow parameters in the direction of the neighboring active regions, we find flows that are consistent with a model of the thermal energy flow being blocked directly below the active region.

NASA's Solar Dynamics Observatory (SDO) was launched 11 February 2010 with 3 instruments on board, including the Helioseismic and Magnetic Imager (HMI). Since beginning normal operations on 1 May 2010, HMI has observed the Sun's entire visible disk almost continuously. HMI collects sequences of polarized filtergrams taken at a fixed cadence with two 4096 x 4096 cameras from which are computed arcsecond-resolution maps of photospheric observables: the line-of-sight (LoS) velocity and magnetic field, continuum intensity, line width, line depth, and the Stokes polarization parameters, I Q U V, at 6 wavelengths. Two processing pipelines implemented at the SDO Joint Science Operations Center (JSOC) at Stanford University compute observables from calibrated Level-1 filtergrams. One generates LoS quantities every 45s, and the other, primarily for the vector magnetic field, computes averages on a 720s cadence. Corrections are made for static and temporally changing CCD characteristics, bad pixels, image alignment and distortion, polarization irregularities, filter-element uncertainty and non-uniformity, as well as Sun-spacecraft velocity. This report explains issues affecting the resulting physical quantities, describes the impact of regular updates to the instrument calibration, and shows how the computations are optimized for actual HMI observations. During the 5 years of the SDO prime mission, regular calibration sequences have been used to regularly improve and update the instrument calibration and to monitor instrument changes. The instrument more than satisfies the original specifications for data quality and continuity. The procedures described here still have significant room for improvement. The most significant remaining systematic errors are associated with the spacecraft orbital velocity.

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.

Ring-diagram analysis of acoustic waves observed at the photosphere can provide a relatively robust determination of the sub-surface flows at a particular time under a particular region. The depth of penetration of the waves is related to the size of the region, hence the depth extent of the measured flows is inversely proportional to the spatial resolution. Most ring-diagram analysis has focused on regions of extent ~15{\deg} (180 Mm) or more in order to provide reasonable mode sets for inversions. HMI data analysis also provides a set of ring fit parameters on a scale three times smaller. These provide flow estimates for the outer 1% (7 Mm) of the Sun only, with very limited depth resolution, but with spatial resolution adequate to map structures potentially associated with the belts and regions of magnetic activity. There are a number of systematic effects affecting the determination of flows from local helioseismic analysis of regions over different parts of the observable disk, not all well understood. In this study we characterize those systematic effects with higher spatial resolution, so that they may more effectively be accounted for in mapping temporal and spatial evolution of the flows. Leaving open the question of the mean structure of the global meridional circulation and the differential rotation, we describe the near-surface flow anomalies in time and latitude corresponding to the torsional oscillation pattern in differential rotation and analogous patterns in the meridional cell structure as observed by SDO/HMI.

Solar subsurface zonal and meridional-flow profiles during the rising phase of solar cycle 24 are studied using time-distance helioseismology technique. The faster zonal bands in the torsional-oscillation pattern show strong hemispheric asymmetries and temporal variations in both width and speed. The faster band in the northern hemisphere is located closer to the equator than the band in the southern hemisphere, and migrates past the equator when the magnetic activity in the southern hemisphere is reaching maximum. The meridional-flow speed decreases substantially with the increase of magnetic activity, and the flow profile shows two zonal structures in each hemisphere. The residual meridional flow, after subtracting a mean meridional-flow profile, converges toward the activity belts and shows faster and slower bands like the torsional-oscillation pattern. More interestingly, the meridional-flow speed above latitude 30 degree shows an anti-correlation with the poleward-transporting magnetic flux, slower when the following-polarity flux is transported and faster when the leading-polarity flux is transported. It is expected that this phenomenon slows the process of magnetic cancellation and polarity reversal in the high-latitude areas.

We report observations of white-light ejecta in the low corona, for two X-class flares on the 2013 May 13, using data from the Helioseismic and Magnetic Imager (HMI) of the Solar Dynamics Observatory. At least two distinct kinds of sources appeared (chromospheric and coronal), in the early and later phases of flare development, in addition to the white-light footpoint sources commonly observed in the lower atmosphere. The gradual emissions have a clear identification with the classical loop-prominence system, but are brighter than expected and possibly seen here in the continuum rather than line emission. We find the HMI flux exceeds the radio/X-ray interpolation of the bremsstrahlung produced in the flare soft X-ray sources by at least one order of magnitude. This implies the participation of cooler sources that can produce free-bound continua and possibly line emission detectable by HMI. One of the early sources dynamically resembles "coronal rain", appearing at a maximum apparent height and moving toward the photosphere at an apparent constant projected speed of 134 $\pm$ 8 $\mathrm{km s^{-1}}$. Not much literature exists on the detection of optical continuum sources above the limb of the Sun by non-coronagraphic instruments, and these observations have potential implications for our basic understanding of flare development, since visible observations can in principle provide high spatial and temporal resolution.

Meridional flow in the solar interior plays an important role in redistributing angular momentum and transporting magnetic flux inside the Sun. Although it has long been recognized that the meridional flow is predominantly poleward at the Sun's surface and in its shallow interior, the location of the equatorward return flow and the meridional flow profile in the deeper interior remain unclear. Using the first two years of continuous helioseismology observations from the Solar Dynamics Observatory / Helioseismic Magnetic Imager, we analyze travel times of acoustic waves that propagate through different depths of the solar interior carrying information about the solar interior dynamics. After removing a systematic center-to-limb effect in the helioseismic measurements and performing inversions for flow speed, we find that the poleward meridional flow of a speed of 15 m/s extends in depth from the photosphere to about 0.91 R_sun. An equatorward flow of a speed of 10 m/s is found between 0.82 to 0.91 R_sun in the middle of the convection zone. Our analysis also shows evidence of that the meridional flow turns poleward again below 0.82 R_sun, indicating an existence of a second meridional circulation cell below the shallower one. This double-cell meridional circulation profile with an equatorward flow shallower than previously thought suggests a rethinking of how magnetic field is generated and redistributed inside the Sun.

We study temporal variations in the amplitudes and widths of high-degree acoustic modes by applying the ring-diagram technique to the GONG+, MDI and HMI Dopplergrams during the declining phase of cycle 23 and rising phase of cycle 24. The mode parameters from all three instruments respond similarly to the varying magnetic activity. The mode amplitudes and widths show consistently lower variation due to smaller magnetic activity in cycle 24 as compared to the previous solar cycle.

We employ ring-diagram analysis to study the sub-surface thermal structure of active regions. We present results using a large number of active regions over the course of Solar Cycle 23. We present both traditional inversions of ring-diagram frequency differences, with a total sample size of 264, and a statistical study using Principal Component Analysis. We confirm earlier results on smaller samples that sound speed and adiabatic index are changed below regions of strong magnetic field. We find that sound speed is decreased in the region between approximately r=0.99R_sun and r=0.995R_sun (depths of 3Mm to 7Mm), and increased in the region between r=0.97R_sun and r=0.985R_sun (depths of 11Mm to 21Mm). The adiabatic index is enhanced in the same deeper layers that sound-speed enhancement is seen. A weak decrease in adiabatic index is seen in the shallower layers in many active regions. We find that the magnitudes of these perturbations depend on the strength of the surface magnetic field, but we find a great deal of scatter in this relation, implying other factors may be relevant.

While the {\it Helioseismic and Magnetic Imager} (HMI) onboard the {\it Solar Dynamics Observatory} (SDO) provides Doppler velocity [$V$], continuum intensity [$I_C$], and line-depth [$Ld$] observations, each of which is sensitive to the five-minute acoustic spectrum, the {\it Atmospheric Imaging Array} (AIA) also observes at wavelengths -- specifically the 1600 and 1700 Angstrom bands -- that are partly formed in the upper photosphere and have good sensitivity to acoustic modes. In this article we consider the characteristics of the spatio--temporal Fourier spectra in AIA and HMI observables for a 15-degree region around NOAA Active Region 11072. We map the spatio--temporal-power distribution for the different observables and the HMI Line Core [$I_L$], or Continuum minus Line Depth, and the phase and coherence functions for selected observable pairs, as a function of position and frequency. Five-minute oscillation power in all observables is suppressed in the sunspot and also in plage areas. Above the acoustic cut-off frequency, the behaviour is more complicated: power in HMI $I_C$ is still suppressed in the presence of surface magnetic fields, while power in HMI $I_L$ and the AIA bands is suppressed in areas of surface field but enhanced in an extended area around the active region, and power in HMI $V$ is enhanced in a narrow zone around strong-field concentrations and suppressed in a wider surrounding area. The relative phase of the observables, and their cross-coherence functions, are also altered around the active region. These effects may help us to understand the interaction of waves and magnetic fields in the different layers of the photosphere, and will need to be taken into account in multi-wavelength local helioseismic analysis of active regions.

We report on a systematic center-to-limb variation in measured helioseismic travel times, which must be taken into account for an accurate determination of solar interior meridional flows. The systematic variation, found in time-distance helioseismology analysis using SDO/HMI and SDO/AIA observations, is different in both travel-time magnitude and variation trend for different observables. It is not clear what causes this systematic effect. Subtracting the longitude-dependent east-west travel times, obtained along the equatorial area, from the latitude-dependent north-south travel times, obtained along the central meridian area, gives remarkably similar results for different observables. We suggest this as an effective procedure for removing the systematic center-to-limb variation. The subsurface meridional flows obtained from inversion of the corrected travel times are approximately 10 m/s slower than those obtained without removing the systematic effect. The detected center-to-limb variation may have important implications in the derivation of meridional flows in the deep interior, and needs a better understanding.

High degree solar mode frequencies as measured by ring diagrams are known to change in the presence of the strong magnetic fields found in active regions. We examine these changes in frequency for a large sample of active regions analyzed with data from the Michelson Doppler Imager (MDI) onboard the SoHO spacecraft, spanning most of solar cycle 23. We confirm that the frequencies increase with increasing magnetic field strength, and that this dependence is generally linear. We find that the dependence is slightly but significantly different for active regions with different sunspot types.

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

We study temporal variations in the amplitudes and widths of high-degree acoustic modes in the quiet and active Sun by applying ring-diagram technique to the GONG+ and MDI Dopplergrams during the declining phase of cycle 23. The increase in amplitudes and decrease in line-widths in the declining phase of the solar activity is in agreement with previous studies. A similar solar cycle trend in the mode parameters is also seen in the quiet-Sun regions but with a reduced magnitude. Moreover, the amplitudes obtained from GONG+ data show long-term variations on top of the solar cycle trend.

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.

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.

In this paper we describe in detail the implementation and main properties of a new inversion code for the polarized radiative transfer equation (VFISV: Very Fast inversion of the Stokes vector). VFISV will routinely analyze pipeline data from the Helioseismic and Magnetic Imager (HMI) on-board of the Solar Dynamics Observatory (SDO). It will provide full-disk maps (4096$\times$4096 pixels) of the magnetic field vector on the Solar Photosphere every 10 minutes. For this reason VFISV is optimized to achieve an inversion speed that will allow it to invert 16 million pixels every 10 minutes with a modest number (approx. 50) of CPUs. Here we focus on describing a number of important details, simplifications and tweaks that have allowed us to significantly speed up the inversion process. We also give details on tests performed with data from the spectropolarimeter on-board of the Hinode spacecraft.

We present results on the structure of the near-surface layers of the Sun obtained by inverting frequencies of high-degree solar modes from "ring diagrams". We have results for eight epochs between June 1996 and October 2003. The frequencies for each epoch were obtained from ring diagrams constructed from MDI Dopplergrams spanning complete Carrington rotations. We find that there is a substantial latitudinal variation of both sound speed and the adiabatic index Gamma_1 in the outer 2% of the Sun. We find that both the sound-speed and Gamma_1 profiles change with changes in the level of solar activity. In addition, we also study differences between the northern and southern hemispheres of the Sun and find a small asymmetry that appears to reflect the difference in magnetic activity between the two hemispheres.