Essays, Posts & Presentations

Fire in the Hole! Exploring the Yellowstone Calderas with GeoGraphics and USGS Data

Jeffrey Bryant

Author

Jeffrey Bryant

Title

Fire in the Hole! Exploring the Yellowstone Calderas with GeoGraphics and USGS Data

Description

A computational exploration of features of the Yellowstone Caldera supervolcano using the GeoGraphics function in the Wolfram Language. Code provided to generate a relief map and mine for details.

Fire in the Hole! Exploring the Yellowstone Calderas with GeoGraphics and USGS Data

January 14, 2020

Jeffrey Bryant, Research Programmer, Wolfram|Alpha Scientific Content

Yellowstone National Park has long been known for its active geysers. These geysers are a surface indication of subterranean volcanic activity in the park. In fact, Yellowstone is actually the location of the Yellowstone Caldera, a supervolcano: a volcano with an exceptionally large magma reservoir. The park has had a history of many explosive eruptions over the last two million years or so.

I’ve found that the United States Geological Survey (USGS) maintains data on the various volcanic calderas and related features, which makes it perfect for computational exploration with the Wolfram Language. This data is in the form of SHP files and related data stored as a ZIP archive. Thanks to the detail of this available data, we can use the Wolfram Language and, in particular,

GeoGraphics

to get a better picture of what this data is telling us.

Generating a Relief Map First, I want to examine the area of interest. The following code generates a relief map of the data within 50 miles of the center of the Yellowstone Caldera:

GeoGraphicsGeoBackground"ReliefMap",

Rule[

]

,GeoCenter

Yellowstone Caldera

VOLCANO

,GeoRangeQuantity[50,"Miles"]

Out[]=

Although there are some interesting features visible in the data, it's not immediately clear where the volcanic calderas are, although some possible candidates are visible as flat areas among the surrounding mountainous terrain.

Mining for Details In order to gain a better understanding of the region, let’s pull in some additional data from the USGS in the form of an SHP file and the corresponding support files stored in the ZIP archive. Ideally, we would be working with the data straight from the Wolfram Data Repository, which makes storing, importing and computing data easy. It may be some time before this data is submitted, due to the complexity of the SHP file format and sheer volume of metadata; so in the meantime, this is a good example of working with data “in the wild,” without the benefit of the standardization found in the Data Repository. In order to use this data, we need to understand the datum used so it can be properly transformed for representation using GeoGraphics:

In[]:=

csi=Import["https://www.wsgs.wyo.gov/gis-files/caldera.zip","CoordinateSystemInformation"]

Out[]=

PROJCSNAD_1983_UTM_Zone_12N,GEOGCSGCS_North_American_1983,DATUM{North_American_Datum_1983,SPHEROID{GRS_1980,6.37814×

,298.257}},PRIMEM{Greenwich,0.},UNIT{Degree,0.0174533},PROJECTION{Transverse_Mercator},PARAMETER{False_Easting,500000.},PARAMETER{False_Northing,0.},PARAMETER{Central_Meridian,-111.},PARAMETER{Scale_Factor,0.9996},PARAMETER{Latitude_Of_Origin,0.},UNIT{Meter,1.}

Next, I extract the data and labels:

In[]:=

shpdata=Import["https://www.wsgs.wyo.gov/gis-files/caldera.zip","Data"];

The data has various “properties” that can be extracted:

In[]:=

shpdata[[1,All,1]]

Out[]=

{LayerName,Geometry,Labels,LabeledData}

For my purposes, I only need the "Geometry" and "LabeledData" properties:

In[]:=

labeleddata="LabeledData"/.shpdata[[1]]

Out[]=

{LTYPE{Outer Ring Fault,Inner Ring Fault,Resurgent Dome,Outer Ring Fault,Outer Ring Fault,Outer Ring Fault,Outer Ring Fault,Outer Ring Fault,Outer Ring Fault,Outer Ring Fault,Caldera Rim,Outer Ring Fault,Outer Ring Fault,Outer Ring Fault,Resurgent Dome,Resurgent Dome,Caldera Rim,Caldera Rim,Caldera Rim,Outer Ring Fault},cycle{West Thumb,III,III,II,II,I,I,I,I,I,III,III,III,III,III,III,III,III,III,III},Shape_Leng{23741.7,175351.,48118.3,60942.8,9980.28,11407.1,725.149,40874.9,35689.,196769.,13871.6,11718.8,86414.3,57881.,37795.4,1197.5,43561.9,44777.,120815.,44777.},Shape_Le_1{23741.7,175351.,48118.3,60942.8,9980.28,11407.1,725.149,40874.9,35689.,196769.,13871.6,11718.8,86414.3,57881.,37795.4,1197.5,43561.9,44777.,120815.,44777.}}

In[]:=

labels="LTYPE"/.labeleddata;

In[]:=

lines="Geometry"/.shpdata[[1]];

With the data extracted, it’s now in suitable shape to add to our map.

Making Mountains out of Data Points Let’s see what the data can reveal to us. I can use the "CoordinateSystemInformation" we extracted to transform the data via

GeoGridPosition

In[]:=

convertedlines=lines/.Line[data_]Line[GeoGridPosition[data,"UTMZone12","NAD831986"]];

Now let’s give this some style! I choose some distinct colors to represent the various labels found in the data:

In[]:=

color["Outer Ring Fault"]=White;color["Inner Ring Fault"]=Red;color["Resurgent Dome"]=Yellow;color["Caldera Rim"]=Black;

Then I combine the colors and labels, and add a tooltip:

In[]:=

coloredtooltip[line_,label_]:={color[label],Tooltip[line,label]}

For a sense of scale, I include the state outlines for Wyoming, Idaho and Montana, with the data drawn using colored lines:

In[]:=

GeoGraphicsEdgeForm[Black],GeoStyling[Opacity[.4]],Red,Polygon

Montana, United States

ADMINISTRATIVE DIVISION

,Purple,Polygon

Idaho, United States

ADMINISTRATIVE DIVISION

,Blue,Polygon

Wyoming, United States

ADMINISTRATIVE DIVISION

,Thickness[.003],coloredtooltip@@@Transpose[{convertedlines,labels}],GeoBackground"ReliefMap",ImageSize600,GeoCenter

Yellowstone Caldera

VOLCANO



Out[]=

The extent of the Yellowstone calderas is quite impressive! Now that we have an idea of the larger scope, I want to try and get a bit more detail with a closer look at the area.

A useful surface marker we can use as a reference point is the location of Yellowstone Lake, which we’ll visually distinguish with a white rectangle:

In[]:=

yellowstonelake=Rectangle[GeoPosition[{Latitude[#[[1]]],Longitude[#[[2]]]}],GeoPosition[{Latitude[#[[3]]],Longitude[#[[4]]]}]]&@

Yellowstone Lake

LAKE

[{"SouthernmostPoint","WesternmostPoint","NorthernmostPoint","EasternmostPoint"}];

Zooming in on Yellowstone with a radius of 50 miles will show more detail:

In[]:=

GeoGraphicsEdgeForm[Black],GeoStyling[Opacity[.4]],Red,Polygon

Montana, United States

ADMINISTRATIVE DIVISION

,Purple,Polygon

Idaho, United States

ADMINISTRATIVE DIVISION

,Blue,Polygon

Wyoming, United States

ADMINISTRATIVE DIVISION

,{EdgeForm[],White,Opacity[.5],yellowstonelake},Thickness[.007],coloredtooltip@@@Transpose[{convertedlines,labels}],GeoBackground"ReliefMap",ImageSize600,GeoCenter

Yellowstone Caldera

VOLCANO

,GeoRangeQuantity[50,"Miles"]

Out[]=

The Feature Presentation As a result of visually consolidating the available data and its properties, we can better understand the interactions and relationships between them. This map represents a flattened timescale of the Yellowstone calderas’ activity and their surrounding area, shaped over many millennia. We can see that some of the faults and rims of the various calderas correspond to visible surface features. A pair of resurgent domes (the yellow circular features) can be found within the data; these domes are areas where uplift occurred due to follow-up eruptions within older calderas. Through its visualization, the data reveals several calderas in the location of the largest—which we know by name as the Yellowstone Caldera.

Analyzing even this one dataset with the Wolfram Language allows us some perspective into geographical features that appear to be immutable from our point of view, yet are actually full of life and transformation. And fortunately for us, there’s a surfeit of data out there, just waiting to be given form—I look forward to seeing what you can create out of it!

Cite this as: Jeffrey Bryant, "Fire in the Hole! Exploring the Yellowstone Calderas with GeoGraphics and USGS Data" from the Notebook Archive (2019), https://notebookarchive.org/2020-01-6h2ohwn