Showing 1–37 of 37 results for author: Schunker, H

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Evolution of the magnetic field and flows of solar active regions with persistent magnetic bipoles before emergence

Authors: Camron S. Alley, Hannah Schunker

Abstract: Magnetic active regions on the Sun are harbingers of space weather events. Understanding the physics of how they form and evolve will improve space weather forecasting. Our aim is to characterise the surface magnetic field and flows for a sample of active regions with persistent magnetic bipoles prior to emergence. We identified 42 emerging active regions (EARs), in the Solar Dynamics Observatory… ▽ More Magnetic active regions on the Sun are harbingers of space weather events. Understanding the physics of how they form and evolve will improve space weather forecasting. Our aim is to characterise the surface magnetic field and flows for a sample of active regions with persistent magnetic bipoles prior to emergence. We identified 42 emerging active regions (EARs), in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey (Schunker et al. 2016), associated with small magnetic bipoles at least one day before the time of emergence. We then identified a contrasting sample of 42 EARs that emerge more abruptly without bipoles before emergence. We computed the supergranulation scale surface flows using helioseismic holography. We averaged the flow maps and magnetic field maps over all active regions in each sample at each time interval from 2 days before emergence to 1 day after. We found that EARs associated with a persistent pre-emergence bipole evolve to be, on average, lower flux active regions than EARs that emerge more abruptly. Further, we found that the EARs that emerge more abruptly do so with a diverging flow of $(3\pm 0.6) \times 10^{-6}$ s$^{-1}$ on the order of 50-100 ms$^{-1}$. Our results suggest that there is a statistical dependence of the surface flow signature throughout the emergence process on the maximum magnetic flux of the active region. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 15 pages, 6 figures, 2 tables

Journal ref: Publications of the Astronomical Society of Australia, 2023, 40

Understanding Active Region Emergence and Origins on the Sun and Other Cool Stars

Authors: Maria A. Weber, Hannah Schunker, Laurène Jouve, Emre Işık

Abstract: The emergence of active regions on the Sun is an integral feature of the solar dynamo mechanism. However, details about the generation of active-region-scale magnetism and the journey of this magnetic flux to the photosphere are still in question. Shifting paradigms are now developing for the source depth of the Sun's large-scale magnetism, the organization of this magnetism into fibril flux tubes… ▽ More The emergence of active regions on the Sun is an integral feature of the solar dynamo mechanism. However, details about the generation of active-region-scale magnetism and the journey of this magnetic flux to the photosphere are still in question. Shifting paradigms are now developing for the source depth of the Sun's large-scale magnetism, the organization of this magnetism into fibril flux tubes, and the role of convection in shaping active-region observables. Here we review the landscape of flux emergence theories and simulations, highlight the role flux emergence plays in the global dynamo process, and make connections between flux emergence on the Sun and other cool stars. As longer-term and higher fidelity observations of both solar active regions and their associated flows are amassed, it is now possible to place new constraints on models of emerging flux. We discuss the outcomes of statistical studies which provide observational evidence that flux emergence may be a more passive process (at least in the upper convection zone); dominated to a greater extent by the influence of convection and to a lesser extent by buoyancy and the Coriolis force acting on rising magnetic flux tubes than previously thought. We also discuss how the relationship between stellar rotation, fractional convection zone depth, and magnetic activity on other stars can help us better understand flux emergence processes. Looking forward, we identify open questions regarding magnetic flux emergence that we anticipate can be addressed in the next decade with further observations and simulations. △ Less

Submitted 10 June, 2023; originally announced June 2023.

Comments: 44 pages, 6 figures, ISSI Workshop "Solar and stellar dynamos"

Impact of spatially correlated fluctuations in sunspots on metrics related to magnetic twist

Authors: C. Baumgartner, A. C. Birch, H. Schunker, R. H. Cameron, L. Gizon

Abstract: The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxi… ▽ More The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxies is based on $α_z$, the ratio of the vertical current to the vertical magnetic field. Another class of proxies is based on the so-called twist density, $q$, which depends on the ratio of the azimuthal field to the vertical field. However, the sensitivity of these proxies to temporal fluctuations of the magnetic field has not yet been well characterized. We aim to determine the sensitivity of twist proxies to temporal fluctuations in the magnetic field as estimated from time-series of SDO/HMI vector magnetic field maps. To this end, we introduce a model of a sunspot with a peak vertical field of 2370 Gauss at the photosphere and a uniform twist density $q= -0.024$ Mm$^{-1}$. We add realizations of the temporal fluctuations of the magnetic field that are consistent with SDO/HMI observations, including the spatial correlations. Using a Monte-Carlo approach, we determine the robustness of the different proxies to the temporal fluctuations. The temporal fluctuations of the three components of the magnetic field are correlated for spatial separations up to 1.4 Mm (more than expected from the point spread function alone). The Monte-Carlo approach enables us to demonstrate that several proxies for the twist of the magnetic field are not biased in each of the individual magnetograms. The associated random errors on the proxies have standard deviations in the range between $0.002$ and $0.006$ Mm$^{-1}$, which is smaller by approximately one order of magnitude than the mean value of $q$. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: accepted in A&A

Journal ref: A&A 664, A183 (2022)

Testing solar surface flux transport models in the first days after active region emergence

Authors: Nils Gottschling, Hannah Schunker, Aaron C. Birch, Robert Cameron, Laurent Gizon

Abstract: Active regions (ARs) play an important role in the magnetic dynamics of the Sun. Solar surface flux transport models (SFTMs) are used to describe the evolution of the radial magnetic field at the solar surface. There is however uncertainty about using these models in the early stage of AR evolution. We aim to test the applicability of SFTMs in the first days after the emergence of ARs by comparing… ▽ More Active regions (ARs) play an important role in the magnetic dynamics of the Sun. Solar surface flux transport models (SFTMs) are used to describe the evolution of the radial magnetic field at the solar surface. There is however uncertainty about using these models in the early stage of AR evolution. We aim to test the applicability of SFTMs in the first days after the emergence of ARs by comparing them with observations. The models we employ range from passive evolution to models where the inflows around ARs are included. We simulate the evolution of the surface magnetic field of 17 emerging active regions using a local surface flux transport simulation. We selected regions that do not form fully-fledged sunspots that exhibit moat flows. The simulation includes diffusion and advection. We use observed flows from local correlation tracking of solar granulation, as well as parametrizations of the inflows around ARs. To evaluate our simulations, we measure the cross correlation between the observed and the simulated magnetic field, as well as the total unsigned flux of the ARs, over time. We also test the validity of our simulations by varying the starting time relative to the emergence of flux. We find that the simulations using observed surface flows can reproduce the evolution of the observed magnetic flux. The effect of buffeting of the field by supergranulation can be described as a diffusion process. The SFTM is applicable after 90% of the peak total unsigned flux of the AR has emerged. Diffusivities in the range between $D=250$ to $720$ km$^2$/s are consistent with the evolution of the AR flux in the first five days after this time. We find that the converging flows around emerging ARs are not important for the evolution of the total flux of the AR in these first five days; their effect of increasing flux cancellation is balanced by the decrease of flux transport away from the AR. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Astronomy & Astrophysics. 13 pages, 11 figures

Journal ref: A&A 660, A6 (2022)

Evolution of solar surface inflows around emerging active regions

Authors: N. Gottschling, H. Schunker, A. C. Birch, B. Löptien, L. Gizon

Abstract: Solar active regions are associated with Evershed outflows in sunspot penumbrae, moat outflows surrounding sunspots, and extended inflows surrounding active regions. The latter have been identified on established active regions by various methods. The evolution of these inflows and their dependence on active region properties as well as their impact on the global magnetic field are not yet underst… ▽ More Solar active regions are associated with Evershed outflows in sunspot penumbrae, moat outflows surrounding sunspots, and extended inflows surrounding active regions. The latter have been identified on established active regions by various methods. The evolution of these inflows and their dependence on active region properties as well as their impact on the global magnetic field are not yet understood. We aim to understand the evolution of the average inflows around emerging active regions and to derive an empirical model for these inflows. We analyze horizontal flows at the surface of the Sun using local correlation tracking of solar granules observed in continuum images of SDO/HMI. We measure average flows of a sample of 182 isolated active regions up to seven days before and after their emergence onto the solar surface with a cadence of 12 hours. We investigate the average inflow properties with respect to active region characteristics of total flux and latitude. We fit a model to these observed inflows for a quantitative analysis. We find that converging flows of around $20$ to $30$ m/s are first visible one day prior to emergence, in agreement with recent results. These converging flows are present independently of active region properties of latitude or flux. We confirm a recently found prograde flow of about $40$ m/s at the leading polarity during emergence. We find that the time after emergence when the latitudinal inflows increase in amplitude depends on the flux of the active region, ranging from one to four days after emergence and increasing with flux. The largest extent of the inflows is up to about $7 \pm 1^\circ$ away from the center of the active region within the first six days after emergence. The inflow velocities have amplitudes of about $50$ m/s. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: Accepted for publication in Astronomy & Astrophysics. 16 pages, 17 figures

Journal ref: A&A 652, A148 (2021)

Application and interpretation of deep learning for identifying pre-emergence magnetic-field patterns

Authors: Dattaraj B. Dhuri, Shravan M. Hanasoge, Aaron C. Birch, Hannah Schunker

Abstract: Magnetic flux generated within the solar interior emerges to the surface, forming active regions (ARs) and sunspots. Flux emergence may trigger explosive events, such as flares and coronal mass ejections and therefore understanding emergence is useful for space-weather forecasting. Evidence of any pre-emergence signatures will also shed light on sub-surface processes responsible for emergence. In… ▽ More Magnetic flux generated within the solar interior emerges to the surface, forming active regions (ARs) and sunspots. Flux emergence may trigger explosive events, such as flares and coronal mass ejections and therefore understanding emergence is useful for space-weather forecasting. Evidence of any pre-emergence signatures will also shed light on sub-surface processes responsible for emergence. In this paper, we present a first analysis of emerging ARs from the Solar Dynamics Observatory/Helioseismic Emerging Active Regions (SDO/HEAR) dataset (Schunker et al. 2016) using deep convolutional neural networks (CNN) to characterize pre-emergence surface magnetic-field properties. The trained CNN classifies between pre-emergence (PE) line-of-sight magnetograms and a control set of non-emergence (NE) magnetograms with a True Skill Statistic (TSS) score of ~85%, 3h prior to emergence and ~40\%, 24h prior to emergence. Our results are better than a baseline classification TSS obtained using discriminant analysis of only the unsigned magnetic flux. We develop a network pruning algorithm to interpret the trained CNN and show that the CNN incorporates filters that respond positively as well as negatively to the unsigned magnetic flux of the magnetograms. Using synthetic magnetograms, we demonstrate that the CNN output is sensitive to the length-scale of the magnetic regions with small-scale and intense fields producing maximum CNN output and possibly a characteristic pre-emergence pattern. Given increasing popularity of deep learning, techniques developed here for interpretation of the trained CNN -- using network pruning and synthetic data -- are relevant for future applications in solar and astrophysical data analysis. △ Less

Submitted 14 September, 2020; originally announced September 2020.

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

The Evolution of Rotation and Magnetic Activity in 94 Aqr Aa from Asteroseismology with TESS

Authors: Travis S. Metcalfe, Jennifer L. van Saders, Sarbani Basu, Derek Buzasi, William J. Chaplin, Ricky Egeland, Rafael A. Garcia, Patrick Gaulme, Daniel Huber, Timo Reinhold, Hannah Schunker, Keivan G. Stassun, Thierry Appourchaux, Warrick H. Ball, Timothy R. Bedding, Sebastien Deheuvels, Lucia Gonzalez-Cuesta, Rasmus Handberg, Antonio Jimenez, Hans Kjeldsen, Tanda Li, Mikkel N. Lund, Savita Mathur, Benoit Mosser, Martin B. Nielsen , et al. (7 additional authors not shown)

Abstract: Most previous efforts to calibrate how rotation and magnetic activity depend on stellar age and mass have relied on observations of clusters, where isochrones from stellar evolution models are used to determine the properties of the ensemble. Asteroseismology employs similar models to measure the properties of an individual star by matching its normal modes of oscillation, yielding the stellar age… ▽ More Most previous efforts to calibrate how rotation and magnetic activity depend on stellar age and mass have relied on observations of clusters, where isochrones from stellar evolution models are used to determine the properties of the ensemble. Asteroseismology employs similar models to measure the properties of an individual star by matching its normal modes of oscillation, yielding the stellar age and mass with high precision. We use 27 days of photometry from the Transiting Exoplanet Survey Satellite to characterize solar-like oscillations in the G8 subgiant of the 94 Aqr triple system. The resulting stellar properties, when combined with a reanalysis of 35 yr of activity measurements from the Mount Wilson HK project, allow us to probe the evolution of rotation and magnetic activity in the system. The asteroseismic age of the subgiant agrees with a stellar isochrone fit, but the rotation period is much shorter than expected from standard models of angular momentum evolution. We conclude that weakened magnetic braking may be needed to reproduce the stellar properties, and that evolved subgiants in the hydrogen shell-burning phase can reinvigorate large-scale dynamo action and briefly sustain magnetic activity cycles before ascending the red giant branch. △ Less

Submitted 25 August, 2020; v1 submitted 24 July, 2020; originally announced July 2020.

Comments: 14 pages including 8 figures and 3 tables (updated Table 3 & Figure 6). ApJ in press

Journal ref: Astrophys. J. 900, 154 (2020)

Average motion of emerging solar active region polarities II: Joy's law

Authors: Hannah Schunker, Christian Baumgartner, Aaron C. Birch, Robert H. Cameron, Douglas C. Braun, Laurent Gizon

Abstract: The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes Joy's law, acting on east-west flows as they rise towards the surface. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we ca… ▽ More The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes Joy's law, acting on east-west flows as they rise towards the surface. Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we can constrain the origins of Joy's law. We measured the tilt angle of the primary bipoles in 153 emerging active regions in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey. We used line-of-sight magnetic field measurements averaged over 6 hours to define the polarities and measure the tilt angle up to four days after emergence. We find that at the time of emergence the polarities are on average aligned east-west, and that neither the separation nor the tilt depends on latitude. We do find, however, that ARs at higher latitudes have a faster north-south separation speed than those closer to the equator at the emergence time. After emergence, the tilt angle increases and Joy's law is evident about two days later. The scatter in the tilt angle is independent of flux until about one day after emergence, when higher-flux regions have a smaller scatter in tilt angle than lower-flux regions. Our finding that active regions emerge with an east-west alignment is surprising since thin flux tube models predict that tilt angles of rising flux tubes are generated below the surface. Previously reported tilt angle relaxation of deeply anchored flux tubes can be largely explained by the change in east-west separation. We conclude that Joy's law is caused by an inherent north-south separation speed present when the flux first reaches the surface, and that the scatter in the tilt angle is consistent with buffeting of the polarities by supergranulation. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: accepted in Astronomy and Astrophysics

Journal ref: A&A 640, A116 (2020)

Average motion of emerging solar active region polarities I: Two phases of emergence

Authors: Hannah Schunker, Aaron C. Birch, Robert H. Cameron, Douglas C. Braun, Laurent Gizon, Raymond B. Burston

Abstract: Our goal is to constrain models of active region formation by tracking the average motion of active region polarity pairs as they emerge onto the surface. We measured the motion of the two main opposite polarities in 153 emerging active regions (EARs) using line-of-sight magnetic field observations from the Solar Dynamics Observatory Helioseismic Emerging Active Region (SDO/HEAR) survey (Schunker… ▽ More Our goal is to constrain models of active region formation by tracking the average motion of active region polarity pairs as they emerge onto the surface. We measured the motion of the two main opposite polarities in 153 emerging active regions (EARs) using line-of-sight magnetic field observations from the Solar Dynamics Observatory Helioseismic Emerging Active Region (SDO/HEAR) survey (Schunker et al. 2016). We first measured the position of each of the polarities eight hours after emergence and tracked their location forwards and backwards in time. We find that, on average, the polarities emerge with an east-west orientation and the separation speed between the polarities increases. At about 0.1 days after emergence, the average separation speed reaches a peak value of 229 +/- 11 m/s, and then starts to decrease, and about 2.5 days after emergence the polarities stop separating. We also find that the separation and the separation speed in the east-west direction are systematically larger for active regions with higher flux. Our results reveal two phases of the emergence process defined by the rate of change of the separation speed as the polarities move apart. Phase 1 begins when the opposite polarity pairs first appear at the surface, with an east-west alignment and an increasing separation speed. We define Phase 2 to begin when the separation speed starts to decrease, and ends when the polarities have stopped separating. This is consistent with the picture of Chen, Rempel, & Fan (2017): the peak of a flux tube breaks through the surface during Phase 1. During Phase 2 the magnetic field lines are straightened by magnetic tension, so that the polarities continue to move apart, until they eventually lie directly above their anchored subsurface footpoints. △ Less

Submitted 30 April, 2019; v1 submitted 28 March, 2019; originally announced March 2019.

Comments: accepted A&A

Journal ref: A&A 625, A53 (2019)

Fragile detection of solar g modes by Fossat et al

Authors: Hannah Schunker, Jesper Schou, Patrick Gaulme, Laurent Gizon

Abstract: The internal gravity modes of the Sun are notoriously difficult to detect, and the claimed detection of gravity modes presented in Fossat et al. 2017 is thus very exciting. Given the importance of these modes for understanding solar structure and dynamics, the results must be robust. While Fossat et al. 2017 described their method and parameter choices in detail, the sensitivity of their results t… ▽ More The internal gravity modes of the Sun are notoriously difficult to detect, and the claimed detection of gravity modes presented in Fossat et al. 2017 is thus very exciting. Given the importance of these modes for understanding solar structure and dynamics, the results must be robust. While Fossat et al. 2017 described their method and parameter choices in detail, the sensitivity of their results to several parameters were not presented. Therefore, we test the sensitivity to a selection of them. The most concerning result is that the detection vanishes when we adjust the start time of the 16.5 year velocity time series by a few hours. We conclude that this reported detection of gravity modes is extremely fragile and should be treated with utmost caution. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: 15 pages, 11 Figures

Observing and modeling the poloidal and toroidal fields of the solar dynamo

Authors: R. H. Cameron, T. L Duvall Jr., M. Schüssler, H. Schunker

Abstract: Context. The solar dynamo consists of a process that converts poloidal field to toroidal field followed by a process which creates new poloidal field from the toroidal field. Aims. Our aim is to observe the poloidal and toroidal fields relevant to the global solar dynamo and see if their evolution is captured by a Babcock-Leighton dynamo. Methods. We use synoptic maps of the surface radial fie… ▽ More Context. The solar dynamo consists of a process that converts poloidal field to toroidal field followed by a process which creates new poloidal field from the toroidal field. Aims. Our aim is to observe the poloidal and toroidal fields relevant to the global solar dynamo and see if their evolution is captured by a Babcock-Leighton dynamo. Methods. We use synoptic maps of the surface radial field from the KPNSO/VT and SOLIS observatories to construct the poloidal field as a function of time and latitude, and Wilcox Solar Observatory and SOHO/MDI full disk images to infer the longitudinally averaged surface azimuthal field. We show that the latter is consistent with an estimate of that due to flux emergence and therefore closely related to the subsurface toroidal field. Results. We present maps of the poloidal and toroidal magnetic field of the global solar dynamo. The longitude-averaged azimuthal field observed at the surface results from flux emergence. At high latitudes this component follows the radial component of the polar fields with a short time lag (1-3 years). The lag increases at lower latitudes. The observed evolution of the poloidal and toroidal magnetic fields is described by the (updated) Babcock-Leighton dynamo model. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: A&A

Journal ref: A&A 609, A56 (2018)

Limits on radial differential rotation in Sun-like stars from parametric fits to oscillation power spectra

Authors: M. B. Nielsen, H. Schunker, L. Gizon, J. Schou, W. H. Ball

Abstract: Rotational shear in Sun-like stars is thought to be an important ingredient in models of stellar dynamos. Thanks to helioseismology, rotation in the Sun is characterized well, but the interior rotation profiles of other Sun-like stars are not so well constrained. Until recently, measurements of rotation in Sun-like stars have focused on the mean rotation, but little progress has been made on measu… ▽ More Rotational shear in Sun-like stars is thought to be an important ingredient in models of stellar dynamos. Thanks to helioseismology, rotation in the Sun is characterized well, but the interior rotation profiles of other Sun-like stars are not so well constrained. Until recently, measurements of rotation in Sun-like stars have focused on the mean rotation, but little progress has been made on measuring or even placing limits on differential rotation. Using asteroseismic measurements of rotation we aim to constrain the radial shear in five Sun-like stars observed by the NASA Kepler mission: KIC004914923, KIC005184732, KIC006116048, KIC006933899, and KIC010963065. We used stellar structure models for these five stars from previous works. These models provide the mass density, mode eigenfunctions, and the convection zone depth, which we used to compute the sensitivity kernels for the rotational frequency splitting of the modes. We used these kernels as weights in a parametric model of the stellar rotation profile of each star, where we allowed different rotation rates for the radiative interior and the convective envelope. This parametric model was incorporated into a fit to the oscillation power spectrum of each of the five Kepler stars. This fit included a prior on the rotation of the envelope, estimated from the rotation of surface magnetic activity measured from the photometric variability. The asteroseismic measurements without the application of priors are unable to place meaningful limits on the radial shear. Using a prior on the envelope rotation enables us to constrain the interior rotation rate and thus the radial shear. In the five cases that we studied, the interior rotation rate does not differ from the envelope by more than approximately +/-30%. Uncertainties in the rotational splittings are too large to unambiguously determine the sign of the radial shear. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: Accepted for publication in Astronomy & Astrophysics. 8 pages, 4 figures, 1 table

Journal ref: A&A 603, A6 (2017)

SDO/HMI survey of emerging active regions for helioseismology

Authors: H. Schunker, D. C. Braun, A. C. Birch, R. B. Burston, L. Gizon

Abstract: Observations from the Solar Dynamics Observatory (SDO) have the potential for allowing the helioseismic study of the formation of hundreds of active regions, which would enable us to perform statistical analyses. Our goal is to collate a uniform data set of emerging active regions observed by the SDO/HMI instrument suitable for helioseismic analysis up to seven days before emergence. We restricted… ▽ More Observations from the Solar Dynamics Observatory (SDO) have the potential for allowing the helioseismic study of the formation of hundreds of active regions, which would enable us to perform statistical analyses. Our goal is to collate a uniform data set of emerging active regions observed by the SDO/HMI instrument suitable for helioseismic analysis up to seven days before emergence. We restricted the sample to active regions that were visible in the continuum and emerged into quiet Sun largely avoiding pre-existing magnetic regions. As a reference data set we paired a control region (CR), with the same latitude and distance from central meridian, with each emerging active region (EAR). We call this data set, which is currently comprised of 105 emerging active regions observed between May 2010 and November 2012, the SDO Helioseismic Emerging Active Region (SDO/HEAR) survey. To demonstrate the utility of a data set of a large number of emerging active regions, we measure the relative east-west velocity of the leading and trailing polarities from the line-of-sight magnetogram maps during the first day after emergence. The latitudinally averaged line-of-sight magnetic field of all the EARs shows that, on average, the leading (trailing) polarity moves in a prograde (retrograde) direction with a speed of 121 +/- 22 m/s (-70 +/- 13 m/s) relative to the Carrington rotation rate in the first day. However, relative to the differential rotation of the surface plasma, the east-west velocity is symmetric, with a mean of 95 +/- 13 m/s. The SDO/HEAR data set will not only be useful for helioseismic studies, but will also be useful to study other features such as the surface magnetic field evolution of a large sample of EARs. △ Less

Submitted 7 September, 2016; v1 submitted 29 August, 2016; originally announced August 2016.

Comments: Accepted by Astronomy and Astrophysics, 11 figures, one longtable; update corrects units in Figure 1

A low upper limit on the subsurface rise speed of solar active regions

Authors: Aaron C. Birch, Hannah Schunker, Douglas C. Braun, Robert Cameron, Laurent Gizon, Björn Löptien, Matthias Rempel

Abstract: Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube, it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around… ▽ More Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube, it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around emerging active regions, in combination with numerical simulations of solar magnetoconvection, can constrain the subsurface rise speed of emerging magnetic flux. The observed flows imply that the rise speed of the magnetic field is no larger than 150 m/s at a depth of 20 Mm, that is, well below the prediction of the (standard) thin flux tube model but in the range expected for convective velocities at this depth. We conclude that convective flows control the dynamics of rising flux tubes in the upper layers of the Sun and cannot be neglected in models of flux emergence. △ Less

Submitted 18 July, 2016; originally announced July 2016.

Comments: Published in Science Advances. DOI: 10.1126/sciadv.1600557

Journal ref: Sci. Adv. 2, e1600557 (2016)

Asteroseismic inversions for radial differential rotation of Sun-like stars: ensemble fits

Authors: Hannah Schunker, Jesper Schou, Warrick Ball, Martin Bo Nielsen, Laurent Gizon

Abstract: Radial differential rotation is an important parameter for stellar dynamo theory and for understanding angular momentum transport. We investigate the potential of using a large number of similar stars simultaneously to constrain their average radial differential rotation gradient: we call this 'ensemble fitting'. We use a range of stellar models along the main sequence, each with a synthetic rotat… ▽ More Radial differential rotation is an important parameter for stellar dynamo theory and for understanding angular momentum transport. We investigate the potential of using a large number of similar stars simultaneously to constrain their average radial differential rotation gradient: we call this 'ensemble fitting'. We use a range of stellar models along the main sequence, each with a synthetic rotation profile. The rotation profiles are step functions with a step of -350 nHz, which is located at the base of the convection zone. These models are used to compute the rotational splittings of the p modes and to model their uncertainties. We then fit an ensemble of stars to infer the average step size. All the uncertainties on the inferred step size for individual stars are of the order 1 micro Hz. Using 15 stellar models in an ensemble fit, we show that the uncertainty on the average step size is reduced to less than the input step size, which allows us to constrain the sign of the radial differential rotation. We show that a solar-like step size (approximately 30 nHz) can be constrained by an ensemble fit of thousands of main-sequence stars. Observing the number of stars required to successfully exploit the ensemble fitting method will be possible with future asteroseismology missions, such as PLATO. We demonstrate the potential of ensemble fitting by showing that any systematic differences in the average step size between F, G, and K-type stars larger than 100 nHz can be detected. △ Less

Submitted 22 December, 2015; originally announced December 2015.

Comments: 9 pages, 1 figure, 6 tables

Journal ref: A&A 586, A79 (2016)

Asteroseismic inversions for radial differential rotation of Sun-like stars: Sensitivity to uncertainties

Authors: Hannah Schunker, Jesper Schou, Warrick H. Ball

Abstract: We quantify the effect of observational spectroscopic and asteroseismic uncertainties on regularised least squares (RLS) inversions for the radial differential rotation of Sun-like and subgiant stars. We first solved the forward problem to model rotational splittings plus the observed uncertainties for models of a Sun-like star, HD 52265, and a subgiant star, KIC 7341231. We randomly perturbed the… ▽ More We quantify the effect of observational spectroscopic and asteroseismic uncertainties on regularised least squares (RLS) inversions for the radial differential rotation of Sun-like and subgiant stars. We first solved the forward problem to model rotational splittings plus the observed uncertainties for models of a Sun-like star, HD 52265, and a subgiant star, KIC 7341231. We randomly perturbed the parameters of the stellar models within the uncertainties of the spectroscopic and asteroseismic constraints and used these perturbed stellar models to compute rotational splittings. We experimented with three rotation profiles: solid body rotation, a step function, and a smooth rotation profile decreasing with radius. We then solved the inverse problem to infer the radial differential rotation profile using a RLS inversion and kernels from the best-fit stellar model.We found that the inversions for Sun-like stars with solar-like radial differential rotation profiles are insensitive to the uncertainties in the stellar models. We found that when the rotation rate below the convection zone is increased to six times that of the surface rotation rate the inferred rotation profile excluded solid body rotation. With the current observational uncertainties, we found that inversions of subgiant stars are sensitive to the uncertainties in the stellar model. Our findings suggest that inversions for the radial differential rotation of subgiant stars would benefit from more tightly constrained stellar models. In Sun-like stars, the insensitivity of the inversions to stellar model uncertainties suggests that it may be possible to perform ensemble inversions for the average radial differential rotation of many stars with a range of stellar types to better constrain the inversions. △ Less

Submitted 26 November, 2015; originally announced November 2015.

Comments: accepted in A&A, 9 pages, 10 figures

Journal ref: A&A 586, A24 (2016)

Constraining differential rotation of Sun-like stars from asteroseismic and starspot rotation periods

Authors: M. B. Nielsen, H. Schunker, L. Gizon, W. H. Ball

Abstract: In previous work we identified six Sun-like stars observed by Kepler with exceptionally clear asteroseismic signatures of rotation. Here, we show that five of these stars exhibit surface variability suitable for measuring rotation. In order to further constrain differential rotation, we compare the rotation periods obtained from light-curve variability with those from asteroseismology. The two rot… ▽ More In previous work we identified six Sun-like stars observed by Kepler with exceptionally clear asteroseismic signatures of rotation. Here, we show that five of these stars exhibit surface variability suitable for measuring rotation. In order to further constrain differential rotation, we compare the rotation periods obtained from light-curve variability with those from asteroseismology. The two rotation measurement methods are found to agree within uncertainties, suggesting that radial differential rotation is weak, as is the case for the Sun. Furthermore, we find significant discrepancies between ages from asteroseismology and from three different gyrochronology relations, implying that stellar age estimation is problematic even for Sun-like stars. △ Less

Submitted 21 August, 2015; v1 submitted 10 August, 2015; originally announced August 2015.

Comments: Accepted for publication in A&A. 5 pages, 4 figures

Measuring stellar rotation periods with Kepler

Authors: M. B. Nielsen, L. Gizon, H. Schunker, C. Karoff

Abstract: We measure rotation periods for 12151 stars in the Kepler field, based on the photometric variability caused by stellar activity. Our analysis returns stable rotation periods over at least six out of eight quarters of Kepler data. This large sample of stars enables us to study the rotation periods as a function of spectral type. We find good agreement with previous studies and vsini measurements f… ▽ More We measure rotation periods for 12151 stars in the Kepler field, based on the photometric variability caused by stellar activity. Our analysis returns stable rotation periods over at least six out of eight quarters of Kepler data. This large sample of stars enables us to study the rotation periods as a function of spectral type. We find good agreement with previous studies and vsini measurements for F, G and K stars. Combining rotation periods, B-V color, and gyrochronology relations, we find that the cool stars in our sample are predominantly younger than ~1Gyr. △ Less

Submitted 31 March, 2015; originally announced March 2015.

Comments: Published 12/2013 in ASP Conference Series, Vol. 479, Progress in Physics of the Sun and Stars, p. 137. http://aspbooks.org/custom/publications/paper/479-0137.html

Journal ref: 2013ASPC..479..137N

Rotational splitting as a function of mode frequency for six Sun-like stars

Authors: M. B. Nielsen, L. Gizon, H. Schunker, J. Schou

Abstract: Asteroseismology offers the prospect of constraining differential rotation in Sun-like stars. Here we have identified six high signal-to-noise main-sequence Sun-like stars in the Kepler field, which all have visible signs of rotational splitting of their p-mode frequencies. For each star, we extract the rotational frequency splitting and inclination angle from separate mode sets (adjacent modes wi… ▽ More Asteroseismology offers the prospect of constraining differential rotation in Sun-like stars. Here we have identified six high signal-to-noise main-sequence Sun-like stars in the Kepler field, which all have visible signs of rotational splitting of their p-mode frequencies. For each star, we extract the rotational frequency splitting and inclination angle from separate mode sets (adjacent modes with l=2, 0, and 1) spanning the p-mode envelope. We use a Markov chain Monte Carlo method to obtain the best fit and errors associated with each parameter. We are able to make independent measurements of rotational splittings of ~8 radial orders for each star. For all six stars, the measured splittings are consistent with uniform rotation, allowing us to exclude large radial differential rotation. This work opens the possibility of constraining internal rotation of Sun-like stars. △ Less

Submitted 9 September, 2014; v1 submitted 19 August, 2014; originally announced August 2014.

Comments: Published in Astronomy and Astrophysics. 4 pages, 3 figures

Journal ref: A&A 568, L12 (2014)

The PLATO 2.0 Mission

Authors: H. Rauer, C. Catala, C. Aerts, T. Appourchaux, W. Benz, A. Brandeker, J. Christensen-Dalsgaard, M. Deleuil, L. Gizon, M. -J. Goupil, M. Güdel, E. Janot-Pacheco, M. Mas-Hesse, I. Pagano, G. Piotto, D. Pollacco, N. C. Santos, A. Smith, J. -C., Suárez, R. Szabó, S. Udry, V. Adibekyan, Y. Alibert, J. -M. Almenara , et al. (137 additional authors not shown)

Abstract: PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small ap… ▽ More PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0. △ Less

Submitted 4 March, 2014; v1 submitted 2 October, 2013; originally announced October 2013.

Comments: 63 pages, 17 figures, submitted to Experimental Astronomy (ExA)

Rotation periods of 12 000 main-sequence Kepler stars: Dependence on stellar spectral type and comparison with v sin i observations

Authors: M. B. Nielsen, L. Gizon, H. Schunker, C. Karoff

Abstract: Aims: We aim to measure the starspot rotation periods of active stars in the Kepler field as a function of spectral type and to extend reliable rotation measurements from F-, G-, and K-type to M-type stars. Methods: Using the Lomb-Scargle periodogram we searched more than 150 000 stellar light curves for periodic brightness variations. We analyzed periods between 1 and 30 days in eight consecuti… ▽ More Aims: We aim to measure the starspot rotation periods of active stars in the Kepler field as a function of spectral type and to extend reliable rotation measurements from F-, G-, and K-type to M-type stars. Methods: Using the Lomb-Scargle periodogram we searched more than 150 000 stellar light curves for periodic brightness variations. We analyzed periods between 1 and 30 days in eight consecutive Kepler quarters, where 30 days is an estimated maximum for the validity of the PDC_MAP data correction pipeline. We selected stable rotation periods, i.e., periods that do not vary from the median by more than one day in at least six of the eight quarters. We averaged the periods for each stellar spectral class according to B - V color and compared the results to archival vsini data, using stellar radii estimates from the Kepler Input Catalog. Results: We report on the stable starspot rotation periods of 12 151 Kepler stars. We find good agreement between starspot velocities and vsini data for all F-, G- and early K-type stars. The 795 M-type stars in our sample have a median rotation period of 15.4 days. We find an excess of M-type stars with periods less than 7.5 days that are potentially fast-rotating and fully convective. Measuring photometric variability in multiple Kepler quarters appears to be a straightforward and reliable way to determine the rotation periods of a large sample of active stars, including late-type stars. △ Less

Submitted 4 September, 2013; v1 submitted 24 May, 2013; originally announced May 2013.

Comments: 5 pages, 3 figures. Online material is available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/557/L10

Helioseismology of sunspots: how sensitive are travel times to the Wilson depression and to the subsurface magnetic field?

Authors: Hannah Schunker, Laurent Gizon, Robert H. Cameron, Aaron C. Birch

Abstract: In order to assess the ability of helioseismology to probe the subsurface structure and magnetic field of sunspots, we need to determine how helioseismic travel times depend on perturbations to sunspot models. Here we numerically simulate the propagation of f, p1, and p2 wave packets through magnetic sunspot models. Among the models we considered, a ~50 km change in the height of the Wilson depres… ▽ More In order to assess the ability of helioseismology to probe the subsurface structure and magnetic field of sunspots, we need to determine how helioseismic travel times depend on perturbations to sunspot models. Here we numerically simulate the propagation of f, p1, and p2 wave packets through magnetic sunspot models. Among the models we considered, a ~50 km change in the height of the Wilson depression and a change in the subsurface magnetic field geometry can both be detected above the observational noise level. We also find that the travel-time shifts due to changes in a sunspot model must be modeled by computing the effects of changing the reference sunspot model, and not by computing the effects of changing the subsurface structure in the quiet-Sun model. For p1 modes the latter is wrong by a factor of four. In conclusion, numerical modeling of MHD wave propagation is an essential tool for the interpretation of the effects of sunspots on seismic waveforms. △ Less

Submitted 25 March, 2013; originally announced March 2013.

Comments: 5 pages, 3 figures: submitted to A&A

Time--distance inversions for horizontal and vertical flows on supergranular scales applied to MDI and HMI data

Authors: Michal Svanda, Hannah Schunker, Raymond Burston

Abstract: We study the possibility of consistent extension of MDI full-disc helioseismic campaigns with the growing data set of HMI observations. To do so, we down-sample and filter the HMI Dopplegrams so that the resulting spatial power spectrum is similar to the spatial power spectrum of MDI full-disc Dopplergrams. The set of co-spatial and co-temporal datacube pairs from both instruments containing no mi… ▽ More We study the possibility of consistent extension of MDI full-disc helioseismic campaigns with the growing data set of HMI observations. To do so, we down-sample and filter the HMI Dopplegrams so that the resulting spatial power spectrum is similar to the spatial power spectrum of MDI full-disc Dopplergrams. The set of co-spatial and co-temporal datacube pairs from both instruments containing no missing and no bad frames were processed using the same codes and inverted independently for all three components of the plasma flow in the near surface layers. The results from the two instruments are highly correlated, however systematically larger (by ~20%) flow magnitudes are derived from HMI. We comment that this may be an effect of the different formation depth of the Doppler signal △ Less

Submitted 4 February, 2013; originally announced February 2013.

Comments: 6 pages, 4 figures, accepted in Journal of Physics: Conference Series

Properties of high-frequency wave power halos around active regions: an analysis of multi-height data from HMI and AIA onboard SDO

Authors: S. P. Rajaguru, S. Couvidat, Xudong Sun, K. Hayashi, H. Schunker

Abstract: We study properties of waves of frequencies above the photospheric acoustic cut-off of $\approx$5.3 mHz, around four active regions, through spatial maps of their power estimated using data from Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The wavelength channels 1600 Å and 1700 Å from AIA are now known to capture clear osc… ▽ More We study properties of waves of frequencies above the photospheric acoustic cut-off of $\approx$5.3 mHz, around four active regions, through spatial maps of their power estimated using data from Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The wavelength channels 1600 Å and 1700 Å from AIA are now known to capture clear oscillation signals due to helioseismic p modes as well as waves propagating up through to the chromosphere. Here we study in detail, in comparison with HMI Doppler data, properties of the power maps, especially the so called 'acoustic halos' seen around active regions, as a function of wave frequencies, inclination and strength of magnetic field (derived from the vector field observations by HMI) and observation height. We infer possible signatures of (magneto-)acoustic wave refraction from the observation height dependent changes, and hence due to changing magnetic strength and geometry, in the dependences of power maps on the photospheric magnetic quantities. We discuss the implications for theories of p mode absorption and mode conversions by the magnetic field. △ Less

Submitted 24 September, 2012; v1 submitted 25 June, 2012; originally announced June 2012.

Comments: 22 pages, 12 figures, Accepted by journal Solar Physics

Constructing and Characterising Solar Structure Models for Computational Helioseismology

Authors: Hannah Schunker, Robert Cameron, Laurent Gizon, Hamed Moradi

Abstract: In this paper, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science, 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose t… ▽ More In this paper, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science, 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose to make additional corrections to the sound speed of Model S before stabilisation. No stabilised model can be perfectly solar-like, so we present three stabilised models with slightly different eigenmodes. The models are appropriate to study the f and p1 to p4 modes with spherical harmonic degrees in the range from 400 to 900. Background model CSM has a modified pressure gradient for stabilisation and has eigenfrequencies within 2% of Model S. Model CSM_A has an additional 10% increase in sound speed in the top 1 Mm resulting in eigenfrequencies within 2% of Model S and eigenfunctions that are, in comparison with CSM, closest to those of Model S. Model CSM_B has a 3% decrease in sound speed in the top 5 Mm resulting in eigenfrequencies within 1% of Model S and eigenfunctions that are only marginally adversely affected. These models are useful to study the interaction of solar waves with embedded three-dimensional heterogeneities, such as convective flows and model sunspots. We have also calculated the response of the stabilised models to excitation by random near-surface sources, using simulations of the propagation of linear waves. We find that the simulated power spectra of wave motion are in good agreement with an observed SOHO/MDI power spectrum. Overall, our convectively stabilised background models provide a good basis for quantitative numerical local helioseismology. The models are available for download from http://www.mps.mpg.de/projects/seismo/NA4/. △ Less

Submitted 6 June, 2011; v1 submitted 1 May, 2011; originally announced May 2011.

Comments: 35 pages, 23 figures Changed title Updated Figure 18

Scattering of helioseismic waves by a sunspot: wavefront healing and folding

Authors: Z. -C. Liang, L. Gizon, H. Schunker

Abstract: We observe and characterize the scattering of acoustic wave packets by a sunspot, in a regime where the wavelength is comparable to the size of the sunspot. Spatial maps of wave traveltimes and amplitudes are measured from the cross-covariance function of the random wave field. The averaging procedure is such that incoming wave packets are plane wave packets. Observations show that the magnitude o… ▽ More We observe and characterize the scattering of acoustic wave packets by a sunspot, in a regime where the wavelength is comparable to the size of the sunspot. Spatial maps of wave traveltimes and amplitudes are measured from the cross-covariance function of the random wave field. The averaging procedure is such that incoming wave packets are plane wave packets. Observations show that the magnitude of the traveltime perturbation caused by the sunspot diminishes as waves propagate away from the sunspot -- a finite-wavelength phenomenon known as wavefront healing. Observations also show a reduction of the amplitude of the waves after their passage through the sunspot. A significant fraction of this amplitude reduction is due to the defocusing of wave energy by the fast wave-speed perturbation introduced by the sunspot. This ``geometrical attenuation'' will contribute to the wave amplitude reduction in addition to the physical absorption of waves. In addition, we observe an enhancement of wave amplitude away from the central path: diffracted rays intersect with unperturbed rays (caustics) and wavefronts fold and triplicate. Thus we find that ray tracing is useful to interpret these phenomena, although it cannot explain wavefront healing. △ Less

Submitted 1 October, 2010; originally announced October 2010.

Comments: 4 pages, 4 figures, submitted to ApJL

Local helioseismology and the active Sun

Authors: Hannah Schunker

Abstract: The goal of local helioseismology is to elicit three-dimensional information about the sub-surface (or far-side) structure and dynamics of the Sun from observations of the helioseismic wave field at the surface. The physical quantities of interest include flows, sound-speed deviations and magnetic fields. However, strong surface magnetic fields induce large perturbations to the waves making inver… ▽ More The goal of local helioseismology is to elicit three-dimensional information about the sub-surface (or far-side) structure and dynamics of the Sun from observations of the helioseismic wave field at the surface. The physical quantities of interest include flows, sound-speed deviations and magnetic fields. However, strong surface magnetic fields induce large perturbations to the waves making inversions difficult to interpret. The purpose of this paper is to outline the methods of analysis used in local helioseismology, review discoveries associated with the magnetic Sun made using local helioseismology from the past three years, and highlight the efforts towards imaging the interior in the presence of strong magnetic fields. △ Less

Submitted 14 May, 2010; originally announced May 2010.

Comments: 6 pages, 4th HELAS International Conference, Lanzarote, Spain, 1-5 February 2010

Acoustic power absorption and enhancement generated by slow and fast MHD waves

Authors: R. Simoniello, W. Finsterle, R. A. Garcia, D. Salabert, A. Jimenez, Y. Elsworth, H. Schunker

Abstract: We used long duration, high quality, unresolved (Sun-as-a star) observations collected by the ground based network BiSON and by the instruments GOLF and VIRGO on board the ESA/NASA SOHO satellite to search for solar-cycle-related changes in mode characteristics in velocity and continuum intensity for the frequency range between 2.5mHz < nu < 6.8mHz. Over the ascending phase of solar cycle 23 we f… ▽ More We used long duration, high quality, unresolved (Sun-as-a star) observations collected by the ground based network BiSON and by the instruments GOLF and VIRGO on board the ESA/NASA SOHO satellite to search for solar-cycle-related changes in mode characteristics in velocity and continuum intensity for the frequency range between 2.5mHz < nu < 6.8mHz. Over the ascending phase of solar cycle 23 we found a suppression in the p-mode amplitudes both in the velocity and intensity data between 2.5mHz <nu< 4.5mHz with a maximum suppression for frequencies in the range between 2.5mHz <nu< 3.5mHz. The size of the amplitude suppression is 13+-2 per cent for the velocity and 9+-2 per cent for the intensity observations. Over the range 4.5mHz <nu< 5.5mHz the findings hint within the errors to a null change both in the velocity and intensity amplitudes. At still higher frequencies, in the so called High-frequency Interference Peaks (HIPs) between 5.8mHz <nu < 6.8mHz, we found an enhancement in the velocity amplitudes with the maximum 36+-7 per cent occurring for 6.3mHz <nu< 6.8mHz. However, in intensity observations we found a rather smaller enhancement of about 5+-2 per cent in the same interval. There is evidence that the frequency dependence of solar-cycle velocity amplitude changes is consistent with the theory behind the mode conversion of acoustic waves in a non-vertical magnetic field, but there are some problems with the intensity data, which may be due to the height in the solar atmosphere at which the VIRGO data are taken. △ Less

Submitted 20 April, 2010; v1 submitted 19 April, 2010; originally announced April 2010.

Comments: Accepted for publication in A&A. 10 pages, 9 figures.

Constructing semi-empirical sunspot models for helioseismology

Authors: R. Cameron, L. Gizon, H. Schunker, A. Pietarila

Abstract: One goal of helioseismology is to determine the subsurface structure of sunspots. In order to do so, it is important to understand first the near-surface effects of sunspots on solar waves, which are dominant. Here we construct simplified, cylindrically-symmetric sunspot models, which are designed to capture the magnetic and thermodynamics effects coming from about 500 km below the quiet-Sun… ▽ More One goal of helioseismology is to determine the subsurface structure of sunspots. In order to do so, it is important to understand first the near-surface effects of sunspots on solar waves, which are dominant. Here we construct simplified, cylindrically-symmetric sunspot models, which are designed to capture the magnetic and thermodynamics effects coming from about 500 km below the quiet-Sun $τ_{5000}=1$ level to the lower chromosphere. We use a combination of existing semi-empirical models of sunspot thermodynamic structure (density, temperature, pressure): the umbral model of Maltby et al. (1986) and the penumbral model of Ding and Fang (1989). The OPAL equation of state tables are used to derive the sound speed profile. We smoothly merge the near-surface properties to the quiet-Sun values about 1mm below the surface. The umbral and penumbral radii are free parameters. The magnetic field is added to the thermodynamic structure, without requiring magnetostatic equilibrium. The vertical component of the magnetic field is assumed to have a Gaussian horizontal profile, with a maximum surface field strength fixed by surface observations. The full magnetic field vector is solenoidal and determined by the on-axis vertical field, which, at the surface, is chosen such that the field inclination is 45$^\circ$ at the umbral-penumbral boundary. We construct a particular sunspot model based on SOHO/MDI observations of the sunspot in active region NOAA 9787. The helioseismic signature of the model sunspot is studied using numerical simulations of the propagation of f, p$_1$, and p$_2$ wave packets. These simulations are compared against cross-covariances of the observed wave field. We find that the sunspot model gives a helioseismic signature that is similar to the observations. △ Less

Submitted 2 March, 2010; originally announced March 2010.

Comments: 18 pages, 9 Figures submitted to Solar Physics

Surface magnetic field effects in local helioseismology

Authors: H. Schunker, D. C. Braun, P. S. Cally

Abstract: Using helioseismic holography strong evidence is presented that the phase (or equivalent travel-time) of helioseismic signatures in Dopplergrams within sunspots depend upon the line-of-sight angle in the plane containing the magnetic field and vertical directions. This is shown for the velocity signal in the penumbrae of two sunspots at 3, 4 and 5 mHz. Phase-sensitive holography demonstrates tha… ▽ More Using helioseismic holography strong evidence is presented that the phase (or equivalent travel-time) of helioseismic signatures in Dopplergrams within sunspots depend upon the line-of-sight angle in the plane containing the magnetic field and vertical directions. This is shown for the velocity signal in the penumbrae of two sunspots at 3, 4 and 5 mHz. Phase-sensitive holography demonstrates that they are significantly affected in a strong, moderately inclined magnetic field. This research indicates that the effects of the surface magnetic field are potentially very significant for local helioseismic analysis of active regions. △ Less

Submitted 11 February, 2010; originally announced February 2010.

Comments: 6 pages, 8 figures

Journal ref: Astronomische Nachrichten, Volume 328, Issue 3, page 292, 2007

Helioseismology of Sunspots: A Case Study of NOAA Region 9787

Authors: L. Gizon, H. Schunker, C. S. Baldner, S. Basu, A. C. Birch, R. S. Bogart, D. C. Braun, R. Cameron, T. L. Duvall Jr., S. M. Hanasoge, J. Jackiewicz, M. Roth, T. Stahn, M. J. Thompson, S. Zharkov

Abstract: Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate… ▽ More Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot. △ Less

Submitted 11 February, 2010; v1 submitted 11 February, 2010; originally announced February 2010.

Comments: 28 pages, 18 figures

Journal ref: Space Science Reviews, Volume 144, Issue 1-4, page 249, 2009

Convectively stabilised background solar models for local helioseismology

Authors: H. Schunker, R. Cameron, L. Gizon

Abstract: In local helioseismology numerical simulations of wave propagation are useful to model the interaction of solar waves with perturbations to a background solar model. However, the solution to the equations of motions include convective modes that can swamp the waves we are interested in. For this reason, we choose to first stabilise the background solar model against convection by altering the ve… ▽ More In local helioseismology numerical simulations of wave propagation are useful to model the interaction of solar waves with perturbations to a background solar model. However, the solution to the equations of motions include convective modes that can swamp the waves we are interested in. For this reason, we choose to first stabilise the background solar model against convection by altering the vertical pressure gradient. Here we compare the eigenmodes of our convectively stabilised model with a standard solar model (Model S) and find a good agreement. △ Less

Submitted 9 February, 2010; originally announced February 2010.

Comments: 3 pages, 3 figures, HELAS NA3, The Acoustic Solar Cycle, Birmingham, 6-8 January 2009

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

Newly identified properties of surface acoustic power

Authors: H. Schunker, D. C. Braun

Abstract: The cause of enhanced acoustic power surrounding active regions, the acoustic halo, is not as yet understood. We explore the properties of the enhanced acoustic power observed near disk center from 21 to 27 January 2002, including AR 9787. We find that (i) there exists a strong correlation of the enhanced high frequency power with magnetic-field inclination, with greater power in more horizontal… ▽ More The cause of enhanced acoustic power surrounding active regions, the acoustic halo, is not as yet understood. We explore the properties of the enhanced acoustic power observed near disk center from 21 to 27 January 2002, including AR 9787. We find that (i) there exists a strong correlation of the enhanced high frequency power with magnetic-field inclination, with greater power in more horizontal fields, (ii) the frequency of the maximum enhancement increases along with magnetic field strength, and (iii) the oscillations contributing to the halos show modal ridges which are shifted to higher wavenumber at constant frequency in comparison to the ridges of modes in the quiet-Sun. △ Less

Submitted 1 March, 2010; v1 submitted 16 November, 2009; originally announced November 2009.

Comments: 16 pages, 10 figures, submitted to solar physics

HELAS Local Helioseismology Activities

Authors: Hannah Schunker, Laurent Gizon

Abstract: The main goals of the HELAS local helioseismology network activity are to consolidate this field of research in Europe, to organise scientific workshops, and to facilitate the distribution of observations and data analysis software. Most of this is currently accomplished via a dedicated website at http://www.mps.mpg.de/projects/seismo/NA4/ . In this paper we list the outreach material, observati… ▽ More The main goals of the HELAS local helioseismology network activity are to consolidate this field of research in Europe, to organise scientific workshops, and to facilitate the distribution of observations and data analysis software. Most of this is currently accomplished via a dedicated website at http://www.mps.mpg.de/projects/seismo/NA4/ . In this paper we list the outreach material, observational data, analysis tools and modelling tools currently available from the website and describe the focus of the scientific workshops and their proceedings. △ Less

Submitted 19 October, 2008; originally announced October 2008.

Comments: 14 pages, 5 figures

Physical Properties of Wave Motion in Inclined Magnetic Fields Within Sunspot Penumbrae

Authors: H. Schunker, D. C. Braun, C. Lindsey, P. S. Cally

Abstract: At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic… ▽ More At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic waves and the local surface motion within two sunspots. These correlations are modeled assuming the surface motion is elliptical, and we explore the properties of the elliptical motion on the magnetic field inclination. We also demonstrate that the phase shift of the outward propagating waves is opposite to the phase shift of the inward propagating waves in stronger, more vertical fields, but similar to the inward phase shifts in weaker, more inclined fields. △ Less

Submitted 22 April, 2008; v1 submitted 29 January, 2008; originally announced January 2008.

Comments: 22 pages, 13 figures