A standard Water and Leith climate diagram for the mean values of precipitation and temperature extracted from the species presence points.
The Walter and Leith diagram assumes that the growing season occurs when rainfall is over 100mm. A more refined method is to extract the values from a bucket model that keeps track of input to the soil profile through precipitation and reductions in soil moisture through evaptranspiration over the course of the year. This can be compared to changes in NDVI at the collection points.
In some cases NDVI will remain fairly constant, even when the balance model shows that soil water constant is lowered for part of the year. Providing SWC is above 50% of maximum levels the vegetation would not experience a great deal of hydric stress.
| Min. | 1st Qu. | Median | Mean | 3rd Qu. | Max. | |
|---|---|---|---|---|---|---|
| x | -110.00 | -99.80 | -96.70 | -95.80 | -90.00 | -84.20 |
| y | 11.90 | 17.60 | 19.40 | 19.20 | 20.60 | 29.20 |
| X | -110.00 | -99.80 | -96.70 | -95.80 | -90.00 | -84.30 |
| Y | 11.90 | 17.60 | 19.40 | 19.20 | 20.60 | 29.20 |
| elev | 0.00 | 128.00 | 909.00 | 911.00 | 1500.00 | 3160.00 |
| anprec | 398.00 | 912.00 | 1140.00 | 1210.00 | 1400.00 | 3580.00 |
| mtemp | 9.35 | 19.20 | 22.30 | 22.00 | 25.40 | 29.10 |
| Trange | 11.00 | 17.00 | 18.50 | 19.50 | 21.30 | 32.40 |
| saet | 399.00 | 862.00 | 994.00 | 984.00 | 1120.00 | 1490.00 |
| msoil | 215.00 | 262.00 | 295.00 | 304.00 | 332.00 | 477.00 |
The following diagrams show kernel densities one two synthestic climate axes (total annual rainfall and mean annual temperature). If their are signs of multimodality this may indicate that the species has not fully explored its climate niche, or that there are disjunct populations with differing characteristics. The method will not show clear results for species with few collection points.
## Warning: weights overwritten by binning
The same analysis can be run to look at spatial clustering. The kernel densities are smoothed, so will only suggest multimodality if the points are very highly clustered.
The significance of any spatial clusters can be checked using the silhouette width method. The width is calculated for values of k between 2 and 5. If any are higher than 0.52 the analysis will produce a diagram showing the clusters.
## [1] "The collections form 5 spatial clusters"
The collections form 5 spatial clusters
If there is evidence that the points fall into at least 2 groups, but fewer than 6 we can look at whether there is significant differences in variability between and within groups in the climatic conditions at the sites using Anosim. This is a sensitive test, as would be MANOVA, so there will often be significant differences. They should only be intepreted as important if R is much larger than 0.3.
##
## Call:
## anosim(dat = dis, grouping = fit$cluster, permutations = 100)
## Dissimilarity: euclidean
##
## ANOSIM statistic R: 0.504
## Significance: 0.0099
##
## Based on 100 permutations
##
## Upper quantiles of permutations (null model):
## 90% 95% 97.5% 99%
## 0.00831 0.01002 0.01399 0.01576
##
## Dissimilarity ranks between and within classes:
## 0% 25% 50% 75% 100% N
## Between 7078 456416 771498 1071292 1339065 990983
## 1 726 67340 139556 303078 1296571 4465
## 2 726 105078 302172 643517 1337727 32131
## 3 726 192250 495050 821956 1338263 207690
## 4 726 140682 335234 593799 1314604 43071
## 5 726 24572 56935 113881 792996 60726
## Analysis of Variance Table
##
## Response: vars$mtemp
## Df Sum Sq Mean Sq F value Pr(>F)
## as.factor(fit$cluster) 4 6821 1705 219 <2e-16 ***
## Residuals 1632 12701 8
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
The collections form 5 spatial clusters
For comparison we can fit a model using mean temperature, temperature range and annual precipitation.
Now fit a gam using temperature range, mean temperature and the annual soil moisture dynamic as input.
