In USGS-R/HASP: Hydrologic AnalySis Package

library(knitr)
opts_chunk$set(echo = TRUE,
               warning = FALSE,
               message = FALSE,
               fig.width = 7,
               fig.height = 6)
library(dplyr)
library(tidyr)
library(ggplot2)
library(HASP)

options(knitr.kable.NA = '')

This vignette shows how to use HASP and other R tools to reproduce:

https://fl.water.usgs.gov/mapper/waterlevel_site_info.php?site=263819081585801

Site Information

wzxhzdk:1

wzxhzdk:2

Get data using dataRetrieval

Using `dataRetrieval` code to get the data for this site: wzxhzdk:3

Using the example data provided in the package: wzxhzdk:4

Water Level Data and Analysis

Weekly frequency

weekly_frequency_plot(gw_level_dv, 
                      gwl_data = gwl_data,
                      parameter_cd = parameterCd,
                      plot_title = site_title,
                      y_axis_label = y_label)

Weekly frequency analysis of daily maximum water level record. Only showing the first 10 rows for this example:

weekly_table <- weekly_frequency_table(gw_level_dv, 
                                       gwl_data,
                                       parameter_cd = parameterCd)
kable(head(weekly_table, 10), digits = 1)

Daily 2-year

daily_gwl_plot(gw_level_dv, 
               gwl_data,
               parameter_cd = parameterCd,
               plot_title = site_title,
               historical_stat = "mean",
               month_breaks = TRUE,
               y_axis_label = y_label)

Statistics of maximum daily water level record (DOY = day of year). Only showing the first 10 rows for this example:

daily_table <- daily_frequency_table(gw_level_dv, 
                                     gwl_data,
                                     parameter_cd = parameterCd) 

kable(head(daily_table, 10))

Daily value trends

gwl_plot_all(gw_level_dv, gwl_data, 
             y_label = y_label,
             parameter_cd = parameterCd,
             plot_title = site_title,
             add_trend = TRUE, flip = FALSE)

Summary statistics for maximum daily water level measurements:

siteDV <- site_data_summary(gw_level_dv,
                            value_col = "X_62610_00001") 

kable(siteDV, digits = 1)

Results of trend analysis on maximum daily water levels:

trend_result <- trend_test(gw_level_dv, 
                          NULL,
                          parameter_cd = parameterCd)
kable(trend_result, digits = 1)

Field GWL values

gwl_plot_all(NULL, gwl_data,
             y_label = y_label,
             parameter_cd = parameterCd,
             plot_title = site_title, 
             flip = FALSE)

Summary statistics for manual water level measurements

quantiles <- gwl_data %>% 
  rename(value = sl_lev_va) %>% 
  site_data_summary() 

kable(quantiles, digits = 1)

USGS-R/HASP documentation built on July 28, 2024, 7:53 a.m.