AC4: Visualize Major Earthquakes around world and its relationship with volcanoes

#install.packages("mapproj") #install.packages("ggmap")
#install.packages("DeducerSpatial") #install.packages("leaflet")
dfv = read.csv("volcano.csv", header=TRUE) # volcano data
dfqb = read.csv("quakes_past_10_year_gt_60.csv", header=TRUE) #quake data

Data

• Volcano data is gathered from naaa website contains data of 1572 volcanoes around the globe
• Eartgquake data gathered from usgs website contains earthquake of > 6.0 mag around the world (past 10 years) [1630 datapoints]
• Please check the reference below for the link to the source of data

Analysis

we analyze following 3 way relation

1. Visualize the relationship between volcanoes and earthquakes.

• From the interactive leaflet map below it is clear that most of the earthquakes are in the vicinity of volcanoes
• Volcanoes are marked with red markers which the earthquake data is Blue or green (Based on intensity)
• Notice that many of the valconoes are situated at continent boundry
• Interactive map at http://rpubs.com/dpant123/quake contains popup feature if a geolocation is clicked

2. Earthquake magnitude as a function of distance from equator (0° latitude) and the prime meridian (0° longitude) and its epicenter depth

• There are 4 seperate cluster based on the haversin distance from the lat,lng =(0,0). Again near the point (0,0) there are less number of quakes

3. Earthquake magnitude as a function of longitude or latitude.

• From the scatter plot below it is evident that very few earthquakes are near the equator.
• Higher magnitude earthquake are shallow in depth, encoded by the marker colors
• Fewer quakes aorund the polar region

library(geosphere)
library(ggplot2)
library(leaflet)
library(gridExtra)
library(data.table)
dfqb$depthColor = cut(dfqb$depth, seq(from = min(dfqb$depth)-1, to = max(dfqb$depth)+1,
                                      length.out= 4))
dfqb$magfact = cut(dfqb$mag, seq(from = min(dfqb$mag)-1, to = max(dfqb$mag)+1, 
                                 length.out= 3))
levels(dfqb$magfact) = c("green", "blue")
p1 = ggplot(dfqb) + geom_point(aes(x=longitude,y=mag,color=dfqb$depthColor))
p2 = ggplot(dfqb) + geom_point(aes(x=latitude,y=mag,color=dfqb$depthColor))
dfqb$distHaversinedistfrom00 = distHaversine(matrix(c(dfqb$longitude,
                                             dfqb$latitude),ncol=2),
                                             matrix(c(rep(0,each=length(dfqb$longitude)),                                                                                        rep(0,each=length(dfqb$longitude))), ncol = 2))
p3 = ggplot(dfqb) + geom_point(aes(x=distHaversinedistfrom00,y=mag))
grid.arrange(p1,p2,p3,ncol=1)

# Leaflet map
leaflet(dfqb) %>%
      addTiles() %>%  # Add default OpenStreetMap map tiles
      addCircleMarkers(lat=dfqb$latitude, lng=dfqb$longitude, popup=dfqb$time,
                 fillOpacity = 0.1,color= dfqb$magfact,weight = 3, radius=2) %>%
      addCircleMarkers(lat=dfv$Latitude, lng=dfv$Longitude,  weight = 3, radius=2,
                       popup=dfv$VolcanoName,color="red", stroke = FALSE, fillOpacity = 0.8) 

References

volcanoes data

http://bit.ly/2mzyXwk

earthquake data

https://earthquake.usgs.gov/earthquakes/search/