Work Package 1: Empirical estimates of fish diversity patterns. I collected occurrence data for all the fishes of interest (434 species) from online databases and the literature. Next, I conceived of a novel, 3-D approach to infer species distributions from occurrences and bundled the resulting software tools into a new R package, voluModel. Species distribution modelling and related analytical methods were first developed in two-dimensional terrestrial systems; many common workflows organize and analyze geographically structured occurrence and environmental data based on horizontal latitude and longitude coordinates. However, marine species also distributed vertically, and water conditions may vary strikingly with depth. voluModel's innovative tools facilitate generation of fine-grained, 3-D distribution models to infer potential species distributions. I then modeled species distributions based on environmental conditions where each species was observed, using both standard methods to model horizontal distributions, and my new method to model 3-D distributions.
Work package 2: Process-based models of marine biodiversity patterns. WP2 has mostly consisted of gathering and processing mapped past and present environmental data for ecological variables. According to theory, response to both oceanographic conditions and how stable these conditions were may have driven the emergence of marine biodiversity patterns. I downloaded decadal mean annual average ocean salinity, temperature, and ideal water age (i.e. time since surface contact) between 22,000 and 400 BCE from climate simulations run using two different climate models and I am converting these data into 100-year climatologies at 1,000 year time steps over the last 20,000 years. I will then calculate Holocene and Late Pleistocene temperature, salinity, and water age stability. I have also obtained modern oceanographic data on water density and nutrients.
Work package 3: Describe effects of various drivers on marine biodiversity patterns. Upon completion of species distribution models (WP1), I estimated Atlantic Ocean species richness for each of the three orders of interest by stacking both 2-D and 3-D species distributions. While 2-D species distributions recreated biodiversity patterns like those of other studies using 2-D distributions, patterns of biodiversity based on 3-D models differed substantially, especially for Scombriformes and Gadiformes, which generally dwell below the surface. Once explanatory variable data processing is complete (WP2), these data will be used to examine which variables are most strongly correlated with observed biodiversity patterns, both in horizontal space and with depth, using a machine learning workflow I am currently developing.
Preliminary results of WP1 were disseminated as a live virtual oral presentation at the Joint Meeting of Ichthyologists and Herpetologists in July 2021. Preliminary results of WP1 and WP3 were presented virtually at the International Biogeography Society (IBS) meeting in Vancouver in June 2021 and received an award for best on-demand talk of the conference. The first MARDIGRA project manuscript, describing the voluModel R software package I created in WP1, is currently in review. The software is publicly archived and indexed in OpenAire. Two additional publications are planned. “The influence of natural history on biodiversity patterns: Lessons from fishes” is in preparation and will present empirical biodiversity patterns uncovered using the 3-D distribution modelling workflow. The final manuscript, presenting WP3 results, “Biodiversity gradients in marine systems: Uniting theory and empirical observation to understand processes behind patterns”, will use explanatory variables generated by WP2 to model WP1 observed empirical biodiversity patterns.