Skip to contents

Single Site Water Level

Single_site_water_level.Rmd

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 

library(HASP)

# siteID <- "263819081585801"
# site_metadata <- site_summary(siteID, markdown = TRUE)
 
# data_info <- data_available(siteID)
# kable(data_info)

Get data using dataRetrieval

Using dataRetrieval code to get the data for this site:

library(dataRetrieval)
siteID <- "USGS-263819081585801"
parameterCd <- "62610"
statCd <- "00001"
# Daily data:
gw_level_dv <- read_waterdata_daily(monitoring_location_id = siteID,
                                    parameter_code = parameterCd,
                                    statistic_id = statCd, 
                                    skipGeometry = TRUE)

# Field GWL measured:
gwl_data <- dataRetrieval::read_waterdata_field_measurements(monitoring_location_id = siteID,
                                                             skipGeometry = TRUE)

site_title <- "L -2701\n263819081585801"

y_label <- dataRetrieval::read_waterdata_parameter_codes(parameter_code = parameterCd)$parameter_name

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)
week p05 p10 p25 p50 p75 p90 p95 nYears minMed maxMed week_start
1 -37.7 -32.9 -28.7 -20.7 -16.2 -10.5 -8.7 45 -40.6 -5.5 01-01
2 -38.0 -32.5 -28.3 -20.2 -16.2 -10.7 -7.5 46 -41.4 -5.5 01-08
3 -37.9 -32.5 -28.5 -20.5 -16.4 -10.8 -6.5 46 -42.1 -5.7 01-15
4 -39.0 -31.8 -28.6 -21.0 -16.8 -10.9 -7.1 46 -42.4 -6.1 01-22
5 -39.9 -31.7 -29.2 -21.2 -17.1 -11.4 -7.2 47 -41.6 -6.8 01-29
6 -40.1 -32.1 -28.9 -21.3 -17.4 -11.9 -7.9 47 -42.1 -7.0 02-05
7 -39.9 -32.5 -28.8 -21.4 -17.7 -12.3 -7.8 47 -43.5 -7.4 02-12
8 -40.7 -33.4 -29.2 -21.9 -18.2 -11.4 -7.9 47 -45.0 -6.1 02-19
9 -41.3 -34.1 -29.4 -21.8 -18.3 -11.3 -7.9 47 -42.6 -6.2 02-26
10 -41.1 -34.3 -29.6 -22.2 -18.3 -11.4 -7.7 47 -43.3 -6.3 03-04

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))
DOY max mean min points
1 -5.29 -21.59489 -41.40 45
2 -5.24 -21.64733 -40.11 45
3 -5.30 -21.42957 -40.18 46
4 -5.55 -21.73378 -41.60 45
5 -5.48 -22.12106 -40.51 47
6 -5.71 -21.67587 -40.64 46
7 -5.64 -21.87533 -42.25 45
8 -5.75 -21.71109 -41.39 46
9 -5.52 -21.65447 -40.51 47
10 -5.30 -21.71739 -40.57 46 
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) 

kable(siteDV, digits = 1)
site min_site max_site mean_site p10 p25 p50 p75 p90 count
USGS-263819081585801 -51.1 -2.8 -21.9 -33.4 -27.5 -21 -15.8 -11.4 16275

Results of trend analysis on maximum daily water levels:

trend_result <- trend_test(gw_level_dv = gw_level_dv, 
                           gwl_data = NULL,
                           parameter_cd = parameterCd)
kable(trend_result, digits = 1)
test tau pValue slope intercept trend
10-year trend 0.6 0 1.2 -2543.1 Up
Period of record -0.2 0 -0.1 260.2 Down

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 |> 
  site_data_summary(site_col = "monitoring_location_id") 

kable(quantiles, digits = 1)
site min_site max_site mean_site p10 p25 p50 p75 p90 count
USGS-263819081585801 -52.4 64.3 -4.6 -34.8 -29 -20.1 32.5 43.2 1461