{"id":677,"date":"2017-01-16T03:15:04","date_gmt":"2017-01-16T03:15:04","guid":{"rendered":"http:\/\/crediblehulk.org\/?p=677"},"modified":"2018-06-30T07:03:09","modified_gmt":"2018-06-30T07:03:09","slug":"mean-field-theory-and-solar-dynamo-modeling","status":"publish","type":"post","link":"https:\/\/www.crediblehulk.org\/index.php\/2017\/01\/16\/mean-field-theory-and-solar-dynamo-modeling\/","title":{"rendered":"Mean Field Theory and Solar Dynamo Modeling"},"content":{"rendered":"

In a recent post<\/a>, I talked about the characteristics of the sun\u2019s 11 and 22 year cycles, the observed laws which describe the behavior of the sunspot cycle, how proxy data is used to reconstruct a record of solar cycles of the past, Grand Solar Maxima and Minima, the relationship between Total Solar Irradiance (TSI<\/a>) and the sunspot cycle, and the relevance of these factors to earth\u2019s climate system. In a follow up post<\/a>, I went over the structure of the sun, and some of the characteristics of each layer, which laid the groundwork for my last post<\/a>, in which I explained the solar dynamo: the physical mechanism underlying solar cycles.<\/p>\n

Before elaborating on the sun\u2019s role in climate change in the installment following this one, I\u2019ll be going over an approach called \u201cMean Field Theory\u201d in this installment, which dynamo theorists and other scientists sometimes use to make the modelling of certain systems more manageable. As was the case with part III, this may be a bit more technical than most of my subscribers are accustomed to, but I think the small subset of readers with the tools to digest it will appreciate it. And to be perfectly blunt, writing this was not just about my subscribers. I wanted to do it. It was an excuse for me to dig more deeply into something that has been going on in modern stellar astrophysics that I thought was interesting. The fact that it happened to be tangentially related to my series on climate science was a mere convenience. Anyone wanting to avoid the math and\/or to cut to the chase with respect to the effects of solar cycles on climate change might want to skip ahead to part V, or perhaps just read only the text portions of this post. However, for those who don’t mind a little bit of math, I present to you the following:<\/p>\n

Mean Field Theory<\/strong><\/p>\n

One approach by which scientists and mathematicians can simplify the models describing large complex systems stochastic effects is called Mean Field Theory (Schrinner 2005<\/a>). This involves subsuming multiple complicated interactions between different parts of a system into a single averaged effect. In this way, multi-body problems, (which are notoriously difficult to solve even with numerical approximation methods with supercomputers), can be reduced to simpler single body problems. For instance, the velocity field u<\/strong> and magnetic field B<\/strong> could each be broken up into two separate terms: A mean term (u0<\/sub> <\/strong>and B0<\/sub><\/strong> respectively), and a fluctuating term (u\u2019<\/strong> and B\u2019 <\/strong>respectively), whereby the mean terms are taken as averages over time and\/or space depending on what is appropriate to the system being modeled.<\/p>\n

In other words, u = u0 <\/sub>+ u\u2019 <\/strong>and B = B0<\/sub> + B\u2019<\/strong>, where (by definition) the average velocity <u> = u0<\/sub><\/strong>, the average magnetic field <B> = B0<\/sub><\/strong>, because u0<\/sub> <\/strong>and B0<\/sub><\/strong> are the mean field terms, and the average of the fluctuating terms <B\u2019><\/strong> and <u\u2019><\/strong> are both equal to zero by definition. The angled brackets simply denote a suitable average of the term they enclose (again taken over time and\/or space as deemed appropriate by the scientist or mathematician). The reason the fluctuating terms average out to zero is because the mean field terms are defined as the average of the entire field, so by definition, the only way that can be true is if the fluctuating terms average out to zero.<\/p>\n

Using the vector calculus identity<\/a> \u2207<\/strong> \u00d7 (<\/strong>\u2207<\/strong> \u00d7 B) = <\/strong>\u2207(\u2207<\/strong>\u00b7<\/strong>B)<\/strong> + <\/strong>\u2207<\/strong>2<\/sup><\/strong>B<\/strong>, and the fact that \u2207<\/strong>\u00b7<\/strong>B = 0<\/strong> by Gauss\u2019s Law<\/a> for magnetic fields, the induction equation, \u2202B\/\u2202t\u00a0 = <\/strong>\u2207<\/strong> \u00d7 (u \u00d7 B \u2212 \u03b7<\/strong>\u2207<\/strong> \u00d7 B)<\/strong> from the previous section can also be expressed as \u2202B\/\u2202t = \u03b7<\/strong>\u2207<\/strong>2<\/sup><\/strong>B + <\/strong>\u2207<\/strong> \u00d7 (u \u00d7 B)<\/strong>, where \u2207<\/strong>2<\/sup><\/strong>B<\/strong> is called the Laplacian<\/a> operator of the magnetic field B<\/strong>.<\/p>\n

Plugging in our mean field equations u<\/strong> and B<\/strong> into this form of the induction equation:<\/p>\n

\u2202B0<\/sub>\/\u2202t + \u2202B\u2019\/\u2202t = \u03b7<\/strong>\u2207<\/strong>2<\/sup><\/strong>B0<\/sub> + \u03b7<\/strong>\u2207<\/strong>2<\/sup><\/strong>B\u2019+<\/strong><\/p>\n

\u2207<\/strong> \u00d7 (u0<\/sub>\u00d7 B0<\/sub>) + <\/strong>\u2207<\/strong> \u00d7 (u0<\/sub> \u00d7 B\u2019) + <\/strong><\/p>\n

\u2207<\/strong> \u00d7 (u\u2019 \u00d7 B0<\/sub>) + <\/strong>\u2207<\/strong> \u00d7 (u\u2019 \u00d7 B\u2019)<\/strong><\/p>\n

Now we take the average of both sides:<\/p>\n

\u2202<B0<\/sub>>\/\u2202t + \u2202<B\u2019>\/\u2202t = \u03b7<\/strong>\u2207<\/strong>2<\/sup><\/strong><B0<\/sub>> + \u03b7<\/strong>\u2207<\/strong>2 <\/sup><\/strong><B\u2019> +\u00a0<\/strong><\/p>\n

\u2207<\/strong> \u00d7 <u0<\/sub>\u00d7 B0<\/sub>> + <\/strong>\u2207<\/strong> \u00d7 <u0<\/sub> \u00d7 B\u2019> + <\/strong><\/p>\n

\u2207<\/strong> \u00d7 <u\u2019 \u00d7 B0<\/sub>> + <\/strong>\u2207<\/strong> \u00d7 <u\u2019 \u00d7 B\u2019><\/strong><\/p>\n

However, we\u2019ve already know that <B0<\/sub>> = B0<\/sub><\/strong>, <u0<\/sub>> = u0<\/sub><\/strong>, and <u\u2019> = < B\u2019> = 0<\/strong>, so this can be simplified:<\/p>\n

\u2202B0<\/sub>\/\u2202t = \u03b7<\/strong>\u2207<\/strong>2<\/sup><\/strong>B0<\/sub> +<\/strong> \u2207<\/strong> \u00d7 <u0<\/sub>\u00d7 B0<\/sub>> + <\/strong>\u2207<\/strong> \u00d7 <u\u2019 \u00d7 B\u2019><\/strong>.<\/p>\n

The \u2207<\/strong> \u00d7 <u\u2019 \u00d7 B\u2019> <\/strong>is typically then replaced with a term \u2207<\/strong> \u00d7 \u03b5<\/strong>, where \u03b5 <\/strong>is called the mean electromotive force (Radler 2007<\/a>). This yields the Mean Field Induction Equation:<\/p>\n

\u2202B0<\/sub>\/\u2202t = \u03b7<\/strong>\u2207<\/strong>2<\/sup><\/strong>B0<\/sub> +<\/strong> \u2207<\/strong> \u00d7 <u0<\/sub>\u00d7 B0 <\/sub>+ \u03b5><\/strong><\/p>\n

Although many details of the theory are still being worked out, models based on the solar dynamo mechanism are consistent with the periodicity of the solar cycle, Hale\u2019s law (the opposing magnetic polarity of sunspots above and below the solar equator and the alternation of polarity in successive 11 year cycles), as well as both Sporer\u2019s and Joy\u2019s laws (the apparent migration of sunspots towards the equator as a cycle progresses, as well as their tilt), which together produce the observed sunspot butterfly diagrams I talked about here<\/a>.<\/p>\n

Some models can even simulate variations in amplitude from one cycle to the next, but the precise manner in which Grand Solar Maxima and Minima emerge is still being worked out. Consequently, our models\u2019 ability to reliably and accurately forecast them is currently still limited. Methods have been developed for estimating sunspot number and solar activity of a cycle\u2019s solar maxima from observations of the poloidal field strength of the preceding cycle\u2019s solar minima (Schatten 1978<\/a>, Svalgaard 2005<\/a>). In addition to only providing information on the solar maxima immediately after the minima being measured, this approach is also limited by the fact that our poloidal field measurements only go back a few cycles, and because the poloidal magnetic fields during solar minima are difficult to measure reliably, because they are weak and have radial as well as meridional components.<\/p>\n

Other researchers have focused on kinematic flux transport solar dynamo models which, in addition to differential rotation, include the effects of meridional flow in the convective envelope, whereby the poloidal magnetic field is regenerated by the decay of the bipolar magnetically active regions subsequent to their emergence at the solar surface (Dikpati 1999<\/a>, Dikpati 2006<\/a>, Choudhuri 2007<\/a>). Active regions are the high magnetic flux regions at which sunspots emerge.<\/em><\/p>\n

\"\"<\/a>

Image <\/a>\u00a0by Andres Munoz-Jaramillo (check out his fantastic presentations on solardynamo.org<\/a>). <\/span><\/p><\/div>\n

This meridional flow sets the period of the cycle, the strength of the poloidal field, and the amplitude of the solar maximum of the subsequent cycle. However, estimates of meridional flow velocities prior to 1996 are highly uncertain (Hathaway 1996<\/a>). All of these models have been criticized by peers of their proponents. A concise summary of the blow by blow can be viewed here<\/a>.<\/p>\n

As for Grand Solar Maxima and Minima, no comprehensive theory has yet emerged on how they arise and decay, let alone a scientific consensus. However, certain constraints have been identified. There is evidence that the dynamo cycle does continue in some modified form during Maunder-type minima periods. The idea is that the dynamo enters Grand Maxima and Minima by way of chaotic and\/or stochastic processes. In the case of Grand Maxima, the dynamo also exits that state via stochastic processes. In the case of Grand Minima, on the other hand, the dynamo then gets \u201ctrapped\u201d in this state, but eventually gets out of it via deterministic internal processes (Usoskin 2007<\/a>). It is also thought that the polarity of the sun\u2019s toroidal magnetic field may lose its equatorial anti-symmetry during such minima, and instead become symmetric (Beer, Tobias and Weiss 1998<\/a>).<\/p>\n

Truly fantastic long term predictive power for solar cycles probably won\u2019t be achieved until poloidal magnetic field generation is better understood, which will likely include improvements in flux transport models, and a more complete characterization of the statistical properties of bipolar magnetic regions (BMRs). For a comprehensive overview of the current state of Solar Dynamo Models and their predictive strengths and limitations, see Charbonneau 2010<\/a>.<\/p>\n

In the next installment, I\u2019ll explain how all of this relates to climate change on earth, and address the elephant in the room: “are solar variations responsible for the current global warming trend?”<\/p>\n

Related Articles:<\/strong><\/p>\n