DRAFT DRAFT DRAFT

Kurt Schwehr
Copyright (C) 2003
kdschwehr _at_ ucsd dot edu



Viz Lecture Series
Lecture 8
Building Eel River
$Revision: 914 $



CONTENTS

INTRODUCTION
ETOPO2
DSDP
ODP
L. MAYER'S CD
SEISMICS
NGDC SAMPLES
WELL LOGS



* INTRODUCTION

This lecture will put together a full presentation of the Eel
River/Humboldt Slide structure for my thesis.  The focus is on
determining if the Humboldt Slide is actually a slump/slide or
sediment waves.


* etopo2/topo8.2

Make etopo2 and Smith/Sandwell topo8.2 maps.

*ODP/DSDP Hole locations

ODP = Ocean Drilling Project
DSDP = Deep Sea Drilling Project

We would like to know if there are any ODP/DSDP drill holes in the
area.  So first we need to get a table of the hole locations.  We can
build one from the available web pages.  The first is the DSDP CD-ROM
Hole list. 

    wget http://www.ngdc.noaa.gov/mgg/geology/dsdp/data/holelist.htm

The file is an html wrapped table in the format of hole, leg, lat,
lon, and ocean basin.  We can write a small script to pull out these
columns from this data:

    hole    leg     latitude        longitude       ocean/sea

    1   1        25.8583         -92.1833       Gulf of Mexico
    2   1        23.0455         -92.0587       Gulf of Mexico
    3   1        23.0300         -92.0433       Gulf of Mexico


    grep index holelist.htm > holelist2.htm
    awk '{print $6,$7,$8,$9}' holelist2.htm | sort -u > basins.txt
    ./dsdpHoles.py holelist2.htm > dsdpHoles.txt

This gives us the whole world.  What about if we just want an area?
We can write another python tool that only passes points within a
long/lat region.

    #!/usr/bin/env python
    import sys
    inFileName=sys.argv[1]
    lonMin=float (sys.argv[2])
    lonMax=float (sys.argv[3])
    latMin=float (sys.argv[4])
    latMax=float (sys.argv[5])
    infile = open(inFileName,"r");
    for line in infile.xreadlines():
        s = line.split()
        lon = float(s[0])
        lat = float (s[1])
        if (lonlonMax): continue
        if (latlatMax): continue
        print line

Now we can find out what is in the neighborhood:

    ./subRegion.py dsdpHoles.txt -128 -124 39 43

Which returns:

-127.4968 39.4747 33 http://www.ngdc.noaa.gov/mgg/geology/dsdp/data/5/33/index.htm

    -127.2757 39.4702 34 http://www.ngdc.noaa.gov/mgg/geology/dsdp/data/5/34/index.htm
    -127.4747 40.6737 35 http://www.ngdc.noaa.gov/mgg/geology/dsdp/data/5/35/index.htm
    -127.5202 40.8392 35A http://www.ngdc.noaa.gov/mgg/geology/dsdp/data/5/35A/index.htm
    -125.4520 39.9618 173 http://www.ngdc.noaa.gov/mgg/geology/dsdp/data/18/173/index.htm

* ODP SITE MAP 

    How to get the Long/Lat?
    http://www-odp.tamu.edu/sitemap/sitemap.html

    http://www-odp.tamu.edu/database/ODP_trivia.xls

    wget http://www.ngdc.noaa.gov/mgg/geology/odp/data/holelist.htm
    mv holelist.htm odp-holelist.htm

    grep index odp-holelist.htm > odp-holelist2.htm
    
    Then edited with emacs.


* MAYER ET AL. CD

    GIS presentation

    
* Seismics


    http://wedge.ig.utexas.edu/Web/cruise_html/NP09/W9605B.htm
    http://www.ig.utexas.edu/research/projects/strataform/ca_mcs/ca_mcs.htm


 Robin.R.Warnken@noaa.gov
Subject: Price of data cd?
Reply-To: kdschweh@ucsd.edu


Hi,

I'm interested int the 24 channel seismic data from northern CA and
oregon.  How much is the CD or is there some place I could just download
the data from an FTP or HTTP site?  I'd like to try to use the data in my
thesis.

On the web page (http://www.ngdc.noaa.gov/mgg/fliers/90-MGG-02) the data
is listed as:

    CD-ROM, Product #129F27001.


* NGDC Samples

    http://www.ngdc.noaa.gov/mgg/curator/cursrch.HTML

* Well logs

    http://oas.ngdc.noaa.gov/mgg/plsql/geolin.set_expand?v_mggid=60945002
    http://www.ngdc.noaa.gov/mgg/fliers/83mgg03.html



* MODIS SAT PHOTOS


http://visibleearth.nasa.gov/cgi-bin/viewrecord?25694
http://visibleearth.nasa.gov/data/ev256/ev25694_Oregon.A2003183.1900.250m.jpg