fishstat (version 2025.1.0.1)

production: Global Aquaculture and Capture Production

Description

Aquaculture and capture production quantity by species, area, and country for the years 1950-2023, compiled and published by FAO (2025).

Usage

production

Format

Data frame containing eight columns:

species	species code
year	year
area	area code
country	country code
value	quantity in tonnes or number of individuals
measure	measure code
status	status code
source	source code

Details

This data frame contains the full set of 1,159,613 data records from the FishStat Production Quantity data table. Column names have been simplified to facilitate quick exploration and plotting in R.

See Also

aquaculture and capture data are also available in a combined production format.

area, country, environment, measure, source, species, and status are lookup tables.

fishstat-package gives an overview of the package.

Examples

head(production)

# Analyze production measured in tonnes
prod <- production[production$measure == "Q_tlw" & production$value > 0,]

# Fast merge yearbook and inlandmarine columns
prod$yearbook <- species$yearbook[match(prod$species, species$species)]
prod$inlandmarine <- area$inlandmarine[match(prod$area, area$area)]

# Select SOFIA species, excluding mammals, reptiles, and plants
prod <- prod[grep("Fish, crustaceans and molluscs", prod$yearbook, fixed=TRUE),]

# Determine origin
prod$origin <- ifelse(prod$source == "CAPTURE", "Capture", "Aquaculture")
prod$w <- ifelse(prod$inlandmarine == "Marine areas", "marine", "inland")
prod$origin <- paste0(prod$origin, " (", prod$w, ")")
cbind(sort(unique(prod$origin)))

# World capture fisheries and aquaculture production
x <- xtabs(value~year+origin, aggregate(value~year+origin, prod, sum))
x <- x[,c(2,1,4,3)] / 1e6
library(areaplot)
areaplot(x, legend=TRUE, args.legend=list(x="topleft"), ylab="million tonnes")