Download FDEP Upper Floridan Aquifer potentiometric surface contour lines

Usage

util_gw_getcontour(
  season = c("dry", "wet"),
  yr,
  max_records = 1000,
  north_dist = 0,
  verbose = TRUE
)

Arguments

season: character, "dry" or "wet"
yr: integer, year for which to retrieve data. Biannual (May/September) observations are available from 2010 through 2022.
max_records: integer, maximum number of records per paginated request. Default is 1000.
north_dist: numeric, distance in CRS units (US Survey Feet for EPSG 6443) to extend the spatial filter and clipping boundary northward beyond the Tampa Bay watershed (tbfullshed). Use a positive value when the potentiometric high point for one or more bay segments lies north of the watershed boundary (e.g., Old Tampa Bay). The same value (or larger) should be passed as north_dist in util_gw_grad so that the returned contours cover the extended search areas. Default 0 (no extension).
verbose: logical, if TRUE (default) progress messages are printed during download.

Value

An sf object of LINESTRING features with columns CONTOUR (integer, feet MSL) and MONTH_YEAR (character), in the same CRS as tbfullshed (EPSG 6443). Returns NULL if no features are found for the requested season/year.

Details

Downloads contour lines representing the potentiometric surface of the Upper Floridan Aquifer from the Florida Department of Environmental Protection / Florida Geological Survey ArcGIS REST service (https://ca.dep.state.fl.us/arcgis/rest/services/OpenData/FGS_PUBLIC/MapServer/8).

Contours are available biannually: "dry" season maps to May of yr and "wet" season maps to September of yr. Results are spatially filtered to the Tampa Bay watershed (tbfullshed), optionally extended northward by north_dist, and clipped to that boundary before return.

The CONTOUR field contains potentiometric surface elevations in feet above mean sea level. These are used to compute the hydraulic gradient driving Floridan Aquifer discharge to Tampa Bay segments (Darcy's Law).

Examples

if (FALSE) { # \dontrun{
dry_contours <- util_gw_getcontour("dry", 2022)
wet_contours <- util_gw_getcontour("wet", 2022)

# extend north by ~28 miles for OTB high-point search
dry_contours <- util_gw_getcontour("dry", 2022, north_dist = 150000)
} # }