Search

Understanding the biotic consequences of Pleistocene range shifts and fragmentation remains a fundamental goal in historical biogeography and evolutionary biology. Here, we combine species distribution models ( SDM) from the present and two late Quaternary time periods with multilocus genetic data (mitochondrial DNA and microsatellites) to evaluate the effect of climate-induced habitat shifts on population genetic structure in the Large-blotched Ensatina ( Ensatina eschscholtzii klauberi), a plethodontid salamander endemic to middle and high-elevation conifer forest in the Transverse and Peninsular Ranges of southern California and northern Baja California. A composite SDM representing the range through time predicts two disjunct refugia, one in southern California encompassing the core of the species range and the other in the Sierra San Pedro Mártir of northern Baja California at the southern limit of the species range. Based on our spatial model, we would expect a pattern of high connectivity among populations within the northern refugium and, conversely, a pattern of isolation due to long-term persistence of the Sierra San Pedro Mártir population. Our genetic results are consistent with these predictions based on the hypothetical refugia in that (i) historical measures of population connectivity among stable areas are correlated with gene flow estimates; and (ii) there is strong geographical structure between separate refugia. These results provide evidence for the role of recent climatic change in shaping patterns of population persistence and connectivity within the Transverse and Peninsular Ranges, an evolutionary hotspot. [ABSTRACT FROM AUTHOR]

We developed interpolated climate surfaces for global land areas (excluding Antarctica) at a spatial resolution of 30 arc s (often referred to as 1‐km spatial resolution). The climate elements considered were monthly precipitation and mean, minimum, and maximum temperature. Input data were gathered from a variety of sources and, where possible, were restricted to records from the 1950–2000 period. We used the thin‐plate smoothing spline algorithm implemented in the ANUSPLIN package for interpolation, using latitude, longitude, and elevation as independent variables. We quantified uncertainty arising from the input data and the interpolation by mapping weather station density, elevation bias in the weather stations, and elevation variation within grid cells and through data partitioning and cross validation. Elevation bias tended to be negative (stations lower than expected) at high latitudes but positive in the tropics. Uncertainty is highest in mountainous and in poorly sampled areas. Data partitioning showed high uncertainty of the surfaces on isolated islands, e.g. in the Pacific. Aggregating the elevation and climate data to 10 arc min resolution showed an enormous variation within grid cells, illustrating the value of high‐resolution surfaces. A comparison with an existing data set at 10 arc min resolution showed overall agreement, but with significant variation in some regions. A comparison with two high‐resolution data sets for the United States also identified areas with large local differences, particularly in mountainous areas. Compared to previous global climatologies, ours has the following advantages: the data are at a higher spatial resolution (400 times greater or more); more weather station records were used; improved elevation data were used; and more information about spatial patterns of uncertainty in the data is available. Owing to the overall low density of available climate stations, our surfaces do not capture of all variation that may occur at a resolution of 1 km, particularly of precipitation in mountainous areas. In future work, such variation might be captured through knowledge‐based methods and inclusion of additional co‐variates, particularly layers obtained through remote sensing. Copyright © 2005 Royal Meteorological Society. [ABSTRACT FROM AUTHOR]