Calculate non-point source (NPS) loads for Tampa Bay

Usage

anlz_nps(
  yrrng = c("2021-01-01", "2023-12-31"),
  tbbase,
  rain,
  mancopth = NULL,
  pincopth = NULL,
  lakemanpth = NULL,
  tampabypth = NULL,
  bellshlpth = NULL,
  vernafl,
  allflo = NULL,
  allwq = NULL,
  usgsflow = NULL,
  summ = c("basin", "segment", "all"),
  summtime = c("month", "year"),
  aslu = FALSE,
  verbose = TRUE
)

Arguments

yrrng: A vector of two dates in 'YYYY-MM-DD' format, specifying the date range to retrieve flow data. Default is from '2021-01-01' to '2023-12-31'.
tbbase: data frame containing polygon areas for the combined data layer of bay segment, basin, jurisdiction, land use data, and soils, see details
rain: data frame of rainfall data, see details
mancopth: character, path to the Manatee County water quality data file, see details
pincopth: character, path to the Pinellas County water quality data file, see details
lakemanpth: character, path to the file containing the Lake Manatee flow data, see details
tampabypth: character, path to the file containing the Tampa Bypass flow data, see details
bellshlpth: character, path to the file containing the Bell shoals data, see details
vernafl: character vector of file path to Verna Wellfield atmospheric concentration data
allflo: data frame of flow data, if already available from util_nps_getflow, otherwise NULL and the function will retrieve the data
allwq: data frame of water quality data, if already available from util_nps_getwq, otherwise NULL and the function will retrieve the data.
usgsflow: data frame of USGS flow data, if already available from util_nps_getusgsflow, otherwise NULL and the function will retrieve the data. Default is NULL.
summ: chr string indicating how the returned data are summarized, see details
summtime: chr string indicating how the returned data are summarized temporally (month or year), see details
aslu: logical indicating whether to summarize by land use type (ungaged loads only), default is FALSE
verbose: logical indicating whether to print verbose output

Value

A data frame of non-point source loads for Tampa Bay, including columns for year, month, bay segment, basin, and loads for total nitrogen (TN), total phosphorus (TP), total suspended solids (TSS), biochemical oxygen demand (BOD), and hydrology using default values for the summ and summtime arguments. TN, TP, TSS, and BOD Loads are tons per month or year depending on the summtime argument. Hydrologic loads are cubic meters per month or year depending on the summtime argument.

Details

The function estimates non-point source (NPS) loads for Tampa Bay by combining ungaged and gaged NPS loads. Ungaged loads are estimated using rainfall, flow, event mean concentration, land use, and soils data, while gaged loads are estimated using water quality data and flow data. The function also incorporates atmospheric concentration data from the Verna Wellfield site.

The data are summarized differently based on the summ and summtime arguments. All loading data are summed based on these arguments, e.g., by bay segment (summ = 'segment') and year (summtime = 'year'). Options for summ are 'basin' to summarize across sub-basins within bay segments, 'segment' to summarize by bay segment, and 'all' to summarize total load. Loads can also be summarized by land use type with the summ and summtime argumets by setting aslu = TRUE. Land use type summaries only apply to ungaged load estimates. Options for summtime are 'month' to summarize by month and 'year' to summarize by year. The default is to summarize by basin and month.

The following functions are used internally and are provided here for reference on the components used in the calculations:

anlz_nps_ungaged: Estimates ungaged NPS loads.
anlz_nps_gaged: Estimates gaged NPS loads.
util_nps_fillmiswq: Fills missing water quality data with linear interpolation.
util_nps_getflow: Gets flow estimates for NPS gaged and ungaged calculations.
util_nps_getusgsflow: Gets USGS flow data for NPS calculations, used in util_nps_getflow.
util_nps_getextflow: Gets external flow data (Lake Manatee, Tampa Bypass, and Bell Shoals), used in util_nps_getflow.
util_nps_getwq: Gets water quality data for NPS gaged calculations.
util_nps_preprain: Prepares and formats rainfall data.
util_nps_preplog: Prepares land use data for logistic regression modeling.
util_nps_segment: Assigns basins to bay segments.
util_prepverna: Prepares and fills missing data with five-year means for the Verna Wellfield site data.

Examples

data(tbbase)
data(rain)
data(allwq)
data(allflo)
vernafl <- system.file('extdata/verna-raw.csv', package = 'tbeploads')

nps <- anlz_nps(
  yrrng = c('2021-01-01', '2023-12-31'), 
  tbbase = tbbase, 
  rain = rain, 
  allwq = allwq,
  allflo = allflo,  
  vernafl = vernafl, 
)
#> Estimating ungaged NPS loads...
#> Estimating gaged NPS loads...
#> Combining atmospheric data with ungaged NPS loads...
#> Combining ungaged and gaged NPS loads, estimating final...

head(nps)
#> # A tibble: 6 × 10
#>    Year Month source segment     basin tn_load tp_load tss_load bod_load hy_load
#>   <dbl> <dbl> <chr>  <chr>       <chr>   <dbl>   <dbl>    <dbl>    <dbl>   <dbl>
#> 1  2021     1 NPS    Boca Ciega… 207-5    2.43   0.398     80.8    14.5   1.22e6
#> 2  2021     2 NPS    Boca Ciega… 207-5    1.65   0.270     54.8     9.85  8.26e5
#> 3  2021     3 NPS    Boca Ciega… 207-5    1.37   0.224     45.5     8.18  6.86e5
#> 4  2021     4 NPS    Boca Ciega… 207-5    1.58   0.258     52.5     9.43  7.91e5
#> 5  2021     5 NPS    Boca Ciega… 207-5    1.21   0.198     40.2     7.23  6.06e5
#> 6  2021     6 NPS    Boca Ciega… 207-5    2.68   0.440     89.4    16.1   1.35e6