In amutaya/easyrbi: Calculates Richard's Baker Index Value for USGS Stream Sites

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.path = "man/figures/README-",
  out.width = "100%"
)

easyrbi R Package

Overview

The easyrbi R package consists of functions that retrieve streamflow data (pCode = 00060) from the USGS website, by using tools in the dataRetrieval package. Another set of functions can be used to calculate the Richards-Baker Flashiness Index (RBI) trends over time for any given gauge stations. The package also contains functions to download time series land use data USGS Gages II collected between 1931-01-01 and 2014-12-31. Some of the functions included in the package:

• sitedata: retrieves hydrology data for given USGS gauge sites for a specified period.
• rbi_df: returns calculated RBI values for given USGS sites for a specified period.
• trends: returns a data frame containing Mann-Kendall and sens.slope summary statistics for a custom number of sites specified by the user.
• basin_id: Check the drainage area, eco-region and classification (reference or non-reference)
• dam_removal: retrieves dam removal data from the USGS Gages II website for any number of given gauge station numbers.

Installation

You can install the development version of easyrbi from Github:

# If you haven't installed remotes yet, do so:
# install.packages("remotes")
remotes::install_github("amutaya/easyrbi")

Examples

library(easyrbi)
library(tidyverse)

Retrieve USGS gauge station data

Site data

• Returns a data.frames with 18 variables from the USGS database which include discharge, drainage area and the waterYear.

sitedata(c("01564500", "01567000"), "1970-10-01", "1980-09-30") %>% 
  head(10)

Calculate the annul Richards Barker Flashiness Index (RBI)

• Returns a data frame with the waterYear and respective RBI values for given sites over a specified time frame

rbi_df(c("01564500", "01567000"), "1970-10-01", "1980-09-30")

Test time series trends

• To test the flashiness trends over time, trends returns the Mann-Kendall and sens.slope summary statistics for any given USGS gauge sites.

data <- rbi_df(c("01564500", "01567000"), "1970-10-01", "1980-09-30")

trends(x = data) 

Basin classification and data

Watershed Identification

• Retrieve site classification (Reference and Non-reference sites) and watershed area from USGS GAGES II.

basin_id(c("01567000", "01490000", "01492500"))

Retrieve land use data

Dam removals

• Retrieves dam removal data for specified sites directly from the USGS Gages II website. This time series data recorded between 1931-01-01 and 2014-12-31.

dam_removal(c("01564500", "01567000"))

Dam removals and flashiness trends

• To check if there is a relationship between trends and dam removals, dam_trends returns a dataframe with dam removals and trends statistics from the Mann-Kendall and sens.slope tests.

dam_trends(c("01092000", "01100000", "01208500"), "1970-10-01", "1990-09-30")

amutaya/easyrbi documentation built on May 9, 2022, 4:04 p.m.