Magnetopause Lab

Jim Crumley and Ari Palczewski

It is contrary to reasoning to say that there is a vacuum or space in which there is absolutely nothing.

1 Introduction

Though outer space is often described as a vacuum, in actuality there are particles everywhere in the solar system, and there is continuous activity in the region between the celestial bodies. The study of this part of the solar system is called Space Physics [Kivelson et al.1995], or sometimes Plasma Space Physics after the particles that fill this region. A plasma is a ionized gas, where some of the electrons have been stripped off of atoms creating a mixtures of free electrons and positive nuclei (ions).

Most of the plasma particles filling the solar system can trace their origin to the solar wind, which comes from the outer region of the sun’s atmosphere continually being blown off of the sun. The solar wind typically has a speed of about 400 km/s, though it can range up to about 1200 km/s.

1.1 The Magnetosphere

When the solar wind reaches the earth, it interacts with the earth’s magnetic field to form a region known as the magnetosphere (See Figure 1). The earth’s magnetosphere is the region of space where the earth’s magnetic field dominates behavior of the plasma, just as the sun’s magnetic field does in the solar wind. The magnetosphere has many regions, though only a few will be of interest to us in this lab.


Figure 1: This figure shows the structure of the magnetosphere. The diagram is from the Southwest Research Institute [2005].

The first region of interest (though it is technically part of the magnetosphere) is known as the bow shock. The bow shock is caused by the solar wind meeting an obstruction its flow in the form of the earth. The solar wind has to flow around the the earth, and in doing a shock formed. The bow shock gets its name from the fact that the physics behind the creation of this shock are similar to the shock formed by the bow of a boat as the boat moves through the water.

Nearer to the earth is the magnetopause. The magnetopause marks the furthest that the earth’s magnetic field extends, which makes the magnetopause the boundary of the magnetosphere. The magnetosheath is the region of space between the bow shock and the magnetopause.

1.2 Solar Activity

The surface of the sun, and consequently the solar wind, experiences a variety of events, including flares and coronal mass ejections (CMEs). These events can change both the speed and the number density (particles per volume) of the plasma in the solar wind. When the solar wind from one of these events reaches the earth, the events are called Solar Storms or Magnetic Storms. These variations in the solar wind cause the location of the both the bow shock and the magnetopause vary. Studying the variation in the distance from the earth to the magnetopause is one of the aims of this lab.

The location of the magnetopause is important because the magnetosphere shields the earth and most satellites that go around the earth from the direct effects of the solar wind. When the magnetopause moves, some satellite that are typically protected by the magnetosphere can find themselves in the solar wind. In particular, many satellites are placed in geosynchronous orbits so that they are always above a fixed location on the earth. Geosynchronous orbits are roughly 42000 km ( or 6.6 RE, where 1 RE equals the radius of the earth) from the center of the earth. During normal solar wind conditions the standoff distance of the magnetopause (the distance from the center of the earth to the point of the magnetopause that is on the line between the magnetopause and the earth) is roughly 10 RE . During extremely strong Solar Storms the magnetopause can move in closer than geosynchronous orbit. This exposes many satellites to direct contact with extreme solar wind conditions, and satellites are often damaged or lost when this happens.

2 Theory

2.1 Magnetopause Location

The location of the magnetopause along the line between the earth and the sun can be derived by balancing the magnetic pressure of the magnetosphere with the dynamic pressure of the solar wind:

r(R  ) = 107.4(n v2 )-16
o  E          sw sw

where ro is the standoff distance from the earth to the sun in RE, nsw is the number density of the plasma in the solar wind in cm-3, and vsw is the speed of the solar wind in km/s [Kivelson et al.1995, pages 173-4]. The entire magnetopause boundary has also been fit to various empirical functions by a variety of researchers. The currently favored model for the earth’s magnetopause position can be predicted during extreme events come from Shue et al. [1997,  1998]. These were both derived from specific Solar Storm events. the expressions are

                               (       )a
                          r = ro---3----                 (2)
                                1+ cosh 1
ro = (10.22+ 1.29 tanh [0.184(Bz + 8.14)])(Dp) 6.6              (3)
         a = (0.58- 0.007Bz)[1 + 0.024ln(Dp)]               (4)
where r is the distance from the earth to the magnetopause boundary, ro is the magnetosphere standoff distance, Bz is the z-component of the solar wind’s magnetic field, Dp is the dynamic pressure of the solar wind, and a is the amount of tail flaring from the night side of the magnetopause.

2.2 Coordinate Systems

A variety of different coordinate systems are used to measure locations in space, and in this lab you will have to use data in several of these coordinate systems. More information on these coordinate systems can be found in books [Kivelson et al.1995] or online [ESA2005].

2.2.1 Geocentric Solar Ecliptic

The geocentric solar ecliptic (GSE) coordinate system has its x-axis point from the earth to the sun and its y-axis is in the ecliptic plane (the plane of the earth’s orbit) pointing in the direction of dusk (opposite the direction of the earth’s orbital motion). Its z-axis is in the direction of the ecliptic pole.

2.2.2 Geocentric Solar Magnetospheric

In geocentric solar magnetospheric (GSM) coordinates, as in GSE coordinates, the x-axis points from the earth to the sun. Its z-axis points in the direction of the north magnetic pole. Its y-axis, is perpendicular to the other two, so that a right-handed coordinate system results. The difference between the this system and the GSE system is a simple rotation about the x-axis.

2.2.3 Geocentric Equatorial Inertial

In geocentric equatorial inertial (GEI) coordinates (also known as geocentric solar inertial (GCI)), the x-axis points from the earth to the the first point of Aries (the location of the sun at vernal equinox). This line is where the equatorial and ecliptic planes intersect. The z-axis is parallel to the earth’s rotation axis, and the y-axis complete the the right-handed coordinate system.

3 Simulation

In this portion of the lab you will use a computer simulation to model the location of the magnetopause. The model we will use is called BATS-R-US [Hansen et al.2002], the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme, which was developed at the University of Michigan. We will be running the simulations on super-computers at NASA-Goddard run by the Community Coordinate Modeling Center (CCMC) [NASA2005].

3.1 BATS-R-US Input File

To use BATS-R-US, you must first create an input file which tells the simulation the type of conditions that you would like to simulate. The input file is a simple text file with the conditions at different times listed. A sample input file is available (batrus_input_template.txt) and a portion of one is included below.

YearMonDayHrMnsmsBx[nT]ByBzV x[km/s]V yV zn [cm-3] T[K]

2000 1 1 0 00 0 5 0 0 -500 0 0 10200000
2000 1 1 0 50 0 5 0 0 -500 0 0 10200000
2000 1 1 0 100 0 5 0 0 -1000 0 0 10200000
2000 1 1 0 150 0 5 0 0 -1000 0 0 10200000

First, note that the first line is a header line which should not be included when you try to run the simulation. Next, note that the first 7 columns of the input file set the time for the conditions that follow. Since we will not be simulating real conditions from a real day, the date we pick is arbitrary, though the usual convention is to pick 1/1/2000 for the date. The time differences are significant, though for our work here we will just use a 5 minute time difference between each of the points that we use. The next six columns give the three components of the solar wind magnetic field in nT, followed by the three components of the solar wind speed in km/s. The final two columns are the solar wind plasma number density in cm-3 and the solar wind plasma temperature in Kelvin.

In this lab you will complete one simulation run where you vary, in turn, two parameter. In one part of the run you will vary V x and in the other part you will vary n. For each run you will compare the location of the subsolar point of the magnetopause to the location predicted by equation 1. Ask your instructor what specific range you should use for each. Note that because of the way the simulation interpolates the parameters between input lines that you will always want to have the same parameters between consecutive input lines. Before submitting your simulation, have your instructor check your input parameters.

3.2 Submitting Your BAT-R-US Simulation Run

As stated explained above, we will be running the BAT-R-US simulation using supercomputers at CCMC. This section describes how to submit your simulation run.

  1. Go to
  2. Click on “Space Weather Models” on left.
  3. Under “Global Magnetosphere Models” click on the “BATS-R-US” link.
  4. Click the “Submit a Run” link.
  5. Under Runs on Request Submission click the “Submit Request” link.
  6. Fill in all required areas and hit the “Submit” button.
  7. Under Step 2, choose “BATSRUS” for the simulation.
  8. Under Step 3, choose “Simulation with Modeled Conditions” and “Time-Dependent.”
  9. Set the time at the end of the run equal to the last time in your input file.
  10. Leave the dipole tilt at 0 for simplified models.
  11. Choose “No Corotation Velocity” and “Semi-Empirical Auroral” and leave the radio flux at 150 since we will not be dealing with the ionosphere, these settings won’t matter.
  12. Load the input file. Making sure that you input data looks like the data above with no units at the top. After you have found your file, hit the “Press” button.
  13. Check over the input data to make sure its right and then hit the “Continue” button.
  14. Do not set anything under “Step 6”, just press the “Continue” button. Hit “Continue” on the next screen as well.
  15. On the “Submit your run” screen, fill in a description of your run in the first text box. In the second text box put in a key word like “test1”. Finally, if you are sure everything is correct, hit the “SUBMIT REQUEST” button.

Once you have submitted the job, you will be sent an email confirmation. After your run is complete, another email will be sent to you.

3.3 Viewing Your Results

After your simulation run has been completed, you can use the online tools to view the results of your run. To get to these tools, hit the “View Simulation Results” link from the CCMC BAT-R-US page. Find your run and click the link. Then click the link to “View Magnetosphere.” This will bring you to a complicated web form which you can use to view your results. The simulation results can be viewed in many different ways, though we will only use a few of them.

3.3.1 Finding Subsolar Magnetopause

First you want want to choose which time point of your simulation to look at under the “Choose data time” option. If you used the data input scheme suggested above where the input parameters are constant for times in a row, then you will want to look at the second time of each pair. In order to find the subsolar point of the magnetopause, it easiest to look at a line plot of various plasma parameters, so choose “Plot Mode” of “Line Plot.” Since the subsolar point in on the x-axis, choose “Y1”,“Y2”, “Z1”, “Z2” all equal to 0 under the “Choose Plot Area” section at the bottom of the screen. For the x-limits you will want to start out with a range that is sure to include your magnetopause, and then zoom in once you have found your subsolar point. In order to determine your magnetopause, you should take a look at several of the following data quantities: V x, B, Bx, By, Bz, and n.

The magnetopause boundary should be visible in your plots as the location where the the plasma parameters are changing rapidly from the solar wind values which you input, to the magnetosphere values. You should be able to determine the same boundary location from looking at more than one plasma parameter. You should assign a value for the uncertainty of the location of the subsolar point based on how precisely you can determine the locations from your plot. Be sure to save the plots that you use for finding the subsolar point to a file. You should include copies of these plots in your lab notebook.

3.4 Making Movie of Results

In order to get a better feel of the global results of your simulation, you should also make a movie of your simulation results. To make this movie you will make plots at a series of times of your simulations results and then combine them to make a simple movie.

To make the movie you will use the same web plotting interface that you used above to make images of your results. You can choose whichever plot type and plasma parameters you would find interesting -- experiment with different choices.

Once you have found a plot type you like, make a plot of that type at each time for your simulation. Save a file copy of each plot (by right clicking on the image).

The program mpeg_creator can be used to create a movie from the image files. This is which asks you a few questions and then calls the program mpeg_encode which actually makes the movie. You can access the program from the class web site or from ~jcrumley/public\_html/332/magnetopause/data/mpeg\_creator. To run the program open up a command shell (xterm), change to the directory your data is stored in, and then run the program by typing
into the command shell. The program will give you instructions as you run it. When you hand in this lab, you should turn in one copy of your movie to your instructor.

3.5 Questions

  1. The distance from the earth to is magnetopause is given roughly by the equation 1. This equation makes some assumption that the solar wind for a typical Solar-Wind momentum flux of 2.6 nPa and with the Earths magnetosphere at about 10RE [Kivelson et al.1995][][page 171]. This is only a general prediction of the location of the magnetopause. From your modeling the solar wind find your empirical constant factor (107.4) for your modeled solar wind. You may use any means necessary to calculate your number from your model. Earth’s dipole field is -31100 nT for BATRUS.
  2. What are some possible reasons why in BATRUS the 107.4 changes? What approximation are being made when dealing with MHD in space?

4 Spacecraft Data

For part of the lab you will take solar wind input data and examine where earth’s magnetopause was located. You will be looking at 3 specific events: October 31 2003, July 15 2000 (Bastille Day storm), and May 4, 1998. For each event the solar wind data was collected from the Advanced Composition Explorer (ACE) satellite from ACE resides at about 210 RE from the earth on the line between the earth and the sun. On board ACE is an array of instruments to measure the solar wind. Solar wind plasma data (SWEPAM) and magnetic field data (MAG) are the two instruments that we will be using data from. The SWEPAM takes the solar wind velocity, ion density, and ion temperature. While MAG take the magnetic field intensity.

Along with ACE data, we will be using data from several other satellites. During the events listed above these satellites spend at least some time inside the magnetosphere. The point of this lab is to check the satellite data for signs of having crossed the magnetopause and compare that to predictions of equation 2.

4.1 Procedure

For each event the earth’s magnetopause location has already been calculated using equation 2 based on ACE data which is also included in the spreadsheet. The position of the magnetopause along with the location of a particular satellite. These results are included in spreadsheet files which have been stored under the data directory on the course web site. Along with the spreadsheet files are data files which you will want to add to the spreadsheet in order to do further analysis. These files contain either magnetic field data or ion number density data.

  1. Magnetic field data usually contains MAG in the file (G0_K0_MAG_65544.txt) the txt is the actual data and the g  if is a graph of the data.
  2. Ion data should be in a MPA (L1_K0_MPA_65778.txt) file from which you want the “low P DENS”‘ info. Or it will be an EP# (G0_K0_EP8_225752.txt) file where a large flux usually means a transfer between barriers. You may also want to use the SPA data for some of the satellites if it is available (L4_K0_SPA_226515.txt)

To interpret the data you must first make it ready for the spreadsheet. After you have the data in the right coordinates and time scale you can then put them into the spreadsheet for that particular event to interpret the data. In the program “gnumeric”, you can add the data to the spreadsheet by opening it as a file of type “Text Import”. Once you have the data imported, you can plot it and compare it to the plot of the position of the magnetopause.

The Geotail from October 2003 has been done for you so you know what you are looking for and how to interpret the data. Some explanations have been saved on the spreadsheet for the event.

For each satellite there is some sort of data to try and confirm the location of the magnetopause. In some cases the satellite never crosses the theoretical magnetopause, but you need to check that the data shows the same thing. When examining the satellite data, pay particular attention to the times at which the theoretical magnetopause location plot shows the satellite crossing the magnetopause. Note any signs of a magnetopause crossing. Based on your examination of the data and the the theoretical magnetopause location, determine if you believe the satellite crossed the magnetopause and if so at what time[s].

Listed below are the events and the data that we have available for them. Your instructor will tell you which events you should look at.

Event May 4, 1998July 15, 2000October 31, 2003

SatellitesGOES 9 GOES 8 Geotail
L4 GOES 10 GOES 10
L7 L1 GOES 12
Polar L4 L0
L9 L1

4.2 Questions

  1. For which satellite and event pairs did you find magnetopause crossings? How strong do you believe the evidence is for a crossing in this case?
  2. In any of your events does it appear that the satellite crossed the magnetopause repeatedly? If so, does it seem as though magnetopause is moving back and forth over the satellite? What could explain this type of behavior?
  3. Which data did you find easier to interpret - ion or magnetic field ? Why?
  4. Which spacecraft data did you find easier to interpret? Why? Do you think that this is a result of the behavior of a specific event? Or of the satellite’s position during that event?

  1. In any of these cases does the data show the magnetopause crossing geosynchronous orbit?


   ESA, Coordinate systems and transformations, 2005,

   Hansen, K., G. Tóth, A. Ridley, and D. DeZeeuw, BATS-R-US User Manual: Code Version 7.5.0, 2002,

   Kivelson, A., M. G. Kivelson, and C. T. Russell, Introduction to Space Physics, Cambridge University Press, London, 1995.

   NASA, 2005,

   Shue, J.-H., J. K. Chao, H. C. Fu, C. T. Russell, P. Song, K. K. Khurana, and H. J. Singer, A new functional form to study the solar wind control of the magnetopause size and shape, J. Geophys. Res., 102, 9497-9512, 1997.

   Shue, J.-H., et al., Magnetopause location under extreme solar wind conditions, J. Geophys. Res., 103, 17,691-17,700, 1998.

   Southwest Research Institute, The magnetosphere, 2005,