Jump to content

Jesterface23

Member
  • Content Count

    30
  • Joined

  • Last visited

  • Days Won

    3

About Jesterface23

  • Rank
    Minor flare

Profile Information

  • Gender
    Male
  • Location
    United States
  • Interests
    Earth and space weather

Recent Profile Visitors

The recent visitors block is disabled and is not being shown to other users.

  1. I am back to update where I am with my research of Coronal Mass Ejections. To start out, I need to go back to the very beginning. Being interested in Earth's severe weather, the Coronal Mass Ejection of June 22nd, 2015 gave a new interest in space weather and CMEs. I was looking up the parameters and tried to learn all of the basics, the solar wind, IMF, K and Kp indices, etc. There was a lack of severe weather at the time and it wasn't too long before the time of Hurricane Patricia that I began my research. I started out with Coronal Holes, but eventually I thought I might as well try to see if I can find a way to predict the time of CME arrivals. My goal for predicting a CME's arrival time was by using one satellite's Coronagraph imagery and the time of launch. Using SOHO's LASCO C3 imagery, SDO's AIA 131 imagery and ACE and Wind solar wind data to try and find some sort of pattern. I made run after run after run of trial and error using different calculations and points in time from the imagery. I did eventually find a pattern and the image below is the final graph where I have points on paper. This begins a chain reaction to several new discoveries. Here are the steps used per CME to put a point on the graph. 1. Get the difference of the CME launch time to arrival time in hours. The launch time that I use is when the CME reaches 2.5Rsun from the surface of the Sun. 2. In the C3 imagery, get the point in time of when the CME first reaches the edge of the FOV and a second point in time when the opposite side of the imagery reaches the edge of the FOV. The difference between these two times in hours is the X-Axis in the graph above. 3. From step 2 get the time between the two imagery points and get the difference from the launch time. Divide the travel time from step 1 by this step's number and this is the Y-Axis on the graph. Briefly moving forward to mid 2018 I was able to discover the exact base formula for what is seen in the image above. The formula is y=cos(pi*x)-(x/sqrt(pi/2)) and with adjustments made for the C3 image Field of View and L1 position it results to the graph below. I will say there was a huge amount of luck because a minimal amount of changes had to be made to the formula for SOHO's FOV and L1 position, compared to changes that would need to be made for the STEREO satellites. So now we have a problem. There are multiple lines and I am now stuck because I was expecting one line. I was even marking out or skipping data points in the first graph thinking I made a mistake somewhere. It can't be due to changes of the distance from the Sun, or the pattern wouldn't exist. The Coronagraph FOV is near constant, so that is ruled out. Randomly shaped CMEs are also ruled out because again, the pattern wouldn't exist. This is where I step away from the project for a while until the CME of September 10th, 2017. Sunspot region 2673 came and went and I am now back to my research project with a sort of cleared mind. The September 10th CME is what brought a major breakthrough. In SOHO's C3 imagery and STEREO A's C2 imagery you can see a bulge ahead of a larger trailing shock. I look around at other CMEs and I now know why there are multiple lines on the graph. A CME can have two faces, a leading and secondary shock. To say if the secondary shock was right behind the leading shock, given the distance from the sun to the secondary shock multiplied by sqrt(pi/2) gives the position of the leading shock. The 2D illusions are gone and now have the leading structure of CMEs. The lines on the graphs have a new name, Multipliers. The Multiplier is only determined by the Coronagraph imagery and not the overall CME. A very, very preliminary image from this time of the research is below to sort of show the shape. Now looking back at both graphs at the top of the post, the first doesn't show all of the lines. It wasn't until I ran calculations on STEREO A's July 23rd, 2012 CME to find the lowest Multiplier. The old Multiplier 0, or M0, then became M2. At this current time I only know of 2 other M0 CME Coronagraph views from SOHO and STEREO B. Like how the graph is created, the same process can be done to see what we are looking at through Coronagraph imagery by using the Multipliers. Each Multiplier's type of Coronagraph view and impacts are below. M0 - These are still partially a mystery to me. What I currently know is that they are narrower CMEs and are closer they appear with a leading shock impact. M1 - These have a leading shock impact and are similar to M0 CMEs, but have a buckled outer shock as the secondary shock with a leading shock impact. M2 - A leading shock impact where the two view points are the leading shock. M3 - A leading shock impact where the two view points are a leading and secondary shock. M4 - A leading shock impact where the two view points are the secondary shock. M5 - A secondary shock impact where the two view points are a leading and secondary shock. M6 - A secondary shock impact where the two view points are the secondary shock. I wanted to keep going and my next goal has changed to figure out the CME structure. I know how important it will be. Going into this I know of sqrt(pi/2) and I know where the leading and secondary shocks are. Knowing those things I did what I know best, messing around with the numbers until I found a pattern. Eventually run after run after run I was able to find the CME structure while comparing imagery and solar wind data. The result shows CMEs have 3 layers, a central core, a middle layer, and an outer layer. A direct impact through the central core will result in hitting 5 internal shocks. Going further, there is a half shock between each layer. The image below shows almost every type of impact possible and is to scale, but does not account for the correct curve. A question would be how does the structure hold together? My guess is has something to do with the magnetic structure of CME's, but at this point it's something I can't answer. Another question would be how do the different shapes come together? My current answer would be that it is determined by the size of the central core and how explosive the CME is. If there is too much push the middle and outer layers will begin to buckle, but the chances are there are a few other factors as well. As for the internal shocks, I call them shocks for a reason. It varies, but CMEs can be shock driven at any shock. It can be almost similar to the arrival of the CME. As far as I currently know, CMEs can be shock driven all the way back to the last internal shock at minimum. As this is most likely new, I gave simple names to the constants and this is basically what I have used on a calculator. This is how the position of how each shock can be calculated. A = sqrt(pi/2) - This is used for the position of a shock B = Second Root of A - This is used for the position of a half shock. Going past beyond the second root is uncertain and will require more research to be able to see if things become more detailed. Examples of how I will refer to this is that the leading and secondary face of a CME is A0. The duration of the far shock is A0-A2. An impact through the edge of the central core would have impacted A0, A1, A2, A4, A5, A6, and not A3. A middle layer impact has half shock impacts at A0B1, A1B1-A4B1, A5B1. Below are examples of the different types of impacts. Inside each graph I have lines marking each shock. Red lines are the first half shock while blue is the second half. A black line ends the CME, unless disrupted by another CME or CH. In the top left corner of each image is Coronagraph imagery. In the top right corner of each image is my best estimated line of impact through the CMEs and they do not account for the correct curve of the CME. There is a graph for the CME's IMF, Velocity, and Proton Density. One of the best CMEs that I know of is a M0 CME that impacted STEREO B on November 8th, 2013. In the image below, the final internal shock is shown as being shock driven with the velocities at A5. All 4 previous internal shocks are seen in the IMF. The line of impact was through the far edge of the A3 shock, where the magnetic field is disrupted, making the central core visible through the IMF data. An example of an impact through the A3 centr al core shock is a M2 CME that Wind on November 24th, 2001. A sign of this type of impact is when there is a close mirrored IMF shift near A2B1.5 to A3B0.5. An example of an impact through the edge of the central core, missing the A3 shock, is the CME that hit STEREO A on July 23rd, 2012. This is a M0 CME. A sign of this type of impact is gradual changes of the IMF through A2-A4. An example of a M2 CME middle layer impact is a CME that impacted ACE on May 23rd, 2002. The IMF remains near constant from A1B1-A4B1 as it runs through the middle layer. An example of a buckled middle core impact is a M3 CME that impacted STEREO A on September 11th, 2010. This CME is closest to the T-intersections of the central core and middle layers. This may result in a sudden change of the IMF like in the image below. The IMF will also gradually change from A1-A3. An example of a buckled outer core impact is a M3 CME that impacted ACE on April 11th, 2001. Due to the leading shock extending from A0-A2, there may be some issues with the CME being shock driven. This CME below may have had a little help from behind from another CME. It became shock driven at A3B1-A4 at the end of the buckled middle layer. An example of a M5 secondary shock impact CME is from Wind on November 10th, 2000. The CME misses the A2 internal shock and becomes shock driven at A2B1-A3. An example of a M5 far shock impact on DSCOVR is the CME of September 12th, 2017. There isn't much to these impacts. Once the CME passes, there may be a little wall of higher proton density from the solar wind trailing it. This has gone well beyond where I thought this would lead to and what I have previously posted about, but this post is to close this chapter of my research. I still need to reach the goal of reliable CME arrival time forecasting using my method, so at some point in the future another chapter may begin.
  2. Hello, I have some new preliminary information. I've already posted about 2 different leading shocks we can be impacted by. After looking into the internal structure of a CME though there may be possibly be 7 shocks(I might need a different name for it later) total with the final being the exit. The positions of each can be calculated. The cores between each shock are made of two halves and further details are still to be determined. The CME degrades from the end to the leading edge over time. While looking into the internal structure of a CME I needed to know when exactly does a CME launch. When backtracking CMEs currently I ended up coming to when the leading edge of a CME is at around 2.5 solar radii from the sun that seems to be the time of launch. This is still preliminary, but after some time I noticed that is around 1.25 solar diameter. That just happens to be similar to the square root of pi/2, so what are the odds that is where that number comes from. This is very preliminary, I'm not sure of all cases of it, but if one CME runs over another CME it seems to be pushed by the other and becomes very small. If it impacts Earth the second CME would arrive a few hours after the first CME impact.
  3. Well, I though I'd give an update on the current progress of this project. It might be a while until I get to working on this again until a get a pile of work done. I've attempted to contact the correct people at NASA that might take interest in this, and so far there's been no luck with that. I guess this will be a solo project along with the help and ideas from others here until I prove what it can do. I have been working on this a little and I think I have the exact base formula, y=cos(pi*x)-(x/sqrt(pi/2)). I have a modified version for the use of my method of CME arrival time prediction, but I know it's has a ways to go and I'll need to do calculations from the STEREO spacecrafts to finish that also. Edit: I did make one more attempt and was finally able to contact the correct person.
  4. I'll post an image on roughly what I think the stages of CME impacts are. There are 3 CMEs, the latest CME, the CME on September 7th, and the CME from the final X-Class flare that arrived on the 12th. There are the velocities, density, and IMF total strength and Z component for each CME. If a CME hits and it is a leading edge, then we can get a decent idea on when the shock will start to come in. Subtracting the launch time from the arrival time, multiply that by the square root of pi/2, add that to the launch time and you have the time of around when the shock will start. With the CME at the bottom we were only hit by the shock just about as far as you can get from the leading edge.
  5. From what I see from the numbers, we took a hit from the leading edge and both time points I had ended up being shock points shock points instead of the first point being a leading edge point. If it's noticeable, I think the shock following the leading edge should hit at some point over the next several hours and shortly after that the solar wind will start to decrease.
  6. Not yet, but I'm thinking that might be what was causing some small differences in some of the first graphs I've made on this project. I want to get through all of the new data this year but so far I've been really busy. I'm waiting for this current CME though to see if we are getting a direct hit or side hit to try and backtrack it.
  7. This latest CME looked clear enough to get the needed points to try predicting the time this morning and now looking at the model from the SWPC, the times look like they are fairly close. Launch: 2/12/2018 00:30Z P1: 9:00Z P2: 15:00Z PDif: 6:00 P3: 12:00Z LaunchDif: 11:30Z Multiplier: 5.1 (Multiplier 3) Travel: 2 Days, 10.65 Hours Arrival: 2/14/2018 11:00Z If I read the CME correctly, the first point was the leading edge on the bottom right of the C3 imagery. What I think happened to the second point is the CME started folding at the shock behind the leading edge of the CME making it a shock point.
  8. Here is the math I have on CME arrival times. The one thing that all of the full halo CMEs have in common is this, 'y=(1/(40x+28500))(x+1765)^2-100' (or at least really really close). First, you need to get three times, 1st the time when the CME exits the sun in a SDO AIA 131 image, 2nd is the time of when the CME edge first hits the outer circle in the LASCO C3 image, and the 3rd from that point where the CME first hit the outer edge of the circle in the C3 image, go around the outer circle 180 degrees and get the time of when the CME hits that point. Then the math, get the difference of time #2 and #3 and write that down as time #4, get the time that is between time #2 and #3 and write that down as time #5, then get the difference of time #1 and time #5 and turn that into all minutes and write it down as number #6. Next, turn time #4 into minutes and put that number as x into the equation y=(1/(40x+28500))(x+1765)^2-100 and you now have #7. Now multiply #7 by #6 and you have the estimated travel time in minutes, just change that to days hours and minutes and add it to time #1. The issues are mainly if the image isn't clear like CME #1 a week ago, or there is a big chunk or missing images, then you are left to guessing the time of when it gets to the edge of the C3 circle edge and the time will likely be way off. Currently what I have is for CMEs where you need to guess one of either time #2 or #3 is likely +-3 hours off the arrival time and perfectly clear CMEs are at +-1.5 hours off the arrival time. For the equation, there is a lack of data before x=80 and after x=580, so I just went with the pattern and guessed and will need to modify the equation once there are some clear events to get that data from. The next thing I want to work on is get the arrival time from inside of the C3 image circle and not need to get times #2 and #3 from faint CMEs all the way at the edge of the C3 circle.
×
×
  • Create New...

Important Information

We have placed cookies on your device to help make this website better. You can adjust your cookie settings, otherwise we'll assume you're okay to continue. By using this site, you also agree to our Terms of Use and our Privacy Policy.