Outputs of current speed and sea otter abundance models in Glacier Bay, Alaska

R scripts to read and process both the current speed and sea otter rasters are available at https://github.com/clint-leach/otter-clam and in the linked Zenodo archive. Code to load the current speed raster and extract its values at a set of lat-long coordinates is provided in the script 'align_data.r' provided in that repository. After extracting the zip, the current speed raster can be loaded into R using the 'raster' package. The zipped archive 'otters.zip' contains R object files ('rds' format) that contain a raster defining the Glacier Bay domain ('Boundaries.rds'), covariate inputs to the sea otter ecological diffusion model ('lambda_covars.rds'), and a list containing MCMC chains for the parameters of the ecological diffusion model ('lambda_mcmc.rds'). The R script `compute_lambda.r` uses these inputs to generate a raster brick containing posterior estimates of sea otter abundance across Glacier Bay through time.