Search NASASearch

Engineering topics

Duvall, Thomas L., Jr.

Publications and source records attributed to Duvall, Thomas L., Jr..

Subsurface Supergranular Vertical Flows as Measured Using Large Distance Separations in Time-Distance Helioseismology

As large-distance rays (say, 10 - 24deg) approach the solar surface approximately vertically, travel times measured from surface pairs for these large separations are mostly sensitive to vertical flows, at least for shallow flows within a few Mm of the solar surface. All previous analyses of supergranulation have used smaller separations and have been hampered by the difficulty of separating the horizontal and vertical flow components. We find that the large-separation travel times associated with supergranulation cannot be studied using the standard phase-speed filters of time-distance helioseismology. These filters, whose use is based upon a refractive model of the perturbations, reduce the resultant travel time signal by at least an order of magnitude at some distances. More effective filters are derived. Modeling suggests that the center-annulus travel-time difference [outward-going time minus inward-going time] in the separation range delta= 10 - 24deg is insensitive to the horizontally diverging flow from the centers of the supergranules and should lead to a constant signal from the vertical flow. Our measurement of this quantity, 5.1+/-0.1 seconds, is constant over the distance range. This magnitude of the signal cannot be caused by the level of upflow at cell centers seen at the photosphere of 10 ms(exp−1) extended in depth. It requires the vertical flow to increase with depth. A simple Gaussian model of the increase with depth implies a peak upward flow of 240 ms(exp−1) at a depth of 2.3 Mm and a peak horizontal flow of 700 ms(exp−1) at a depth of 1.6 Mm.

convection

The Rapidly Rotating Sun

Convection in the solar interior is thought to comprise structures at a continuum of scales, from large to small. This conclusion emerges from phenomenological studies and numerical simulations though neither covers the proper range of dynamical parameters of solar convection. In the present work, imaging techniques of time-distance helioseismology applied to observational data reveal no long-range order in the convective motion. We conservatively bound the associated velocity magnitudes, as a function of depth and the spherical-harmonic degree l to be 20-100 times weaker than prevailing estimates within the wavenumber band l < 60. The observationally constrained kinetic energy is approximately a thousandth of the theoretical prediction, suggesting the prevalence of an intrinsically different paradigm of turbulence. A fundamental question arises: what mechanism of turbulence transports the heat ux of a solar luminosity outwards? The Sun is seemingly a much faster rotator than previously thought, with advection dominated by Coriolis forces at scales l < 60.

Hanasoge, Shravan M.

Time-Distance Helioseismology Data-Analysis Pipeline for Helioseismic and Magnetic Imager Onboard Solar Dynamics Observatory (SDO-HMI) and Its Initial Results

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.

Sun: helioseismology

Time-Distance Helioseismology

Time-distance helioseismology is a method of ambient noise imaging using the solar oscillations. The basic realization that led to time-distance helioseismology was that the temporal cross correlation of the signals at two 'surface' (or photospheric) locations should show a feature at the time lag corresponding to the subsurface travel time between the locations. The temporal cross correlation, as a function of the location separation, is the Fourier transform of the spatio-temporal power spectrum of the solar oscillations, a commonly used function in helioseismology. It is therefore likely the characteristic ridge structure of the correlation function had been seen before without appreciation of its significance. Travel times are measured from the cross correlations. The times are sensitive to a number of important subsurface solar phenomena. These include sound speed variations, flows, and magnetic fields. There has been much interesting progress in the 17 years since the first paper on this subject (Duvall et al., Nature, 1993, 362, 430-432). This progress will be reviewed in this paper.

Duvall, Thomas L., Jr.

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory-Helioseismic and Magnetic Imager Time-Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a time - distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time - distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference traveltime perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.

HMI

Time-Distance Helioseismology with the HMI Instrument

We expect considerable improvement of time-distance results from the Helioseismic and Magnetic Imager (HMI) instrument as opposed to the earlier MDI and GONG data. The higher data rate makes possible several improvements, including faster temporal sampling (45 sec), smaller spatial pixels (0.5 arc sec), better wavelength coverage (6 samples across the line all transmitted to the ground), and year-round coverage of the full disk. The higher spatial resolution makes possible better longitude coverage of active regions and supergranulation and also better latitude coverage. Doppler, continuum, and line depth images have a strong granulation signal. Line core images show little granulation. Analyses to test the limits of these new capabilities will be presented.

Duvall, Thomas L., Jr.

A New Component of Solar Dynamics: North-South Diverging Flows Migrating toward the Equator with an 11 Year Period

Time-distance helioseismology analysis of dopplergrams provides maps of torsional oscillations and meridional flows. Meridional flow maps show a time-varying component that has a banded structure which matches the torsional oscillations with an equatorward migration over the solar cycle. The time-varying component of meridional flow consists of a flow diverging from the dominant latitude of magnetic activity. These maps are compared with other torsional oscillation maps and with magnetic flux maps, showing a strong correlation with active latitudes. These results demonstrate a strong link between the time-varying component of the meridional flow and the torsional oscillations.

Beck, J. G.

Advances in Time-Distance Helioseismology

Time-distance helioseismology is a way to measure travel times between surface locations for waves traversing the solar interior. Coupling the travel with an extensive modeling effort has proven to be a powerful tool for measuring flows and other wave speed inhomogeneities in the solar interior. Problems receiving current attention include studying the time variation of the meridional circulation and torsional oscillation and active region emergence and evolution, current results on these topics will be presented.

Duvall, Thomas L., Jr.

Time-Distance Helioseismology with f Modes as a Method for Measurement of Near-Surface Flows

Travel times measured for the f mode have been used to study flows near the solar surface in conjunction with simultaneous measurements of the magnetic field. Previous flow measurements of doppler surface rotation, small magnetic feature rotation, supergranular pattern rotation, and surface meridional circulation have been confirmed. In addition, the flow in supergranules due to Coriolis forces has been measured. The spatial and temporal power spectra for a six-day observing sequence has been measured.

Duvall, Thomas L., Jr.

Plane-wave analysis of solar acoustic-gravity waves: A (slightly) new approach

The plane-wave decomposition of the acoustic-gravity wave effects observed in the photosphere provides a computationally efficient technique that probes the structure of the upper convective zone and boundary. In this region, the flat sun approximation is considered as being reasonably accurate. A technique to be used for the systematic plane-wave analysis of Michelson Doppler imager data, as part of the solar oscillations investigation, is described. Estimates of sensitivity are presented, and the effects of using different planar mappings are discussed. The technique is compared with previous approaches to the three dimensional plane-wave problem.

Bogart, Richard S.

Other groundbased observations

Ground-based solar observational programs are discussed. The Birmingham (U.K.) solar oscillation network (BiSON) program, the Crimean (Ukraine) low degree program, the Stanford (California) low degree program, the Tenerife (Spain) group program, the South Pole program, the Mount Wilson (California) program, the international research on the interior of the sun (IRIS) network program, the high-l helioseismometer instrument, the ground luminosity oscillations imager (LOI) program, and the Taiwanese oscillation network (TON) program, are outlined.

Duvall, Thomas L., Jr.

Time-distance helioseismology in the vicinity of sunspots

We use the ray description of acoustic-gravity modes to calculate time-distance diagrams for the quiet Sun and for regions in the vicinity of a sunspot with a monolithic flux-tube structure. Time-distance curves for the quiet Sun match the observations of Duvall et al. In the vicinity of a sunspot these quiet Sun curves split into a family of closely spaced curves. The structure of this bandlike feature is found to be sensitive to the sunspot model and can be a diagnostic of the subsurface geometry of the sunspot flux tube.

D'Silva, Sydney

The NASA/NSO spectromagnetograph

The Kitt Peak Vacuum Telescope's spectromagnetograph features real-time digital analysis of long-slit spectra formed on a 2D CCD detector, and is located at an exit port of a Littrow spectrograph. Commercial video processing boards are used to digitize the spectral images at video rates and to separate, accumulate, and buffer the spectra in the two polarization states. A processor removes fixed-pattern bias and gain from the spectra in cadence with spatial scanning of the image across the entrance slit; a data-control computer conducts position and width analysis of the line profiles and records line-of-sight magnetic field, Doppler shift, and other computed parameters, as demonstrated by the early observations presented.

Jones, Harrison P.

Modeling of solar oscillation power spectra

To produce accurate estimates of the line-profile parameters of a model used to represent the spectral features in a solar oscillation power spectrum, it is necessary to (1) select the appropriate probability density function when deriving the maximum-likelihood function to be employed for the parameter estimation and (2) allow for the redistribution of spectral power caused by gaps in the data string. This paper describes a maximum-likelihood method for estimating the model parameters (based on the observed power spectrum statistics) that accounts for redistribution of spectral power caused by gaps in the data string, by convolving the model with the power spectrum of the observed window function. The accuracy and reliability of the method were tested using both artificial and authentic solar oscillation power spectrum data. A comparison of this method with various least-squares techniques is also presented.

Anderson, Edwin R.

A test of a modified algorithm for computing spherical harmonic coefficients using an FFT

The Dilts (1985) algorithm for computing the spherical harmonic expansion coefficients for a function on a sphere, on the basis of a two-dimensional FFT, is presently modified, tested, and found to eliminate problems of overflow and large storage requirements associated with the encounter of harmonic degree values greater than 16. Results from timing tests show the Dilts program to be impractical, however, for the computation of spherical harmonic expansion coefficients for large harmonic degree values.

Elowitz, Mark