azores.seamounts

The goal of {azores.seamounts} is to provide easy access in R to seamount details provided in Table A1 in Morato et al. (2008).

Installation

# install.packages("pak")
pak::pak("patterninstitute/azores.seamounts")

Data

The exported data set seamounts provides easy access to the data in Table A1 in Morato et al. (2008).

library(azores.seamounts)
library(sf)
#> Linking to GEOS 3.13.0, GDAL 3.10.0, PROJ 9.5.0; sf_use_s2() is TRUE

# `seamounts` is a tibble.
# Coordinates longitude and latitude in WGS84
seamounts
#> # A tibble: 461 × 13
#>    seamount_index category longitude latitude depth height basal_area
#>             <dbl> <chr>        <dbl>    <dbl> <dbl>  <dbl>      <dbl>
#>  1              1 Large        -24.5     37.4   821   1252       1291
#>  2              2 Large        -25.9     37.6   265   1050        974
#>  3              3 Large        -29.6     38.0   364   1161        895
#>  4              4 Large        -26.1     38.6  1154   1038        964
#>  5              5 Large        -29.9     39.7   886   1020        796
#>  6              6 Large        -26.9     39.7   756   1125       1159
#>  7              7 Large        -25.7     36.6  1957   1144       1180
#>  8              8 Large        -26.1     37.1  1048   1721       1004
#>  9              9 Large        -28.2     37.2  1646   1141        654
#> 10             10 Large        -31.4     37.2   413   1305       1295
#> # ℹ 451 more rows
#> # ℹ 6 more variables: height_to_radius_ratio <dbl>, avg_slope <dbl>,
#> #   distance_to_nearest_largest_seamount <dbl>,
#> #   nearest_largest_seamount_index <dbl>, distance_to_nearest_seamount <dbl>,
#> #   nearest_seamount_index <dbl>

# `seamounts_sf` is a sf object.
# Coordinates in a custom Lambert Azimuthal Equal-Area (LAEA) projection,
# centered on the Azores.
seamounts_sf
#> Simple feature collection with 461 features and 11 fields
#> Geometry type: POINT
#> Dimension:     XY
#> Bounding box:  xmin: -640040.3 ymin: -522187 xmax: 630274.9 ymax: 493607.8
#> Projected CRS: +proj=laea +lat_0=38.5 +lon_0=-28 +datum=WGS84 +units=m +no_defs
#> # A tibble: 461 × 12
#>    seamount_index category depth height basal_area height_to_radius_ratio
#>  *          <dbl> <chr>    <dbl>  <dbl>      <dbl>                  <dbl>
#>  1              1 Large      821   1252       1291                   0.1 
#>  2              2 Large      265   1050        974                   0.09
#>  3              3 Large      364   1161        895                   0.11
#>  4              4 Large     1154   1038        964                   0.08
#>  5              5 Large      886   1020        796                   0.1 
#>  6              6 Large      756   1125       1159                   0.09
#>  7              7 Large     1957   1144       1180                   0.11
#>  8              8 Large     1048   1721       1004                   0.15
#>  9              9 Large     1646   1141        654                   0.1 
#> 10             10 Large      413   1305       1295                   0.09
#> # ℹ 451 more rows
#> # ℹ 6 more variables: avg_slope <dbl>,
#> #   distance_to_nearest_largest_seamount <dbl>,
#> #   nearest_largest_seamount_index <dbl>, distance_to_nearest_seamount <dbl>,
#> #   nearest_seamount_index <dbl>, geometry <POINT [m]>

Create a SpatRaster with distances to seamounts

Distances to all (large and small) seamounts:

(rast <- seamounts_dist(resolution = 20000, seamount_size = "All"))
#> |---------|---------|---------|---------|=========================================                                          
#> class       : SpatRaster 
#> dimensions  : 58, 69, 1  (nrow, ncol, nlyr)
#> resolution  : 20000, 20000  (x, y)
#> extent      : -7e+05, 680000, -590000, 570000  (xmin, xmax, ymin, ymax)
#> coord. ref. : +proj=laea +lat_0=38.5 +lon_0=-28 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source(s)   : memory
#> name        :       lyr.1 
#> min value   :    567.6739 
#> max value   : 429754.2835

tidyterra::autoplot(rast)

Distances to small seamounts:

(rast <- seamounts_dist(resolution = 20000, seamount_size = "Small"))
#> |---------|---------|---------|---------|=========================================                                          
#> class       : SpatRaster 
#> dimensions  : 58, 69, 1  (nrow, ncol, nlyr)
#> resolution  : 20000, 20000  (x, y)
#> extent      : -7e+05, 680000, -590000, 570000  (xmin, xmax, ymin, ymax)
#> coord. ref. : +proj=laea +lat_0=38.5 +lon_0=-28 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source(s)   : memory
#> name        :       lyr.1 
#> min value   :    567.6739 
#> max value   : 429754.2835

tidyterra::autoplot(rast)

Distances to large seamounts:

(rast <- seamounts_dist(resolution = 20000, seamount_size = "Large"))
#> |---------|---------|---------|---------|=========================================                                          
#> class       : SpatRaster 
#> dimensions  : 58, 69, 1  (nrow, ncol, nlyr)
#> resolution  : 20000, 20000  (x, y)
#> extent      : -7e+05, 680000, -590000, 570000  (xmin, xmax, ymin, ymax)
#> coord. ref. : +proj=laea +lat_0=38.5 +lon_0=-28 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source(s)   : memory
#> name        :     lyr.1 
#> min value   :   1050.98 
#> max value   : 665153.80

tidyterra::autoplot(rast)

Distances to large seamounts at a higher resolution (500 meters):

(rast <- seamounts_dist(resolution = 500, seamount_size = "Large"))
#> |---------|---------|---------|---------|=========================================                                          
#> class       : SpatRaster 
#> dimensions  : 2300, 2770, 1  (nrow, ncol, nlyr)
#> resolution  : 500, 500  (x, y)
#> extent      : -7e+05, 685000, -590000, 560000  (xmin, xmax, ymin, ymax)
#> coord. ref. : +proj=laea +lat_0=38.5 +lon_0=-28 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs 
#> source(s)   : memory
#> name        :        lyr.1 
#> min value   :     15.08156 
#> max value   : 677453.60236

tidyterra::autoplot(rast)
#> <SpatRaster> resampled to 500520 cells.