The map below shows earthquakes from the NGDC/WDS with economic damage exceeding $1bn. Circles and hover functionality for magnitude.
Source: National Geophysical Data Center / World Data Service (NGDC/WDS): Significant Earthquake Database. National Geophysical Data Center, NOAA.
library(leaflet)
df<-read.csv("signif.csv",stringsAsFactors=FALSE)
df2<-subset(df,df$LONGITUDE!="" & df$DAMAGE_MILLIONS_DOLLARS>1000)
df3<-data.frame(lat=df2$LATITUDE,lng=df2$LONGITUDE,magnitude=df2$EQ_PRIMARY)
quakeIcons<-makeIcon(iconUrl="https://cdn4.iconfinder.com/data/icons/san-francisco-cityset/32/earthquake-512.png",iconWidth=15,iconHeight=15)
leaflet(data=df3) %>% addTiles() %>% addCircles(weight=1,radius=5^df3$magnitude) %>% addMarkers(~lng,~lat,icon=quakeIcons,popup=df3$intensity,label=as.character(df3$magnitude))## Assuming "lng" and "lat" are longitude and latitude, respectively