We present a comprehensive study of one method for measuring various parameters of global modes of oscillation of the Sun. Using velocity data taken by the Michelson Doppler Imager (MDI), we analyze spherical harmonic degrees l <= 300. Both current and historical methodologies are explained, and the various differences between the two are investigated to determine their effects on global-mode parameters and systematic errors in the analysis. These differences include a number of geometric corrections made during spherical harmonic decomposition; updated routines for generating window functions, detrending timeseries, and filling gaps; and consideration of physical effects such as mode profile asymmetry, horizontal displacement at the solar surface, and distortion of eigenfunctions by differential rotation. We apply these changes one by one to three years of data, and then reanalyze the entire MDI mission applying all of them, using both the original 72-day long timeseries and 360-day long timeseries. We find significant changes in mode parameters, both as a result of the various changes to the processing, as well as between the 72-day and 360-day analyses. We find reduced residuals of inversions for internal rotation, but seeming artifacts remain, such as the peak in the rotation rate near the surface at high latitudes. An annual periodicity in the f-mode frequencies is also investigated.

A new fitting methodology is presented which is equally well suited for the estimation of low-, medium-, and high-degree mode parameters from $m$-averaged solar oscillation power spectra of widely differing spectral resolution. This method, which we call the "Windowed, MuLTiple-Peak, averaged spectrum", or WMLTP Method, constructs a theoretical profile by convolving the weighted sum of the profiles of the modes appearing in the fitting box with the power spectrum of the window function of the observing run using weights from a leakage matrix that takes into account both observational and physical effects, such as the distortion of modes by solar latitudinal differential rotation. We demonstrate that the WMLTP Method makes substantial improvements in the inferences of the properties of the solar oscillations in comparison with a previous method that employed a single profile to represent each spectral peak. We also present an inversion for the internal solar structure which is based upon 6,366 modes that we have computed using the WMLTP method on the 66-day long 2010 SOHO/MDI Dynamics Run. To improve both the numerical stability and reliability of the inversion we developed a new procedure for the identification and correction of outliers in a frequency data set. We present evidence for a pronounced departure of the sound speed in the outer half of the solar convection zone and in the subsurface shear layer from the radial sound speed profile contained in Model~S of Christensen-Dalsgaard and his collaborators that existed in the rising phase of Solar Cycle~24 during mid-2010.

We present the first accurate characterization of high-degree modes, derived using the best MDI full-disk full-resolution data set available. A ninety day long time series of full-disk two arc-second per pixel resolution dopplergrams was acquired in 2001. These dopplergrams were spatially decomposed using our best estimate of the image scale and the known components of MDI's image distortion. A multi-taper power spectrum estimator was used to generate power spectra up to l = 1000, with a large number of tapers to reduce the realization noise, while the blending at high degrees negates the need for high spectral resolution. These power spectra were fitted for all degrees and all azimuthal orders, between l = 100 and l = 1000. This fitting generated in excess of 6x10^6 individual estimates of ridge frequencies, line-widths, amplitudes and asymmetries (singlets: l,n,m), corresponding to some 6,000 multiplets (l,n). Fitting at high degrees generates ridge characteristics, characteristics that do not correspond to the underlying mode characteristics. We used a sophisticated forward modeling to recover the best possible estimate of the underlying mode characteristics (mode frequencies, as well as line-widths, amplitudes and asymmetries). We describe in detail this modeling and its validation. The modeling has been extensively reviewed and refined, by including an iterative process to improve its input parameters to better match the observations. Also, the contribution of the leakage matrix on the accuracy of the procedure has been carefully assessed. We present the derived set of corrected mode characteristics, discuss their uncertainties and the precision of the ridge to mode correction schemes. In our conclusions, we address how to further improve these estimates, and the implications for other data sets, like GONG+ and HMI.

Normal modes of oscillation of the Sun are useful probes of the solar interior. In this work, we use the even-order splitting coefficients to study the evolution of magnetic fields in the convection zone over solar cycle 23, assuming that the frequency splitting is only due to rotation and a large scale magnetic field. We find that the data are best fit by a combination of a poloidal field and a double-peaked near-surface toroidal field. The toroidal fields are centered at r=0.999R_solar and r=0.996R_solar and are confined to the near-surface layers. The poloidal field is a dipole field. The peak strength of the poloidal field is 124 +/- 17G. The toroidal field peaks at 380 +/- 30G and 1.4 +/- 0.2kG for the shallower and deeper fields respectively. The field strengths are highly correlated with surface activity. The toroidal field strength shows a hysteresis-like effect when compared to the global 10.7 cm radio flux. The poloidal field strength shows evidence of saturation at high activity.

We study the effect of localized sound-speed perturbations on global mode frequencies by applying techniques of global helioseismology on numerical simulations of the solar acoustic wave field. Extending the method of realization noise subtraction (e.g. Hanasoge et al. 2007) to global modes and exploiting the luxury of full spherical coverage, we are able to achieve very highly resolved frequency differences that are used to study sensitivities and the signatures of the thermal asphericities. We find that (1) global modes are almost twice as sensitive to sound-speed perturbations at the bottom of the convection zone as in comparison to anomalies well in the radiative interior ($r\lesssim0.55 R_\odot$), (2) the $m$-degeneracy is lifted ever so slightly, as seen in the $a$ coefficients, and (3) modes that propagate in the vicinity of the perturbations show small amplitude shifts ($\sim 0.5%$).