Alexander G. Kosovichev

Helioseismic Studies of Differential Rotation in the Solar Envelope by the Solar Oscillations Investigation Using the Michelson Doppler Imager

The splitting of the frequencies of the global resonant acoustic modes of the Sun by large-scale flows and rotation permits study of the variation of angular velocity Ω with both radius and latitude within the turbulent convection zone and the deeper radiative interior. The nearly uninterrupted Doppler imaging observations, provided by the Solar Oscillations Investigation (SOI) using the Michelson Doppler Imager (MDI) on the Solar and Heliospheric Observatory (SOHO) spacecraft positioned at the L1 Lagrangian point in continuous sunlight, yield oscillation power spectra with very high signal-to-noise ratios that allow frequency splittings to be determined with exceptional accuracy. This paper reports on joint helioseismic analyses of solar rotation in the convection zone and in the outer part of the radiative core. Inversions have been obtained for a medium-l mode set (involving modes of angular degree l extending to about 250) obtained from the first 144 day interval of SOI-MDI observations in 1996. Drawing inferences about the solar internal rotation from the splitting data is a subtle process. By applying more than one inversion technique to the data, we get some indication of what are the more robust and less robust features of our inversion solutions. Here we have used seven different inversion methods. To test the reliability and sensitivity of these methods, we have performed a set of controlled experiments utilizing artificial data. This gives us some confidence in the inferences we can draw from the real solar data. The inversions of SOI-MDI data have confirmed that the decrease of Ω with latitude seen at the surface extends with little radial variation through much of the convection zone, at the base of which is an adjustment layer, called the tachocline, leading to nearly uniform rotation deeper in the radiative interior. A prominent rotational shearing layer in which Ω increases just below the surface is discernible at low to mid latitudes. Using the new data, we have also been able to study the solar rotation closer to the poles than has been achieved in previous investigations. The data have revealed that the angular velocity is distinctly lower at high latitudes than the values previously extrapolated from measurements at lower latitudes based on surface Doppler observations and helioseismology. Furthermore, we have found some evidence near latitudes of 75° of a submerged polar jet which is rotating more rapidly than its immediate surroundings. Superposed on the relatively smooth latitudinal variation in Ω are alternating zonal bands of slightly faster and slower rotation, each extending some 10° to 15° in latitude. These relatively weak banded flows have been followed by inversion to a depth of about 5% of the solar radius and appear to coincide with the evolving pattern of torsional oscillations reported from earlier surface Doppler studies.

The current state of solar modeling

Data from the Global Oscillation Network Group (GONG) project and other helioseismic experiments provide a test for models of stellar interiors and for the thermodynamic and radiative properties, on which the models depend, of matter under the extreme conditions found in the sun. Current models are in agreement with the helioseismic inferences, which suggests, for example, that the disagreement between the predicted and observed fluxes of neutrinos from the sun is not caused by errors in the models. However, the GONG data reveal subtle errors in the models, such as an excess in sound speed just beneath the convection zone. These discrepancies indicate effects that have so far not been correctly accounted for; for example, it is plausible that the sound-speed differences reflect weak mixing in stellar interiors, of potential importance to the overall evolution of stars and ultimately to estimates of the age of the galaxy based on stellar evolution calculations.

Torsional Oscillation, Meridional Flows, and Vorticity Inferred in the Upper Convection Zone of the Sun by Time-Distance Helioseismology

By applying time-distance helioseismology measurements and inversions to Solar and Heliospheric Observatory (SOHO) Michelson Doppler Imager (MDI) dynamics campaign data, we obtain synoptic maps of subsurface plasma-flow fields at a depth of 0-12 Mm for seven solar Carrington rotations, covering the years 1996-2002, from solar-activity minimum to maximum. Vorticity distribution and both zonal and meridional flows are derived from such synoptic flow maps, which contain an enormous amount of information about solar dynamics. The results for the zonal flows agree well with previous results. The meridional flows of an order of 20 m s-1 are found to remain poleward during the whole period of observations. In addition to the poleward meridional flows observed at the solar minimum, extra meridional circulation cells of flow converging toward the activity belts are found in both hemispheres, which may imply plasma downdrafts in the activity belts. These converging flow cells migrate toward the solar equator together with the activity belts as the solar cycle evolves. The vorticity distributions are largely linear with latitude, and the deviations from the vorticity caused by the mean differential rotation are presented. Patterns of large-scale flows are investigated for a large active region at different depths. Converging flows toward the center of the active region are found near the solar surface, and divergent flows in this large active region are found to be rooted much deeper than similar flows observed in individual sunspots. We conclude that the extremely rich and complicated dynamics of the upper convection zone reveal remarkable organization on the large scale, which can be correlated with the magnetic activity zones.

Differential rotation and dynamics of the solar interior

Splitting of the suns global oscillation frequencies by large-scale flows can be used to investigate how rotation varies with radius and latitude within the solar interior. The nearly uninterrupted observations by the Global Oscillation Network Group (GONG) yield oscillation power spectra with high duty cycles and high signal-to-noise ratios. Frequency splittings derived from GONG observations confirm that the variation of rotation rate with latitude seen at the surface carries through much of the convection zone, at the base of which is an adjustment layer leading to latitudinally independent rotation at greater depths. A distinctive shear layer just below the surface is discernible at low to mid-latitudes.

Investigation of Mass Flows beneath a Sunspot by Time-Distance Helioseismology

A time-distance helioseismic technique is employed to analyze a set of high-resolution Dopplergram observations of a large sunspot by SOHO/MDI on 1998 June 18. A regularized, damped least-squares inversion is applied to the measurements of travel times to infer mass flows around the sunspot below the solar surface. Powerful converging and downward directed flows are detected at depths of 1.5-5 Mm, which may provide observational evidence for the downdrafts and vortex flows that were suggested by Parker for a cluster model of sunspots. Strong outflows extending more than 30 Mm are found below the downward and converging flows. It is suggested that the sunspot might be a relatively shallow phenomenon, with a depth of 5-6 Mm, as defined by its thermal and hydrodynamic properties. A strong mass flow across the sunspot is found at depths of 9-12 Mm, which may provide more evidence in support of the cluster model, as opposed to the monolithic sunspot model. We suggest that a new magnetic emergence that was found 5 hr after our analysis period is related to this mass flow.

The Seismic Structure of the Sun

Global Oscillation Network Group data reveal that the internal structure of the sun can be well represented by a calibrated standard model. However, immediately beneath the convection zone and at the edge of the energy-generating core, the sound-speed variation is somewhat smoother in the sun than it is in the model. This could be a consequence of chemical inhomogeneity that is too severe in the model, perhaps owing to inaccurate modeling of gravitational settling or to neglected macroscopic motion that may be present in the sun. Accurate knowledge of the suns structure enables inferences to be made about the physics that controls the sun; for example, through the opacity, the equation of state, or wave motion. Those inferences can then be used elsewhere in astrophysics.

MAGNETIC ENERGY RELEASE AND TRANSIENTS IN THE SOLAR FLARE OF 2000 JULY 14

High-resolution observations of a large solar flare on 2000 July 14 (“Bastille Day Flare”) from the Michelson Doppler Imager instrument on the SOHO spacecraft reveal rapid variations of the magnetic field in the lower solar atmosphere during the flare. Some of these variations were irreversible, occurred in the vicinity of magnetic neutral lines, and likely were related to magnetic energy release in the flare. A surprising result is that these variations happened very rapidly on the scale of 10‐15 minutes in a large area of »50 Mm 2 at the beginning of the flare. Other, more localized and impulsive magnetic field variations somewhat similar to “magnetic transients” observed by Zirin and coworkers were accompanied by impulses in continuum intensity and Doppler velocity. These impulses have dynamic characteristics similar to Ellerman’s “bombs” and Severny’s “mustaches” and were probably caused by high-energy particles bombarding the solar surface.

A PRECISE ASTEROSEISMIC AGE AND RADIUS FOR THE EVOLVED SUN-LIKE STAR KIC 11026764

The primary science goal of the Kepler Mission is to provide a census of exoplanets in the solar neighborhood, including the identification and characterization of habitable Earth-like planets. The asteroseismic capabilities of the mission are being used to determine precise radii and ages for the target stars from their solar-like oscillations. Chaplin et al. published observations of three bright G-type stars, which were monitored during the first 33.5 days of science operations. One of these stars, the subgiant KIC 11026764, exhibits a characteristic pattern of oscillation frequencies suggesting that it has evolved significantly. We have derived asteroseismic estimates of the properties of KIC 11026764 from Kepler photometry combined with ground-based spectroscopic data. We present the results of detailed modeling for this star, employing a variety of independent codes and analyses that attempt to match the asteroseismic and spectroscopic constraints simultaneously. We determine both the radius and the age of KIC 11026764 with a precision near 1%, and an accuracy near 2% for the radius and 15% for the age. Continued observations of this star promise to reveal additional oscillation frequencies that will further improve the determination of its fundamental properties.

Helioseismic Constraints on the Gradient of Angular Velocity at the Base of the Solar Convection Zone

The layer of transition from the nearly rigid rotation of the radiative interior to the latitudinal differential rotation of the convection zone plays a significant role in the internal dynamics of the Sun. Using rotational splitting coefficients of the p-mode frequencies, obtained during 1986-1990 at the Big Bear Solar Observatory, we have found that the thickness of the transitional layer is 0.09 ± 0.04 solar radii (63 ± 28 Mm), and that most of the transition occurs beneath the adiabatically stratified part of the convection zone, as suggested by the dynamo theories of the 22 yr solar activity cycle.

STRUCTURE AND ROTATION OF THE SOLAR INTERIOR: INITIAL RESULTS FROM THE MDI MEDIUM-L PROGRAM

The medium-l program of the Michelson Doppler Imager instrument on board SOHO provides continuous observations of oscillation modes of angular degree, l, from 0 to ∼ 300. The data for the program are partly processed on board because only about 3% of MDI observations can be transmitted continuously to the ground. The on-board data processing, the main component of which is Gaussian-weighted binning, has been optimized to reduce the negative influence of spatial aliasing of the high-degree oscillation modes. The data processing is completed in a data analysis pipeline at the SOI Stanford Support Center to determine the mean multiplet frequencies and splitting coefficients.