It has been a while but here is what I have been up to!
I met up with the family for a reunion and to celebrate my Dad's birthday in Las Vegas this past weekend. It was really great to see the family and take a couple days off from the computer near the midpoint of my time at USC. I cannot say that I caught up on my sleep while I was there but oh well.
Last week I started messing with the Matlab script that I will be using to compute S receiver functions. My lack of experience with Matlab coupled with my nonexistent knowledge of linear algebra made it really challenging to make sense of the whole thing and actually make the script work. Luckily Rob helped me out and we got it to start working earlier this week. In the end we changed a "b" to an "o". Woe is me. But now that the script is working my life has basically been consumed by looking through thousands of traces and either pressing 1 or 0 to keep the trace or get rid of it. I guess this is my first opportunity to contribute my own human error to the proccess. cheeky
Here is a little bit of background info on exactly what I am doing. Seismometers measure ground motion in three directions: E-W, N-S and vertical. S waves, or shear waves, are registered mostly on the horizontal plane, i.e., E-W and N-W. When I say that I am looking at traces, I mean that I am looking at the what each component of the siesmometer has recorded. I am focusing on the E-W and N-W traces and looking for distinct blips that would indicate an S wave arrival. Here are some examples of what I am looking for and what I never want to see.
The blips I am looking for, then some ugly noisy garbage, then the classic flatline, and then the art you can get IF you're lucky. Unfortunately, the pictures on my blog aren't working right now so the first and middle picture are missing... maybe I will fix that later. Maybe. Either way these are the kinds of squiggles that I am sifting through.
I have run into two problems while I have been working with the S wave data. The first one is not so bad. Basically some of the seismometers were either too close or too far away from the event that they were recording. We are looking for a distance between 55-86 radians. Any event where all/most of the distances are outside of this range goes into a "notgoodevents" folder and therefore is taken out of our dataset. They are also joined by ugly traces. The second problem that I have been having is more annoying but happens less often. Some of the event folders are missing station information either because some components were not recording or because certain component data was not retrieved by our script. Basically, some of the events are missing data and making the script crash. In order to fix this I have been going into the event folders and looking through 100 or so SAC files to try and figure out which station is missing information. It's basically a big matching game. In the end both of these problems mean getting rid of data which is bittersweet because we want as much data as we can get so that our results are reliable but at the same time getting rid of data will make calculations go a little big faster in the future.
In conclusion, the P receiver funtion script is still running... BUT there is finally an end in sight.
Even though I am a day early, HAPPY FOURTH!!
I have successfully downloaded all of my data and I am now in the thick of processing said data. It turns out that this step takes even longer! I started a P recevier function script last Tuesday and I am not very optimistic that it will be done by the end of this week. In the mean time I have been working with surface wave dispersion data and making preliminary maps with GMT to view the results. I have also started setting up directory structures and scripts to compute S receiver functions.
This weekend I went on two day trips to explore the San Andreas Fault with the SCEC program. Between seeing some really interesting geologic features, getting free food, waking up at 5 am, getting the bus stuck for over an hour in 102 degree weather, watching some of the best and worst geophysics-related movies ever, and hanging out with friends, the trip definitely had its ups and downs. Fun fact that I learned from this weekend: dates grow on palm trees. Who knew????!!? Definitely not me.
Today I started gathering USArray seismic data from the IRIS database. I am basically aquiring data from every seismometer in my study region that has recorded an event, i.e. an earthquake, in the last 8 years. Needless to say, it is going to take a while to download all of these files. I am going to be manipulating this data with scripts that have already been written by Dr. Robert Herrmann from Saint Louis University or by my advisor Rob Porritt. I will write/use some of my own scripts, but they are going to be extremely simplistic in comparison and will basically just rename/redirect files or run other scripts. Nevertheless, I will definitely get some experience with reading and writing bash scripts. I will be running some of these scripts through a program call gsac which was also written by Dr. Robert Herrmann (apparently he is a really smart guy). Later on I will probably use GMT to create some figures and maps. I am going to be working with the raw data, so I am basically starting from the ground up. I am pretty happy about this because I think that it will enable me to get a pretty sound understanding of whats happening to the data during each step of the process.
Enough with the science, I made it to the beach! My sister's friend lives in Santa Monica so I went and visited her this past Saturday. Very long story short, it was really weird and too fun. I also visited my Dad's childhood friend and his family in Encino on Sunday, so not only did I get beach time, but I also got pool time(!) along with some delicious food, a crazy dog that will play fetch until she dies, and some diehard Nascar fans. I took the bus to get to both of my weekend destinations, so all in all I probably spent a solid 4 or 5 hours on the bus. As awful as that may sound, it was actually really nice to get familiar with LA's transportation system and get a bus card for future adventures. Speaking of which, I will be using my bus card to go to a Dodgers game... TONIGHT! Maybe I will start the wave...
I am not much of a blogger, but here we go!
It's my fourth day at USC and everything has gone really well so far. The house that I am staying at for the first week of my internship is very pleasant and conveniently close to campus. My first night there I met two Bulgarians, Vera and Alexander, who introduced me to a very tasty Bulgarian soup, so what beats that?? At USC I have my very own desk amongst the graduate student cubicles, so I am feeling pretty official. Everyone in the department has been very nice and welcoming, which has made the transition a lot easier. My work so far has consisted of reading articles in order to familiarize myself with the study region and with current theories concerning the geologic processes going on over there under our feet.
My goals for this internship go as follows:
During the first third of my internship I want to get comfortable with the basic concepts behinds receiver functions, inversion theory, and common conversion point stacking. I do not think that I will be writing scripts this summer, but rather I will be using ones that Rob has already written. With that said, I want to learn general scripting conventions so that I can understand what Rob's scripts are doing on a fundmental level. After learning about the scripts I want to USE them and start doing data collection and processing. I want to make a skeleton of my poster for AGU and start thinking about the abstart that I will submit. Last but not least, I want to meet as many people as I can inside and outside of the department so that I can make some connections and (fingers crossed!) make some new friends.
During the second third of my internship I want to complete data processing and start anaylzing the resulting models. Having done some analysis, I hope to get my abstract a little more finalized. I also hope that I will have gotten to the beach by now!
During the final third of my internship I want to submit my abstrat to AGU and begin filling in my poster. Hopefully at this point I will know enough about what I am doing to start drafting a paper of my findings, eek!
Well, that is all that I have for now...
I am sunburned