Validated helioseismic inversions for 3-D vector flows

According to time-distance helioseismology, information about internal fluid motions is encoded in the travel times of solar waves. The inverse problem consists of inferring 3-D vector flows from a set of travel-time measurements. Here we investigate the potential of time-distance helioseismology to infer 3-D convective velocities in the near-surface layers of the Sun. We developed a new Subtractive Optimally Localised Averaging (SOLA) code suitable for pipeline pseudo-automatic processing. Compared to its predecessor, the code was improved by accounting for additional constraints in order to get the right answer within a given noise level. The main aim of this study is to validate results obtained by our inversion code. We simulate travel-time maps using a snapshot from a numerical simulation of solar convective flows, realistic Born travel-time sensitivity kernels, and a realistic model of travel-time noise. These synthetic travel times are inverted for flows and the results compared with the known input flow field. Additional constraints are implemented in the inversion: cross-talk minimization between flow components and spatial localization of inversion coefficients. Using modes f, p1 through p4, we show that horizontal convective flow velocities can be inferred without bias, at a signal-to-noise ratio greater than one in the top 3.5 Mm, provided that observations span at least four days. The vertical component of velocity (v_z), if it were to be weak, is more difficult to infer and is seriously affected by cross-talk from horizontal velocity components. We emphasise that this cross-talk must be explicitly minimised in order to retrieve v_z in the top 1 Mm. We also show that statistical averaging over many different areas of the Sun allows for reliably measuring of average properties of all three flow components in the top 5.5 Mm of the convection zone.Comment: 14 pages main paper, 9 pages electronic supplement, 28 figures. Accepted for publication in Astronomy & Astrophysic

Probing sunspots with two-skip time-distance helioseismology

Previous helioseismology of sunspots has been sensitive to both the structural and magnetic aspects of sunspot structure. We aim to develop a technique that is insensitive to the magnetic component so the two aspects can be more readily separated. We study waves reflected almost vertically from the underside of a sunspot. Time-distance helioseismology was used to measure travel times for the waves. Ray theory and a detailed sunspot model were used to calculate travel times for comparison. It is shown that these large distance waves are insensitive to the magnetic field in the sunspot. The largest travel time differences for any solar phenomena are observed. With sufficient modeling effort, these should lead to better understanding of sunspot structure

Balltracking: an highly efficient method for tracking flow fields

We present a method for tracking solar photospheric flows that is highly efficient, and demonstrate it using high resolution MDI continuum images. The method involves making a surface from the photospheric granulation data, and allowing many small floating tracers or balls to be moved around by the evolving granulation pattern. The results are tested against synthesised granulation with known flow fields and compared to the results produced by Local Correlation tracking (LCT). The results from this new method have similar accuracy to those produced by LCT. We also investigate the maximum spatial and temporal resolution of the velocity field that it is possible to extract, based on the statistical properties of the granulation data. We conclude that both methods produce results that are close to the maximum resolution possible from granulation data. The code runs very significantly faster than our similarly optimised LCT code, making real time applications on large data sets possible. The tracking method is not limited to photospheric flows, and will also work on any velocity field where there are visible moving features of known scale length

Greater Forearm Blood Flow is Associated With Higher Physical Activity in Older Individuals

Please refer to the pdf version of the abstract located adjacent to the title

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.Comment: Accepted by Solar Physics topical issue 'Solar Dynamics Observatory

Scattering of the f-mode by small magnetic flux elements from observations and numerical simulations

The scattering of f-modes by magnetic tubes is analyzed using three-dimensional numerical simulations. An f-mode wave packet is propagated through a solar atmosphere embedded with three different flux tube models which differ in radius and total magnetic flux. A quiet Sun simulation without a tube present is also performed as a reference. Waves are excited inside the flux tube and propagate along the field lines, and jacket modes are generated in the surroundings of the flux tube, carrying 40% as much energy as the tube modes. The resulting scattered wave is mainly an f-mode composed of a mixture of m=0 and m=+/-1 modes. The amplitude of the scattered wave approximately scales with the magnetic flux. A small amount of power is scattered into the p_1-mode. We have evaluated the absorption and phase shift from a Fourier-Hankel decomposition of the photospheric vertical velocities. They are compared with the results obtained from the emsemble average of 3400 small magnetic elements observed in high-resolution MDI Doppler datacubes. The comparison shows that the observed dependence of the phase shift with wavenumber can be matched reasonably well with the simulated flux tube model. The observed variation of the phase-shifts with the azimuthal order $m$ appears to depend on details of the ensemble averaging, including possible motions of the magnetic elements and asymmetrically shaped elements.Comment: Accepted for publication in The Astrophysical Journa

Phylogeny of the Grasses (Poaceae) Revisited

The most robust previously published phylogeny for the overall structure of the grass family (Poaceae) shows three early diverging lineages and two major derived clades, the BEP clade and the PACCAD clade (Grass Phylogeny Working Group 2001). A few key taxa were incompletely sampled, however, and support for the BEP clade was moderate at best and relationships among the major lineages within the PACCAD clade remained unresolved. In addition, recent studies indicated that the sister group to Poaceae may be Joinvilleaceae and/or Ecdeiocoleaceae, the latter of which were not previously sampled. In this study, missing structural data were determined and analyzed as well as sequence data for ndhF and rbcL, the two most complete plastid sequence data sets. Sampling was increased with a particular focus on key taxa such as Danthoniopsis, Eriachne, Micraira, and Streptogyna and a representative of the outgroup, Ecdeiocoleaceae. A total of 61 ingroup and two outgroup taxa were analyzed using maximum parsimony for total data, and maximum parsimony, Bayesian inference, and neighbor joining for the molecular data. A strongly supported clade of ((Eriachneae, Isachne) Micraira) was recovered as a sister subfamily to Arundinoideae and excluded from Panicoideae. Arundinaria was strongly united with Bambusoideae. The position of Streptogyna was weakly supported among Ehrhartoideae, and is still unresolved. An outgroup effect on ingroup topology was observed, demonstrating that highly divergent outgroups may unpredictably alter ingroup relationships

Test-Anchored Vibration Response Predictions for an Acoustically Energized Curved Orthogrid Panel with Mounted Components

A rich body of vibroacoustic test data was recently generated at Marshall Space Flight Center for a curved orthogrid panel typical of launch vehicle skin structures. Several test article configurations were produced by adding component equipment of differing weights to the flight-like vehicle panel. The test data were used to anchor computational predictions of a variety of spatially distributed responses including acceleration, strain and component interface force. Transfer functions relating the responses to the input pressure field were generated from finite element based modal solutions and test-derived damping estimates. A diffuse acoustic field model was employed to describe the assumed correlation of phased input sound pressures across the energized panel. This application demonstrates the ability to quickly and accurately predict a variety of responses to acoustically energized skin panels with mounted components. Favorable comparisons between the measured and predicted responses were established. The validated models were used to examine vibration response sensitivities to relevant modeling parameters such as pressure patch density, mesh density, weight of the mounted component and model form. Convergence metrics include spectral densities and cumulative root-mean squared (RMS) functions for acceleration, velocity, displacement, strain and interface force. Minimum frequencies for response convergence were established as well as recommendations for modeling techniques, particularly in the early stages of a component design when accurate structural vibration requirements are needed relatively quickly. The results were compared with long-established guidelines for modeling accuracy of component-loaded panels. A theoretical basis for the Response/Pressure Transfer Function (RPTF) approach provides insight into trends observed in the response predictions and confirmed in the test data. The software modules developed for the RPTF method can be easily adapted for quick replacement of the diffuse acoustic field with other pressure field models; for example a turbulent boundary layer (TBL) model suitable for vehicle ascent. Wind tunnel tests have been proposed to anchor the predictions and provide new insight into modeling approaches for this type of environment. Finally, component vibration environments for design were developed from the measured and predicted responses and compared with those derived from traditional techniques such as Barrett scaling methods for unloaded and component-loaded panels

Helioseismology of Sunspots: A Case Study of NOAA Region 9787

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.Comment: 28 pages, 18 figure