Showing 1–50 of 87 results for author: Yokoyama, T

The Chemical Composition of Ryugu: Prospects as a Reference Material for Solar System Composition

Authors: Tetsuya Yokoyama, Nicolas Dauphas, Ryota Fukai, Tomohiro Usui, Shogo Tachibana, Maria Schönbächler, Henner Busemann, Masanao Abe, Toru Yada

Abstract: The Hayabusa 2 spacecraft sampled approximately 5.4 g of asteroid material from the Cb-type asteroid Ryugu. Initial analysis of the Ryugu materials revealed a mineralogical, chemical, and isotopic kinship to the CI chondrites. The pristine nature of Ryugu makes the returned samples ideal for constraining the composition of the Solar System. However, some elements (e.g., P, Ca, Mn, and rare earth e… ▽ More The Hayabusa 2 spacecraft sampled approximately 5.4 g of asteroid material from the Cb-type asteroid Ryugu. Initial analysis of the Ryugu materials revealed a mineralogical, chemical, and isotopic kinship to the CI chondrites. The pristine nature of Ryugu makes the returned samples ideal for constraining the composition of the Solar System. However, some elements (e.g., P, Ca, Mn, and rare earth elements) show large relative dispersions compared to the other elements in the returned materials studied so far, most likely due to the presence of aqueously formed secondary minerals (e.g., carbonates, phosphates) in Ryugu. Therefore, the estimation of the Solar System composition using currently available Ryugu data is challenging due to the so-called nugget effect of carbonates, phosphates, and possibly other accessory minerals. The nugget effect can be mitigated by analyzing a homogenized, relatively large amount of sample. We estimate that for approximately 0.1 g of Ryugu sample, the dispersion (2SD) of the bulk Mn/Cr and Rb/Sr ratios are +/-13% and +/-15%, respectively, while they will be improved to be better than +/-5% for approximately 1 g of homogenized Ryugu sample. To further constrain the Solar System composition and to evaluate if previous estimates based on CI chondrites stored in museums for decades to centuries are reliable, it is strongly recommended to determine the chemical and isotopic compositions of Ryugu using a homogenized sample prepared from relatively large (approx. 1 g) returned material. Determining Ryugu reference compositions will be used by multidisciplinary communities, including Earth and planetary sciences, astronomy, physics, and chemistry. △ Less

Submitted 7 May, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 13 pages, 5 figures, 2 Tables

Comprehensive Synthesis of Magnetic Tornado: Co-spatial Incidence of Chromospheric Swirls and EUV Brightening

Authors: Hidetaka Kuniyoshi, Souvik Bose, Takaaki Yokoyama

Abstract: Magnetic tornadoes, characterized as impulsive Alfven waves initiated by photospheric vortices in intergranular lanes, are considered efficient energy channels to the corona. Despite their acknowledged importance for solar coronal heating, their observational counterparts from the corona have not been well understood. To address this issue, we use a radiative MHD simulation of a coronal loop with… ▽ More Magnetic tornadoes, characterized as impulsive Alfven waves initiated by photospheric vortices in intergranular lanes, are considered efficient energy channels to the corona. Despite their acknowledged importance for solar coronal heating, their observational counterparts from the corona have not been well understood. To address this issue, we use a radiative MHD simulation of a coronal loop with footpoints rooted in the upper convection zone, and synthesize the chromospheric and coronal emissions corresponding to a magnetic tornado. Considering SDO/AIA 171 A and Solar Orbiter/EUI 174 A channels, our synthesis reveals that the coronal response to magnetic tornadoes can be observed as an EUV brightening of which width is ~2 Mm. This brightening is located above the synthesized chromospheric swirl observed in Ca II 8542 A, Ca II K, and Mg II k lines, which can be detected by instruments such as SST/CRISP, GST/FISS, and IRIS. Considering the height correspondence of the synthesized brightening, magnetic tornadoes can be an alternative mechanism for the small-scale EUV brightenings such as the solar "campfires''. Our findings indicate that coordinated observations encompassing the chromosphere to the corona are indispensable for comprehending the origin of coronal EUV brightenings. △ Less

Submitted 2 July, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJL

Can the solar p-modes contribute to the high-frequency transverse oscillations of spicules?

Authors: Hidetaka Kuniyoshi, Munehito Shoda, Richard J. Morton, Takaaki Yokoyama

Abstract: Lateral motions of spicules serve as vital indicators of transverse waves in the solar atmosphere, and their study is crucial for understanding the wave heating process of the corona. Recent observations have focused on "high-frequency" transverse waves (periods < 100 s), which have the potential to transport sufficient energy for coronal heating. These high-frequency spicule oscillations are dist… ▽ More Lateral motions of spicules serve as vital indicators of transverse waves in the solar atmosphere, and their study is crucial for understanding the wave heating process of the corona. Recent observations have focused on "high-frequency" transverse waves (periods < 100 s), which have the potential to transport sufficient energy for coronal heating. These high-frequency spicule oscillations are distinct from granular motions, which have much longer time scales of 5-10 min. Instead, it is proposed that they are generated through the mode conversion from high-frequency longitudinal waves that arise from a shock steepening process. Therefore, these oscillations may not solely be produced by the horizontal buffeting motions of granulation but also by the leakage of p-mode oscillations. To investigate the contribution of p-modes, our study employs a two-dimensional magneto-convection simulation spanning from the upper convection zone to the corona. During the course of the simulation, we introduce a p-mode-like driver at the bottom boundary. We reveal a notable increase in the mean velocity amplitude of the transverse oscillations in spicules, ranging from 10% to 30%, and attribute this to the energy transfer from longitudinal to transverse waves. This effect results in an enhancement of the estimated energy flux by 30-80%. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ, 14 pages, 6 figures, 2 tables, Comments are welcome

Magnetic Tornado Properties: A Substantial Contribution to the Solar Coronal Heating via Efficient Energy Transfer

Authors: Hidetaka Kuniyoshi, Munehito Shoda, Haruhisa Iijima, Takaaki Yokoyama

Abstract: In solving the solar coronal heating problem, it is crucial to comprehend the mechanisms by which energy is conveyed from the photosphere to the corona. Recently, magnetic tornadoes, characterized as coherent, rotating magnetic field structures extending from the photosphere to the corona, have drawn growing interest as a possible means of efficient energy transfer. Despite its acknowledged import… ▽ More In solving the solar coronal heating problem, it is crucial to comprehend the mechanisms by which energy is conveyed from the photosphere to the corona. Recently, magnetic tornadoes, characterized as coherent, rotating magnetic field structures extending from the photosphere to the corona, have drawn growing interest as a possible means of efficient energy transfer. Despite its acknowledged importance, the underlying physics of magnetic tornadoes remains still elusive. In this study, we conduct a three-dimensional radiative magnetohydrodynamic simulation that encompasses the upper convective layer and extends into the corona, with a view to investigating how magnetic tornadoes are generated and efficiently transfer energy into the corona. We find that a single event of magnetic flux concentration merger on the photosphere gives rise to the formation of a single magnetic tornado. The Poynting flux transferred into the corona is found to be four times greater in the presence of the magnetic tornado, as compared to its absence. This increase is attributed to a reduction in energy loss in the chromosphere, resulting from the weakened magnetic energy cascade. Based on an evaluation of the fraction of the merging events, our results suggest that magnetic tornadoes contribute approximately 50% of the Poynting flux into the corona in regions where the coronal magnetic field strength is 10 G. Potentially, the contribution could be even greater in areas with a stronger coronal magnetic field. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: submitted to the Astrophysical Journal. first referee report received. comments welcome

Presolar stardust in asteroid Ryugu

Authors: Jens Barosch, Larry R. Nittler, Jianhua Wang, Conel M. O'D. Alexander, Bradley T. De Gregorio, Cécile Engrand, Yoko Kebukawa, Kazuhide Nagashima, Rhonda M. Stroud, Hikaru Yabuta, Yoshinari Abe, Jérôme Aléon, Sachiko Amari, Yuri Amelin, Ken-ichi Bajo, Laure Bejach, Martin Bizzarro, Lydie Bonal, Audrey Bouvier, Richard W. Carlson, Marc Chaussidon, Byeon-Gak Choi, George D. Cody, Emmanuel Dartois, Nicolas Dauphas , et al. (99 additional authors not shown)

Abstract: We have conducted a NanoSIMS-based search for presolar material in samples recently returned from C-type asteroid Ryugu as part of JAXA's Hayabusa2 mission. We report the detection of all major presolar grain types with O- and C-anomalous isotopic compositions typically identified in carbonaceous chondrite meteorites: 1 silicate, 1 oxide, 1 O-anomalous supernova grain of ambiguous phase, 38 SiC, a… ▽ More We have conducted a NanoSIMS-based search for presolar material in samples recently returned from C-type asteroid Ryugu as part of JAXA's Hayabusa2 mission. We report the detection of all major presolar grain types with O- and C-anomalous isotopic compositions typically identified in carbonaceous chondrite meteorites: 1 silicate, 1 oxide, 1 O-anomalous supernova grain of ambiguous phase, 38 SiC, and 16 carbonaceous grains. At least two of the carbonaceous grains are presolar graphites, whereas several grains with moderate C isotopic anomalies are probably organics. The presolar silicate was located in a clast with a less altered lithology than the typical extensively aqueously altered Ryugu matrix. The matrix-normalized presolar grain abundances in Ryugu are 4.8$^{+4.7}_{-2.6}$ ppm for O-anomalous grains, 25$^{+6}_{-5}$ ppm for SiC grains and 11$^{+5}_{-3}$ ppm for carbonaceous grains. Ryugu is isotopically and petrologically similar to carbonaceous Ivuna-type (CI) chondrites. To compare the in situ presolar grain abundances of Ryugu with CI chondrites, we also mapped Ivuna and Orgueil samples and found a total of SiC grains and 6 carbonaceous grains. No O-anomalous grains were detected. The matrix-normalized presolar grain abundances in the CI chondrites are similar to those in Ryugu: 23 $^{+7}_{-6}$ ppm SiC and 9.0$^{+5.3}_{-4.6}$ ppm carbonaceous grains. Thus, our results provide further evidence in support of the Ryugu-CI connection. They also reveal intriguing hints of small-scale heterogeneities in the Ryugu samples, such as locally distinct degrees of alteration that allowed the preservation of delicate presolar material. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 12 pages, 3 figures, 2 tables. Published in ApJL

Journal ref: 2022, The Astrophysical Journal Letters, 935, L3 (12pp)

Mean-field analysis on large-scale magnetic fields at high Reynolds numbers

Authors: Ryota Shimada, Hideyuki Hotta, Takaaki Yokoyama

Abstract: Solar magnetic fields comprise an 11-year activity cycle, represented by the number of sunspots. The maintenance of such a solar magnetic field can be attributed to fluid motion in the convection zone, i.e. a dynamo. This study conducts the mean-field analyses of the global solar dynamo simulation presented by Hotta et al. (2016). Although the study succeeds in producing coherent large-scale magne… ▽ More Solar magnetic fields comprise an 11-year activity cycle, represented by the number of sunspots. The maintenance of such a solar magnetic field can be attributed to fluid motion in the convection zone, i.e. a dynamo. This study conducts the mean-field analyses of the global solar dynamo simulation presented by Hotta et al. (2016). Although the study succeeds in producing coherent large-scale magnetic fields at high Reynolds numbers, the detailed physics of the maintenance of this field have not been fully understood. This study extracts the alpha-tensor and the turbulent magnetic diffusivity tensor through mean-field analyses. The turbulent magnetic diffusivity exhibits a significant decrease towards high Reynolds numbers. The decrease in the turbulent magnetic diffusivity suppresses the energy conversion of large-scale field to small-scale field. This implies that the decrease in the turbulent magnetic diffusivity contributes to the maintenance of a large-scale magnetic field at high Reynolds numbers. A significant downward turbulent pumping is observed; it is enhanced in the weak phase of the large-scale field. This study proposes a cyclic reversal process of a large-scale field which is dominantly driven by the alpha-effect and is possibly triggered by downward pumping. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 17 pages, 20 figures, Accepted for publication in ApJ

Fast magnetic wave could heat the solar low-beta chromosphere

Authors: Yikang Wang, Takaaki Yokoyama, Haruhisa Iijima

Abstract: Magnetohydrodynamic (MHD) waves are candidates for heating the solar chromosphere, although it is still unclear which mode of the wave is dominant in heating. We perform two-dimensional radiative MHD simulation to investigate the propagation of MHD waves in the quiet region of the solar chromosphere. We identify the mode of the shock waves by using the relationship between gas pressure and magneti… ▽ More Magnetohydrodynamic (MHD) waves are candidates for heating the solar chromosphere, although it is still unclear which mode of the wave is dominant in heating. We perform two-dimensional radiative MHD simulation to investigate the propagation of MHD waves in the quiet region of the solar chromosphere. We identify the mode of the shock waves by using the relationship between gas pressure and magnetic pressure across the shock front and calculate their corresponding heating rate through the entropy jump to obtain a quantitative understanding of the wave heating process in the chromosphere. Our result shows that the fast magnetic wave is significant in heating the low-beta chromosphere. The low-beta fast magnetic waves are generated from high-beta fast acoustic waves via mode conversion crossing the equipartition layer. Efficient mode conversion is achieved by large attacking angles between the propagation direction of the shock waves and the chromospheric magnetic field. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: Accepted by ApJL

Variable refractory lithophile element compositions of planetary building blocks: insights from components of enstatite chondrites

Authors: Takashi Yoshizaki, Richard D. Ash, Marc D. Lipella, Tetsuya Yokoyama, William F. McDonough

Abstract: Chondrites are sediments of materials left over from the earliest stage of the solar system history. Based on their undifferentiated nature and less fractionated chemical compositions, chondrites are widely considered to represent the unprocessed building blocks of the terrestrial planets and their embryos. Models of chemical composition of the terrestrial planets generally find chondritic relativ… ▽ More Chondrites are sediments of materials left over from the earliest stage of the solar system history. Based on their undifferentiated nature and less fractionated chemical compositions, chondrites are widely considered to represent the unprocessed building blocks of the terrestrial planets and their embryos. Models of chemical composition of the terrestrial planets generally find chondritic relative abundances of refractory lithophile elements (RLE) in the bulk bodies ("constant RLE ratio rule"), based on limited variations of RLE ratios among chondritic meteorites and the solar photosphere. Here, we show that ratios of RLE, such as Nb/Ta, Zr/Hf, Sm/Nd and Al/Ti, are fractionated from the solar value in chondrules from enstatite chondrites (EC). The fractionated RLE ratios of individual EC chondrules document different chalcophile affinities of RLE under highly reducing environments and a separation of RLE-bearing sulfides from silicates before and/or during chondrule formation. In contrast, the bulk EC have solar-like RLE ratios, indicating that a physical sorting of silicates and sulfides was negligible before and during the accretion of EC parent bodies. Likewise, if the Earth's accretion were dominated by EC-like materials, as supported by multiple isotope systematics, the physical sorting of silicates and sulfides in the accretionary disk should not have occurred among the Earth's building blocks. Alternatively, the Earth's precursors might have been high-temperature nebular materials that condensed before the RLE fractionation due to precipitation of the RLE-bearing sulfides...Highly reduced planets that have experienced selective removal or accretion of silicates or metal/sulfide phases, such as Mercury, might have fractionated, non-solar bulk RLE ratios. △ Less

Submitted 30 May, 2021; v1 submitted 26 November, 2020; originally announced November 2020.

Comments: 73 pages, 16 figures, 3 tables

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena throughout the Universe

Authors: H. Ji, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, L. -J. Chen, Y. Chen, A. Chien, D. Craig, J. Dahlin, W. Daughton, E. DeLuca, C. F. Dong, S. Dorfman, J. Drake, F. Ebrahimi , et al. (75 additional authors not shown)

Abstract: This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys

Simulation of Alfven wave propagation in magnetic chromosphere with radiative loss: effects of non-linear mode coupling on chromospheric heating

Authors: Yikang Wang, Takaaki Yokoyama

Abstract: We perform magnetohydrodynamic (MHD) simulations to investigate the propagation of Alfven wave in magnetic chromosphere. We use the 1.5-dimensional expanding flux tube geometry setting and transverse perturbation at the bottom to generate the Alfven wave. Compared with previous studies, our expansion is that we include the radiative loss term introduced by Carlsson & Leenaarts (2012). We find that… ▽ More We perform magnetohydrodynamic (MHD) simulations to investigate the propagation of Alfven wave in magnetic chromosphere. We use the 1.5-dimensional expanding flux tube geometry setting and transverse perturbation at the bottom to generate the Alfven wave. Compared with previous studies, our expansion is that we include the radiative loss term introduced by Carlsson & Leenaarts (2012). We find that when an observation based transverse wave generator is applied, the spatial distribution of the time-averaged radiative loss profile in our simulation is consistent with that in the classic atmospheric model. In addition, the energy flux in the corona is larger than the required value for coronal heating in the quiet region. Our study shows that the Alfven wave driven model has the potential to explain chromospheric heating and transport enough energy to the corona simultaneously. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: Accepted by ApJ

Inverse FIP Effects in Giant Solar Flares Found from Earth X-Ray Albedo with Suzaku/XIS

Authors: Satoru Katsuda, Masanori Ohno, Koji Mori, Tatsuhiko Beppu, Yoshiaki Kanemaru, Makoto S. Tashiro, Yukikatsu Terada, Kosuke Sato, Kae Morita, Hikari Sagara, Futa Ogawa, Haruya Takahashi, Hiroshi Murakami, Masayoshi Nobukawa, Hiroshi Tsunemi, Kiyoshi Hayashida, Hironori Matsumoto, Hirofumi Noda, Hiroshi Nakajima, Yuichiro Ezoe, Yohko Tsuboi, Yoshitomo Maeda, Takaaki Yokoyama, Noriyuki Narukage

Abstract: We report X-ray spectroscopic results for four giant solar flares occurred on 2005 September 7 (X17.0), 2005 September 8 (X5.4), 2005 September 9 (X6.2), and 2006 December 5 (X9.0), obtained from Earth albedo data with the X-ray imaging spectrometer (XIS) onboard Suzaku. The good energy resolution of the XIS (FWHM$\sim$100 eV) enables us to separate a number of line-like features and detect the un… ▽ More We report X-ray spectroscopic results for four giant solar flares occurred on 2005 September 7 (X17.0), 2005 September 8 (X5.4), 2005 September 9 (X6.2), and 2006 December 5 (X9.0), obtained from Earth albedo data with the X-ray imaging spectrometer (XIS) onboard Suzaku. The good energy resolution of the XIS (FWHM$\sim$100 eV) enables us to separate a number of line-like features and detect the underlying continuum emission. These features include Si He$α$, Si Ly$α$, S He$α$, S Ly$α$, Ar He$α$, and Ca He$α$ originating from solar flares as well as fluorescent Ar K$α$ and Ar K$β$ from the Earth atmosphere. Absolute elemental abundances (X/H) averaged over the four flares are obtained to be $\sim$2.0 (Ca), $\sim$0.7 (Si), $\sim$0.3 (S), and $\sim$0.9 (Ar) at around flare peaks. This abundance pattern is similar to those of active stars' coronae showing inverse first ionization potential (i-FIP) effects, i.e., elemental abundances decrease with decreasing FIP with a turnover at the low end of FIP. The abundances are almost constant during the flares, with an exception of Si which increases by a factor of $\sim$2 in the decay phase. The evolution of the Si abundance is consistent with a picture that the i-FIP plasma originates from the chromosphere evaporation and then mixes with the surrounding low-FIP biased materials. Flare-to-flare abundance varied by a factor of 2, agreeing with past observations of solar flares. Finally, we emphasize that Earth albedo data acquired by X-ray astronomy satellites like Suzaku and XRISM can significantly contribute to studies of solar physics. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: 19 pages, 10 figures, accepted for publication in The Astrophysical Journal

Estimating the temperature and density of a spicule from 100 GHz data obtained with ALMA

Authors: Masumi Shimojo, Tomoko Kawate, Takenori J. Okamoto, Takaaki Yokoyama, Noriyuki Narukage, Taro Sakao, Kazumasa Iwai, Gregory D. Fleishman, Kazunari Shibata

Abstract: We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304Å images (MgII/304A spicule). The other one is a spicule seen in the 100GHz im… ▽ More We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304Å images (MgII/304A spicule). The other one is a spicule seen in the 100GHz images obtained with ALMA (100GHz spicule). Although the 100GHz spicule overlapped with the MgII/304A spicule in the early phase, it did not show any corresponding structures in the IRIS Mg II and AIA 304A images after the early phase. It suggests that the spicules are individual events and do not have a physical relationship. To obtain the physical parameters of the 100GHz spicule, we estimate the optical depths as a function of temperature and density using two different methods. One is using the observed brightness temperature by assuming a filling factor, and the other is using an emission model for the optical depth. As a result of comparing them, the kinetic temperature of the plasma and the number density of ionized hydrogens in the 100GHz spicule are ~6800 K and 2.2 x 10^10 cm^-3. The estimated values can explain the absorbing structure in the 193A image, which appear as a counterpart of the 100GHz spicule. These results suggest that the 100GHz spicule presented in this paper is classified to a macrospicule without a hot sheath in former terminology. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 8 pages, 4 figures, Accepted to ApJL

Journal ref: Astrophysical Journal Letters, Vol. 888, L28 (2020)

Diffusion regions and 3D energy mode development in spontaneous reconnection

Authors: Shuoyang Wang, Takaaki Yokoyama

Abstract: The understanding of magnetic reconnection in three-dimension (3D) is far shallower than its counterpart in two-dimension (2D) due to its potential complication, not to mention the evolving of the spontaneously growing turbulence. We investigate the reason for reconnection acceleration on the characters and development of diffusion regions and sheared 3D energy modes (energy modes that are not par… ▽ More The understanding of magnetic reconnection in three-dimension (3D) is far shallower than its counterpart in two-dimension (2D) due to its potential complication, not to mention the evolving of the spontaneously growing turbulence. We investigate the reason for reconnection acceleration on the characters and development of diffusion regions and sheared 3D energy modes (energy modes that are not parallel to the anti-parallel magnetic fields) during the turbulence building stage. We found that multiple reconnection layers emerge due to the growth of 3D sheared tearing instability. Diffusion regions on adjacent reconnection layers form an inflow-outflow coupling that enhances the local reconnection. Further coupling of the existing energy modes breeds new energy modes near the current sheet edge. As reconnection layers span and interact with each other across the whole current sheet, global magnetic energy consumption accelerates. The significant contribution of 3D energy modes and their interaction to the reconnection rate enhancement seems to be independent on magnetic diffusivity. On the other hand, the global guide field changes the layout of 3D reconnection layer thus determines whether the system is fast-reconnection-preferable. △ Less

Submitted 10 June, 2019; originally announced June 2019.

Comments: 13 pages, 14 figures

Three-dimensional simulation of the fast solar wind driven by compressible magnetohydrodynamic turbulence

Authors: Munehito Shoda, Takeru Ken Suzuki, Mahboubeh Asgari-Targhi, Takaaki Yokoyama

Abstract: Using a three-dimensional compressible magnetohydrodynamic (MHD) simulation, we have reproduced the fast solar wind in a direct and self-consistent manner, based on the wave/turbulence driven scenario. As a natural consequence of Alfvenic perturbations at its base, highly compressional and turbulent fluctuations are generated, leading to heating and acceleration of the solar wind. The analysis of… ▽ More Using a three-dimensional compressible magnetohydrodynamic (MHD) simulation, we have reproduced the fast solar wind in a direct and self-consistent manner, based on the wave/turbulence driven scenario. As a natural consequence of Alfvenic perturbations at its base, highly compressional and turbulent fluctuations are generated, leading to heating and acceleration of the solar wind. The analysis of power spectra and structure functions reveals that the turbulence is characterized by its imbalanced (in the sense of outward Alfvenic fluctuations) and anisotropic nature. The density fluctuation originates from the parametric decay instability of outwardly propagating Alfven waves and plays a significant role in the Alfven wave reflection that triggers turbulence. Our conclusion is that the fast solar wind is heated and accelerated by compressible MHD turbulence driven by parametric decay instability and resultant Alfven wave reflection. △ Less

Submitted 19 July, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: published from The Astrophysical Journal Letters

Observations of Photospheric Magnetic Structure below a Dark Filament using the Hinode Spectro-Polarimeter

Authors: Takaaki Yokoyama, Yukio Katsukawa, Masumi Shimojo

Abstract: The structure of the photospheric vector magnetic field below a dark filament on the Sun is studied using the observations of the Spectro-Polarimeter attached to the Solar Optical Telescope onboard Hinode. Special attention is paid to discriminate the two suggested models, a flux rope or a bent arcade. "Inverse-polarity" orientation is possible below the filament in a flux rope, whereas "normal-po… ▽ More The structure of the photospheric vector magnetic field below a dark filament on the Sun is studied using the observations of the Spectro-Polarimeter attached to the Solar Optical Telescope onboard Hinode. Special attention is paid to discriminate the two suggested models, a flux rope or a bent arcade. "Inverse-polarity" orientation is possible below the filament in a flux rope, whereas "normal-polarity" can appear in both models. We study a filament in active region NOAA 10930, which appeared on the solar disk during 2006 December. The transverse field perpendicular to the line of sight has a direction almost parallel to the filament spine with a shear angle of 30 deg, whose orientation includes the 180-degree ambiguity. To know whether it is in the normal orientation or in the inverse one, the center-to-limb variation is used for the solution under the assumption that the filament does not drastically change its magnetic structure during the passage. When the filament is near the east limb, we found that the line-of-site magnetic component below is positive, while it is negative near the west limb. This change of sign indicates that the horizontal photospheric field perpendicular to the polarity inversion line beneath the filament has an "inverse-polarity", which indicates a flux-rope structure of the filament supporting field. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 21 pages, 8 figures, accepted for publication in PASJ

Chemically defining the building blocks of the Earth

Authors: Takashi Yoshizaki, Richard D. Ash, Tetsuya Yokoyama, Marc D. Lipella, William F. McDonough

Abstract: Chondrites are undifferentiated sediments of material left over from the earliest solar system and are widely considered as representatives of the unprocessed building blocks of the terrestrial planets. The chondrites, along with processed igneous meteorites, have been divided into two broad categories based upon their isotopic signatures; these have been termed the CC and NC groups and have been… ▽ More Chondrites are undifferentiated sediments of material left over from the earliest solar system and are widely considered as representatives of the unprocessed building blocks of the terrestrial planets. The chondrites, along with processed igneous meteorites, have been divided into two broad categories based upon their isotopic signatures; these have been termed the CC and NC groups and have been interpreted as reflecting their distinctive birth places within the solar system. The isotopic distinctiveness of NC and CC meteorites document limited radial-mixing in the accretionary disk. The enstatite and ordinary chondrites are NC-type and likely represent samples from inner solar system (likely $<$4 AU). Measurement and modeling of ratios of refractory lithophile elements (RLE) in enstatite chondrites establish these meteorites as the closest starting materials for the bulk of the silicate Earth and the core. Comparing chondritic and terrestrial RLE ratios demonstrate that the Bulk Silicate Earth, not the core, host the Earth's inventory of Ti, Zr, Nb, and Ta, but not the full complement of V. △ Less

Submitted 13 January, 2019; v1 submitted 31 December, 2018; originally announced December 2018.

Comments: 28 pages, 7 figures

Impact of Dynamic State on the Mass Condensation Rate of Solar Prominences

Authors: Takafumi Kaneko, Takaaki Yokoyama

Abstract: The interiors of quiescent prominences are filled with turbulent flows. The evolution of upflow plumes, descending pillars, and vortex motions has been clearly detected in high-resolution observations. The Rayleigh-Taylor instability is thought to be a driver of such internal flows. Descending pillars are related to the mass budgets of prominences. There is a hypothesis of dynamic equilibrium wher… ▽ More The interiors of quiescent prominences are filled with turbulent flows. The evolution of upflow plumes, descending pillars, and vortex motions has been clearly detected in high-resolution observations. The Rayleigh-Taylor instability is thought to be a driver of such internal flows. Descending pillars are related to the mass budgets of prominences. There is a hypothesis of dynamic equilibrium where the mass drainage via descending pillars and the mass supply via radiative condensation are balanced to maintain the prominence mass; however, the background physics connecting the two different processes is poorly understood. In this study, we reproduced the dynamic interior of a prominence via radiative condensation and the mechanism similar to the Rayleigh-Taylor instability using a three-dimensional magnetohydrodynamic simulation including optically thin radiative cooling and nonlinear anisotropic thermal conduction. The process to prominence formation in the simulation follows the reconnection-condensation model, where topological change in the magnetic field caused by reconnection leads to radiative condensation. Reconnection is driven by converging motion at the footpoints of the coronal arcade fields. In contrast to the previous model, by randomly changing the speed of the footpoint motion along a polarity inversion line, the dynamic interior of prominence is successfully reproduced. We find that the mass condensation rate of the prominence is enhanced in the case with dynamic state. Our results support the observational hypothesis that the condensation rate is balanced with the mass drainage rate and suggest that a self-induced mass maintenance mechanism exists. △ Less

Submitted 2 November, 2018; originally announced November 2018.

Comments: Accepted for publication in ApJ, 26 pages, 12 figures

ALMA Observations of the Solar Chromosphere on the Polar Limb

Authors: Takaaki Yokoyama, Masumi Shimojo, Takenori J. Okamoto, Haruhisa Iijima

Abstract: We report the results of the Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the solar chromosphere on the southern polar limb. Coordinated observations with the Interface Region Imaging Spectrograph (IRIS) are also conducted. ALMA provided unprecedented high spatial resolution in the millimeter band ($\approx$ 2.0 arcsec) at 100 GHz frequency with a moderate cadence (20 s). T… ▽ More We report the results of the Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the solar chromosphere on the southern polar limb. Coordinated observations with the Interface Region Imaging Spectrograph (IRIS) are also conducted. ALMA provided unprecedented high spatial resolution in the millimeter band ($\approx$ 2.0 arcsec) at 100 GHz frequency with a moderate cadence (20 s). The results are as follows: (1) The ALMA 100 GHz images show saw-tooth patterns on the limb, and a comparison with SDO/AIA 171Å images shows a good correspondence of the limbs with each other. (2) The ALMA 100 GHz movie shows a dynamic thorn-like structure elongating from the saw-tooth patterns on the limb, with lengths reaching at least 8 arcsec, thus suggesting jet-like activity in the ALMA microwave range. These ALMA jets are in good correspondence with IRIS jet clusters. (3) A blob ejection event is observed. By comparing with the IRIS Mg II slit-jaw images, the trajectory of the blob is located along the spicular patterns. △ Less

Submitted 3 July, 2018; originally announced July 2018.

Comments: Accepted for publication in The Astrophysical Journal

Anisotropic magnetohydrodynamic turbulence driven by parametric decay instability: the onset of phase mixing and Alfvén wave turbulence

Authors: Munehito Shoda, Takaaki Yokoyama

Abstract: We conduct a three-dimensional magnetohydrodynamic (MHD) simulation of the parametric decay instability of Alfvén waves and resultant compressible MHD turbulence, which is likely to develop in the solar wind acceleration region. Because of the presence of the mean magnetic field, the nonlinear stage is characterized by filament-like structuring and anisotropic cascading. By calculating the timesca… ▽ More We conduct a three-dimensional magnetohydrodynamic (MHD) simulation of the parametric decay instability of Alfvén waves and resultant compressible MHD turbulence, which is likely to develop in the solar wind acceleration region. Because of the presence of the mean magnetic field, the nonlinear stage is characterized by filament-like structuring and anisotropic cascading. By calculating the timescales of phase mixing and the evolution of Alfvén wave turbulence, we have found that the early nonlinear stage is dominated by phase mixing, while the later phase is dominated by imbalanced Alfvén wave turbulence. Our results indicate that the regions in the solar atmosphere with large density fluctuation, such as the coronal bottom and wind acceleration region, are heated by phase-mixed Alfvén waves, while the other regions are heated by Alfvén wave turbulence. △ Less

Submitted 25 May, 2018; v1 submitted 1 May, 2018; originally announced May 2018.

Comments: Accepted for publication in The Astrophysical Journal Letters

Frequency-dependent Alfvén-wave propagation in the solar wind: Onset and suppression of parametric decay instability

Authors: Munehito Shoda, Takaaki Yokoyama, Takeru K. Suzuki

Abstract: Using numerical simulations, we investigate the onset and suppression of parametric decay instability (PDI) in the solar wind, focusing on the suppression effect by the wind acceleration and expansion. We also discuss the role of PDI in generating the density fluctuation in the solar wind and the cross-helicity evolution up to $1 {\rm \ au}$. Numerical simulations from the coronal base to 1 au are… ▽ More Using numerical simulations, we investigate the onset and suppression of parametric decay instability (PDI) in the solar wind, focusing on the suppression effect by the wind acceleration and expansion. We also discuss the role of PDI in generating the density fluctuation in the solar wind and the cross-helicity evolution up to $1 {\rm \ au}$. Numerical simulations from the coronal base to 1 au are conducted in a self-consistent manner; we take into account the feedback of wave energy and pressure in the background. Monochromatic waves with various injection frequencies $f_0$ are injected to discuss the suppression of PDI, while broadband waves are applied to compare the numerical results with observation. We find, via the simulations with monochromatic waves, that high-frequency ($f_0 \gtrsim10^{-3} {\rm \ Hz}$) Alfvén waves are subject to PDI. Meanwhile, the maximum growth rate of the PDI of low-frequency ($f_0 \lesssim10^{-4} {\rm \ Hz}$) Alfvén waves becomes negative due to acceleration and expansion effects. Medium-frequency ($f_0 \approx 10^{-3.5} {\rm \ Hz}$) Alfvén waves have a positive growth rate but do not show the signature of PDI up to 1 au because the growth rate is too small. The medium-frequency waves experience neither PDI nor reflection so they propagate through the solar wind most efficiently. This is called the frequency-filtering mechanism of the solar wind. The simulations with broadband waves indicate that the observed trend of the density fluctuation is well explained by the evolution of PDI while the observed cross-helicity evolution is in agreement with low-frequency wave propagation. △ Less

Submitted 1 May, 2018; v1 submitted 7 March, 2018; originally announced March 2018.

Comments: Accepted for publication in The Astrophysical Journal

High-frequency spicule oscillations generated via mode conversion

Authors: Munehito Shoda, Takaaki Yokoyama

Abstract: Spicule oscillations involve high-frequency components with a typical period approximately corresponding to $40-50$ s. The typical time scale of the photospheric oscillation is a few minutes, and thus, the origin of this high-frequency component is not trivial. In this study, a one-dimensional numerical simulation is performed to demonstrate that the observed spicule oscillations originate from lo… ▽ More Spicule oscillations involve high-frequency components with a typical period approximately corresponding to $40-50$ s. The typical time scale of the photospheric oscillation is a few minutes, and thus, the origin of this high-frequency component is not trivial. In this study, a one-dimensional numerical simulation is performed to demonstrate that the observed spicule oscillations originate from longitudinal-to-transverse mode conversion that occurs around the equipartition layer in the chromosphere. Calculations are conducted in a self-consistent manner with the exception of additional heating to maintain coronal temperature. The analyses indicate the following features: (1) mode conversion efficiently excites high-frequency transverse waves; (2) the typical period of the high-frequency waves corresponds to the sound-crossing time of the mode conversion region; and (3) simulated root-mean-square velocity of the high-frequency component is consistent with the observed value. These results indicate that the observation of spicule oscillation provides direct evidence of mode conversion in the chromosphere. △ Less

Submitted 4 January, 2018; originally announced January 2018.

Comments: accepted for publication in The Astrophysical Journal

A self-consistent model of the coronal heating and solar wind acceleration including compressible and incompressible heating processes

Authors: Munehito Shoda, Takaaki Yokoyama, Takeru K. Suzuki

Abstract: We propose a novel one-dimensional model that includes both shock and turbulence heating and qualify how these processes contribute to heating the corona and driving the solar wind. Compressible MHD simulations allow us to automatically consider shock formation and dissipation, while turbulent dissipation is modeled via a one-point closure based on Alfvén wave turbulence. Numerical simulations wer… ▽ More We propose a novel one-dimensional model that includes both shock and turbulence heating and qualify how these processes contribute to heating the corona and driving the solar wind. Compressible MHD simulations allow us to automatically consider shock formation and dissipation, while turbulent dissipation is modeled via a one-point closure based on Alfvén wave turbulence. Numerical simulations were conducted with different photospheric perpendicular correlation lengths $λ_0$, which is a critical parameter of Alfvén wave turbulence, and different root-mean-square photospheric transverse-wave amplitudes $δv_0$. For the various $λ_0$, we obtain a low-temperature chromosphere, high-temperature corona, and supersonic solar wind. Our analysis shows that turbulence heating is always dominant when $λ_0 \lesssim 1{\rm \ Mm}$. This result does not mean that we can ignore the compressibility because the analysis indicates that the compressible waves and their associated density fluctuations enhance the Alfvén wave reflection and therefore the turbulence heating. The density fluctuation and the cross helicity are strongly affected by $λ_0$, while the coronal temperature and mass loss rate depend weakly on $λ_0$. △ Less

Submitted 14 February, 2018; v1 submitted 20 December, 2017; originally announced December 2017.

Comments: Accepted for publicaiton in The Astrophysical Journal

Convective velocity suppression via the enhancement of subadiabatic layer: Role of the effective Prandtl number

Authors: Y. Bekki, H. Hotta, T. Yokoyama

Abstract: It has recently been recognized that the convective velocities achieved in the current solar convection simulations might be over-estimated. The newly-revealed effects of the prevailing small-scale magnetic field within the convection zone may offer possible solutions to this problem. The small-scale magnetic fields can reduce the convective amplitude of small-scale motions through the Lorentz-for… ▽ More It has recently been recognized that the convective velocities achieved in the current solar convection simulations might be over-estimated. The newly-revealed effects of the prevailing small-scale magnetic field within the convection zone may offer possible solutions to this problem. The small-scale magnetic fields can reduce the convective amplitude of small-scale motions through the Lorentz-force feedback, which concurrently inhibits the turbulent mixing of entropy between upflows and downflows. As a result, the effective Prandtl number may exceed unity inside the solar convection zone. In this paper, we propose and numerically confirm a possible suppression mechanism of convective velocity in the effectively high-Prandtl number regime. If the effective horizontal thermal diffusivity decreases (the Prandtl number accordingly increases), the subadiabatic layer which is formed near the base of the convection zone by continuous depositions of low entropy transported by adiabatically downflowing plumes is enhanced and extended. The global convective amplitude in the high-Prandtl thermal convection is thus reduced especially in the lower part of the convection zone via the change in the mean entropy profile which becomes more subadiabatic near the base and less superadiabatic in the bulk. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 17 pages, 13 figures, accepted for publication in ApJ

Three-dimensional Magnetohydrodynamic Simulation of the Formation of Solar Chromospheric Jets with Twisted Magnetic Field Lines

Authors: H. Iijima, T. Yokoyama

Abstract: This paper presents a three-dimensional simulation of chromospheric jets with twisted magnetic field lines. Detailed treatments of the photospheric radiative transfer and the equation of states allow us to model realistic thermal convection near the solar surface, which excites various MHD waves and produces chromospheric jets in the simulation. A tall chromospheric jet with a maximum height of 10… ▽ More This paper presents a three-dimensional simulation of chromospheric jets with twisted magnetic field lines. Detailed treatments of the photospheric radiative transfer and the equation of states allow us to model realistic thermal convection near the solar surface, which excites various MHD waves and produces chromospheric jets in the simulation. A tall chromospheric jet with a maximum height of 10--11 Mm and lifetime of 8--10 min is formed above a strong magnetic field concentration. The magnetic field lines are strongly entangled in the chromosphere, which helps the chromospheric jet to be driven by the Lorentz force. The jet exhibits oscillatory motion as a natural consequence of its generation mechanism. We also find that the produced chromospheric jet forms a cluster with a diameter of several Mm with finer strands. These results imply a close relationship between the simulated jet and solar spicules. △ Less

Submitted 5 September, 2017; originally announced September 2017.

Comments: 33 pages, 15 figures, accepted for publication in ApJ

Reconnection-Condensation Model for Solar Prominence Formation

Authors: Takafumi Kaneko, Takaaki Yokoyama

Abstract: We propose a reconnection-condensation model in which topological change in a coronal magnetic field via reconnection triggers radiative condensation, thereby resulting in prominence formation. Previous observational studies have suggested that reconnection at a polarity inversion line of a coronal arcade field creates a flux rope that can sustain a prominence; however, they did not explain the or… ▽ More We propose a reconnection-condensation model in which topological change in a coronal magnetic field via reconnection triggers radiative condensation, thereby resulting in prominence formation. Previous observational studies have suggested that reconnection at a polarity inversion line of a coronal arcade field creates a flux rope that can sustain a prominence; however, they did not explain the origin of cool dense plasmas of prominences. Using three-dimensional magnetohydrodynamic simulations including anisotropic nonlinear thermal conduction and optically thin radiative cooling, we demonstrate that reconnection can lead not only to flux rope formation but also to radiative condensation under a certain condition. In our model, this condition is described by the Field length, which is defined as the scale length for thermal balance between radiative cooling and thermal conduction. This critical condition depends weakly on the artificial background heating. The extreme ultraviolet emissions synthesized with our simulation results have good agreement with observational signatures reported in previous studies. △ Less

Submitted 3 July, 2017; v1 submitted 29 June, 2017; originally announced June 2017.

Comments: 22 pages, 11 figures, accepted to ApJ

Observational signatures of transverse MHD waves and associated dynamic instabilities

Authors: Patrick Antolin, Ineke De Moortel, Tom Van Doorsselaere, Takaaki Yokoyama

Abstract: MHD waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave power. Detection is limited by instrumental constraints, but also by wave processes that localise the wave power in undetectable spatial scales. In this study we conduct 3D MHD simulations and forward modelling of standing trans… ▽ More MHD waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave power. Detection is limited by instrumental constraints, but also by wave processes that localise the wave power in undetectable spatial scales. In this study we conduct 3D MHD simulations and forward modelling of standing transverse MHD waves in coronal loops with uniform and non-uniform temperature variation in the perpendicular cross-section. The observed signatures are largely dominated by the combination of the Kelvin-Helmholtz instability (KHI), resonant absorption and phase mixing. In the presence of a cross-loop temperature gradient we find that emission lines sensitive to the loop core catch different signatures than those more sensitive to the loop boundary and the surrounding corona, leading to an out-of-phase intensity modulation produced by the KHI mixing. Common signatures to all considered models include an intensity and loop width modulation at half the kink period, fine strand-like structure, a characteristic arrow-shaped structure in the Doppler maps, overall line broadening in time but particularly at the loop edges. For our model, most of these features can be captured with a spatial resolution of $0.33\arcsec$ and spectral resolution of 25~km~s$^{-1}$, although severe over-estimation of the line width is obtained. Resonant absorption leads to a significant decrease of the observed kinetic energy from Doppler motions over time, which is not recovered by a corresponding increase in the line width from phase mixing and the KHI motions. We estimate this hidden wave energy to be a factor of $5-10$ of the observed value. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 24 pages, 19 figures. Accepted for publication in ApJ. In this ArXiv version, some of the figures have reduced quality due to size limits

Chemical evolution of 244Pu in the solar vicinity and its implication for the properties of r-process production

Authors: Takuji Tsujimoto, Tetsuya Yokoyama, Kenji Bekki

Abstract: Meteoritic abundances of r-process elements are analyzed to deduce the history of chemical enrichment by r-process from the beginning of disk formation to the present time in the solar vicinity, by combining the abundance information from short-lived radioactive nuclei such as 244Pu with that from stable r-process nuclei such as Eu. These two types of nuclei can be associated with one r-process ev… ▽ More Meteoritic abundances of r-process elements are analyzed to deduce the history of chemical enrichment by r-process from the beginning of disk formation to the present time in the solar vicinity, by combining the abundance information from short-lived radioactive nuclei such as 244Pu with that from stable r-process nuclei such as Eu. These two types of nuclei can be associated with one r-process event and cumulation of events till formation of the solar system, respectively. With help of the observed local star formation history, we deduce the chemical evolution of 244Pu and obtain three main results: (i) the last r-process event occurred 130-140 Myr before formation of the solar system, (ii) the present-day low 244Pu abundance as measured in deep sea reservoirs results from the low recent star formation rate compared to ~4.5 - 5 Gyr ago, and (iii) there were ~15 r-process events in the solar vicinity from formation of the Galaxy to the time of solar system formation and ~30 r-process events to the present time. Then, adopting a reasonable hypothesis that a neutron star merger is the r-process production site, we find that the ejected r-process elements are extensively spread out and mixed with interstellar matter with a mass of ~3.5 million solar masses, which is about 100 times larger than that for supernova ejecta. In addition, the event frequency of r-process production is estimated to be one per about 1400 core-collapse supernovae, which is identical to the frequency of neutron star mergers estimated from the analysis of stellar abundances. △ Less

Submitted 9 January, 2017; originally announced January 2017.

Comments: 5 pages including 2 figures, accepted for publication in ApJL

Double-Cell Type Solar Meridional Circulation Based on Mean-Field Hydrodynamic Model

Authors: Yuto Bekki, Takaaki Yokoyama

Abstract: The main object of the paper is to present the condition of the non-diffusive part of the Reynolds stress for driving the double-cell structure of the solar meridional circulation, which has been revealed by recent helioseismic observations. By conducting a set of mean-field hydrodynamic simulations, we confirm for the first time that the double-cell meridional circulation can be achieved along wi… ▽ More The main object of the paper is to present the condition of the non-diffusive part of the Reynolds stress for driving the double-cell structure of the solar meridional circulation, which has been revealed by recent helioseismic observations. By conducting a set of mean-field hydrodynamic simulations, we confirm for the first time that the double-cell meridional circulation can be achieved along with the solar-like differential rotation when the Reynolds stress transports the angular momentum upward in the lower part and downward in the upper part of the convection zone. It is concluded that, in a stationary state, the accumulated angular momentum via the Reynolds stress in the middle layer is advected to both the upper and lower parts of the convection zone by each of the two meridional circulation cells, respectively. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: 16 pages, 9 figures, accepted for publication in ApJ

Modelling observed decay-less oscillations as resonantly enhanced Kelvin-Helmholtz vortices from transverse MHD waves and their seismological application

Authors: Patrick Antolin, Ineke De Moortel, Tom Van Doorsselaere, Takaaki Yokoyama

Abstract: In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfvén waves. Due to its localised nature, a direct detection of this mechanism is extremely difficult. Yet, it is the leading theory explaining the observed fast damping of the global transverse waves. However, at odds with this theoretical… ▽ More In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfvén waves. Due to its localised nature, a direct detection of this mechanism is extremely difficult. Yet, it is the leading theory explaining the observed fast damping of the global transverse waves. However, at odds with this theoretical prediction, recent observations indicate that in the low amplitude regime such transverse MHD waves can also appear decay-less, a yet unsolved phenomenon. Recent numerical work has shown that Kelvin-Helmholtz instabilities (KHI) often accompany transverse MHD waves. In this work, we combine 3D MHD simulations and forward modelling to show that for currently achieved spatial resolution and observed small amplitudes, an apparent decay-less oscillation is obtained. This effect results from the combination of periodic brightenings produced by the KHI and the coherent motion of the KHI vortices amplified by resonant absorption. Such effect is especially clear in emission lines forming at temperatures that capture the boundary dynamics rather than the core, and reflects the low damping character of the local azimuthal Alfvén waves resonantly coupled to the kink mode. Due to phase mixing, the detected period can vary depending on the emission line, with those sensitive to the boundary having shorter periods than those sensitive to the loop core. This allows to estimate the density contrast at the boundary. △ Less

Submitted 30 September, 2016; originally announced September 2016.

Comments: 8 pages, 4 figures; Accepted for publication in The Astrophysical Journal Letters

Nonlinear reflection process of linearly-polarized, broadband Alfvén waves in the fast solar wind

Authors: Munehito Shoda, Takaaki Yokoyama

Abstract: Using one-dimensional numerical simulations, we study the elementary process of Alfvén wave reflection in a uniform medium, including nonlinear effects. In the linear regime, Alfvén wave reflection is triggered only by the inhomogeneity of the medium, whereas in the nonlinear regime, it can occur via nonlinear wave-wave interactions. Such nonlinear reflection (backscattering) is typified by decay… ▽ More Using one-dimensional numerical simulations, we study the elementary process of Alfvén wave reflection in a uniform medium, including nonlinear effects. In the linear regime, Alfvén wave reflection is triggered only by the inhomogeneity of the medium, whereas in the nonlinear regime, it can occur via nonlinear wave-wave interactions. Such nonlinear reflection (backscattering) is typified by decay instability. In most studies of decay instabilities, the initial condition has been a circularly polarized Alfvén wave. In this study we consider a linearly polarized Alfvén wave, which drives density fluctuations by its magnetic pressure force. For generality, we also assume a broadband wave with a red-noise spectrum. In the data analysis, we decompose the fluctuations into characteristic variables using local eigenvectors, thus revealing the behaviors of the individual modes. Different from circular-polarization case, we find that the wave steepening produces a new energy channel from the parent Alfvén wave to the backscattered one. Such nonlinear reflection explains the observed increasing energy ratio of the sunward to the anti-sunward Alfvénic fluctuations in the solar wind with distance against the dynamical alignment effect. △ Less

Submitted 11 February, 2016; originally announced February 2016.

Comments: 31 pages, 12 figures, accepted by The Astrophysical Journal

Doppler shift of the quiet region measured by meridional scans with the EUV Imaging Spectrometer onboard Hinode

Authors: N. Kitagawa, H. Hara, T. Yokoyama

Abstract: Spatially averaged (> 50'') EUV spectral lines in the transition region of solar quiet regions are known to be redshifted. Because the mechanism underlying this phenomenon is unclear, we require additional physical information on the lower corona for limiting the theoretical models. To acquire this information, we measured the Doppler shifts over a wide coronal temperature range (log T[K]=5.7--6.3… ▽ More Spatially averaged (> 50'') EUV spectral lines in the transition region of solar quiet regions are known to be redshifted. Because the mechanism underlying this phenomenon is unclear, we require additional physical information on the lower corona for limiting the theoretical models. To acquire this information, we measured the Doppler shifts over a wide coronal temperature range (log T[K]=5.7--6.3) using the spectroscopic data taken by the Hinode EUV Imaging Spectrometer. By analyzing the data over the center-to-limb variations covering the meridian from the south to the north pole, we successfully measured the velocity to an accuracy of 3 km/s. Below log T[K] = 6.0, the Doppler shifts of the emission lines were almost zero with an error of 1--3 km/s; above this temperature, they were blueshifted with a gradually increasing magnitude, reaching - 6.3 +/- 2.1 km/s at log T[K]=6.25. △ Less

Submitted 16 November, 2015; originally announced November 2015.

Comments: accepted for publication in ApJ

Detection of Flux Emergence, Splitting, Merging, and Cancellation of Network Fields. II Apparent Unipolar Flux Change and Cancellation

Authors: Y. Iida, H. Hagenaar, T. Yokoyama

Abstract: In this second paper in the series, we investigate occurrence frequencies of apparent unipolar processes, cancellation, and emergence of patch structures in quiet regions. Apparent unipolar events are considerably more frequent than cancellation and emergence as per our definition, which is consistent with Lamb et al. (2013). Furthermore, we investigate the frequency distributions of changes in fl… ▽ More In this second paper in the series, we investigate occurrence frequencies of apparent unipolar processes, cancellation, and emergence of patch structures in quiet regions. Apparent unipolar events are considerably more frequent than cancellation and emergence as per our definition, which is consistent with Lamb et al. (2013). Furthermore, we investigate the frequency distributions of changes in flux during apparent unipolar processes are and found that they concentrate around the detection limit of the analysis. Combining these findings with the results of our previous paper, Iida et al. (2012), that merging and splitting are more dominant than emergence and cancellation, these results support the understanding that apparent unipolar processes are actually interactions with and among patches below the detection limit and that there still are numerous flux interactions between the flux range in this analysis and below the detection limit. We also investigate occurrence frequency distributions of flux decrease during cancellation. We found a relatively strong dependence, 2.48$\pm$0:26 as a power-law index. This strong dependence on flux is consistent with the model, which is suggested in the previous paper. △ Less

Submitted 15 October, 2015; originally announced October 2015.

Comments: 38 pages, 11 figures, accepted for ApJ

Effect of coronal temperature on the scale of solar chromospheric jets

Authors: H. Iijima, T. Yokoyama

Abstract: We investigate the effect of coronal temperature on the formation process of solar chromospheric jets using two-dimensional magnetohydrodynamic simulations of the region from the upper convection zone to the lower corona. We develop a new radiative magnetohydrodynamic code for the dynamic modeling of the solar atmosphere, employing a LTE equation of state, optically thick radiative loss in the pho… ▽ More We investigate the effect of coronal temperature on the formation process of solar chromospheric jets using two-dimensional magnetohydrodynamic simulations of the region from the upper convection zone to the lower corona. We develop a new radiative magnetohydrodynamic code for the dynamic modeling of the solar atmosphere, employing a LTE equation of state, optically thick radiative loss in the photosphere, optically thin radiative loss in the chromosphere and the corona, and thermal conduction along the magnetic field lines. Many chromospheric jets are produced in the simulations by shock waves passing through the transition region. We find that these jets are projected farther outward when the coronal temperature is lower (similar to that in coronal holes) and shorter when the coronal temperature is higher (similar to that in active regions). When the coronal temperature is high, the deceleration of the chromospheric jets is consistent with the model in which deceleration is determined by the periodic chromospheric shock waves. However, when the coronal temperature is low, the gravitational deceleration becomes more important and the chromospheric jets approach ballistic motion. △ Less

Submitted 22 September, 2015; originally announced September 2015.

Comments: 19 pages, 4 figures, accepted for publication in ApJL

Apparent cross-field superslow propagation of magnetohydrodynamic waves in solar plasmas

Authors: T. Kaneko, M. Goossens, R. Soler, J. Terradas, T. Van Doorsselaere, T. Yokoyama, A. N. Wright

Abstract: In this paper we show that the phase mixing of continuum Alfvén waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent b… ▽ More In this paper we show that the phase mixing of continuum Alfvén waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfvén and slow spectra in 2D Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfvén waves and continuum slow waves that naturally live on those structures and phase mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfvén/slow frequency across the magnetic surfaces and is slower than the Alfvén/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations. △ Less

Submitted 10 September, 2015; originally announced September 2015.

Comments: 41 pages, 7 figures, accepted for publication in ApJ

Three-Dimensional MHD Magnetic Reconnection Simulations with Finite Guide Field: Proposal of the Shock-Evoking Positive-Feedback Model

Authors: Shuoyang Wang, Takaaki Yokoyama, Hiroaki Isobe

Abstract: Using a three-dimensional magnetohydrodynamic model, we simulate the magnetic reconnection in a single current sheet. We assume a finite guide field, a random perturbation on the velocity field and uniform resistivity. Our model enhances the reconnection rate relative to the classical Sweet-Parker model in the same configuration. The efficiency of magnetic energy conversion is increased by interac… ▽ More Using a three-dimensional magnetohydrodynamic model, we simulate the magnetic reconnection in a single current sheet. We assume a finite guide field, a random perturbation on the velocity field and uniform resistivity. Our model enhances the reconnection rate relative to the classical Sweet-Parker model in the same configuration. The efficiency of magnetic energy conversion is increased by interactions between the multiple tearing layers coexisting in the global current sheet. This interaction, which forms a positive-feedback system, arises from coupling of the inflow and outflow regions in different layers across the current sheet. The coupling accelerates the elementary reconnection events, thereby enhancing the global reconnection rate. The reconnection establishes flux tubes along each tearing layer. Slow-mode shocks gradually form along the outer boundaries of these tubes, further accelerating the magnetic energy conversion.Such positive-feedback system is absent in two-dimensional simulation, three-dimensional reconnection without a guide field and a reconnection under a single perturbation mode. We refer to our model as the "shock-evoking positive-feedback" model. △ Less

Submitted 25 August, 2015; v1 submitted 13 August, 2015; originally announced August 2015.

Comments: 36 pages, 18 figures, accepted for publication in ApJ

Resonant Absorption of Transverse Oscillations and Associated Heating in a Solar Prominence. II- Numerical aspects

Authors: Patrick Antolin, Takenori J. Okamoto, Bart De Pontieu, Han Uitenbroek, Tom Van Doorsselaere, Takaaki Yokoyama

Abstract: Transverse magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere and may be responsible for generating the Sun's million-degree outer atmosphere. However, direct evidence of the dissipation process and heating from these waves remains elusive. Through advanced numerical simulations combined with appropriate forward modeling of a prominence flux tube, we provide the observational s… ▽ More Transverse magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere and may be responsible for generating the Sun's million-degree outer atmosphere. However, direct evidence of the dissipation process and heating from these waves remains elusive. Through advanced numerical simulations combined with appropriate forward modeling of a prominence flux tube, we provide the observational signatures of transverse MHD waves in prominence plasmas. We show that these signatures are characterized by thread-like substructure, strong transverse dynamical coherence, an out-of-phase difference between plane-of-the-sky motions and LOS velocities, and enhanced line broadening and heating around most of the flux tube. A complex combination between resonant absorption and Kelvin-Helmholtz instabilities (KHI) takes place in which the KHI extracts the energy from the resonant layer and dissipates it through vortices and current sheets, which rapidly degenerate into turbulence. An inward enlargement of the boundary is produced in which the turbulent flows conserve the characteristic dynamics from the resonance, therefore guaranteeing detectability of the resonance imprints. We show that the features described in the accompanying paper (Okamoto et al. 2015) through coordinated Hinode and IRIS observations match well the numerical results. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Comments: This is part 2 of a series of 2 papers. Part 1 corresponds to Okamoto et al. (2015, accepted). 36 Pages (single column), 10 figures. Accepted for publication in The Astrophysical Journal

Resonant Absorption of Transverse Oscillations and Associated Heating in a Solar Prominence. I- Observational aspects

Authors: Takenori J. Okamoto, Patrick Antolin, Bart De Pontieu, Han Uitenbroek, Tom Van Doorsselaere, Takaaki Yokoyama

Abstract: Transverse magnetohydrodynamic (MHD) waves have been shown to be ubiquitous in the solar atmosphere and can in principle carry sufficient energy to generate and maintain the Sun's million-degree outer atmosphere or corona. However, direct evidence of the dissipation process of these waves and subsequent heating has not yet been directly observed. Here we report on high spatial, temporal, and spect… ▽ More Transverse magnetohydrodynamic (MHD) waves have been shown to be ubiquitous in the solar atmosphere and can in principle carry sufficient energy to generate and maintain the Sun's million-degree outer atmosphere or corona. However, direct evidence of the dissipation process of these waves and subsequent heating has not yet been directly observed. Here we report on high spatial, temporal, and spectral resolution observations of a solar prominence that show a compelling signature of so-called resonant absorption, a long hypothesized mechanism to efficiently convert and dissipate transverse wave energy into heat. Aside from coherence in the transverse direction, our observations show telltale phase differences around 180 degrees between transverse motions in the plane-of-sky and line-of-sight velocities of the oscillating fine structures or threads, and also suggest significant heating from chromospheric to higher temperatures. Comparison with advanced numerical simulations support a scenario in which transverse oscillations trigger a Kelvin-Helmholtz instability (KHI) at the boundaries of oscillating threads via resonant absorption. This instability leads to numerous thin current sheets in which wave energy is dissipated and plasma is heated. Our results provide direct evidence for wave-related heating in action, one of the candidate coronal heating mechanisms. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Comments: 28 pages, 9 figures, accepted for publication in ApJ. Part II by Patrick Antolin et al. will appear soon

Numerical Study on In-Situ Prominence Formation by Radiative Condensation in the Solar Corona

Authors: Takafumi Kaneko, Takaaki Yokoyama

Abstract: We propose an in-situ formation model for inverse-polarity solar prominence and demonstrate it using self-consistent 2.5-dimensional magnetohydrodynamics simulations, including thermal conduction along magnetic fields and optically thin radiative cooling. The model enables us to form cool dense plasma clouds inside a flux rope by radiative condensation, which is regarded as an inverse-polarity pro… ▽ More We propose an in-situ formation model for inverse-polarity solar prominence and demonstrate it using self-consistent 2.5-dimensional magnetohydrodynamics simulations, including thermal conduction along magnetic fields and optically thin radiative cooling. The model enables us to form cool dense plasma clouds inside a flux rope by radiative condensation, which is regarded as an inverse-polarity prominence. Radiative condensation is triggered by changes in the magnetic topology, i.e., formation of the flux rope from the sheared arcade field, and by thermal imbalance due to the dense plasma trapped inside the flux rope. The flux rope is created by imposing converging and shearing motion on the arcade field. Either when the footpoint motion is in the anti-shearing direction or when heating is proportional to local density, the thermal state inside the flux rope becomes cooling-dominant, leading to radiative condensation. By controlling the temperature of condensation, we investigate the relationship between the temperature and density of prominence and derive a scaling formula for this relationship. This formula suggests that the proposed model reproduce the observed density of prominence, which is 10-100 times larger than the coronal density. Moreover, the time evolution of the extreme ultraviolet emission synthesized by combining our simulation results with the response function of the Solar Dynamics Observatory Atmospheric Imaging Assembly filters agrees with the observed temporal and spatial intensity shift among multi-wavelength during in-situ condensation. △ Less

Submitted 30 April, 2015; v1 submitted 30 April, 2015; originally announced April 2015.

Comments: 31 pages, 12 figures; accepted for publicaion in The Astrophysical Journal

Electron density of active region outflows measured by the EUV Imaging Spectrometer onboard Hinode

Authors: Naomasa Kitagawa, Takaaki Yokoyama

Abstract: In order to better understand the nature of active region outflows, the electron density was measured by using a density-sensitive line pair Fe xiv 264.78A/274.20A.Since coronal line profiles of the outflow region are composed of a major component with a Doppler shift of < 10 km s^-1 and a minor component (enhanced blue wing: EBW) blueshifted by up to 100 km s^-1, we extracted EBW from the line pr… ▽ More In order to better understand the nature of active region outflows, the electron density was measured by using a density-sensitive line pair Fe xiv 264.78A/274.20A.Since coronal line profiles of the outflow region are composed of a major component with a Doppler shift of < 10 km s^-1 and a minor component (enhanced blue wing: EBW) blueshifted by up to 100 km s^-1, we extracted EBW from the line profiles through double-Gaussian fitting. We tried applying the simultaneous fitting to those two Fe xiv lines with several physical restrictions. Electron density for both components (n_Major and n_EBW, respectively) was calculated by referring to the theoretical intensity ratio as a function of electron density as per the CHIANTI database. We studied six locations in the outflow regions around NOAA AR10978. The average electron density was n_Major = 10^(9.16 +- 0.16) cm^-3 and n_EBW = 10^(8.74 +- 0.29) cm^-3. The magnitude relationship between n_Major and n_EBW was opposite in the eastern and western outflow regions. The column depth was also calculated for each component, which leads to the result that the outflows possess only a small fraction (~ 0.1) in the eastern region, while they dominate over the major component in the line profiles by a factor of five in the western region. When taking into account the extending coronal structures, the western region can be thought to represent the mass leakage. In contrast, we suggest a possibility that the eastern region actually contributes to the mass supply to coronal loops. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: 31 pages, 15 figures

Efficient small-scale dynamo in solar convection zone

Authors: H. Hotta, M. Rempel, T. Yokoyama

Abstract: We investigate small-scale dynamo action in the solar convection zone through a series of high resolution MHD simulations in a local Cartesian domain with 1$R_\odot$ (solar radius) of horizontal extent and a radial extent from 0.715 to 0.96$R_\odot$. The dependence of the solution on resolution and diffusivity is studied. For a grid spacing of less than 350 km, the root mean square magnetic field… ▽ More We investigate small-scale dynamo action in the solar convection zone through a series of high resolution MHD simulations in a local Cartesian domain with 1$R_\odot$ (solar radius) of horizontal extent and a radial extent from 0.715 to 0.96$R_\odot$. The dependence of the solution on resolution and diffusivity is studied. For a grid spacing of less than 350 km, the root mean square magnetic field strength near the base of the convection zone reaches 95% of the equipartition field strength (i.e. magnetic and kinetic energy are comparable). For these solutions the Lorentz force feedback on the convection velocity is found to be significant. The velocity near the base of the convection zone is reduced to 50% of the hydrodynamic one. In spite of a significant decrease of the convection velocity, the reduction in the enthalpy flux is relatively small, since the magnetic field also suppresses the horizontal mixing of the entropy between up- and downflow regions. This effect increases the amplitude of the entropy perturbation and makes convective energy transport more efficient. We discuss potential implications of these results for solar global convection and dynamo simulations. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: 46 pages, 25 figures, 1 table, accepted by ApJ

High-resolution calculation of the solar global convection with the reduced speed of sound technique: II. Near surface shear layer with the rotation

Authors: H. Hotta, M. Rempel, T. Yokoyama

Abstract: We present a high-resolution, highly stratified numerical simulation of rotating thermal convection in a spherical shell. Our aim is to study in detail the processes that can maintain a near surface shear layer (NSSL) as inferred from helioseismology. Using the reduced speed of sound technique we can extend our global convection simulation to $0.99\,R_{\odot}$ and include near the top of our domai… ▽ More We present a high-resolution, highly stratified numerical simulation of rotating thermal convection in a spherical shell. Our aim is to study in detail the processes that can maintain a near surface shear layer (NSSL) as inferred from helioseismology. Using the reduced speed of sound technique we can extend our global convection simulation to $0.99\,R_{\odot}$ and include near the top of our domain small-scale convection with short time scales that is only weakly influenced by rotation. We find the formation of a NSSL preferentially in high latitudes in the depth range $r=0.95-0.975R_\odot$. The maintenance mechanisms are summarized as follows. Convection under weak influence of rotation leads to Reynolds stresses that transport angular momentum radially inward in all latitudes. This leads to the formation of a strong poleward directed meridional flow and a NSSL, which is balanced in the meridional plane by forces resulting from the $\langle v'_r v'_θ\rangle$ correlation of turbulent velocities. The origin of the required correlations depends to some degree on latitude. In high latitudes a positive correlation $\langle v'_rv'_θ\rangle$ is induced in the NSSL by the poleward meridional flow whose amplitude increases with the radius, while a negative correlation is generated by the Coriolis force in bulk of the convection zone. In low latitudes a positive correlation $\langle v'_rv'_θ\rangle$ results from rotationally aligned convection cells ("banana cells"). The force caused by these Reynolds stresses is in balance with the Coriolis force in the NSSL. △ Less

Submitted 26 October, 2014; originally announced October 2014.

Comments: 49 pages, 20 figures, accepted by ApJ

Simulation Study of Solar Plasma Eruptions Caused by Interactions between Emerging Flux and Coronal Arcade Fields

Authors: T. Kaneko, T. Yokoyama

Abstract: We investigate the triggering mechanisms of plasma eruptions in the solar atmosphere due to interactions between emerging flux and coronal arcade fields by using two-dimensional MHD simulations. We perform parameter surveys with respect to arcade field height, magnetic field strength, and emerging flux location. Our results show that two possible mechanisms exist, and which mechanism is dominant d… ▽ More We investigate the triggering mechanisms of plasma eruptions in the solar atmosphere due to interactions between emerging flux and coronal arcade fields by using two-dimensional MHD simulations. We perform parameter surveys with respect to arcade field height, magnetic field strength, and emerging flux location. Our results show that two possible mechanisms exist, and which mechanism is dominant depends mostly on emerging flux location. One mechanism appears when the location of emerging flux is close to the polarity inversion line (PIL) of an arcade field. This mechanism requires reconnection between the emerging flux and the arcade field, as pointed out by previous studies. The other mechanism appears when the location of emerging flux is around the edge of an arcade field. This mechanism does not require reconnection between the emerging flux and the arcade field but does demand reconnection in the arcade field above the PIL. Furthermore, we found that the eruptive condition for this mechanism can be represented by a simple formula. △ Less

Submitted 1 October, 2014; originally announced October 2014.

Comments: 40 pages, 13 figures, accepted for publication in ApJ

Statistical Analysis of the Horizontal Divergent Flow in Emerging Solar Active Regions

Authors: Shin Toriumi, Keiji Hayashi, Takaaki Yokoyama

Abstract: Solar active regions (ARs) are thought to be formed by magnetic fields from the convection zone. Our flux emergence simulations revealed that a strong horizontal divergent flow (HDF) of unmagnetized plasma appears at the photosphere before the flux begins to emerge. In our earlier study, we analyzed HMI data for a single AR and confirmed presence of this precursor plasma flow in the actual Sun. In… ▽ More Solar active regions (ARs) are thought to be formed by magnetic fields from the convection zone. Our flux emergence simulations revealed that a strong horizontal divergent flow (HDF) of unmagnetized plasma appears at the photosphere before the flux begins to emerge. In our earlier study, we analyzed HMI data for a single AR and confirmed presence of this precursor plasma flow in the actual Sun. In this paper, as an extension of our earlier study, we conducted a statistical analysis of the HDFs to further investigate their characteristics and better determine the properties. From SDO/HMI data, we picked up 23 flux emergence events over a period of 14 months, the total flux of which ranges from 10^{20} to 10^{22} Mx. Out of 23 selected events, 6 clear HDFs were detected by the method we developed in our earlier study, and 7 HDFs detected by visual inspection were added to this statistic analysis. We found that the duration of the HDF is on average 61 minutes and the maximum HDF speed is on average 3.1 km s^{-1}. We also estimated the rising speed of the subsurface magnetic flux to be 0.6-1.4 km s^{-1}. These values are highly consistent with our previous one-event analysis as well as our simulation results. The observation results lead us to the conclusion that the HDF is rather a common feature in the earliest phase of AR emergence. Moreover, our HDF analysis has capability of determining the subsurface properties of emerging fields that cannot be directly measured. △ Less

Submitted 11 August, 2014; originally announced August 2014.

Comments: 40 pages, 8 figures, 3 tables, accepted for publication in ApJ. A high-quality version can be found at http://www-space.eps.s.u-tokyo.ac.jp/~toriumi/toriumi2014apj.pdf

Orbital Perturbations of the Galilean Satellites During Planetary Encounters

Authors: R. Deienno, D. Nesvorny, D. Vokrouhlicky, T. Yokoyama

Abstract: The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present Solar System, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here we study the surviva… ▽ More The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present Solar System, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here we study the survival of the Galilean satellites in the jumping-Jupiter model. This is an important concern because the ice-giant encounters, if deep enough, could dynamically perturb the orbits of the Galilean satellites, and lead to implausible results. We performed numerical integrations where we tracked the effect of planetary encounters on the Galilean moons. We considered three instability cases from Nesvorny & Morbidelli (2012) that differed in the number and distribution of encounters. We found that in one case, where the number of close encounters was relatively small, the Galilean satellite orbits were not significantly affected. In the other two, the orbital eccentricities of all moons were excited by encounters, Callisto's semimajor axis changed, and, in a large fraction of trials, the Laplace resonance of the inner three moons was disrupted. The subsequent evolution by tides damps eccentricities and can recapture the moons in the Laplace resonance. A more important constraint is represented by the orbital inclinations of the moons, which can be excited during the encounters and not appreciably damped by tides. We find that one instability case taken from Nesvorny & Morbidelli (2012) clearly fails this constraint. This shows how the regular satellites of Jupiter can be used to set limits on the properties of encounters in the jumping-Jupiter model, and help us to better understand how the early Solar System evolved. △ Less

Submitted 8 May, 2014; originally announced May 2014.

Comments: The Astronomical Journal, in press

Fine strand-like structure in the solar corona from MHD transverse oscillations

Authors: P. Antolin, T. Yokoyama, T. Van Doorsselaere

Abstract: Current analytical and numerical modelling suggest the existence of ubiquitous thin current sheets in the corona that could explain the observed heating requirements. On the other hand, new high resolution observations of the corona indicate that its magnetic field may tend to organise itself in fine strand-like structures of few hundred kilometres widths. The link between small structure in model… ▽ More Current analytical and numerical modelling suggest the existence of ubiquitous thin current sheets in the corona that could explain the observed heating requirements. On the other hand, new high resolution observations of the corona indicate that its magnetic field may tend to organise itself in fine strand-like structures of few hundred kilometres widths. The link between small structure in models and the observed widths of strand-like structure several orders of magnitude larger is still not clear. A popular theoretical scenario is the nanoflare model, in which each strand is the product of an ensemble of heating events. Here, we suggest an alternative mechanism for strand generation. Through forward modelling of 3D MHD simulations we show that small amplitude transverse MHD waves can lead in a few periods time to strand-like structure in loops in EUV intensity images. Our model is based on previous numerical work showing that transverse MHD oscillations can lead to Kelvin-Helmholtz instabilities that deform the cross-sectional area of loops. While previous work has focused on large amplitude oscillations, here we show that the instability can occur even for low wave amplitudes for long and thin loops, matching those presently observed in the corona. We show that the vortices generated from the instability are velocity sheared regions with enhanced emissivity hosting current sheets. Strands result as a complex combination of the vortices and the line-of-sight angle, last for timescales of a period and can be observed for spatial resolutions of a tenth of loop radius. △ Less

Submitted 30 April, 2014; originally announced May 2014.

Comments: Accepted for publication in ApJL

High-resolution calculations of the solar global convection with the reduced speed of sound technique: I. The structure of the convection and the magnetic field without the rotation

Authors: H. Hotta, M. Rempel, T. Yokoyama

Abstract: We carry out non-rotating high-resolution calculations of the solar global convection, which resolve convective scales of less than 10 Mm. To cope with the low Mach number conditions in the lower convection zone, we use the reduced speed of sound technique (RSST), which is simple to implement and requires only local communication in the parallel computation. In addition, the RSST allows us to expa… ▽ More We carry out non-rotating high-resolution calculations of the solar global convection, which resolve convective scales of less than 10 Mm. To cope with the low Mach number conditions in the lower convection zone, we use the reduced speed of sound technique (RSST), which is simple to implement and requires only local communication in the parallel computation. In addition, the RSST allows us to expand the computational domain upward to about $0.99 R_{\odot}$ as it can also handle compressible flows. Using this approach, we study the solar convection zone on the global scale, including small-scale near-surface convection. In particular, we investigate the influence of the top boundary condition on the convective structure throughout the convection zone as well as on small-scale dynamo action. Our main conclusions are: 1. The small-scale downflows generated in the near-surface layer penetrate into deeper layers to some extent and excite small-scale turbulence in the region $>0.9R_\odot$, where $R_\odot$ is the solar radius. 2. In the deeper convection zone ($<0.9R_\odot$), the convection is not influenced by the location of the upper boundary. 3. Using an LES approach, we can achieve small-scale dynamo action and maintain a field of about $0.15-0.25 B_\mathrm{eq}$ throughout the convection zone, where $B_\mathrm{eq}$ is the equipartition magnetic field to the kinetic energy. 4. The overall dynamo efficiency varies significantly in the convection zone as a consequence of the downward directed Poynting flux and the depth variation of the intrinsic convective scales. △ Less

Submitted 20 February, 2014; originally announced February 2014.

Comments: 60 pages, 24 figures, accepted for publication in ApJ

Probing the Shallow Convection Zone: Rising Motion of Subsurface Magnetic Fields in the Solar Active Region

Authors: Shin Toriumi, Stathis Ilonidis, Takashi Sekii, Takaaki Yokoyama

Abstract: In this Letter we present a seismological detection of a rising motion of magnetic flux in the shallow convection zone of the Sun, and show estimates of the emerging speed and its decelerating nature. In order to evaluate the speed of subsurface flux that creates an active region, we apply six Fourier filters to the Doppler data of NOAA AR 10488, observed with SOHO/MDI, to detect the reduction of… ▽ More In this Letter we present a seismological detection of a rising motion of magnetic flux in the shallow convection zone of the Sun, and show estimates of the emerging speed and its decelerating nature. In order to evaluate the speed of subsurface flux that creates an active region, we apply six Fourier filters to the Doppler data of NOAA AR 10488, observed with SOHO/MDI, to detect the reduction of acoustic power at six different depths from -15 to -2 Mm. All the filtered acoustic powers show reductions, up to 2 hours before the magnetic flux first appears at the visible surface. The start times of these reductions show a rising trend with a gradual deceleration. The obtained velocity is first several km s^-1 in a depth range of 15--10 Mm, then ~1.5 km s^-1 at 10-5 Mm, finally ~0.5 km s^-1 at 5-2 Mm. If we assume that the power reduction is actually caused by the magnetic field, the velocity of order of 1 km s^-1 is well in accordance with previous observations and numerical studies. Moreover, the gradual deceleration strongly supports the theoretical model that the emerging flux slows down in the uppermost convection zone before it expands into the atmosphere to build an active region. △ Less

Submitted 14 May, 2013; originally announced May 2013.

Comments: 19 pages, 4 figures, accepted for publication in ApJL

Three-dimensional magnetohydrodynamic simulation of the solar magnetic flux emergence: Parametric study on the horizontal divergent flow

Authors: Shin Toriumi, Takaaki Yokoyama

Abstract: Solar active regions are formed through the emergence of magnetic flux from the deeper convection zone. Recent satellite observations have shown that a horizontal divergent flow (HDF) stretches out over the solar surface just before the magnetic flux appearance. The aims of this study are to investigate the driver of the HDF and to see the dependency of the HDF on the parameters of the magnetic fl… ▽ More Solar active regions are formed through the emergence of magnetic flux from the deeper convection zone. Recent satellite observations have shown that a horizontal divergent flow (HDF) stretches out over the solar surface just before the magnetic flux appearance. The aims of this study are to investigate the driver of the HDF and to see the dependency of the HDF on the parameters of the magnetic flux in the convection zone. We conduct three-dimensional magnetohydrodynamic (3D MHD) numerical simulations of the magnetic flux emergence and vary the parameters in the initial conditions. An analytical approach is also taken to explain the dependency. The horizontal gas pressure gradient is found to be the main driver of the HDF. The maximum HDF speed shows positive correlations with the field strength and twist intensity. The HDF duration has a weak relation with the twist, while it shows negative dependency on the field strength only in the case of the stronger field regime. Parametric dependencies analyzed in this study may allow us to probe the structure of the subsurface magnetic flux by observing properties of the HDF. △ Less

Submitted 19 March, 2013; originally announced March 2013.

Comments: 7 pages, 5 figures, 1 table, accepted for publication in Astronomy & Astrophysics

Temporal and Spatial Analyses of Spectral Indices of Nonthermal Emissions Derived from Hard X-Rays and Microwaves

Authors: Ayumi Asai, Junko Kiyohara, Hiroyuki Takasaki, Noriyuki Narukage, Takaaki Yokoyama, Satoshi Masuda, Masumi Shimojo, Hiroshi Nakajima

Abstract: We studied electron spectral indices of nonthermal emissions seen in hard X-rays (HXRs) and in microwaves. We analyzed 12 flares observed by the Hard X-ray Telescope aboard {\it Yohkoh}, Nobeyama Radio Polarimeters (NoRP), and the Nobeyama Radioheliograph (NoRH), and compared the spectral indices derived from total fluxes of hard X-rays and microwaves. Except for four events, which have very soft… ▽ More We studied electron spectral indices of nonthermal emissions seen in hard X-rays (HXRs) and in microwaves. We analyzed 12 flares observed by the Hard X-ray Telescope aboard {\it Yohkoh}, Nobeyama Radio Polarimeters (NoRP), and the Nobeyama Radioheliograph (NoRH), and compared the spectral indices derived from total fluxes of hard X-rays and microwaves. Except for four events, which have very soft HXR spectra suffering from the thermal component, these flares show a gap $Δδ$ between the electron spectral indices derived from hard X-rays $δ_{X}$ and those from microwaves $δ_μ$ ($Δδ= δ_{X} - δ_μ$) of about 1.6. Furthermore, from the start to the peak times of the HXR bursts, the time profiles of the HXR spectral index $δ_{X}$ evolve synchronously with those of the microwave spectral index $δ_μ$, keeping the constant gap. We also examined the spatially resolved distribution of the microwave spectral index by using NoRH data. The microwave spectral index $δ_μ$ tends to be larger, which means a softer spectrum, at HXR footpoint sources with stronger magnetic field than that at the loop tops. These results suggest that the electron spectra are bent at around several hundreds of keV, and become harder at the higher energy range that contributes the microwave gyrosynchrotron emission. △ Less

Submitted 8 December, 2012; originally announced December 2012.

Comments: 24 pages, 6 figures, accepted for publication in ApJ