Published October 23, 2020 | Version v1

Climatic and Vegetational Drivers of Insect Beta Diversity at the Continental Scale

  • 1. Chinese Academy of Sciences
  • 2. Forest Ecology and Ecophysiology
  • 3. University of Guelph
  • 4. Hong Kong Baptist University

Description

Aim

We construct a framework for mapping pattern and drivers of insect diversity at the continental scale, and use it to test whether and which environmental gradients drive insect beta diversity.

Location

Global; North and Central America; Western Europe

Time period

21st century

Major taxa studied

Insects

Methods

An informatics system was developed to integrate terrestrial data on insects with environmental parameters. We mined repositories of data for distribution, climatic data was retrieved (WorldClim), and vegetation parameters inferred from remote sensing analysis (MODIS Vegetation Continuous Fields). Beta diversity between sites was calculated and then modelled with two methods, Mantel test with Multiple Regression and Generalized Dissimilarity Modelling.

Results

Geographic distance was usually the main driver of insect beta diversity. Independent of geographic distance, bioclimate variables explained more variance in dissimilarity than vegetation variables, although the particular variables found to be significant were more consistent in the latter, particularly, tree cover. Tree cover gradients drove compositional dissimilarity at denser coverages, in both continental case studies. For climate, gradients in temperature parameters were significant in driving beta diversity more so than gradients in precipitation parameters.

Main conclusions

Although environmental gradients drive insect beta diversity independently of geography, the relative contribution of different climatic and vegetational parameters is not expected to be consistent in different study systems. With further incorporation of additional temporal information and variables, this approach will enable development of a predictive framework for conserving insect biodiversity at the global-scale.

 

Notes

 

# will need to make equivelent rows using the intersect R command. example:

overlapping_samples <- intersect(row.names(params_MCL), row.names(comm_MCL));length(overlapping_samples)
comm_MCL <- comm_MCL[overlapping_samples,]
params_MCL <- params_MCL[overlapping_samples,]

 

# consider removing species which are found at very few sites:
column_names <- colnames(comm_MCL)
region_counts_per_species <- seq(from=0, to=0, length=length(comm_MCL[1 , ]))
for(j in 1:length(region_counts_per_species))
    {
    current_column <- comm_MCL[ , j ];current_column_ID <- column_names[j]
    indexs <- (1:length(current_column))[current_column >= 1]; # print(c( j , current_column_ID, length(region_counts_per_species) , length(indexs) ))
    if(length(indexs) >= 3){region_counts_per_species[j] <- 1}
    }
rm_cols2 <- (1:length(region_counts_per_species))[region_counts_per_species == 0];print (c("removing species:" , length(rm_cols2) )); print (c( "proportion:" , length(rm_cols2) /  length(comm_MCL[ 1 ,  ])))

 

# and if removing sites with few observations:

species_counts_per_region <- seq(from=0, to=0, length=length(comm_MCL[ , 1]))
for(j in 1:length(species_counts_per_region))
    {
    current_row <- comm_MCL[j , ]
    indexs <- (1:length(current_row))[current_row >= 1]
    if(length(indexs) >= 30)
        {species_counts_per_region[j] <- 1}
    }

rm_rows<-(1:length(species_counts_per_region))[species_counts_per_region == 0];

if(length(rm_rows) >= 1)
    {
    comm_MCL<-comm_MCL[-rm_rows , ];params_MCL <- params_MCL[-rm_rows , ];
    }

 

Files

comm_table_NCAmerica.zip

Files (4.5 MB)

Name Size Download all
md5:d66ce8cdfec1e58dc31c3244e158971e
334.1 kB Preview Download
md5:a7a035bb6ad269f8b7f0e648342de992
1.9 MB Preview Download
md5:a6bae708432da87553632a9b8bc2294b
554.8 kB Download
md5:8ef534b2c4f618652970a54c0e145822
1.7 MB Download