Showing 1–37 of 37 results for author: Hindman, B W

Inertial Waves in a Nonlinear Simulation of the Sun's Convection Zone and Radiative Interior

Authors: Catherine C. Blume, Bradley W. Hindman, Loren I. Matilsky

Abstract: Recent observations of Rossby waves and other more exotic forms of inertial oscillations in the Sun's convection zone have kindled the hope that such waves might be used as a seismic probe of the Sun's interior. Here we present a 3D numerical simulation in spherical geometry that models the Sun's convection zone and upper radiative interior. This model features a wide variety of inertial oscillati… ▽ More Recent observations of Rossby waves and other more exotic forms of inertial oscillations in the Sun's convection zone have kindled the hope that such waves might be used as a seismic probe of the Sun's interior. Here we present a 3D numerical simulation in spherical geometry that models the Sun's convection zone and upper radiative interior. This model features a wide variety of inertial oscillations, including both sectoral and tesseral equatorial Rossby waves, retrograde mixed inertial modes, prograde thermal Rossby waves, the recently observed high-frequency retrograde (HFR) vorticity modes, and what may be latitudinal overtones of these HFR modes. With this model, we demonstrate that sectoral and tesseral Rossby waves are ubiquitous within the radiative interior as well as within the convection zone. We suggest that there are two different Rossby-wave families in this simulation that live in different wave cavities: one in the radiative interior and one in the convection zone. Finally, we suggest that many of the retrograde inertial waves that appear in the convection zone, including the HFR modes, are in fact all related, being latitudinal overtones that are mixed modes with the prograde thermal Rossby waves. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 23 pages, 14 figures, submitted to The Astrophysical Journal

Confinement of the Solar Tachocline by a Non-Axisymmetric Dynamo

Authors: Loren I. Matilsky, Nicholas H. Brummell, Bradley W. Hindman, Juri Toomre

Abstract: We recently presented the first 3D numerical simulation of the solar interior for which tachocline confinement was achieved by a dynamo-generated magnetic field. In this followup study, we analyze the degree of confinement as the magnetic field strength changes (controlled by varying the magnetic Prandtl number) in a coupled radiative zone (RZ) and convection zone (CZ) system. We broadly find thre… ▽ More We recently presented the first 3D numerical simulation of the solar interior for which tachocline confinement was achieved by a dynamo-generated magnetic field. In this followup study, we analyze the degree of confinement as the magnetic field strength changes (controlled by varying the magnetic Prandtl number) in a coupled radiative zone (RZ) and convection zone (CZ) system. We broadly find three solution regimes, corresponding to weak, medium, and strong dynamo magnetic field strengths. In the weak-field regime, the large-scale magnetic field is mostly axisymmetric with regular, periodic polarity reversals (reminiscent of the observed solar cycle), but fails to create a confined tachocline. In the strong-field regime, the large-scale field is mostly non-axisymmetric with irregular, quasi-periodic polarity reversals, and creates a confined tachocline. In the medium-field regime, the large-scale field resembles a strong-field dynamo for extended intervals, but intermittently weakens to allow temporary epochs of strong differential rotation. In all regimes, the amplitude of poloidal field strength in the RZ is very well explained by skin-depth arguments, wherein the oscillating field that gives rise to the skin depth (in the medium- and strong-field cases) is a non-axisymmetric field structure rotating with respect to the RZ. These new simulations reaffirm that tachocline confinement by the solar dynamo (the so-called fast magnetic confinement scenario) is possible, but suggest a new picture in which non-axisymmetric field components rotating with respect to the RZ play the primary role, instead of the regularly reversing axisymmetic field associated with the 22-year cycle. △ Less

Submitted 8 January, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: accepted for publication in ApJ; 29 pages, 9 figures, 4 tables, 3 appendices

Dwindling Surface Cooling of a Rotating Jovian Planet Leads to a Convection Zone that Grows to a Finite Depth

Authors: Bradley W. Hindman, J. R. Fuentes

Abstract: Recent measurements of Jupiter's gravitational field (by Juno) and seismology of Saturn's rings (by Cassini) strongly suggest that both planets have a stably-stratified core that still possesses a primordial gradient in the concentration of heavy elements. The existence of such a "diffusely" stratified core has been a surprise as it was long expected that the Jovian planets should be fully convect… ▽ More Recent measurements of Jupiter's gravitational field (by Juno) and seismology of Saturn's rings (by Cassini) strongly suggest that both planets have a stably-stratified core that still possesses a primordial gradient in the concentration of heavy elements. The existence of such a "diffusely" stratified core has been a surprise as it was long expected that the Jovian planets should be fully convective and hence fully mixed. A vigorous zone of convection, driven by surface cooling, forms at the surface and deepens through entrainment of fluid from underneath. In fact, it was believed that this convection zone should grow so rapidly that the entire planet would be consumed in less than a million years. Here we suggest that two processes, acting in concert, present a solution to this puzzle. All of the giant planets are rapidly rotating and have a cooling rate that declines with time. Both of these effects reduce the rate of fluid entrainment into the convection zone. Through the use of an analytic prescription of entrainment in giant planets, we demonstrate that these two effects, rotation and dwindling surface cooling, result in a convection zone which initially grows but eventually stalls. The depth to which the convective interface asymptotes depends on the rotation rate and on the stratification of the stable interior. Conversely, in a nonrotating planet, or in a planet that maintains a higher level of cooling than current models suggest, the convection zone deepens forever, eventually spanning the entire planet. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 7 pages, 2 figures, accepted for publication by Astrophysical Journal Letters

Latitudinal Propagation of Thermal Rossby Waves in Stellar Convection Zones

Authors: Rekha Jain, Bradley W. Hindman

Abstract: Using an analytic model, we derive the eigenfrequencies for thermal Rossby waves that are trapped radially and latitudinally in an isentropically stratified atmosphere. We ignore the star's curvature and work in an equatorial f-plane geometry. The propagation of inertial waves is found to be sensitive to the relative direction of the wave vector to the zonal direction. Prograde propagating thermal… ▽ More Using an analytic model, we derive the eigenfrequencies for thermal Rossby waves that are trapped radially and latitudinally in an isentropically stratified atmosphere. We ignore the star's curvature and work in an equatorial f-plane geometry. The propagation of inertial waves is found to be sensitive to the relative direction of the wave vector to the zonal direction. Prograde propagating thermal Rossby waves are naturally trapped in the radial direction for frequencies above a critical threshold, which depends on the angle of propagation. Below the threshold frequency, there exists a continuous spectrum of prograde and retrograde inertial waves that are untrapped in an isentropic atmosphere, but can be trapped by gradients in the specific entropy density such as occurs in a stellar convection zone. Finally, we discuss the implications of these waves on recent observations of inertial oscillations in the Sun, as well as in numerical simulations. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: 14 pages, 5 figures, accepted by ApJ

Low-Frequency Internal Gravity Waves are Pseudo-incompressible

Authors: Bradley W. Hindman, Keith Julien

Abstract: Starting from the fully compressible fluid equations in a plane-parallel atmosphere, we demonstrate that linear internal gravity waves are naturally pseudo-incompressible in the limit that the wave frequency $ω$ is much less than that of surface gravity waves, i.e., $ω\ll \sqrt{g k_h}$ where $g$ is the gravitational acceleration and $k_h$ is the horizontal wavenumber. We accomplish this by perform… ▽ More Starting from the fully compressible fluid equations in a plane-parallel atmosphere, we demonstrate that linear internal gravity waves are naturally pseudo-incompressible in the limit that the wave frequency $ω$ is much less than that of surface gravity waves, i.e., $ω\ll \sqrt{g k_h}$ where $g$ is the gravitational acceleration and $k_h$ is the horizontal wavenumber. We accomplish this by performing a formal expansion of the wave functions and the local dispersion relation in terms of a dimensionless frequency $\varepsilon = ω/ \sqrt{g k_h}$. Further, we show that in this same low-frequency limit, several forms of the anelastic approximation, including the Lantz-Braginsky-Roberts (LBR) formulation, poorly reproduce the correct behavior of internal gravity waves. The pseudo-incompressible approximation is achieved by assuming that Eulerian fluctuations of the pressure are small in the continuity equation. Whereas, in the anelastic approximation Eulerian density fluctuations are ignored. In an adiabatic stratification, such as occurs in a convection zone, the two approximations become identical. But, in a stable stratification, the differences between the two approximations are stark and only the pseudo-incompressible approximation remains valid. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 18 pages, 2 figures, submitted to the Astrophysical Journal

Journal ref: ApJ, 960, 64 (2024)

Assessing the Observability of Deep Meridional Flow Cells in the Solar Interior

Authors: J. R. Fuentes, Bradley W. Hindman, Junwei Zhao, Catherine C. Blume, Maria E. Camisassa, Nicholas A. Featherstone, Thomas Hartlep, Lydia Korre, Loren I. Matilsky

Abstract: Meridional circulation regulates the Sun's interior dynamics and magnetism. While it is well accepted that meridional flows are poleward at the Sun's surface, helioseismic observations have yet to provide a definitive answer for the depth at which those flows return to the equator, or the number of circulation cells in depth. Here, we explore the observability of multiple circulation cells stacked… ▽ More Meridional circulation regulates the Sun's interior dynamics and magnetism. While it is well accepted that meridional flows are poleward at the Sun's surface, helioseismic observations have yet to provide a definitive answer for the depth at which those flows return to the equator, or the number of circulation cells in depth. Here, we explore the observability of multiple circulation cells stacked in radius. Specifically, we examine the seismic signature of several meridional flow profiles by convolving time-distance averaging kernels with mean flows obtained from a suite of 3D hydrodynamic simulations. At mid and high latitudes, we find that weak flow structures in the deep convection zone can be obscured by signals from the much stronger surface flows. This contamination of 1--2 m s$^{-1}$ is caused by extended side lobes in the averaging kernels, which produce a spurious equatorward signal with flow speeds that are one order of magnitude stronger than the original flow speeds in the simulations. At low latitudes, the flows in the deep layers of the simulations are stronger ($> 2$ m s$^{-1}$) and multiple cells across the convection zone can produce a sufficiently strong signal to survive the convolution process. Now that meridional flows can be measured over two decades of data, the uncertainties arising from convective noise have fallen to a level where they are comparable in magnitude to the systematic biases caused by non-local features in the averaging kernels. Hence, these systematic errors are beginning to influence current helioseismic deductions and need broader consideration. △ Less

Submitted 4 December, 2023; v1 submitted 14 August, 2023; originally announced August 2023.

Comments: New version including suggestions by the reviewers

Report number: ApJ 961 78

Rotation reduces convective mixing in Jupiter and other gas giants

Authors: J. R. Fuentes, Evan H. Anders, Andrew Cumming, Bradley W. Hindman

Abstract: Recent measurements of Jupiter's gravitational moments by the Juno spacecraft and seismology of Saturn's rings suggest that the primordial composition gradients in the deep interior of these planets have persisted since their formation. One possible explanation is the presence of a double-diffusive staircase below the planet's outer convection zone, which inhibits mixing across the deeper layers.… ▽ More Recent measurements of Jupiter's gravitational moments by the Juno spacecraft and seismology of Saturn's rings suggest that the primordial composition gradients in the deep interior of these planets have persisted since their formation. One possible explanation is the presence of a double-diffusive staircase below the planet's outer convection zone, which inhibits mixing across the deeper layers. However, hydrodynamic simulations have shown that these staircases are not long-lasting and can be disrupted by overshooting convection. In this paper we suggests that planetary rotation could be another factor for the longevity of primordial composition gradients. Using rotational mixing-length theory and 3D hydrodynamic simulations, we demonstrate that rotation significantly reduces both the convective velocity and the mixing of primordial composition gradients. In particular, for Jovian conditions at $t\sim10^{8}~\mathrm{yrs}$ after formation, rotation reduces the convective velocity by a factor of 6, and in turn, the kinetic energy flux available for mixing gets reduced by a factor of $6^3\sim 200$. This leads to an entrainment timescale that is more than two orders of magnitude longer than without rotation. We encourage future hydrodynamic models of Jupiter and other gas giants to include rapid rotation, because the decrease in the mixing efficiency could explain why Jupiter and Saturn are not fully mixed. △ Less

Submitted 22 May, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Astrophysical Journal Letters

Journal ref: ApJL 950 L4 (2023)

The Puzzling Structure of Solar Convection: Window into the Dynamo

Authors: Nicholas A. Featherstone, Kyle C. Augustson, Jonathan M. Aurnou, Catherine Blume, Benjamin P. Brown, Nicholas Brummell, Keaton J. Burns, Michael A. Calkins, Maria Camisassa, Mausumi Dikpati, Yuhong Fan, J. R. Fuentes, Gustavo Guerrero, Bradley W. Hindman, Keith Julien, Irina N. Kitiashvili, Lydia Korre, Daniel Lecoanet, Bhishek Manek, Loren Matilsky, Mark Miesch, Nicholas J. Nelson, Jeffrey S. Oishi, Whitney T. Powers, Matthias Rempel , et al. (3 additional authors not shown)

Abstract: The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective… ▽ More The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics 2024-2033. 7 pages; 4 figures

Overstable Convective Modes in a Polytropic Stellar Atmosphere

Authors: Bradley W. Hindman, Rekha Jain

Abstract: Within the convection zone of a rotating star, the presence of the Coriolis force stabilizes long-wavelength convective modes. These modes, which would have been unstable if the star lacked rotation, are called overstable convective modes or thermal Rossby waves. We demonstrate that the Sun's rotation rate is sufficiently rapid that the lower half of its convection zone could possess overstable mo… ▽ More Within the convection zone of a rotating star, the presence of the Coriolis force stabilizes long-wavelength convective modes. These modes, which would have been unstable if the star lacked rotation, are called overstable convective modes or thermal Rossby waves. We demonstrate that the Sun's rotation rate is sufficiently rapid that the lower half of its convection zone could possess overstable modes. Further, we present an analytic solution for atmospheric waves that reside within a polytropic stratification. We explore in detail the properties of the overstable and unstable wave modes that exist when the polytrope is weakly unstable to convective overturning. Finally, we discuss how the thermal Rossby waves that reside within the convection zone of a star might couple with the prograde branch of the $g$ modes that are trapped within the star's radiative zone. We suggest that such coupling might enhance the photospheric visibility of a subset of the Sun's $g$ modes. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: 11 pages, 6 figures, published in the Astrophysical Journal (https://iopscience.iop.org/article/10.3847/1538-4357/acaec4/pdf)

Journal ref: ApJ, 943, 127 (2023)

Confinement of the Solar Tachocline by Dynamo Action in the Radiative Interior

Authors: Loren I. Matilsky, Bradley W. Hindman, Nicholas A. Featherstone, Catherine C. Blume, Juri Toomre

Abstract: A major outstanding problem in solar physics is the confinement of the solar tachocline, the thin shear layer that separates nearly solid-body rotation in the radiative interior from strong differential rotation in the convection zone. Here, we present the first 3-D, global solar simulation in which a tachocline is confined by a self-excited dynamo. The non-axisymmetric magnetism is initially buil… ▽ More A major outstanding problem in solar physics is the confinement of the solar tachocline, the thin shear layer that separates nearly solid-body rotation in the radiative interior from strong differential rotation in the convection zone. Here, we present the first 3-D, global solar simulation in which a tachocline is confined by a self-excited dynamo. The non-axisymmetric magnetism is initially built in the convection zone and then diffusively imprints downward. Additionally, the field is locally amplified throughout the radiative interior by vigorous horizontal motions that arise from equatorial Rossby waves and possibly shear instabilities. Our work thus challenges the long-held notion that the Sun's dynamo magnetic field is amplified only as deep as the tachocline and stored in a quiescent radiative interior. △ Less

Submitted 27 September, 2022; v1 submitted 26 June, 2022; originally announced June 2022.

Comments: 13 pages, 5 figures, 1 table, 2 appendices, accepted for publication in ApJL

Radial Trapping of Thermal Rossby Waves within the Convection Zones of Low-Mass Stars

Authors: Bradley W. Hindman, Rekha Jain

Abstract: We explore how thermal Rossby waves propagate within the gravitationally stratified atmosphere of a low-mass star with an outer convective envelope. Under the conditions of slow, rotationally constrained dynamics, we derive a local dispersion relation for atmospheric waves in a fully compressible stratified fluid. This dispersion relation describes the zonal and radial propagation of acoustic wave… ▽ More We explore how thermal Rossby waves propagate within the gravitationally stratified atmosphere of a low-mass star with an outer convective envelope. Under the conditions of slow, rotationally constrained dynamics, we derive a local dispersion relation for atmospheric waves in a fully compressible stratified fluid. This dispersion relation describes the zonal and radial propagation of acoustic waves and gravito-inertial waves. Thermal Rossby waves are just one class of prograde-propagating gravito-inertial wave that manifests when the buoyancy frequency is small compared to the rotation rate of the star. From this dispersion relation, we identify the radii at which waves naturally reflect and demonstrate how thermal Rossby waves can be trapped radially in a waveguide that permits free propagation in the longitudinal direction. We explore this trapping further by presenting analytic solutions for thermal Rossby waves within an isentropically stratified atmosphere that models a zone of efficient convective heat transport. We find that within such an atmosphere, waves of short zonal wavelength have a wave cavity that is radially thin and confined within the outer reaches of the convection zone near the star's equator. The same behavior is evinced by the thermal Rossby waves that appear at convective onset in numerical simulations of convection within rotating spheres. Finally, we suggest that stable thermal Rossby waves could exist in the lower portion of the Sun's convection zone, despite that region's unstable stratification. For long wavelengths, the Sun's rotation rate is sufficiently rapid to stabilize convective motions and the resulting overstable convective modes are identical to thermal Rossby waves. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 38 pages, 9 figures, accepted for publication in the Astrophysical Journal

Constraining Global Solar Models through Helioseismic Analysis

Authors: Andrey M. Stejko, Alexander G. Kosovichev, Nicholas A. Featherstone, Gustavo Guerrero, Bradley W. Hindman, Loren I. Matilsky, Jörn Warnecke

Abstract: Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this… ▽ More Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this controversy and develop a more robust understanding of global flow regimes in the solar interior, we apply a "forward-modeling" approach to the analysis of helioseismic signatures of meridional circulation profiles obtained from numerical simulations. We employ the global acoustic modeling code GALE to simulate the propagation of acoustic waves through regimes of mean mass flows generated by global hydrodynamic and magnetohydrodynamic models: EULAG, the Pencil Code, and the Rayleigh code. These models are used to create synthetic dopplergram data products, used as inputs for local time-distance helioseismology techniques. Helioseismic travel-time signals from solutions obtained through global numerical simulations are compared directly with inferences from solar observations, in order to set additional constraints on global model parameters in a direct way. We show that even though these models are able to replicate solar-like differential rotation, the resulting rotationally-constrained convection develops a multi-cell global meridional circulation profile that is measurably inconsistent with local time-distance inferences of solar observations. However, we find that the development of rotationally-unconstrained convection close to the model surface is able to maintain solar-like differential rotation, while having a significant impact on the helioseismic travel-time signal, replicating solar observations within one standard deviation of the error due to noise. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 13 pages, 7 figures

Do Coronal Loops Oscillate in Isolation?

Authors: Bradley W. Hindman, Rekha Jain

Abstract: Images of the solar corona by extreme-ultraviolet (EUV) telescopes reveal elegant arches of glowing plasma that trace the corona's magnetic field. Typically, these loops are preferentially illuminated segments of an arcade of vaulted field lines and such loops are often observed to sway in response to nearby solar flares. A flurry of observational and theoretical effort has been devoted to the exp… ▽ More Images of the solar corona by extreme-ultraviolet (EUV) telescopes reveal elegant arches of glowing plasma that trace the corona's magnetic field. Typically, these loops are preferentially illuminated segments of an arcade of vaulted field lines and such loops are often observed to sway in response to nearby solar flares. A flurry of observational and theoretical effort has been devoted to the exploitation of these oscillations with the grand hope that seismic techniques might be used as probes of the strength and structure of the corona's magnetic field. The commonly accepted viewpoint is that each visible loop oscillates as an independent entity and acts as a one-dimensional (1D) wave cavity for magnetohydrodynamic (MHD) kink waves. We argue that for many events, this generally accepted model for the wave cavity is fundamentally flawed. In particular, the 3D magnetic arcade in which the bright loop resides participates in the oscillation. Thus, the true wave cavity is larger than the individual loop and inherently multidimensional. We derive the skin depth of the near-field response for an oscillating loop and demonstrate that most loops are too close to other magnetic structures to oscillate in isolation. Further, we present a simple model of a loop embedded within an arcade and explore how the eigenmodes of the arcade and the eigenmodes of the loop become coupled. In particular, we discuss how distinguishing between these two types of modes can be difficult when the motions within the arcade are often invisible. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: 14 pages, 7 figures, an mp4 movie has been included as an ancillary file as a supplement to figure 1, Accepted for publication in the Astrophysical Journal

Revisiting the Sun's Strong Differential Rotation along Radial Lines

Authors: Loren I. Matilsky, Bradley W. Hindman, Juri Toomre

Abstract: Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the… ▽ More Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the resulting differential rotation is markedly different between the two sets. The fixed-entropy-gradient simulations have strong differential rotation contrast and isocontours tilted along radial lines (in good agreement with the Sun's interior rotation revealed by helioseismology), whereas the fixed-entropy simulations have weaker contrast and contours tilted in the opposite sense. We examine in detail the force balances in our models and find that the poleward transport of heat by Busse columns drives a thermal wind responsible for the different rotation profiles. We conclude that the Sun's strong differential rotation along radial lines may result from the solar emissivity being invariant with latitude (which is similar to the fixed-entropy-gradient condition in our models) and the poleward transport of heat by Busse columns. In future work on convection in the solar context, we strongly advise modelers to use a fixed-gradient outer boundary condition. △ Less

Submitted 16 June, 2020; v1 submitted 31 March, 2020; originally announced April 2020.

Comments: 13 pages, 15 figures, accepted for publication in the Astrophysical Journal

An improved multi-ridge fitting method for ring-diagram helioseismic analysis

Authors: Kaori Nagashima, Aaron C. Birch, Jesper Schou, Bradley W. Hindman, Laurent Gizon

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

Submitted 18 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 633, A109 (2020)

A New Analysis Procedure for Detecting Periodicities within Complex Solar Coronal Arcades

Authors: Farhad Allian, Rekha Jain, B. W Hindman

Abstract: We study intensity variations, as measured by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), in a solar coronal arcade using a newly developed analysis procedure that employs spatio-temporal auto-correlations. We test our new procedure by studying large-amplitude oscillations excited by nearby flaring activity within a complex arcade and detect a dominant per… ▽ More We study intensity variations, as measured by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), in a solar coronal arcade using a newly developed analysis procedure that employs spatio-temporal auto-correlations. We test our new procedure by studying large-amplitude oscillations excited by nearby flaring activity within a complex arcade and detect a dominant periodicity of 12.5 minutes. We compute this period in two ways: from the traditional time-distance fitting method and using our new auto-correlation procedure. The two analyses yield consistent results. The auto-correlation procedure is then implemented on time series for which the traditional method would fail due to the complexity of overlapping loops and a poor contrast between the loops and the background. Using this new procedure, we discover the presence of small-amplitude oscillations within the same arcade with 8-minute and 10-minute periods prior and subsequent to the large-amplitude oscillations, respectively. Consequently, we identify these as "decayless" oscillations that have only been previously observed in non-flaring loop systems. △ Less

Submitted 18 February, 2019; originally announced February 2019.

Comments: Submitted to APJ

Exploring the Influence of Density Contrast on Solar Near-Surface Shear

Authors: Loren I. Matilsky, Bradley W. Hindman, Juri Toomre

Abstract: The advent of helioseismology has determined in detail the average rotation rate of the Sun as a function of radius and latitude. These data immediately reveal two striking boundary layers of shear in the solar convection zone (CZ): a tachocline at the base, where the differential rotation of the CZ transitions to solid-body rotation in the radiative zone, and a 35-Mm-thick near-surface shear laye… ▽ More The advent of helioseismology has determined in detail the average rotation rate of the Sun as a function of radius and latitude. These data immediately reveal two striking boundary layers of shear in the solar convection zone (CZ): a tachocline at the base, where the differential rotation of the CZ transitions to solid-body rotation in the radiative zone, and a 35-Mm-thick near-surface shear layer (NSSL) at the top, where the rotation rate slows by about 5% with increasing radius. Though asteroseismology cannot probe the differential rotation of distant stars to the same level of detail that helioseismology can achieve for the Sun, it is possible that many cool stars with outer convective envelopes possess similar differential rotation characteristics, including both a tachocline and a NSSL. Here we present the results of 3D global hydrodynamic simulations of spherical-shell convection for a Sun-like star at different levels of density contrast across the shell. The simulations with high stratification possess characteristics of near-surface shear, especially at low latitudes. We discuss in detail the dynamical balance of torques giving rise to the NSSL in our models and interpret what these balances imply for the real Sun. We further discuss the dynamical causes that may serve to wipe out near-surface shear at high latitudes, and conclude by offering some theories as to how this problem might be tackled in future work. △ Less

Submitted 1 November, 2018; originally announced November 2018.

Comments: 5 pages, 6 figures; In the 20th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, edited by Scott J. Wolk

The Role of Downflows in Establishing Solar Near-Surface Shear

Authors: Loren I. Matilsky, Bradley W. Hindman, J. Toomre

Abstract: The dynamical origins of the Sun's tachocline and near-surface shear layer (NSSL) are still not well understood. We have attempted to self-consistently reproduce a NSSL in numerical simulations of a solar-like convection zone by increasing the density contrast across rotating, 3D spherical shells. We explore the hypothesis that high density contrast leads to near-surface shear by creating a rotati… ▽ More The dynamical origins of the Sun's tachocline and near-surface shear layer (NSSL) are still not well understood. We have attempted to self-consistently reproduce a NSSL in numerical simulations of a solar-like convection zone by increasing the density contrast across rotating, 3D spherical shells. We explore the hypothesis that high density contrast leads to near-surface shear by creating a rotationally unconstrained layer of fast flows near the outer surface. Although our high-contrast models do have near-surface shear, it is confined primarily to low latitudes (between $\pm15^\circ$). Two distinct types of flow structures maintain the shear dynamically: rotationally $\textit{constrained}$ Busse columns aligned with the rotation axis and fast, rotationally $\textit{unconstrained}$ downflow plumes that deplete angular momentum from the outer fluid layers. The plumes form at all latitudes, and in fact are more efficient at transporting angular momentum inward at high latitudes. The presence of Busse columns at low latitudes thus appears essential to creating near-surface shear in our models. We conclude that a solar-like NSSL is unobtainable from a rotationally unconstrained outer fluid layer alone. In numerical models, the shear is eliminated through the advection of angular momentum by the meridional circulation. Therefore, a detailed understanding how the solar meridional circulation is dynamically achieved will be necessary to elucidate the origin of the Sun's NSSL. △ Less

Submitted 28 September, 2018; originally announced October 2018.

Comments: 13 pages, 16 figures

A Novel Approach to Resonant Absorption of the Fast MHD Eigenmodes of a Coronal Arcade

Authors: Bradley W. Hindman, Rekha Jain

Abstract: The arched field lines forming coronal arcades are often observed to undulate as magnetohydrodynamic (MHD) waves propagate both across and along the magnetic field. These waves are most likely a combination of resonantly coupled fast magnetoacoustic waves and Alfvén waves. The coupling results in resonant absorption of the fast waves, converting fast wave energy into Alfvén waves. The fast eigenmo… ▽ More The arched field lines forming coronal arcades are often observed to undulate as magnetohydrodynamic (MHD) waves propagate both across and along the magnetic field. These waves are most likely a combination of resonantly coupled fast magnetoacoustic waves and Alfvén waves. The coupling results in resonant absorption of the fast waves, converting fast wave energy into Alfvén waves. The fast eigenmodes of the arcade have proven difficult to compute or derive analytically, largely because of the mathematical complexity that the coupling introduces. When a traditional spectral decomposition is employed, the discrete spectrum associated with the fast eigenmodes is often subsumed into the continuous Alfvén spectrum. Thus fast eigenmodes, become collective modes or quasi-modes. Here we present a spectral decomposition that treats the eigenmodes as having real frequencies but complex wavenumbers. Using this procedure we derive dispersion relations, spatial damping rates, and eigenfunctions for the resonant, fast eigenmodes of the arcade. We demonstrate that resonant absorption introduces a fast mode that would not exist otherwise. This new mode is heavily damped by resonant absorption, only travelling a few wavelengths before losing most of its energy. △ Less

Submitted 23 March, 2018; originally announced March 2018.

Comments: 43 pages, 12 figures, accepted by the Astrophysical Journal

Prandtl-Number Effects in High-Rayleigh-Number Spherical Convection

Authors: Ryan J. Orvedahl, Michael A. Calkins, Nicholas A. Featherstone, Bradley W. Hindman

Abstract: Convection is the predominant mechanism by which energy and angular momentum are transported in the outer portion of the Sun. The resulting overturning motions are also the primary energy source for the solar magnetic field. An accurate solar dynamo model therefore requires a complete description of the convective motions, but these motions remain poorly understood. Studying stellar convection num… ▽ More Convection is the predominant mechanism by which energy and angular momentum are transported in the outer portion of the Sun. The resulting overturning motions are also the primary energy source for the solar magnetic field. An accurate solar dynamo model therefore requires a complete description of the convective motions, but these motions remain poorly understood. Studying stellar convection numerically remains challenging; it occurs within a parameter regime that is extreme by computational standards. The fluid properties of the convection zone are characterized in part by the Prandtl number $\mathrm{Pr}=ν/κ$, where $ν$ is the kinematic viscosity and $κ$ is the thermal diffusion; in stars, $\mathrm{Pr}$ is extremely low, $\mathrm{Pr}\approx 10^{-7}$. The influence of $\mathrm{Pr}$ on the convective motions at the heart of the dynamo is not well understood since most numerical studies are limited to using $\mathrm{Pr}\approx1$. We systematically vary $\mathrm{Pr}$ and the degree of thermal forcing, characterized through a Rayleigh number, to explore its influence on the convective dynamics. For sufficiently large thermal driving, the simulations reach a so-called convective free-fall state where diffusion no longer plays an important role in the interior dynamics. Simulations with a lower $\mathrm{Pr}$ generate faster convective flows and broader ranges of scales for equivalent levels of thermal forcing. Characteristics of the spectral distribution of the velocity remain largely insensitive to changes in $\mathrm{Pr}$. Importantly, we find that $\mathrm{Pr}$ plays a key role in determining when the free-fall regime is reached by controlling the thickness of the thermal boundary layer. △ Less

Submitted 19 March, 2018; originally announced March 2018.

Journal ref: ApJ, Vol 856, Page 13, Number 1, 2018

The Emergence of Solar Supergranulation as a Natural Consequence of Rotationally-Constrained Interior Convection

Authors: Nicholas A. Featherstone, Bradley W. Hindman

Abstract: We investigate how rotationally-constrained, deep convection might give rise to supergranulation, the largest distinct spatial scale of convection observed in the solar photosphere. While supergranulation is only weakly influenced by rotation, larger spatial scales of convection sample the deep convection zone and are presumably rotationally influenced. We present numerical results from a series o… ▽ More We investigate how rotationally-constrained, deep convection might give rise to supergranulation, the largest distinct spatial scale of convection observed in the solar photosphere. While supergranulation is only weakly influenced by rotation, larger spatial scales of convection sample the deep convection zone and are presumably rotationally influenced. We present numerical results from a series of nonlinear, 3-D simulations of rotating convection and examine the velocity power distribution realized under a range of Rossby numbers. When rotation is present, the convective power distribution possesses a pronounced peak, at characteristic wavenumber $\ell_\mathrm{peak}$, whose value increases as the Rossby number is decreased. This distribution of power contrasts with that realized in non-rotating convection, where power increases monotonically from high to low wavenumbers. We find that spatial scales smaller than $\ell_\mathrm{peak}$ behave in analogy to non-rotating convection. Spatial scales larger than $\ell_\mathrm{peak}$ are rotationally-constrained and possess substantially reduced power relative to the non-rotating system. We argue that the supergranular scale emerges due to a suppression of power on spatial scales larger than $\ell\approx100$ owing to the presence of deep, rotationally-constrained convection. Supergranulation thus represents the largest non-rotationally-constrained mode of solar convection. We conclude that the characteristic spatial scale of supergranulation bounds that of the deep convective motions from above, making supergranulation an indirect measure of the deep-seated dynamics at work in the solar dynamo. Using the spatial scale of supergranulation in conjunction with our numerical results, we estimate an upper bound of 10 m s$^{-1}$ for the Sun's bulk rms convective velocity. △ Less

Submitted 16 September, 2016; originally announced September 2016.

Comments: Accepted for publication in ApJ Letters; Sep. 12, 2016

Eigenmodes of three-dimensional magnetic arcades in the Sun's corona

Authors: Bradley W. Hindman, Rekha Jain

Abstract: We develop a model of coronal-loop oscillations that treats the observed bright loops as an integral part of a larger 3-D magnetic structure comprised of the entire magnetic arcade. We demonstrate that magnetic arcades within the solar corona can trap MHD fast waves in a 3-D waveguide. This is accomplished through the construction of a cylindrically symmetric model of a magnetic arcade with a pote… ▽ More We develop a model of coronal-loop oscillations that treats the observed bright loops as an integral part of a larger 3-D magnetic structure comprised of the entire magnetic arcade. We demonstrate that magnetic arcades within the solar corona can trap MHD fast waves in a 3-D waveguide. This is accomplished through the construction of a cylindrically symmetric model of a magnetic arcade with a potential magnetic field. For a magnetically dominated plasma, we derive a governing equation for MHD fast waves and from this equation we show that the magnetic arcade forms a 3-D waveguide if the Alfvén speed increases monotonically beyond a fiducial radius. Both magnetic pressure and tension act as restoring forces, instead of just tension as is generally assumed in 1-D models. Since magnetic pressure plays an important role, the eigenmodes involve propagation both parallel and transverse to the magnetic field. Using an analytic solution, we derive the specific eigenfrequencies and eigenfunctions for an arcade possessing a discontinuous density profile. The discontinuity separates a diffuse cylindrical cavity and an overlying shell of denser plasma that corresponds to the bright loops. We emphasize that all of the eigenfunctions have a discontinuous axial velocity at the density interface; hence, the interface can give rise to the Kelvin-Helmholtz instability. Further, we find that all modes have elliptical polarization with the degree of polarization changing with height. However, depending on the line of sight, only one polarization may be clearly visible. △ Less

Submitted 30 November, 2015; originally announced November 2015.

Comments: 24 pages, 6 figures, 1 ancillary animation file

Journal ref: ApJ, 814, 105 (2015)

The Spectral Amplitude of Stellar Convection and its Scaling in the High-Rayleigh-Number Regime

Authors: Nicholas A. Featherstone, Bradley W. Hindman

Abstract: Convection plays a central role in the dynamics of any stellar interior, and yet its operation remains largely-hidden from direct observation. As a result, much of our understanding concerning stellar convection necessarily derives from theoretical and computational models. The Sun is, however, exceptional in that regard. The wealth of observational data afforded by its proximity provides a unique… ▽ More Convection plays a central role in the dynamics of any stellar interior, and yet its operation remains largely-hidden from direct observation. As a result, much of our understanding concerning stellar convection necessarily derives from theoretical and computational models. The Sun is, however, exceptional in that regard. The wealth of observational data afforded by its proximity provides a unique testbed for comparing convection models against observations. When such comparisons are carried out, surprising inconsistencies between those models and observations become apparent. Both photospheric and helioseismic measurements suggest that convection simulations may overestimate convective flow speeds on large spatial scales. Moreover, many solar convection simulations have difficulty reproducing the observed solar differential rotation due to this apparent overestimation. We present a series of 3-dimensional (3-D) stellar convection simulations designed to examine how the amplitude and spectral distribution of convective flows are established within a star's interior. While these simulations are non-magnetic and non-rotating in nature, they demonstrate two robust phenomena. When run with sufficiently high Rayleigh number, the integrated kinetic energy of the convection becomes effectively independent of thermal diffusion, but the spectral distribution of that kinetic energy remains sensitive to both of these quantities. A simulation that has converged to a diffusion-independent value of kinetic energy will divide that energy between spatial scales such that low-wavenumber power is overestimated, and high-wavenumber power is underestimated relative to a comparable system possessing higher Rayleigh number. We discuss the implications of these results in light of the current inconsistencies between models and observations. △ Less

Submitted 7 November, 2015; originally announced November 2015.

Comments: Submitted to The Astrophysical Journal

Fundamental-Mode Oscillations of Two Coronal Loops within a Solar Magnetic Arcade

Authors: Rekha Jain, Ram A. Maurya, B. W. Hindman

Abstract: We analyse intensity variations, as measured by the Atmospheric Imaging Assembly (AIA) in the 171 Å passband, in two coronal loops embedded within a single coronal magnetic arcade. We detect oscillations in the fundamental mode with periods of roughly 2 minutes and decay times of 5 minutes. The oscillations were initiated by interaction of the arcade with a large wavefront issuing from a flare sit… ▽ More We analyse intensity variations, as measured by the Atmospheric Imaging Assembly (AIA) in the 171 Å passband, in two coronal loops embedded within a single coronal magnetic arcade. We detect oscillations in the fundamental mode with periods of roughly 2 minutes and decay times of 5 minutes. The oscillations were initiated by interaction of the arcade with a large wavefront issuing from a flare site. Further, the power spectra of the oscillations evince signatures consistent with oblique propagation to the field lines and for the existence of a 2-D waveguide instead of a 1-D one. △ Less

Submitted 29 April, 2015; originally announced April 2015.

Helioseismic Imaging of Fast Convective Flows Throughout the Near-Surface Shear Layer

Authors: Benjamin J. Greer, Bradley W. Hindman, Nicholas A. Featherstone, Juri Toomre

Abstract: Using a new implementation of ring-diagram helioseismology, we ascertain the strength and spatial scale of convective flows throughout the near-surface shear layer. Our ring-diagram technique employs highly overlapped analysis regions and an efficient method of 3D inversion to measure convective motions with a resolution that ranges from $3 \ \mathrm{Mm}$ at the surface to $80 \ \mathrm{Mm}$ at th… ▽ More Using a new implementation of ring-diagram helioseismology, we ascertain the strength and spatial scale of convective flows throughout the near-surface shear layer. Our ring-diagram technique employs highly overlapped analysis regions and an efficient method of 3D inversion to measure convective motions with a resolution that ranges from $3 \ \mathrm{Mm}$ at the surface to $80 \ \mathrm{Mm}$ at the base of the layer. We find the rms horizontal flow speed to peak at $427 \ \mathrm{m \ s^{-1}}$ at the photosphere and fall to a minimum of $124 \ \mathrm{m \ s^{-1}}$ between $20 \ \mathrm{Mm}$ and $30 \ \mathrm{Mm}$. From the velocity amplitude and the dominant horizontal scales seen at each depth, we infer the level of rotational influence on convection to be low near the surface, but transition to a significant level at the base of the near-surface shear layer with a Rossby number varying between 2.2 to as low as 0.1. △ Less

Submitted 2 April, 2015; originally announced April 2015.

Comments: 12 pages, 5 figures

Five-minute oscillation power within magnetic elements in the solar atmosphere

Authors: Rekha Jain, Andrew Gascoyne, Bradley W. Hindman, Benjamin Greer

Abstract: It has long been known that magnetic plage and sunspots are regions in which the power of acoustic waves is reduced within the photospheric layers. Recent observations now suggest that this suppression of power extends into the low chromosphere and is also present in small magnetic elements far from active regions. In this paper we investigate the observed power supression in plage and magnetic el… ▽ More It has long been known that magnetic plage and sunspots are regions in which the power of acoustic waves is reduced within the photospheric layers. Recent observations now suggest that this suppression of power extends into the low chromosphere and is also present in small magnetic elements far from active regions. In this paper we investigate the observed power supression in plage and magnetic elements, by modelling each as a collection of vertically aligned magnetic fibrils and presuming that the velocity within each fibril is the response to buffeting by incident $p$ modes in the surrounding field-free atmosphere. We restrict our attention to modeling observations made near solar disk center, where the line-of-sight velocity is nearly vertical and hence, only the longitudinal component of the motion within the fibril contributes. Therefore, we only consider the excitation of axisymmetric sausage waves and ignore kink oscillations as their motions are primarily horizontal. We compare the vertical motion within the fibril with the vertical motion of the incident $p$ mode by constructing the ratio of their powers. In agreement with observational measurements we find that the total power is suppressed within strong magnetic elements for frequencies below the acoustic cut-off frequency. We also find that the magnitude of the power deficit increases with the height above the photosphere at which the measurement is made. Further, we argue that the area of the solar disk over which the power suppression extends increases as a function of height. △ Less

Submitted 6 May, 2014; v1 submitted 4 May, 2014; originally announced May 2014.

Multi-Ridge Fitting for Ring-Diagram Helioseismology

Authors: Benjamin J. Greer, Bradley W. Hindman, Juri Toomre

Abstract: Inferences of sub-surface flow velocities using local domain ring-diagram helioseismology depend on measuring the frequency splittings of oscillation modes seen in acoustic power spectra. Current methods for making these measurements utilize maximum-likelihood fitting techniques to match a model of modal power to the spectra. The model typically describes a single oscillation mode, and each mode i… ▽ More Inferences of sub-surface flow velocities using local domain ring-diagram helioseismology depend on measuring the frequency splittings of oscillation modes seen in acoustic power spectra. Current methods for making these measurements utilize maximum-likelihood fitting techniques to match a model of modal power to the spectra. The model typically describes a single oscillation mode, and each mode in a given power spectrum is fit independently. We present a new method that produces measurements with greater reliability and accuracy by fitting multiple modes simultaneously. We demonstrate how this method permits measurements of sub-surface flows deeper into the Sun while providing higher uniformity in data coverage and velocity response closer to the limb of the solar disk. While the previous fitting method performs better for some measurements of low-phase-speed modes, we find this new method to be particularly useful for high phase-speed modes and small spatial areas. △ Less

Submitted 20 February, 2014; originally announced February 2014.

Comments: 26 pages, 11 figures

An Interpretation of Flare-Induced and Decayless Coronal-Loop Oscillations as Interference Patterns

Authors: Bradley W. Hindman, Rekha Jain

Abstract: We present an alternative model of coronal-loop oscillations which considers that the waves are trapped in a 2D waveguide formed by the entire arcade of field lines. This differs from the standard 1D model which treats the waves as the resonant oscillations of just the visible bundle of field lines. Within the framework of our 2D model, the two types of oscillations that have been observationally… ▽ More We present an alternative model of coronal-loop oscillations which considers that the waves are trapped in a 2D waveguide formed by the entire arcade of field lines. This differs from the standard 1D model which treats the waves as the resonant oscillations of just the visible bundle of field lines. Within the framework of our 2D model, the two types of oscillations that have been observationally identified, flare-induced waves and "decayless" oscillations, can both be attributed to MHD fast waves. The two components of the signal differ only because of the duration and spatial extent of the source that creates them. The flare-induced waves are generated by strong localized sources of short duration, while the decayless background can be excited by a continuous, stochastic source. Further, the oscillatory signal arising from a localized, short-duration source can be interpreted as a pattern of interference fringes produced by waves that have traveled diverse routes of various pathlengths through the waveguide. The resulting amplitude of the fringes slowly decays in time with an inverse square root dependence. The details of the interference pattern depend on the shape of the arcade and the spatial variation of the Alfvén speed. The rapid decay of this wave component, which has previously been attributed to physical damping mechanisms that remove energy from resonant oscillations, occurs as a natural consequence of the interference process without the need for local dissipation. △ Less

Submitted 6 December, 2013; originally announced December 2013.

Comments: 26 pages, 8 figures, submitted to the Astrophysical Journal

Equilibrium models of coronal loops that involve curvature and buoyancy

Authors: Bradley W. Hindman, Rekha Jain

Abstract: We construct magnetostatic models of coronal loops in which the thermodynamics of the loop is fully consistent with the shape and geometry of the loop. This is achieved by treating the loop as a thin, compact, magnetic fibril that is a small departure from a force-free state. The density along the loop is related to the loop's curvature by requiring that the Lorentz force arising from this deviati… ▽ More We construct magnetostatic models of coronal loops in which the thermodynamics of the loop is fully consistent with the shape and geometry of the loop. This is achieved by treating the loop as a thin, compact, magnetic fibril that is a small departure from a force-free state. The density along the loop is related to the loop's curvature by requiring that the Lorentz force arising from this deviation is balanced by buoyancy. This equilibrium, coupled with hydrostatic balance and the ideal gas law, then connects the temperature of the loop with the curvature of the loop without resorting to a detailed treatment of heating and cooling. We present two example solutions: one with a spatially invariant magnetic Bond number (the dimensionless ratio of buoyancy to Lorentz forces) and the other with a constant radius of curvature of the loop's axis. We find that the density and temperature profiles are quite sensitive to curvature variations along the loop, even for loops with similar aspect ratios. △ Less

Submitted 2 August, 2013; originally announced August 2013.

Comments: 31 pages, 6 figures

Kink Oscillations of a Curved, Gravitationally Stratified Coronal Loop

Authors: Bradley W. Hindman, Rekha Jain

Abstract: Loops of magnetic field in the corona are observed to oscillate and these oscillations have been posited to be the superposition of resonant kink waves. To date, most analyses of these oscillations have concentrated on calculating the frequency shifts that result from spatial variation in the kink wave speed. Further, most have ignored gravity and treated the loop as a straight tube. Here we ignor… ▽ More Loops of magnetic field in the corona are observed to oscillate and these oscillations have been posited to be the superposition of resonant kink waves. To date, most analyses of these oscillations have concentrated on calculating the frequency shifts that result from spatial variation in the kink wave speed. Further, most have ignored gravity and treated the loop as a straight tube. Here we ignore spatial variation in the wave speed, but self-consistently include the effects of gravity and loop curvature in both the equilibrium loop model and in the wave equation. We model a coronal loop as an isolated, thin, magnetic fibril that is anchored at two points in the photosphere. The equilibrium shape of the loop is determined by a balance between magnetic buoyancy and magnetic tension, which is characterized by a Magnetic Bond Number ε, that is typically small |ε| << 1. This balance produces a loop that has a variable radius of curvature. The resonant kink waves of such a loop come in two polarizations that are decoupled from each other: waves with motion completely within the plane of the loop (normal oscillations) and waves with motions that are completely horizontal, perpendicular to the plane of the loop (binormal oscillations). We solve for the eigensolutions of both polarizations using perturbation theory for small Magnetic Bond Number. For modes of the same order, normal oscillations have smaller eigenfrequencies than binormal oscillations. The additional forces of buoyancy and magnetic tension from the curvature of the loop increase and decrease the mode frequencies, respectively. The ratio of the frequencies of the first overtone to the fundamental mode is modified by the inclusion of buoyancy and curvature. We find that the normal polarization possesses a frequency ratio that exceeds the canonical value of 2, whereas the binormal polarization has a ratio less than 2. △ Less

Submitted 25 September, 2012; originally announced September 2012.

Comments: Submitted to ApJ 8/21/2012, 28 pages, 4 figures

An Improved Method for Fitting p-Mode Profile Asymmetries

Authors: Bradley W. Hindman

Abstract: In a power spectrum of the Sun's acoustic waves, the p modes have distinctly skewed frequency profiles. Furthermore, the asymmetry is observed to have the opposite sign in power spectra made from line-of-sight velocity and continuum intensity. The asymmetry and its reversal in sign has previously been explained using a combination of mechanisms that involve the acoustic source. A localized acousti… ▽ More In a power spectrum of the Sun's acoustic waves, the p modes have distinctly skewed frequency profiles. Furthermore, the asymmetry is observed to have the opposite sign in power spectra made from line-of-sight velocity and continuum intensity. The asymmetry and its reversal in sign has previously been explained using a combination of mechanisms that involve the acoustic source. A localized acoustic source within an acoustic cavity naturally generates asymmetric profiles through wave interference; however, the sign of the asymmetry due to this mechanism is identical for all observables. The reversal of the asymmetry between velocity and intensity observations has been attributed to the visibility of the source itself (i.e., "correlated noise"). In this paper, I will show that asymmetry generated by a localized source can be interpreted as either a wave interference effect in physical space, or a mode interference effect in spectral space. I advocate a new mode-fitting procedure based on this new interpretation, whereby the complex phases of all the modes determine the mode asymmetries. Further, I suggest that information about the acoustic source function is encapsulated in the amplitude of each mode, and present a scheme by which the source function can be obtained from measured mode amplitudes by standard helioseismic inversion techniques. △ Less

Submitted 20 December, 2011; originally announced December 2011.

Comments: 27 pages, 5 figures, to be submitted to ApJ

Axisymmetric Scattering of p Modes by Thin Magnetic Tubes

Authors: Bradley W. Hindman, Rekha Jain

Abstract: We examine the scattering of acoustic p-mode waves from a thin magnetic fibril embedded in a gravitationally stratified atmosphere. The scattering is mediated through the excitation of slow sausage waves on the magnetic tube, and only the scattering of the monopole component of the wavefield is considered. Since such tube waves are not confined by the acoustic cavity and may freely propagate along… ▽ More We examine the scattering of acoustic p-mode waves from a thin magnetic fibril embedded in a gravitationally stratified atmosphere. The scattering is mediated through the excitation of slow sausage waves on the magnetic tube, and only the scattering of the monopole component of the wavefield is considered. Since such tube waves are not confined by the acoustic cavity and may freely propagate along the field lines removing energy from the acoustic wavefield, the excitation of fibril oscillations is a source of acoustic wave absorption as well as scattering. We compute the mode mixing that is achieved and the absorption coefficients and phase shifts. We find that for thin tubes the mode mixing is weak and the absorption coefficient is small and is a smooth function of frequency over the physically relevant band of observed frequencies. The prominent absorption resonances seen in previous studies of unstratified tubes are absent. Despite the relatively small absorption, the phase shift induced can be surprisingly large, reaching values as high as 15 degree for f modes. Further, the phase shift can be positive or negative depending on the incident mode order and the frequency. △ Less

Submitted 18 November, 2011; v1 submitted 24 June, 2011; originally announced June 2011.

Comments: 39 pages, 6 figures

Gyroscopic Pumping in the Solar Near-Surface Shear Layer

Authors: Mark S. Miesch, Bradley W. Hindman

Abstract: We use global and local helioseismic inversions to explore the prevailing dynamical balances in the solar Near-Surface Shear Layer (NSSL). The differential rotation and meridional circulation are intimately linked, with a common origin in the turbulent stresses of the upper solar convection zone. The existence and structure of the NSSL cannot be attributed to the conservation of angular momentum b… ▽ More We use global and local helioseismic inversions to explore the prevailing dynamical balances in the solar Near-Surface Shear Layer (NSSL). The differential rotation and meridional circulation are intimately linked, with a common origin in the turbulent stresses of the upper solar convection zone. The existence and structure of the NSSL cannot be attributed to the conservation of angular momentum by solar surface convection, as is often supposed. Rather, the turbulent angular momentum transport accounts for the poleward meridional flow while the often overlooked meridional force balance maintains the mid-latitude rotational shear. We suggest that the base of the NSSL is marked by a transition from baroclinic to turbulent stresses in the meridional plane which suppress Coriolis-induced circulations that would otherwise establish a cylindrical rotation profile. The turbulent angular momentum transport must be non-diffusive and directed radially inward. Inferred mean flows are consistent with the idea that turbulent convection tends to mix angular momentum but only if the mixing efficiency is inhomogeneous and/or anisotropic. The latitudinal and longitudinal components of the estimated turbulent transport are comparable in amplitude and about an order of magnitude larger than the vertical component. We estimate that it requires 2--4% of the solar luminosity to maintain the solar NSSL against the inertia of the mean flow. Most of this energy is associated with the turbulent transport of angular momentum out of the layer, with a spin-down time scale of $\sim$ 600 days. We also address implications of these results for numerical modeling of the NSSL. △ Less

Submitted 26 August, 2011; v1 submitted 21 June, 2011; originally announced June 2011.

Comments: 26 pages, 11 Figures, Accepted on August 25, 2011 for publication in the Astrophysical Journal

Journal ref: The Astrophysical Journal, 743: 79 (25pp), 2011 Dec 10

Modeling the Subsurface Structure of Sunspots

Authors: H. Moradi, C. Baldner, A. C. Birch, D. Braun, R. Cameron, T. L. Duvall Jr., L. Gizon, D. Haber, S. Hanasoge, B. W. Hindman, J. Jackiewicz, E. Khomenko, R. Komm, P. Rajaguru, M. Rempel, M. Roth, R. Schlichenmaier, H. Schunker, H. Spruit, K. Strassmeier, M. J. Thompson, S. Zharkov

Abstract: While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismol… ▽ More While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat. △ Less

Submitted 23 August, 2010; v1 submitted 25 December, 2009; originally announced December 2009.

Comments: 73 pages, 19 figures, accepted by Solar Physics

Helioseismic detection of deep meridional flow

Authors: Douglas Gough, Bradley W. Hindman

Abstract: Steady meridional flow makes no first-order perturbation to the frequencies of helioseismic normal modes. It does, however, Doppler shift the local wavenumber, thereby distorting the eigenfunctions. For high-degree modes, whose peaks in a power spectrum are blended into continuous ridges, the effect of the distortion is to shift the locations of those ridges. From this blended superposition of m… ▽ More Steady meridional flow makes no first-order perturbation to the frequencies of helioseismic normal modes. It does, however, Doppler shift the local wavenumber, thereby distorting the eigenfunctions. For high-degree modes, whose peaks in a power spectrum are blended into continuous ridges, the effect of the distortion is to shift the locations of those ridges. From this blended superposition of modes, one can isolate oppositely directed wave components with the same local horizontal wavenumber and measure a frequency difference which can be safely used to infer the subsurface background flow. But such a procedure fails for the components of the more-deeply-penetrating low-degree modes that are not blended into ridges. Instead, one must analyze the spatial distortions explicitly. With a simple toy model, we illustrate one method by which that might be accomplished by measuring the spatial variation of the oscillation phase. We estimate that by this procedure it might be possible to infer meridional flow deep in the solar convection zone. △ Less

Submitted 10 November, 2009; originally announced November 2009.

Comments: 23 pages, 9 color figures, submitted to the Astrophysical Journal

Journal ref: Astrophys.J.714:960-970,2010

Subsurface Circulations within Active Regions

Authors: Bradley W. Hindman, Deborah A. Haber, Juri Toomre

Abstract: Using high-resolution ring analysis (HRRA) we deduce subsurface flows within magnetic active regions and within quiet sun. With this procedure we are capable of measuring flows with a horizontal spatial resolution of 2 degrees in heliographic angle (or roughly 20 Mm). From the resulting flow fields we deduce mean inflow rates into active regions, mean circulation speeds around active regions, an… ▽ More Using high-resolution ring analysis (HRRA) we deduce subsurface flows within magnetic active regions and within quiet sun. With this procedure we are capable of measuring flows with a horizontal spatial resolution of 2 degrees in heliographic angle (or roughly 20 Mm). From the resulting flow fields we deduce mean inflow rates into active regions, mean circulation speeds around active regions, and probability density functions (PDFs) of properties of the flow field. These analyses indicate that active regions have a zonal velocity that exceeds that of quiet sun at the same latitude by 20 m/s, yet active regions advect poleward at the same rate as quiet sun. We also find that almost all active regions possess a mean inflow (20-30 m/s) and a cyclonic circulation (~5 m/s) at their peripheries, whereas their cores, where the sunspots are located, are zones of strong anticyclonic outflow (~50 m/s). From the PDFs, we find that active regions modify the structure of convection with a scale greater than that of supergranulation. Instead of possessing an asymmetry between inflows and outflows (with a larger percentage of the surface occupied by outflows), as is seen in quiet sun, active regions possess symmetric distributions. △ Less

Submitted 9 April, 2009; originally announced April 2009.

Comments: 34 pages, 12 figures, submitted to the Astrophysical Journal

Journal ref: Astrophys.J.698:1749-1760,2009

The Generation of Coronal Loop Waves below the Photosphere by p-Mode Forcing

Authors: Bradley W. Hindman, Rekha Jain

Abstract: Recent observations of coronal-loop waves by TRACE and within the corona as a whole by CoMP clearly indicate that the dominant oscillation period is 5 minutes, thus implicating the solar p modes as a possible source. We investigate the generation of tube waves within the solar convection zone by the buffeting of p modes. The tube waves--in the form of longitudinal sausage waves and transverse ki… ▽ More Recent observations of coronal-loop waves by TRACE and within the corona as a whole by CoMP clearly indicate that the dominant oscillation period is 5 minutes, thus implicating the solar p modes as a possible source. We investigate the generation of tube waves within the solar convection zone by the buffeting of p modes. The tube waves--in the form of longitudinal sausage waves and transverse kink waves--are generated on the many magnetic fibrils that lace the convection zone and pierce the solar photosphere. Once generated by p-mode forcing, the tube waves freely propagate up and down the tubes, since the tubes act like light fibers and form a waveguide for these magnetosonic waves. Those waves that propagate upward pass through the photosphere and enter the upper atmosphere where they can be measured as loop oscillations and other forms of propagating coronal waves. We treat the magnetic fibrils as vertically aligned, thin flux tubes and compute the energy flux of tube waves that can generated and driven into the upper atmosphere. We find that a flux in excess of 10^5 ergs/cm^2/s can be produced, easily supplying enough wave energy to explain the observations. Furthermore, we compute the associated damping rate of the driving p modes and find that the damping is significant compared to observed line widths only for the lowest order p modes. △ Less

Submitted 13 May, 2008; originally announced May 2008.

Comments: AAS LaTex v5.0, 23 pages, 7 figures, previously published in ApJ

Journal ref: Astrophys.J.677:769-780,2008