% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/anthropogenic.R
\name{score_anthropogenic_disturbance}
\alias{score_anthropogenic_disturbance}
\title{Score anthropogenic disturbance at survey locations}
\usage{
score_anthropogenic_disturbance(
  result,
  trawling_data = NULL,
  anchor_data = NULL,
  dredge_areas = NULL,
  fetch_live = FALSE,
  match_radius_m = 5000,
  verbose = TRUE
)
}
\arguments{
\item{result}{Dataframe from \code{\link[=predict_oyster]{predict_oyster()}}. Must contain \code{lat}, \code{lon}.}

\item{trawling_data}{Dataframe of ICES VMS SAR data (see Details), or NULL.}

\item{anchor_data}{Dataframe of anchoring density data, or NULL.}

\item{dredge_areas}{Dataframe of dredging/extraction area centroids, or NULL.}

\item{fetch_live}{Logical. Query ICES VMS GeoServer (default FALSE).}

\item{match_radius_m}{Numeric. Spatial matching radius in metres (default
5000 m, consistent with ICES c-square ~0.05degrees resolution).}

\item{verbose}{Logical. Default TRUE.}
}
\value{
\code{result} with additional columns:
\itemize{
\item \code{disturbance_sar}: swept-area ratio at or near each site (NA if no data)
\item \code{disturbance_score} [0,1]: composite anthropogenic disturbance index
\item \code{disturbance_class}: "Negligible" / "Low" / "Moderate" / "High" / "Active"
\item \code{disturbance_gear_gate}: logical -- TRUE if SAR exceeds gear deployment
threshold (used by \code{assess_gear_feasibility()})
\item \code{disturbance_note}: management flag
}
}
\description{
Calculates a disturbance index [0,1] from bottom trawling intensity,
dredging activity, and anchor damage. The primary data input is the
ICES swept-area ratio (SAR), which can be supplied manually or fetched
live from the ICES VMS data portal.

This output is primarily relevant to the gear feasibility and economic
viability modules. Very high trawling intensity (SAR > 0.5/yr) constitutes
a hard gate for bottom culture and restoration reef deployment -- physical
gear deployment is not viable on actively trawled ground regardless of
ecological suitability.

\strong{Input options:}
\enumerate{
\item \code{trawling_data} -- ICES VMS SAR data, manually downloaded as CSV/dataframe.
\item \code{fetch_live = TRUE} -- queries ICES VMS GeoServer endpoint.
\item Neither -- returns zero disturbance with a warning.
}
}
\section{trawling_data format}{

A dataframe with columns:
\itemize{
\item \code{lat}, \code{lon} -- centroid of c-square or fishing location
\item \code{sar} -- swept-area ratio (numeric; dimensionless ratio per year)
\item \code{gear_type} (optional) -- "otter", "beam", "dredge", "seine" etc.
Dredge gear has higher impact per SAR unit.
}
}

\section{Additional disturbance inputs}{

\itemize{
\item \code{anchor_data}: dataframe with \code{lat}, \code{lon}, \code{density} (vessel anchoring
density index) for recreational/commercial anchorage areas.
\item \code{dredge_areas}: dataframe with \code{lat}, \code{lon} identifying licensed
aggregate extraction areas (all sites within the polygon receive score 1.0).
}
}

\examples{
\donttest{
# Manual ICES VMS data (downloaded from ICES data portal)
vms <- read.csv("ices_vms_sar_2022.csv")  # columns: lat, lon, sar, gear_type
result <- score_anthropogenic_disturbance(result, trawling_data = vms)

# Live ICES VMS fetch
result <- score_anthropogenic_disturbance(result, fetch_live = TRUE)
}
}
