Chapter 38 L6 Meteorite Analysis

38.1 Distribution of L6 Meteorite Landings

The following plot shows the distribution of the L6 meteorite landings all over the world.

MetLandings_L6 = MetLandings %>%
  filter(recclass == 'L6')

center_lon = median(MetLandings_L6$reclong,na.rm = TRUE)
center_lat = median(MetLandings_L6$reclat,na.rm = TRUE)


leaflet(MetLandings_L6) %>% addProviderTiles("Esri.NatGeoWorldMap") %>%
  addCircles(lng = ~reclong, lat = ~reclat, 
             color = c("red"))  %>%
  # controls
  setView(lng=center_lon, lat=center_lat,zoom = 1)

38.2 Summary of the L6 Meteorites Mass

summary(MetLandings_L6$mass)

##     Min.  1st Qu.   Median     Mean  3rd Qu.     Max.     NA's 
##      0.0      7.6     34.3   1694.6    202.0 564000.0        6

38.3 Time Series from Year 1970 for L6

The time series shows the Count of the L6 Meteorites and their corresponding Years.

MetLandings_L6 %>%
    filter(year >= 1970) %>%
    group_by(year) %>%
    summarise(Count = n()) %>%
    arrange(year) %>%

    ggplot(aes(x = year,y = Count)) +
    geom_bar(stat='identity',colour="white", fill = fillColor) +
    
    labs(x = 'Meteorites Year', 
         y = 'Count', 
         title = 'Meteorites Year and Count') +
        theme_bw()

38.4 Distribution of L6 US Meteorite Landings with Meteorite Mass

The mass of the Meteorites are indicated by the Radius of the Circles.

top = 49.3457868 # north lat
left = -124.7844079 # west long
right = -66.9513812 # east long
bottom =  24.7433195 # south lat


USMetLandings = MetLandings_L6 %>%
  filter(reclat >= bottom) %>%
  filter ( reclat <= top) %>%
  filter( reclong >= left ) %>%
  filter(reclong <= right)

center_lon = median(USMetLandings$reclong,na.rm = TRUE)
center_lat = median(USMetLandings$reclat,na.rm = TRUE)

factpal <- colorFactor(c("red","blue"), 
                       USMetLandings$fall)

leaflet(USMetLandings) %>% addProviderTiles("Esri.NatGeoWorldMap") %>%
  addCircles(lng = ~reclong, lat = ~reclat,radius = ~(mass) ,
             color = ~factpal(fall))  %>%
  # controls
  setView(lng=center_lon, lat=center_lat,zoom = 4) %>%
  
  addLegend("bottomright", pal = factpal, values = ~fall,
            title = "Meteorites landings and fall",
            opacity = 1)

38.5 Distribution of Indian Meteorite L6 Landings with Meteorite Mass

The mass of the Meteorites are indicated by the Radius of the Circles.

IndiaMetLandings = IndiaMetLandings %>% filter(recclass == 'L6')

center_lon = median(IndiaMetLandings$reclong,na.rm = TRUE)
center_lat = median(IndiaMetLandings$reclat,na.rm = TRUE)


factpal <- colorFactor(c("red"), 
                       IndiaMetLandings$fall)

leaflet(IndiaMetLandings) %>% addProviderTiles("Esri.NatGeoWorldMap") %>%
  addCircles(lng = ~reclong, lat = ~reclat,radius = ~(mass) ,
             color = ~factpal(fall))  %>%
  # controls
  setView(lng=center_lon, lat=center_lat,zoom = 5) %>%
  
  addLegend("bottomright", pal = factpal, values = ~fall,
            title = "Meteorites landings and fall",
            opacity = 1)