AcousticDensity: Calculate number density from NASC in length intervals

View source: R/Density.R

AcousticDensityR Documentation

Calculate number density from NASC in length intervals

Description

This function converts NASC to area number density for each species category based on the acoustic target strength as a function of length for each acoustic category.

Usage

AcousticDensity(
  MeanNASCData,
  AssignmentLengthDistributionData,
  AcousticTargetStrength,
  SpeciesLink = data.table::data.table()
)

Arguments

MeanNASCData

The MeanNASCData data.

AssignmentLengthDistributionData

The AssignmentLengthDistributionData data.

AcousticTargetStrength

The AcousticTargetStrength process data.

SpeciesLink

A table of the two columns AcousticCategory and SpeciesCategory.

Details

The AcousticDensity function calculates length distributed densities as number of fish per square nautical mile by vertical layer. Length based density distributions are calculated for each NASC value of the input acoustic data set (MacLennan et al, 2002) . Usually, the NASC values have a horizontal resolution at PSU level (e.g. transect). By combining a NASC value with a length distribution (usually a total combined length distribution derived from more than one biotic station) and applying a TS vslength relationship, a corresponding density length distribution can be calculated.

The horizontal resolution (EDSU, PSU or Stratum) of the NASC object determines the horizontal resolution of the densities. The same principal applies for the vertical layer resolution.

To combine a NASC value and a length distribution into a density distribution, the TS vs length relationship for the actual species and acoustic frequency, needs to be known.Constants in the TS vs length formula has to be given. If the vertical layer resolution is channels referred from surface and with a fixed thickness, the mean depth of each channel can be calculated and a depth dependent TS vs length formula may be applied. The calculation of densities by length from a NASC value and a corresponding length distribution (as percentage or proportion) is performed as follows (Ona, 2003, Simmond and MacLennan, 2005, Johnsen et al, 2019):

TS_l = m \log_{10}{(l)}+a + d \log_{10}{(1 + \frac{r_y}{10})}

where:

TS_l = target strength (dB re 1 m^{2}) of a fish with length l (cm)

m = constant in the TS vs length relationship for the given species

a = constant in the TS vs length relationship for the given species

d = constant in the TS vs length relationship related to depth dependent TS

l = length of the fish (cm). Typically, the center length of a length group

r_y = average depth (m) of the NASC channel y

\sigma_{bs,l} = 10^{ ( \frac{TS_l}{10} ) }

where:

\sigma_{bs,l} = acoustic backscattering cross-section (m^{2}) for a fish of length l

NASC_l = NASC \frac{\sigma_{bs,l} p_l}{\sum_l{( \sigma_{bs,l} p_l )} }

where:

NASC = the total NASC which is used to calculate densities by length

NASC_l = the proportion of the total NASC which can be attributed to length group l. The sum of NASC_l for all length groups in the total length distribution is equal to NASC

p_l = proportion of fish of length l in the input length distribution. Sum of all p_l is 1.

\rho_l = \frac{NASC_l}{(4 \pi \sigma_{bs,l})}

where:

\rho_l = area density of fish (ind. per sqare nautical mile) in length group l

References:

Johnsen, E,.Totland, A.,Skaalevik, A., et al., 2019, StoX: An open source software for marine survey analyses. Methods Ecol Evol. 2019;10:1523_1528. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1111/2041-210X.13250")}

MacLennan, D. N., Fernandes, P. G., and Dalen, J. 2002. A consistent approach to definitions and symbols in fisheries acoustics. ICES Journal of Marine Science, 59: 365_369.

Ona, E. 2003, An expanded target strength relationship for herring, ICES Journal of Marine Science, Volume 60, Issue 3, 2003, Pages 493_499, \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1016/S1054-3139(03)00031-6")}

Simmonds, J., and MacLennan, D. 2005. Fisheries Acoustics. Theory and Practice, Blackwell Science, Oxford. 437 pp.

Value

An object of StoX data type DensityData. Note that the Density column of the Data table has unit nmi^-2 and DensityType "AreaNumberDensity".

See Also

See SweptAreaDensity for swept-area density. To define the acoustic target strength, see AcousticTargetStrength. The DensityData leads to the QuantityData in a typical survey estimation model.


StoXProject/RstoxBase documentation built on Nov. 6, 2024, noon