
Sunspots are of basic interest in the study of the Sun. Their relevance
ranges from them being an activity indicator of magnetic fields to being the
place where coronal mass ejections and flares erupt. They are therefore also an
important ingredient of space weather. Their formation, however, is still an
unresolved problem in solar physics. Observations utilize just 2D surface
information near the spot, but it is debatable how to infer deep structures and
properties from local helioseismology. For a long time, it was believed that
flux tubes rising from the bottom of the convection zone are the origin of the
bipolar sunspot structure seen on the solar surface. However, this theory has
been challenged, in particular recently by new surface observation,
helioseismic inversions, and numerical models of convective dynamos. In this
article we discuss another theoretical approach to the formation of sunspots:
the negative effective magnetic pressure instability. This is a largescale
instability, in which the total (kinetic plus magnetic) turbulent pressure can
be suppressed in the presence of a weak largescale magnetic field, leading to
a converging downflow, which eventually concentrates the magnetic field within
it. Numerical simulations of forced stratified turbulence have been able to
produce strong superequipartition flux concentrations, similar to sunspots at
the solar surface. In this framework, sunspots would only form close to the
surface due to the instability constraints on stratification and rotation.
Additionally, we present some ideas from local helioseismology, where we plan
to use the Hankel analysis to study the preemergence phase of a sunspot and to
constrain its deep structure and formation mechanism.

The meanfield theory of magnetized stellar convection gives rise to the two
possibility of distinct instabilities: the largescale dynamo instability,
operating in the bulk of the convection zone, and a negative effective magnetic
pressure instability (NEMPI) operating in the strongly stratified surface
layers. The latter might be important in connection with magnetic spot
formation, but the growth rate of NEMPI is suppressed with increasing rotation
rates, although recent direct numerical simulations (DNS) have shown a
subsequent increase in the growth rate. We examine quantitatively whether this
increase in the growth rate of NEMPI can be explained by an alpha squared
meanfield dynamo, and whether both NEMPI and the dynamo instability can
operate at the same time. We use both DNS and meanfield simulations (MFS) to
solve the underlying equations numerically either with or without an imposed
horizontal field. We use the testfield method to compute relevant dynamo
coefficients. DNS show that magnetic flux concentrations are still possible up
to rotation rates above which the largescale dynamo effect produces mean
magnetic fields. The resulting DNS growth rates are quantitatively well
reproduced with MFS. As expected, for weak or vanishing rotation, the growth
rate of NEMPI increases with increasing gravity, but there is a correction term
for strong gravity and large turbulent magnetic diffusivity. Magnetic flux
concentrations are still possible for rotation rates above which dynamo action
takes over. For the solar rotation rate, the corresponding turbulent turnover
time is about 5 hours, with dynamo action commencing in the layers beneath.

The surface layers of the Sun are strongly stratified. In the presence of
turbulence with a weak mean magnetic field, a largescale instability resulting
in the formation of nonuniform magnetic structures, can be excited over the
scale of many turbulent eddies or convection cells. This instability is caused
by a negative contribution of turbulence to the effective (meanfield) magnetic
pressure and has previously been discussed in connection with the formation of
active regions and perhaps sunspots. We want to understand the effects of
rotation on this instability in both two and three dimensions. We use
meanfield magnetohydrodynamics in a parameter regime in which the properties
of the negative effective magnetic pressure instability have previously been
found to be in agreement with those of direct numerical simulations. We find
that the instability is suppressed already for relatively slow rotation with
Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is
strongest at the equator. In the nonlinear regime, we find traveling wave
solutions with propagation in the prograde direction at the equator with
additional poleward migration away from the equator. The prograde rotation of
the magnetic pattern near the equator is argued to be a possible explanation
for the faster rotation speed of magnetic tracers found on the Sun. In the bulk
of the domain, kinetic and current helicities are negative in the northern
hemisphere and positive in the southern.