Context: There is a wide discrepancy in current estimates of the strength of convection flows in the solar interior obtained using different helioseismic methods applied to observations from SDO/HMI. The cause for these disparities is not known. Aims: As one step in the effort to resolve this discrepancy, we aim to characterize the multi-ridge fitting code for ring-diagram helioseismic analysis that is used to obtain flow estimates from local power spectra of solar oscillations. Methods: We updated the multi-ridge fitting code developed by Greer et al.(2014) to solve several problems we identified through our inspection of the code. In particular, we changed the merit function to account for the smoothing of the power spectra, model for the power spectrum, and noise estimates. We used Monte Carlo simulations to generate synthetic data and to characterize the noise and bias of the updated code by fitting these synthetic data. Results: The bias in the output fit parameters, apart from the parameter describing the amplitude of the p-mode resonances in the power spectrum, is below what can be measured from the Monte-Carlo simulations. The amplitude parameters are underestimated; this is a consequence of choosing to fit the logarithm of the averaged power. We defer fixing this problem as it is well understood and not significant for measuring flows in the solar interior. The scatter in the fit parameters from the Monte-Carlo simulations is well-modeled by the formal error estimates from the code. Conclusions: We document and demonstrate a reliable multi-ridge fitting method for ring-diagram analysis. The differences between the updated fitting results and the original results are less than one order of magnitude and therefore we suspect that the changes will not eliminate the aforementioned orders-of-magnitude discrepancy in the amplitude of convective flows in the solar interior.

We developed an experimental platform for studying magnetic reconnection in an external magnetic field with simultaneous measurements of plasma imaging, flow velocity, and magnetic-field variation. Here, we investigate the stagnation and acceleration in counter-streaming plasmas generated by high-power laser beams. A plasma flow perpendicular to the initial flow directions is measured with laser Thomson scattering. The flow is, interestingly, accelerated toward the high-density region, which is opposite to the direction of the acceleration by pressure gradients. This acceleration is possibly interpreted by the interaction of two magnetic field loops initially generated by Biermann battery effect, resulting in a magnetic reconnection forming a single field loop and additional acceleration by a magnetic tension force.

Filamentary enstatite crystals are found in interplanetary dust particles of likely cometary origin but are very rare or absent in meteorites. Crystallographic characteristics of filamentary enstatites indicate that they condensed directly from vapor. We measured the O isotopic composition of an enstatite ribbon from a giant cluster interplanetary dust particle to be $\delta^{18}\rm{O}{=25{\pm}55}$, $\delta^{17}\rm{O}{=-19{\pm}129}$, $\Delta^{17}\rm{O}{=-32{\pm}134}$ (2$\sigma$ errors), which is inconsistent at the 2$\sigma$ level with the composition of the Sun inferred from the Genesis solar wind measurements. The particle's O isotopic composition, consistent with the terrestrial composition, implies that it condensed from a gas of non-solar O isotopic composition, possibly as a result of vaporization of disk region enriched in $^{16}$O-depleted solids. The relative scarcity of filamentary enstatite in asteroids compared to comets implies either that this crystal condensed from dust vaporized \textit{in-situ} in the outer Solar System where comets formed, or it condensed in the inner Solar System and was subsequently transported outward to the comet-forming region.

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.

Context. In time-distance helioseismology, wave travel times are measured from the two-point cross-covariance function of solar oscillations and are used to image the solar convection zone in three dimensions. There is, however, also information in the amplitude of the cross-covariance function, for example about seismic wave attenuation. Aims. Here we develop a convenient procedure to measure the amplitude of the cross-covariance function of solar oscillations. Methods. In this procedure, the amplitude of the cross-covariance function is linearly related to the cross-covariance function and can be measured even for high levels of noise. Results. As an example application, we measure the amplitude perturbations of the seismic waves that propagate through the sunspot in active region NOAA 9787. We can recover the amplitude variations due to the scattering and attenuation of the waves by the sunspot and associated finite-wavelength effects. Conclusions. The proposed definition of cross-covariance amplitude is robust to noise, can be used to relate measured amplitudes to 3D perturbations in the solar interior under the Born approximation, and will provide independent information from the travel times.

Solar flares and coronal mass ejections (CMEs), especially the larger ones, emanate from active regions (ARs). With the aim to understand the magnetic properties that govern such flares and eruptions, we systematically survey all flare events with GOES levels of >=M5.0 within 45 deg from disk center between May 2010 and April 2016. These criteria lead to a total of 51 flares from 29 ARs, for which we analyze the observational data obtained by the Solar Dynamics Observatory. More than 80% of the 29 ARs are found to exhibit delta-sunspots and at least three ARs violate Hale's polarity rule. The flare durations are approximately proportional to the distance between the two flare ribbons, to the total magnetic flux inside the ribbons, and to the ribbon area. From our study, one of the parameters that clearly determine whether a given flare event is CME-eruptive or not is the ribbon area normalized by the sunspot area, which may indicate that the structural relationship between the flaring region and the entire AR controls CME productivity. AR characterization show that even X-class events do not require delta-sunspots or strong-field, high-gradient polarity inversion lines. An investigation of historical observational data suggests the possibility that the largest solar ARs, with magnetic flux of 2x10^23 Mx, might be able to produce "superflares" with energies of order of 10^34 erg. The proportionality between the flare durations and magnetic energies is consistent with stellar flare observations, suggesting a common physical background for solar and stellar flares.

Context: The cross-covariance of solar oscillations observed at pairs of points on the solar surface is a fundamental ingredient in time-distance helioseismology. Wave travel times are extracted from the cross-covariance function and are used to infer the physical conditions in the solar interior. Aims: Understanding the statistics of the two-point cross-covariance function is a necessary step towards optimizing the measurement of travel times. Methods: By modeling stochastic solar oscillations, we evaluate the variance of the cross-covariance function as function of time-lag and distance between the two points. Results: We show that the variance of the cross-covariance is independent of both time-lag and distance in the far field, i.e., when they are large compared to the coherence scales of the solar oscillations. Conclusions: The constant noise level for the cross-covariance means that the signal-to-noise ratio for the cross-covariance is proportional to the amplitude of the expectation value of the cross-covariance. This observation is important for planning data analysis efforts.

Individual particles from comet 81P/Wild 2 collected by NASA's Stardust mission vary in size from small sub-$\mu$m fragments found in the walls of the aerogel tracks, to large fragments up to tens of $\mu$m in size found towards the termini of tracks. The comet, in an orbit beyond Neptune since its formation, retains an intact a record of early-Solar-System processes that was compromised in asteroidal samples by heating and aqueous alteration. We measured the O isotopic composition of seven Stardust fragments larger than $\sim$2 $\mu$m extracted from five different Stardust aerogel tracks, and 63 particles smaller than $\sim$2 $\mu$m from the wall of a Stardust track. The larger particles show a relatively narrow range of O isotopic compositions that is consistent with $^{16}$O-poor phases commonly seen in meteorites. Many of the larger Stardust fragments studied so far have chondrule-like mineralogy which is consistent with formation in the inner Solar System. The fine-grained material shows a very broad range of O isotopic compositions ($-70<\Delta^{17}$O$<+60$) suggesting that Wild 2 fines are either primitive outer-nebula dust or a very diverse sampling of inner Solar System compositional reservoirs that accreted along with a large number of inner-Solar-System rocks to form comet Wild 2.

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe I 6173.3 {\AA} line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multi-height velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams. We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30 - 40 km above (line-center) and 30 - 40 km below (average-wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic cutoff frequency.

The determination of an isotope ratio by secondary ion mass spectrometry (SIMS) traditionally involves averaging a number of ratios collected over the course of a measurement. We show that this method leads to an additive positive bias in the expectation value of the estimated ratio that is approximately equal to the true ratio divided by the counts of the denominator isotope of an individual ratio. This bias does not decrease as the number of ratios used in the average increases. By summing all counts in the numerator isotope, then dividing by the sum of counts in the denominator isotope, the estimated ratio is less biased: the bias is approximately equal to the ratio divided by the summed counts of the denominator isotope over the entire measurement. We propose a third ratio estimator (Beale's estimator) that can be used when the bias from the summed counts is unacceptably large for the hypothesis being tested. We derive expressions for the variance of these ratio estimators as well as the conditions under which they are normally distributed. Finally, we investigate a SIMS dataset showing the effects of ratio bias, and discuss proper ratio estimation for SIMS analysis.

We study Doppler velocity measurements at multiple heights in the solar atmosphere using a set of six filtergrams obtained by the Helioseismic magnetic Imager on board the Solar Dynamics Observatory. There are clear and significant phase differences between core and wing Dopplergrams in the frequency range above the photospheric acoustic cutoff frequency, which indicates that these are really "multi-height" datasets.

We report on a systematic center-to-limb variation in measured helioseismic travel times, which must be taken into account for an accurate determination of solar interior meridional flows. The systematic variation, found in time-distance helioseismology analysis using SDO/HMI and SDO/AIA observations, is different in both travel-time magnitude and variation trend for different observables. It is not clear what causes this systematic effect. Subtracting the longitude-dependent east-west travel times, obtained along the equatorial area, from the latitude-dependent north-south travel times, obtained along the central meridian area, gives remarkably similar results for different observables. We suggest this as an effective procedure for removing the systematic center-to-limb variation. The subsurface meridional flows obtained from inversion of the corrected travel times are approximately 10 m/s slower than those obtained without removing the systematic effect. The detected center-to-limb variation may have important implications in the derivation of meridional flows in the deep interior, and needs a better understanding.

We report the petrology, O isotopic composition, and Al-Mg isotope systematics of a chondrule fragment from the Jupiter-family comet Wild 2, returned to Earth by NASA's Stardust mission. This object shows characteristics of a type II chondrule that formed from an evolved oxygen isotopic reservoir. No evidence for extinct 26Al was found, with (26Al/ 27Al)0 < 3.0 x 10^-6. Assuming homogenous distribution of 26Al in the solar nebula, this particle crystallized at least 3 Myr after the earliest solar system objects--relatively late compared to most chondrules in meteorites. We interpret the presence of this object in a Kuiper Belt body as evidence of late, large-scale transport of small objects between the inner and outer solar nebula. Our observations constrain the formation of Jupiter (a barrier to outward transport if it formed further from the Sun than this cometary chondrule) to be more than 3 Myr after calcium-aluminum-rich inclusionss.

It is believed that Al-26, a short-lived (t1/2 = 0.73 Ma) and now extinct radionuclide, was uniformly distributed in the nascent Solar System with the initial Al-26/Al-27 ratio of ~5.2\times10-5, suggesting its external stellar origin. However, the stellar source of Al-26 and the manner in which it was injected into the solar system remain controversial: the Al-26 could have been produced by an asymptotic giant branch star, a supernova, or a Wolf-Rayet star and injected either into the protosolar molecular cloud or protoplanetary disk. Corundum (Al2O3) is thermodynamically predicted to be the first condensate from a cooling gas of solar composition. Here we show that micron-sized corundum condensates from O-16-rich gas (Big Delta O-17 ~ -25%) of solar composition recorded heterogeneous distribution of Al-26 at the birth of the solar system: the inferred initial Al-26/Al-27 ratio ranges from ~6.5x10-5 to <2x10-6; ~50% of the corundum grains measured are Al-26-poor. Other Al-26-poor, O-16-rich refractory objects include grossite (CaAl4O7)- and hibonite(CaAl12O19)-rich calcium-aluminum-rich inclusions (CAIs) in CH chondrites, platy hibonite crystals in CM chondrites, and FUN (fractionation and unidentified nuclear isotopic anomalies) CAIs in CV, CO, and CR chondrites. Considering the apparently early and short duration (<0.3 Ma) of condensation of refractory O-16-rich solids in the solar system, we infer that Al-26 was injected into the collapsing protosolar molecular cloud and later homogenized in the protoplanetary disk. The apparent lack of correlation between Al-26 abundance and O-isotope compositions of corundum grains put important constraints on the stellar source of Al-26 in the solar system.

We report on a new phenomenon of `alignment' of supergranulation cells in the polar regions of the Sun. Recent high-resolution datasets obtained by the Solar Optical Telescope onboard the Hinode satellite enabled us to investigate supergranular structures in high-latitude regions of the Sun. We have carried out a local helioseismology time-distance analysis of the data, and detected acoustic travel-time variations due to the supergranular flows. The supergranulation cells in both the north and south polar regions show systematic alignment patterns in the north-south direction. The south-pole datasets obtained in a month-long Hinode campaign indicate that the supergranulation alignment property may be quite common in the polar regions. We also discuss the latitudinal dependence of the supergranulation cell sizes; the data show that the east-west cell size decreases towards higher latitudes.

We report on a signature of chromospheric downflows in two emerging-flux regions detected by time-distance helioseismology analysis. We use both chromospheric intensity oscillation data in the Ca II H line and photospheric Dopplergrams in the Fe I 557.6nm line obtained by Hinode/SOT for our analyses. By cross-correlating the Ca II oscillation signals, we have detected a travel-time anomaly in the plage regions; outward travel times are shorter than inward travel times by 0.5-1 minute. However, such an anomaly is absent in the Fe I data. These results can be interpreted as evidence of downflows in the lower chromosphere. The downflow speed is estimated to be below 10 km/s. This result demonstrates a new possibility of studying chromospheric flows by time-distance analysis.

Exploiting high-resolution observations made by the Solar Optical Telescope onboard Hinode, we investigate the spatial distribution of power spectral density of oscillatory signal in and around NOAA active region 10935. The G-band data show that in the umbra the oscillatory power is suppressed in all frequency ranges. On the other hand, in Ca II H intensity maps oscillations in the umbra, so-called umbral flashes, are clearly seen with the power peaking around 5.5 mHz. The Ca II H power distribution shows the enhanced elements with the spatial scale of the umbral flashes over most of the umbra but there is a region with suppressed power at the center of the umbra. The origin and property of this node-like feature remain unexplained.

On 2005 September 13 a filament eruption accompanied by a halo CME occurred in the most flare-productive active region NOAA 10808 in Solar Cycle 23. Using multi-wavelength observations before the filament eruption on Sep. 13th, we investigate the processes leading to the catastrophic eruption. We find that the filament slowly ascended at a speed of 0.1km/s over two days before the eruption. During slow ascending, many small flares were observed close to the footpoints of the filament, where new magnetic elements were emerging. On the basis of the observational facts we discuss the triggering mechanism leading to the filament eruption. We suggest the process toward the eruption as follows: First, a series of small flares played a role in changing the topology of the loops overlying the filament. Second, the small flares gradually changed the equilibrium state of the filament and caused the filament to ascend slowly over two days. Finally, a C2.9 flare that occurred when the filament was close to the critical point for loss of equilibrium directly led to the catastrophic filament eruption right after itself.

We report a statistical study of flares observed with the Soft X-ray Telescope (SXT) onboard Yohkoh in the year of 2000. We measure physical parameters of 77 flares, such as the temporal scale, the size, and the magnetic flux density and find that the sizes of flares tend to be distributed more broadly as the GOES class becomes weaker and that there is a lower limit of magnetic flux density that depends on the GOES class. We also examine the relationship between these parameters and find weak correlation between temporal and spatial scales of flares. We estimate reconnection inflow velocity, coronal Alfven velocity, and reconnection rate using above observed values. The inflow velocities are distributed from a few km/s to several tens km/s and the Alfven velocities in the corona are in the range from 10^3 to 10^4 km/s. Hence the reconnection rate is 10^-3 - 10^-2. We find that the reconnection rate in a flare tends to decrease as the GOES class of the flare increases. This value is within one order of magnitude from the theoretical maximum value predicted by the Petschek model, although the dependence of the reconnection rate on the magnetic Reynolds number tends to be stronger than that in the Petschek model.