GD
Geodynamics

Geodynamics

The world’s largest magnet

The world’s largest magnet

The Geodynamics 101 series serves to show the diversity of topics and methods in the geodynamics community in an understandable manner for every geodynamicist. PhD’s, postdocs, full professors, and everyone in between can introduce their field of expertise in a lighthearted, entertaining manner and touch on some of the outstanding questions and problems related to their method of choice.
This week Maurits Metman, PhD student at the Deep Earth Research group at the University of Leeds in the United Kingdom, explains the dynamics of the core. Do you want to talk about your research area? Contact us!

Rock bottom

Approximately 3,000 km below our relatively minuscule feet lies the Earth’s core. It is our planet’s innermost and therefore most secluded region. It is also the primary source of Earth’s magnetic field that we observe here at the surface. With its dynamics, composition, magnetic field generation, and thermal history not yet completely understood, the core remains amongst the most enigmatic parts of the Earth. It has been established that the core can be partitioned in an inner and outer region, which have distinct physical and chemical properties. For example, the two regions are in a different state of matter: the inner core being solid and the outer core liquid. Therefore, it is the outer core that is of particular geodynamical interest – here we will touch upon some important aspects of the dynamics that take place within the outer core.

The outer core consists of an electrically conducting iron alloy liquid, which circulates throughout the outer core volume. In terms of the forces that drive these motions, there are similarities to the dynamics of other terrestrial systems such as the mantle, oceans, and atmosphere. For example, in all cases gravitational forces are overcome by the process of convection, through which relatively hot and buoyant material at the base of the system rises towards the surface, while elsewhere cold material sinks. Additionally, the flows in these systems are subject to forces due to pressure differences and those associated with the deformation of the material.

Figure 1: An impression of convection in the outer core (not to scale), which is aligned along columnar rolls, and flow in- and outside the tangent cylinder is separated (Credit: United States Geological Survey).

Nevertheless, the dynamics in the outer core are certainly different to other geophysical flows. For one, it is estimated that a typical velocity for outer core flow U ∼ 10-1 mm s-1, which is relatively high for the solid Earth. In fact, recent work has shown that these velocities may locally be as high as roughly 1 mm s-1 (Livermore et al., 2016). Additionally, rotational effects (e.g. centrifugal force, Coriolis effect) have a tremendous impact on the style of convection. This is for example not the case in the mantle, due to the fact that flow velocity is comparatively low there and viscosity is large. In this respect, the so-called Taylor-Proudman theorem provides an important constraint on the style of core motions, and states that for rapidly rotating systems flow is two-dimensional: it can not change parallel to the axis of rotation. More generally, the style of convection inside of the outer core is strongly cylindrical, in the sense that flow is aligned in ‘columnar rolls’ aligned with the axis of rotation (Fig. 1).

Magnetic soup

With its ability to generate a magnetic field, the outer core further distinguishes itself from other parts of the Earth. That this field must indeed be generated somewhere inside the Earth was already demonstrated by Gilbert (1600), but the fact that it is linked to core fluid flow remained unknown for centuries. We now know that the convective motion of the electrically conducting outer core liquid generates such a magnetic field. This conversion of kinetic to magnetic energy is a process that has fittingly been coined the geodynamo.

What clues do we have that this field must be generated internally? A relatively simple argument can be made from the age of the magnetic field, which paleomagnetic observations have shown to be over 109 year. However, if we were to assume there would be no field generation, the present-day field would decay through simple diffusion (or equivalently due to Joule heating of the fluid) on a timescale of 105 year, inconsistent with these observations. Therefore, it is required that some field generation in the outer core acts to sustain the magnetic field against diffusion, which can be accomplished with a specific core fluid flow there.

Initially, some rejected the existence of such a flow. One well known so-called anti-dynamo theorem is Cowling’s (1933), who showed that a steady and axisymmetric flow field can never maintain a magnetic field indefinitely, which led to the general consensus that sustained dynamo action through fluid flow would not be possible.

Figure 2: A schematic of magnetic field generation through the α-effect at different timesteps ti. Here, u, B and j represent fluid velocity, magnetic field and electric current density.

The development of the mean-field theory, which describes how small-scale flow perturbations can on average create a large-scale magnetic field, changed this. An example of such field generation is through the α-effect (Parker, 1955). In this case there is a rising and rotating flow (imagine a corkscrew-shaped motion) moving and twisting a magnetic field line (Fig. 2). The magnetic loops created this way induce an electric current parallel to the the field, which in turn generate a secondary magnetic field that is perpendicular to the initial field. A similar conversion the other way around is also possible, and a planetary dynamo that relies on these two processes is considered to be of the α2-type. Another source of field generation that follows from mean field theory is the ω-effect. This process is the bending of magnetic field lines, due to the differences in rotation rate, also creating magnetic field in a direction opposite to the initial direction (Fig. 3). A dynamo that generates a magnetic field through the α- and ω-effect is referred to as an αω-dynamo.

 

Therefore, it is quite clear that core flow and Earth’s internal magnetic field are deeply intertwined. As mentioned earlier, field generation through fluid motion and field diffusion are two competing processes that control time variations in the field. The magnetic Reynolds number is the ratio of the magnitudes of these contributions, and is defined as:

Figure 3: A schematic of the ω-effect which converts the magnetic field from the initial direction (aligned South-North) to a secondary direction (West-East and vice versa), at different timesteps ti. The solid and dashed curve represent the magnetic field and rotation axis, respectively.

Rm = UL / η

where η is the magnetic diffusivity and L is a length scale for the magnetic field (Roberts and Scott, 1965). For the outer core it is estimated that Rm ∼ 102, and therefore the diffusion term is often considered negligible (at least for relatively large length scales). This is referred to as the frozen-flux approximation. As the name suggests, magnetic field lines are then dynamically ‘frozen’ into the liquid, so that they evolve as though they were material line elements. How realistic is this particular scenario? From the above equation it should be clear that the frozen-flux approximation can break down if the typical length scale decreases. This may for example be the case for flux expulsion, i.e. when a radially expelled field is concentrated below the core-mantle boundary (Bloxham, 1986). This concentration increases the gradient of the field locally, enhancing radial diffusion. However, to what extent this process is realistic remains a subject for debate.

Forecast: reversals?

For the last two decades, advances in computing power have allowed numerical models to reproduce certain properties of the Earth’s magnetic field. For example, such models have been shown to exhibit magnetic polarity reversals (Glatzmaier & Roberts, 1995) and Rm that are similar to the outer core’s. Despite this numerical success and despite the fact that reversals have been documented extensively within the field of paleomagnetism, it remains unknown what physical process underlies these phenomena. This is particularly interesting as the most recent reversal occurred around 0.78 Myr years ago, which has led to speculation that a future reversal is imminent. Future numerical work, increases in computing power, better theoretical understanding of the internal dynamics of the core, and more geomagnetic observations may in time provide a physical explanation for these events.

References
Bloxham, J. (1986). The expulsion of magnetic flux from the Earth’s core. Geophysical Journal International, 87(2):669-678.
Cowling T. G. (1933) The magnetic field of sunspots. Monthly Notices of the Royal Astronomical Society 94: 39-48.
Gilbert W. (1600) De Magnete. London: P. Short.
Glatzmaier G. and Roberts P. (1995) A three-dimensional self-consistent computer simulation of a geomagnetic field reversal. Nature 337: 203-209.
Livermore, P. W., Hollerbach, R. and Finlay C. C. (2016). An accelerating high-latitude jet in Earth's core. Nature Geoscience 10: 62-68.
Parker E. N. (1955) Hydrodynamic dynamo models. Astrophysical Journal 122: 293-314.
Roberts, P. H. and Scott, S. (1965). On Analysis of the Secular Variation. 1. A Hydromagnetic Constraint: Theory. Journal of Geomagnetism and Geoelectricity, 17(2):137-151.

A Geodynamicist and an Early Career Scientist

A Geodynamicist and an Early Career Scientist

This week Adina Pusok, postdoctoral fellow at the Institute of Geophysics and Planetary Physics (IGPP) at Scripps Institution of Oceanography, University of California San Diego, USA, discusses what it is like to be an Early Career Scientist within the EGU Geodynamics division.

The terms “Early Career Scientist” (ECS) or “Young Scientist” (YS) are now so widely used in the scientific community, that certain meetings, sessions, awards, and social events are entirely dedicated to promote this group of scientists. But what are ECS and YS? These terms were created to describe scientists in the early stages of their career. Because the term Young Scientist might have an age connotation and make some scientists feel excluded, the EGU has recently adopted the ECS term instead. The definition of ECS on the EGU page says “an Early Career Scientists (ECS) is an undergraduate or postgraduate (Masters/PhD) student or a scientist who has received his or her highest degree (BSc, MSc, or PhD) within the past seven years*. (*with additional year(s) of parental leave time per child, where appropriate)”.

Yet, it is still not clear why there is a need to create a subgroup and organize specific activities for this group, when the research community is considered open, where everyone can bring contributions in an equal manner? To answer this, I look back at the first meetings I attended as a PhD student. I remember I felt intimidated by the experience and scientific eloquence of the more established scientists (especially the “big names” in the field), only because I knew there was still a lot for me to learn and experience to gain. And I am sure I was not the only PhD student experiencing this. The truth is, ECS have different needs compared to established scientists. The list of challenges for ECS is long (see this special issue in Nature). Besides doing research, ECS need to promote their work, socialize more to expand their network, enter the very competitive job market, and not to mention present and defend their scientific achievements to an unknown audience. For many established scientists, these are stages past gone. Everyone knows them, everyone talks about their work.

As an ECS, I value the interactions with established scientists, but I also welcome events organized for ECS (i.e. career development or social events, such as grant writing and academic presentation courses). However, I feel the best outcome is achieved when both established and early career scientists support and participate in these events. Which is why I agreed to be the ECS representative for the EGU Geodynamics division (GD) – to bring together and improve visibility of ECS in the GD community. And for the occasional beer and karaoke special sessions (see last week’s blog post) with other ECS Geodynamicists.

What does EGU do to support ECS?

A large proportion of the Union’s members consists of scientists in the early stages of their career. For this reason, EGU wants to offer support to this group, by providing reduced conference fees, recognizing outstanding students, awarding travel grants, organizing short courses, arranging networking possibilities and more. Besides, the EGU encourages further support at the division level, such as outreach activities (e.g. blogs and social media), social events, mentoring programs or short courses at the General Assembly. All ECS need to do is pay attention to the opportunities provided.

Probably a lot more could be done, but small contributions can already make a huge difference. For example, the ability to attend a conference due to an awarded travel grant can be very important to meet other scientists and create exposure for your research, create future collaborations or even sign up for that future job.

What can ECS do for themselves?

Other than help create the environment they want to be in? It might sound idealistic or as too much work involved, but again, small steps can go a long way. Anything from organizing to participating, from being informed to inform others, from taking part in outreach to supporting outreach activities, from mentoring to being mentored, the ECS can contribute in various ways to their own and their community’s development.

In Geodynamics, there are many enthusiastic and fun ECS that get together at meetings and workshops, create friendships and collaborations across national and academic borders. I’ve met enough of them to believe that GD ECS can create a more coherent structure under the EGU umbrella. For example, within less than 1 year, the GD ECS have managed to launch the EGU GD Facebook page, organize social events at major conferences (AGU Fall Meeting 2016, EGU General Assembly 2017), and recently, launch the EGU Geodynamics Blog.

Adapted from www.xkcd.com and www.soest.hawaii.edu

More things can be done though. There is currently a strong need for scientists to become more actively involved in science public outreach worldwide, and this is one of the directions where GD ECS can contribute easily. Geodynamics is a field known for “beautiful pictures” that show the complexity and dynamics of natural processes. Thus, highlighting the diversity of Geodynamics studies (numerical simulations, laboratory experiments, or data compilations) to a wider audience can benefit the research community at large.

Because all these initiatives are run by members, there is always a need for motivated people with refreshing ideas. This is why, I encourage other ECS to bring forward new ideas that we can develop within the GD division (see how to get involved). Along with the rest of GD ECS volunteers, I look forward to working with you!

By Adina Pusok
As the ECS GD representative, I am the link between the EGU and the ECS GD community. I provide EGU with feedback from students and early career researchers, so that the union can take action to improve the ECS activities at the EGU General Assembly and maintain the support for early career scientists throughout the year. I am also involved in public and community outreach of the Geodynamics division.

Karaoke, geodynamics, and a bit of history

Karaoke, geodynamics, and a bit of history

Let me just talk to you about what I have been doing with my free time recently: I discovered a feature from Google Books named Ngram viewer which allows you to make graphs that show how words or phrases have occurred in a selection of books (e.g., English) over the selected years. I have of course been playing with this thing! You can imagine how exciting my weekends are. In all seriousness, though: I highly recommend you check it out sometime, as you can do some pretty fun things with it.

First of all, by looking at the number of mentions in books, you can quickly determine when a certain research field was first established as illustrated in the graph below. (Structural) Geology has always been an older branch of Earth sciences, as it originated in field observations rather than instruments or computers. James Hutton (1726-1797) was the Father of Modern Geology, as he developed the theory of uniformitarianism (= the processes on Earth that we see today also occurred in the past). His work was popularised in the 1830s by Charles Lyell, who also coined the phrase ‘The present is the key to the past’. Indeed, mentions of structural geology start to occur somewhere in the 1850s when Hutton’s and Lyell’s ideas have been firmly established in the scientific community.
According to our graph seismology originated in the 1850s. Indeed, a quick google search will tell you that the word ‘seismology’ was coined in 1857 by Robert Mallet, who also laid the foundation of instrumental seismology.
Both geodynamics and tectonophysics (fields that only really originated after the general acceptance of plate tectonics in the 1960s and benefited greatly from advances in computer science) are only starting to flourish in the late 1960s. Note that all graphs show a decline around the 1990s…

Occurrence of seismology, tectonophysics, structural geology, and geodynamics in English books from 1800-2008 with smoothing factor 4.

Apart from discovering when a certain field was established, it is also possible to see the direct effect of certain global events on the publishing history of a particular field. One of the most convincing cases stems from the tsunami research area. After the devastating 2004 Sumatra Boxing day earthquake and tsunami, interest in tsunamis surged in 2005 and has since remained much more popular than before (at least up until 2008).

Occurrence of tsunami in English books from 1800-2008 without a smoothing factor. There is an increase in publishing after the 2004 Sumatra earthquake and tsunami.

Now on to the really fun part: karaoke. A phenomemon in the geodynamics community that is hard to get around. We have all been working and networking very seriously at conferences before we suddenly got swept away to the nearest karaoke bar. If you look at the graph comparing the amount of mentions of seismology, geodynamics etc., you will notice that there is a decline in mentions starting roughly in the 1990s. There could of course be many reasons for this: maybe the Google database is, as of yet, still incomplete, so the ongoing upward trend cannot be seen; maybe too much has been published, so that the relative percentages of words mentioned has declined, even though the absolute amount of published works containing these works has increased; but maybe.. just maybe, the reason is karaoke, as it started to spread around the world in the 1990s.

Now, I am not saying I did any fancy statistics on these results. I am also not saying that there is any causality involved (because we all know how hard it is to determine causality). I am just saying: look at the graph and draw your own conclusions…

Occurrence of seismology, tectonophysics, structural geology, geodynamics, and karaoke in English books from 1800-2008 with smoothing factor 4.

Don’t be a hero – unless you have to

Don’t be a hero – unless you have to

The Geodynamics 101 series serves to show the diversity of topics and methods in the geodynamics community in an understandable manner for every geodynamicist. PhD’s, postdocs, full professors, and everyone in between can introduce their field of expertise in a lighthearted, entertaining manner and touch on some of the outstanding questions and problems related to their method of choice.
This week Dr. Cedric Thieulot, assistant professor at the Mantle dynamics & theoretical geophysics group at Utrecht University in The Netherlands, discusses the advantages and disadvantages of writing your own numerical code. Do you want to talk about your research area? Contact us!

 
In December 2013, I was invited to give a talk about the ASPECT code [1] at the American Geological Union conference in San Francisco. Right after my talk, Prof. Louis Moresi took the stage and gave a talk entitled: Underworld: What we set out to do, How far did we get, What did we Learn?

The abstract went as follows:
 

Underworld was conceived as a tool for modelling 3D lithospheric deformation coupled with the underlying / surrounding mantle flow. The challenges involved were to find a method capable of representing the complicated, non-linear, history dependent rheology of the near surface as well as being able to model mantle convection, and, simultaneously, to be able to solve the numerical system efficiently. […] The elegance of the method is that it can be completely described in a couple of sentences. However, there are some limitations: it is not obvious how to retain this elegance for unstructured or adaptive meshes, arbitrary element types are not sufficiently well integrated by the simple quadrature approach, and swarms of particles representing volumes are usually an inefficient representation of surfaces.

Aside from the standard numerical modelling jargon, Louis used a term during his talk which I thought at the time had a nice ring to it: hero codes. In short, I believe he meant the codes written essentially by one or two people who at some point in time spent great effort into writing a code (usually choosing a range of applications, a geometry, a number of dimensions, a particular numerical method to solve the relevant PDEs(1), and a tracking method for the various fields of interest).

In the long list of Hero codes, one could cite (in alphabetical order) CITCOM [1], DOUAR [8], FANTOM [2], IELVIS [5], LaMEM [3], pTatin [4], SLIM3D [10], SOPALE [7], StaggYY [6], SULEC [11], Underworld [9], and I apologise to all other heroes out there whom I may have overlooked. And who does not want to be a hero? The Spiderman of geodynamics, the Superwoman of modelling?

Louis’ talk echoed my thoughts on two key choices we (computational geodynamicists) are facing: Hero or not, and if yes, what type?
 

Hero or not?

Speaking from experience, it is an intense source of satisfaction when peer-reviewed published results are obtained with the very code one has painstakingly put together over months, if not years. But is it worth it?

On the one hand, writing one own’s code is a source of deep learning, a way to ensure that one understands the tool and knows its limitations, and a way to ensure that the code has the appropriate combination of features which are necessary to answer the research question at hand. On the other hand, it is akin to a journey; a rather long term commitment; a sometimes frustrating endeavour, with no guarantee of success. Let us not deny it – many a student has started with one code only to switch to plan B sooner or later. Ultimately, this yiels a satisfactory tool with often little to no perennial survival over the 5 year mark, a scarce if at all existent documentation, and almost always not compliant with the growing trend of long term repeatability. Furthermore, the resulting code will probably bear the marks of its not-all-knowing creator in its DNA and is likely not to be optimal nor efficient by modern computational standards.

This brings me to the second choice: elegance & modularity or taylored code & raw performance? Should one develop a code in a very broad framework using as much external libraries as possible or is there still space for true heroism?

It is my opinion that the answer to this question is: both. The current form of heroism no more lies in writing one’s own FEM(2)/FDM(3) packages, meshers, or solvers from scratch, but in cleverly taking advantage of state-of-the-art packages such as for example p4est [15] for Adaptive Mesh Refinement, PetSc [13] or Trilinos [14] for solvers, Saint Germain [17] for particle tracking, deal.ii [12] or Fenics [16] for FEM, and sharing their codes through platforms such as svn, bitbucket or github.

In reality, the many different ways of approaching the development or usage of a (new) code is linked to the diversity of individual projects, but ultimately anyone who dares to touch a code (let alone write one) is a hero in his/her own right: although (super-)heroes can be awesome on their own, they often complete each other, team up and join forces for maximum efficiency. Let us all be heroes, then, and join efforts to serve Science to the best of our abilities.

Abbreviations
(1) PDE: Partial Differential Equation (2) FEM: Finite Element Method (3) FDM: Finite Difference Method

References
[1] Zhong et al., JGR 105, 2000; [2] Thieulot, PEPI 188, 2011; [3] Kaus et al., NIC Symposium proceedings, 2016; [4] May et al, CMAME 290, 2015 [5] Gerya and Yuen, PEPI 163, 2007 [6] Tackley, PEPI 171, 2008 [7] Fullsack, GJI 120, 1995 [8] Braun et al., PEPI 171, 2008 [9] http://www.underworldcode.org/ [10] Popov and Sobolev, PEPI 171, 2008 [11] http://www.geodynamics.no/buiter/sulec.html [12] Bangerth et al., J. Numer. Math., 2016; http://www.dealii.org/ [13] http://www.mcs.anl.gov/petsc/ [14] https://trilinos.org/ [15] Burstedde et al., SIAM journal on Scientific Computing, 2011; http://www.p4est.org/ [16] https://fenicsproject.org/ [17] Quenette et al., Proceedings 19th IEEE, 2007