Morphological differences are predicted by geographic location for Maastrichtian Gunnarites from the James Ross Basin, Antarctica
Data files
Apr 01, 2026 version files 465.96 KB
-
Location_Pairwise_Effects.csv
16.24 KB
-
Mohr_et_al_Gunnarites_Supplementary_Data.xlsx
238.42 KB
-
Morphometric_Data.csv
184.16 KB
-
README.md
13.81 KB
-
Specimen_Information.csv
13.32 KB
Abstract
The early Maastrichtian kossmaticeratid ammonite Gunnarites serves as an important biostratigraphic index fossil for correlating between localities within the James Ross Basin (JRB), Antarctica. The JRB is the highest southern paleolatitude marine sedimentary outcrop record of the Upper Cretaceous. Gunnarites is abundant within a relatively narrow stratigraphic interval in the JRB, and its presence in outcrops is used to recognize this specific interval, which is coeval with environmental disturbances including cooling temperatures and local sea level regression. Early taxonomic work on the genus (e.g., Spath, 1953) qualitatively described morphological differences between Gunnarites specimens from different locations within the JRB; if real, these location-based differences could have important implications for the utility of Gunnarites for biostratigraphy in the JRB and could contribute to our understanding of controls (e.g., environmental or temporal) on ammonite morphology.
To test whether JRB Gunnarites specimens exhibit location-based morphological differences, we collected morphometric measurements (conch morphology and ribbing parameters) from multiple ontogenetic stages on 118 specimens from seven distinct localities within the JRB. We used linear mixed models and generalized additive mixed models to quantitatively evaluate mean differences in morphological variables by location, while accounting for variation in specimen size (i.e., characterizing scaling relationships between size and shape). Every morphological variable except whorl radius expansion rate exhibited significant location-based differences, although patterns of variation across sites differed by morphological variable. These location effects on Gunnarites morphology may be a result of phenotypic variation along an environmental gradient (e.g., depth), a result of temporal effects (e.g., evolution), or a combination of both. Large magnitude morphological differences between nearby locations may also be due, in part, to the influence of stratigraphic repetitions produced by previously hypothesized structural features in the JRB such as faults or folds, which could disrupt spatial or temporal gradients in the basin. Untangling the exact mechanisms behind location-based morphological differences in Gunnarites must await more precise, independent, age-dating of each outcrop.
Dataset DOI: 10.5061/dryad.d7wm37qfx
Description of the data and file structure
Supplemental Information for:
"Morphological Differences are Predicted by Geographic Location for Maastrichtian Gunnarites from the James Ross Basin, Antarctica"
Rachel C. Mohr
Emily M. Tompkins
Eduardo B. Olivero
Caroline G. Scott
Stephanie J. White
Thomas S. Tobin
Files and variables
File: "Mohr_et_al_Gunnarites_Supplementary_Data.xlsx"
Contents:
Sheet 1: "Metadata": overview of contents, metadata about all variables
Sheet 2: "Specimen Information": basic information about n=118 Gunnarites specimens (accession numbers, previous taxonomic identifications, locality information, etc.). This sheet is also provided in an alternative csv format: "Specimen_Information.csv"
Sheet 3: "Morphometric Data": morphometric data collected from n=118 Gunnarites specimens. This sheet is also provided in an alternative csv format: "Morphometric_Data.csv"
Sheet 4: "Location Pairwise Effects": coefficient estimates and associated p-values for pairwise effects of location from our "Location Models" (linear mixed models including both size and specimen location as fixed effects). This sheet is also provided in an alternative csv format: "Location_Pairwise_Effects.csv"
See Sheet 1 ("Metadata") for descriptions of all variables (column names), including units (if applicable), and explanations of NA values for each variable. Descriptions of all variables (column names) are also provided below:
Variables / Columns:
Specimen Info
specimen_ID: unique identifier for each Gunnarites specimen
shell: unique identifier for each Gunnarites specimen (same as 'specimen_ID')
collection: institution where Gunnarites specimen is reposited (ALMNH: Alabama Museum of Natural History [at time of publishing, now known as "The University of Alabama Museums Research & Collections" Paleontology Collection]; NHM: Natural History Museum, London; CADIC: Centro Austral de Investigaciones Científicas (CADIC-CONICET); MACN: Museo Argentino de Ciencias Naturales “Bernardino Rivadavia”)
catalog_number: official accession number for each specimen at the institution where it is reposited.
publication: indicates the publication in which the specimen has previously been figured, if applicable. ("Olivero1992": Olivero, E. B. 1992: Asociaciones de Amonites de la Formacion Santa Marta (Cretacico Tardio), Isla James Ross, Antartida. Pp.47–75 in Geologia de la Isla James Ross. Instituto Antártico Argentino, Buenos Aires; "Spath1953": Spath, L. F. 1953: The Upper Cretaceous cephalopod fauna of Graham Land. Falkland Islands Dependencies Survey, London, 1–60 p.; "Tobinetal2018": Tobin, T. S., D. Flannery, and F. J. Sousa. 2018: Stratigraphy, sedimentology and paleontology of Upper Cretaceous deposits of Day Nunatak, Snow Hill Island, Antarctica. Cretaceous Research 84:407–419.; "NA": specimen has not previously been figured in published literature, to our knowledge)
pub_fig: indicates the plate and/or figure number in which the specimen is depicted (in the publication listed under "publication" column). "NA" indicates that the specimen has not previously been figured in the literature, to our knowledge.
genus: genus identification (as indicated in publication where specimen was previously figured, if applicable; otherwise, identified by the authors of this study)
species: species identification, as reported by Spath (1953); "unknown" is listed for all specimens not figured in Spath (1953), for which species identification was not attempted.
species_variety: subspecific variety identification, as reported by Spath (1953); "unknown" is listed for all specimens not figured in Spath (1953), for which subspecific identification was not attempted.
species_author: taxonomic author citation for the species name listed in the "species" column, as reported by Spath (1953)
variety_author: taxonomic author citation for the subspecific variety name listed in the "species_variety" column, as reported by Spath (1953)
loc: location in the James Ross Basin, Antarctica where the specimen was collected. See Manuscript Fig. 1 for map.
island: island in the James Ross Basin, Antarctica where the specimen was collected from. See "loc" column for more precise locality information.
elev: relative elevation information (in meters) from direct on-site measurements during the collection of the specimen. "NA" if no stratigraphic information is available for a specimen.
Morphometrics
unique_pos: unique code for each measured position on an individual Gunnarites specimen
angle: distance (in degrees) between a given 'unique_pos' and the most adapertural (latest ontogeny) measured position on a specimen
WH: whorl height (mm) measured at "unique_pos"; if both ImageJ ('WH_ij') and caliper ('WH_cal') measurements were collected at a given 'unique_pos', the ImageJ measurement is reported here for 'WH'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
WH_opp: whorl height (mm) measured at the position 180 degrees prior (ontogenetically) to the given 'unique_pos' on the specimen. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
WW: whorl width (mm) measured at "unique_pos"; if both ImageJ ('WW_ij') and caliper ('WW_cal') measurements were collected at a given 'unique_pos', the ImageJ measurement is reported here for 'WW'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
R: radius (mm) measured at "unique_pos". "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
R_opp: radius (mm) measured at the position 180 degrees prior (ontogenetically) to the given 'unique_pos' on the specimen. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
D: diameter (mm), measured adapically from the 'unique_pos'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
U: umbilical width (mm), measured adapically from the 'unique_pos'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
X: the straight distance (mm) between both ends of the rib, measured for a single rib nearest the 'unique_pos'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
Y: the curved length (mm) of the rib; measured on the same rib as 'X' for a given 'unique_pos'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
Z: the angle (in degrees) between both ends of the rib, with the vertex placed at the estimated position of the protoconch for a specimen; value is positive when the rib is proiradiate and negative when the rib is rursiradiate; measured on the same rib as 'X' and 'Y' for a given 'unique_pos'. "NA" if Z could not be calculated due to preservation of the specimen.
OR: a count of the number of outer ribs (at the ventral margin of the specimen), in the 45 degrees of whorl adapical to the position of 'unique_pos'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
IR: a count of the number of inner ribs (at the umbilical margin of the specimen), in the 45 degrees of whorl adapical to the position of 'unique_pos'. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
WH_ij: whorl height (mm) collected at position 'unique_pos' using ImageJ. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
WH_cal: whorl height (mm) collected at position 'unique_pos' using digital calipers on the physical specimen; "NA" if a measurement could not be taken at 'unique_pos', either due to preservation of the specimen or due to lack of physical access to the specimen
WW_ij: whorl width (mm) collected at position 'unique_pos' using ImageJ. "NA" if a measurement could not be taken at 'unique_pos', due to preservation of the specimen.
WW_cal: whorl width (mm) collected at position 'unique_pos' using digital calipers on the physical specimen; "NA" if a measurement could not be taken at 'unique_pos', either due to preservation of the specimen or due to lack of physical access to the specimen
Calculated Variables
WWI: whorl width index (unitless). "NA" if 'WWI' could not be calculated due to an NA value for at least one of the required measurements, 'WW' and 'WH' at 'unique_pos'. See manuscript Fig. 3 for variable equation.
UWI: UWI, umbilical width index (unitless). "NA" if 'UWI' could not be calculated due to an NA value for at least one of the required measurements, 'U' and 'D' at 'unique_pos'. See manuscript Fig. 3 for variable equation.
RUWI: RUWI, radial umbilical width index (unitless). "NA" if 'RUWI' could not be calculated due to an NA value for at least one of the required measurements, 'WH' and 'R' at 'unique_pos'. See manuscript Fig. 3 for variable equation.
WHER: whorl height expansion rate (unitless). "NA" if 'WHER' could not be calculated due to an NA value for at least one of the required measurements, 'WH' and 'WH_opp' at 'unique_pos". See manuscript Fig. 3 for variable equation.
WRER: whorl radius expansion rate (unitless). "NA" if 'WRER' could not be calculated due to an NA value for at least one of the required measurements, 'R' and 'R_opp' at 'unique_pos'. See manuscript Fig. 3 for variable equation.
rib_sin: Ribsin, rib sinuosity (unitless); for the same rib measured by 'X' and 'Y'. "NA" if rib_sin could not be calculated due to lack of 'X' or 'Y' measurement at 'unique_pos'. See manuscript Fig. 3 for variable equation.
Log-Transformed Variables
log_WH: log-transformed WH. "NA" if no observation for 'WH' at 'unique_pos'.
log_D: log-transformed D. "NA" if no observation for 'D' at 'unique_pos'.
log_WWI: log-transformed WWI. "NA" if no observation for 'WWI' at 'unique_pos'.
log_UWI: log-transformed UWI. "NA" if no observation for 'UWI' at 'unique_pos'.
log_RUWI: log-transformed RUWI. "NA" if no observation for 'RUWI' at 'unique_pos'.
log_WHER: log-transformed WHER. "NA" if no observation for 'WHER' at 'unique_pos'.
log_WRER: log-transformed WRER. "NA" if no observation for 'WRER' at 'unique_pos'.
Models Information
response_var: the response variable of the model, for results in any given row of sheet "Location Pairwise Effects"
model_type: the type of model used to generate results in any given row of sheet "Location Pairwise Effects"; "single-slope Location LMM": Location Model fit for response variables best fit by a single-slope scaling relationship, "threshold Location LMM": Location Model fit for response variables best fit by a threshold scaling relationship
baseline: for each pairwise comparison of two locations, "baseline" indicates the location set as the base level for the multi-level factor of location ("loc") in the model
loc_effect: for each pairwise comparison of two locations, "loc_effect" indicates the location being compared to the location set as the model's "baseline"
loc_effect_coefficient: coefficient estimate for "loc_effect", expressing the difference in the mean value of "response_var" relative to "baseline"
p.value: the p-value associated with "loc_effect_coefficient"
adjustedP: the p-value associated with "loc_effect_coefficient", after an adjustment for multiple comparisons using the Bonferroni correction method
adjP_sig: indicates if "adjustedP" is a value <0.05; "not_sig": adjustedP is greater than or equal to 0.05, "SIG": adjustedP is less than 0.05
Uploaded to Zenodo
File: "Mohr_et_al Gunnarites Supplement_Text.pdf"
Contents:
Supplement Text
Supplementary Tables 1-8
Supplementary Figures 1-3
Supplementary References
Code/software
Contents of "MOHR_et_al Gunnarites R code scripts" ZIP folder:
This folder includes:
1 R script file ("00_Gunnarites-Manuscript-Code.R") used to analyze the data for the manuscript and produce the figures
1 "renv.lock" file to manage versions of packages used by the R script (renv.lock file is produced by renv package). R version 4.3.2 (2023-10-31) was used for all analyses.
1 csv data file ("SUPP_DATA_morphometric.csv")
"SUPP_DATA_morphometric.csv" file:
Sheet 3 ("Morphometric Data") from "MOHR_et_al Gunnarites Supplementary Data.xlsx", saved as a .csv file.
Data for use in associated R code scripts (below)
"renv.lock" file:
produced by renv package in R; used to maintain a reproducible R environment
"00_Gunnarites-Manuscript-Code.R" script file includes the following 9 sub-sections:
Script 01: Load and Format Data
loads and formats data
Script 02: Histograms
produces main text Figure 5
Script 03: Range Plot
produces main text Figure 4
Script 04: Correlation Matrix
produces Supplementary Figure 3
Script 05: Naze Scaling Models
produces Supplementary Figures 1-2
produces Supplementary Tables 3-4
Script 06: Naze Elevation Models
produces Supplementary Tables 1-2
Script 07: Location LMMs
produces Supplementary Tables 5-6
produces "Location Pairwise Effects" Supplementary Data sheet
Script 08: Location LMMs Plots
produces main text Figure 6
Script 09: GAMMs
produces main text Figures 7-9
produces Supplementary Table 7
Morphometric measurements were collected using ImageJ; some measurements were also collected using digital calipers. See associated manuscript for details.
