Mon. Not. R. Astron. Soc. 000, 1–18 (2010) Printed 27 August 2018 (MN LATEX style file v2.2) A more realistic representation of overshoot at the base of the solar convective envelope as seen by helioseismology J. Christensen-Dalsgaard1,2⋆, M.J.P.F.G. Monteiro3,4, M. Rempel2, M.J. Thompson2,5 1Danish AsteroSeismology Centre, and Department of Physics and Astronomy, Aarhus University, Ny Munkegade, DK-8000 Aarhus C, Denmark 2High Altitude Observatory, NCAR†, Boulder CO 80307-3000, USA 3Centro de Astrof´ısica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 4Departamento de F´ısica e Astronomia, Faculdade de Ciˆencias, Universidade do Porto, Portugal 5School of Mathematics & Statistics, University of Sheffield, Hounsfield Road, Sheffield S3 7RH, United Kingdom Accepted ; Received ; in original form 2009 September 11 ABSTRACT The stratification near the base of the Sun’s convective envelope is governed by processes of convective overshooting and element diffusion, and the region is widely believed to play a key role in the solar dynamo. The stratification in that region gives rise to a characteristic signal in the frequencies of solar p modes, which has been used to determine the depth of the solar convection zone and to investigate the extent of convective overshoot. Previous helioseismic investigations have shown that the Sun’s spherically symmetric stratification in this region is smoother than that in a standard solar model without overshooting, and have ruled out simple models incorporating overshooting, which extend the region of adiabatic stratification and have a more-or-less abrupt transition to subadiabatic stratification at the edge of the overshoot region. In this paper we consider physically motivated models which have a smooth transition in stratification bridging the region from the lower convection zone to the radiative interior beneath. We find that such a model is in better agreement with the helioseismic data than a standard solar model. Key words: convection – Sun: helioseismology – Sun: interior – stars: interior 1 INTRODUCTION conditions inside the Sun, because the frequencies of resonant modes set up by acoustic waves propagating in the solar interior An understanding of the overshoot region at the bottom of the Sun’s depend in particular on the adiabatic sound speed c which is given convective envelope is important for a number of reasons. The over- by shoot region approximately coincides with the solar tachocline, a Γ p Γ k T region of rotational shear which is generally believed to play a key c2 = 1 ≃ 1 B . (1) role in the solar dynamo: overshooting is likely to be important for ρ muµ helping to store the magnetic flux below the convection zone during Here Γ1 is the logarithmic derivative of pressure p with respect to the solar cycle. Bulk motion in the overshoot region also affects the density ρ at constant specific entropy, T is temperature, µ is the thermal stratification and it may contribute to significant mixing of mean molecular weight, kB is Boltzmann’s constant and mu is the arXiv:1102.0235v2 [astro-ph.SR] 20 Feb 2011 chemical elements, for example transporting fragile elements such atomic mass unit. Hence the sound-speed gradient with respect to as lithium to hotter regions where they are destroyed more easily depth depends on the temperature gradient, which itself depends than in the convection zone. More generally, convective overshoot on the mechanism by which heat is transported. The transition be- in stars (particularly those with convective cores) is likely to be an tween fully radiative heat transport beneath the convection zone important and as yet imperfectly understood process affecting age and convective heat transport within the convection zone is mani- estimates of stars, and so improved constraints on theories of over- fest in the temperature gradient and hence too in the sound-speed shooting obtained from a study of how overshooting works in the gradient. If the transition in sound-speed gradient takes place over solar case can be important for understanding stars more widely. a distance that is small compared with the vertical wavelength of Helioseismology provides a means of probing directly the the acoustic waves near the base of the convection zone, then the transition appears to the waves to be more-or-less sharp and this gives rise to an oscillatory signal in the mode frequencies ω: the ⋆ E-mail:[email protected]; [email protected]; [email protected]; [email protected] form of the signal gives information about the location and nature † The National Center for Atmospheric Research is operated by the Uni- of this “acoustic glitch” (Gough 2002a). versity Corporation for Atmospheric Research under sponsorship of the Na- Monteiro et al. (1994) represented the effect of the base of the tional Science Foundation convection zone in terms of an additive contribution δωp to the c 2010 RAS 2 J. Christensen-Dalsgaard et al. frequencies, relative to those of a corresponding model in which one tenth of a pressure scale-height. A similar limit was found by the transition had been smoothed out. For low-degree modes (see Basu et al. (1994), while Roxburgh & Vorontsov (1994) obtained a Monteiro et al. 2000) a break in the first derivative of the sound somewhat weaker limit. Basu & Antia (1994) noted that the com- speed gives rise to a contribution of the form position gradient caused by the inclusion of helium settling pro- duced a sharper transition in the sound-speed gradient, even in δω = A(ω) cos (2ωτ¯ + 2φ ) , (2) p d 0 models without overshoot, and hence a larger oscillatory signal. while a break in the second derivative gives rise to a similar signal From these analyses it would appear that the transition in sound- but with a sine term instead of a cosine, and with a different fre- speed gradient at the base of the solar convection zone is if any- quency dependence of A(ω). Hereτ ¯d is essentially the value of τ at thing smoother than in the non-overshoot models. A caveat is that the location of the acoustic glitch, where what we purport to measure in the above studies is the spherically symmetric component of the structure: departures from sphericity, R dr′ τ = (3) such as a latitudinal dependence to the shape of the base of the Zr c convection zone, could make the transition appear smoother than is acoustic depth beneath the surface, r being the corresponding dis- it is locally. The helioseismic evidence, however, is that the loca- tance to the centre and R the surface radius of the Sun. Also, φ0 is a tion of the base of the convection zone is independent of latitude phase introduced by the reflection of the mode at the turning points, (Monteiro & Thompson 1998; Basu & Antia 2001). Changes of the depending in particular on the near-surface structure. The function base of the convection zone on time scales shorter than the obser- A(ω) is an amplitude which depends on the sharpness and nature of vation interval could also have a similar effect by introducing a the convection-zone base: the smoother the transition, the smaller time-averaging effect on the mode frequencies that would mimic in general will be the amplitude. However, if moderate-degree data a smoother transition. The importance of this effect is difficult to are used, as in the case of Sun where we have accurate data for estimate, however, and depends strongly on the 3D nature of con- modes whose degree is above 3, the above expression needs to in- vection at the base of the envelope. clude additional terms, both in the amplitude (Monteiro et al. 1994) Overshoot has been addressed in the last two decades by a va- and in the argument of the signal (Christensen-Dalsgaard et al. riety of 2D and 3D numerical simulations (Roxburgh & Simmons 1995), to account for the first-order effect of the mode degree on 1993; Hurlburt et al. 1994; Singh et al. 1995, 1998; Saikia et al. the signal. 2000; Brummell et al. 2002; Rogers & Glatzmaier 2005b,a). Overshoot at the base of the solar convection has traditionally Whilst the non-local mixing-length models have clearly predicted been modelled using non-local mixing-length theory (e.g., Zahn an adiabatic overshoot region of a sizeable fraction of a pressure 1991). Such models mostly predict an overshoot region that is scale height and a rather sharp transition to the radiative zone be- nearly adiabatically stratified; in terms of the logarithmic tempera- neath, the numerical simulations show a greater variety of possi- −6 ture gradient ∇ = d ln T/d ln p one finds that δ ≡∇−∇ad ∼−10 , ble behaviours. The work by Brummell et al. (2002) is currently where ∇ad = (∂ ln T/∂ ln p)s is the adiabatic temperature gradient, one of the best resolved and most turbulent investigations: it shows s being specific entropy. The depth of overshoot region is typi- strongly subadiabatic overshoot with very smooth transition to- cally between 0.2Hp and 0.4Hp, where Hp is the pressure scale wards the radiative temperature gradient. Most of the earlier 2D height at the base of the convection zone, with a very steep tran- and 3D simulation were more in the laminar regime and found, sition towards the radiative temperature gradient. The results have depending on their parameters (mainly the stiffness of the subadia- seemed rather robust, since the above behaviour is found in models batic layer), both nearly adiabatic overshoot and extended subadia- incorporating quite different large-scale flow structures: for exam- batic overshoot. ple van Ballegooijen (1982) assumed overturning convective rolls Rempel (2004) tried to understand the cause of the discrepan- whilst Schmitt et al.