Implementation of Helioseismic Data Reduction and Diagnostic Techniques on Massively Parallel Architectures

Under the direction of Dr. Rhodes, and the technical supervision of Dr. Korzennik, the data assimilation of high spatial resolution solar dopplergrams has been carried out throughout the program on the Intel Delta Touchstone supercomputer. With the help of a research assistant, partially supported by this grant, and under the supervision of Dr. Korzennik, code development was carried out at SAO, using various available resources. To ensure cross-platform portability, PVM was selected as the message passing library. A parallel implementation of power spectra computation for helioseismology data reduction, using PVM was successfully completed. It was successfully ported to SMP architectures (i.e. SUN), and to some MPP architectures (i.e. the CM5). Due to limitation of the implementation of PVM on the Cray T3D, the port to that architecture was not completed at the time

Seismic Study of the Subsurface Structure and Dynamics of the Solar Interior from High Spatial Resolution Observations

We have carried out the data reduction and analysis of Mt. Wilson 60' solar tower high spatial resolution observations. The reduction of the 100-day-long summer of 1990 observation campaign in terms of rotational splittings was completed leading to an excess of 600,000 splittings. The analysis of these splittings lead to a new inference of the solar internal rotation rate as a function of depth and latitude

Sensitivity analysis of the solar rotation to helioseismic data from GONG, GOLF and MDI observations

Accurate determination of the rotation rate in the radiative zone of the sun from helioseismic observations requires rotational frequency splittings of exceptional quality as well as reliable inversion techniques. We present here inferences based on mode parameters calculated from 2088-days long MDI, GONG and GOLF time series that were fitted to estimate very low frequency rotational splittings (nu < 1.7 mHz). These low frequency modes provide data of exceptional quality, since the width of the mode peaks is much smaller than the rotational splitting and hence it is much easier to separate the rotational splittings from the effects caused by the finite lifetime and the stochastic excitation of the modes. We also have implemented a new inversion methodology that allows us to infer the rotation rate of the radiative interior from mode sets that span l=1 to 25. Our results are compatible with the sun rotating like a rigid solid in most of the radiative zone and slowing down in the core (R_sun < 0.2). A resolution analysis of the inversion was carried out for the solar rotation inverse problem. This analysis effectively establishes a direct relationship between the mode set included in the inversion and the sensitivity and information content of the resulting inferences. We show that such an approach allows us to determine the effect of adding low frequency and low degree p-modes, high frequency and low degree p-modes, as well as some g-modes on the derived rotation rate in the solar radiative zone, and in particular the solar core. We conclude that the level of uncertainties that is needed to infer the dynamical conditions in the core when only p-modes are included is unlikely to be reached in the near future, and hence sustained efforts are needed towards the detection and characterization of g-modes.Comment: Accepted for publication in Astrophysical journal. 15 pages, 19 figure

On The Determination of MDI High-Degree Mode Frequencies

The characteristic of the solar acoustic spectrum is such that mode lifetimes get shorter and spatial leaks get closer in frequency as the degree of a mode increases for a given order. A direct consequence of this property is that individual p-modes are only resolved at low and intermediate degrees, and that at high degrees, individual modes blend into ridges. Once modes have blended into ridges, the power distribution of the ridge defines the ridge central frequency and it will mask the true underlying mode frequency. An accurate model of the amplitude of the peaks that contribute to the ridge power distribution is needed to recover the underlying mode frequency from fitting the ridge. We present the results of fitting high degree power ridges (up to l = 900) computed from several two to three-month-long time-series of full-disk observations taken with the Michelson Doppler Imager (MDI) on-board the Solar and Heliospheric Observatory between 1996 and 1999. We also present a detailed discussion of the modeling of the ridge power distribution, and the contribution of the various observational and instrumental effects on the spatial leakage, in the context of the MDI instrument. We have constructed a physically motivated model (rather than some ad hoc correction scheme) resulting in a methodology that can produce an unbiased determination of high-degree modes, once the instrumental characteristics are well understood. Finally, we present changes in high degree mode parameters with epoch and thus solar activity level and discuss their significance.Comment: 59 pages, 38 figures -- High-resolution version at http://www-sgk.harvard.edu:1080/~sylvain/preprints/ -- Manuscript submitted to Ap

Precision and systematic errors in global helioseismology mode fitting and inversions: Leveraging some 25 years of nearly uninterrupted observations

We have on hand some 25 years of nearly uninterrupted high-quality and high-cadence global helioseismic data. The Global Oscillations Network Group (GONG) project has been producing science quality data since 1995, the Michelson Doppler Imager (MDI) started in 1996, and the Helioseismic and Magnetic Imager (HMI) took over in 2010. Fundamental new constraints have been imposed by helioseismic inferences, yet global helioseismology data processing seems somewhat frozen in time for some of its methodologies. I review and discuss some specific aspects of global helioseismology data analysis, with an emphasis on the issues and challenges presented by mode fitting and inversion techniques. I compare and contrast results derived by different fitting methods, whether using different techniques, different lengths of time series, or different fitting parameters, like leakage matrices or the inclusion or omission of the mode profile asymmetry, leading to our current best handle on the residual systematic errors

An Upper Limit on the Reflected Light from the Planet Orbiting the Star tau Bootis

The planet orbiting tau Boo at a separation of 0.046 AU could produce a reflected light flux as bright as 1e-4 relative to that of the star. A spectrum of the system will contain a reflected light component which varies in amplitude and Doppler-shift as the planet orbits the star. Assuming the secondary spectrum is primarily the reflected stellar spectrum, we can limit the relative reflected light flux to be less than 5e-5. This implies an upper limit of 0.3 for the planetary geometric albedo near 480 nm, assuming a planetary radius of 1.2 R_Jup. This albedo is significantly less than that of any of the giant planets of the solar system, and is not consistent with certain published theoretical predictions.Comment: 5 pages, 1 figure, accepted by ApJ Letter

The solar differential rotation in the 18th century

The sunspot drawings of Johann Staudacher of 1749--1799 were used to determine the solar differential rotation in that period. These drawings of the full disk lack any indication of their orientation. We used a Bayesian estimator to obtain the position angles of the drawings, the corresponding heliographic spot positions, a time offset between the drawings and the differential rotation parameter \delta\Omega, assuming the equatorial rotation period is the same as today. The drawings are grouped in pairs, and the resulting marginal distributions for \delta\Omega were multiplied. We obtain \delta\Omega=-0.048 \pm 0.025 d^-1 (-2.75^o/d) for the entire period. There is no significant difference to the value of the present Sun. We find an (insignificant) indication for a change of \delta\Omega throughout the observing period from strong differential rotation, \delta\Omega\approx -0.07 d^-1, to weaker differential rotation, \delta\Omega\approx-0.04 d^-1.Comment: 6 pages, 6 figures, accepted for Astronomy and Astrophysic