Modeling monthly mean temperatures for the mountain regions of Taiwan by generalized additive models

Abstract Taiwan is a mountainous island with a mountain flora that has a high degree of endemism. In order to facilitate the assessment of Taiwan’s mountain flora integrity under various climate change scenarios, the objective of this study was to construct the monthly mean temperature models for the island’s mountain regions (elevation ≥ 1000 m a.s.l.) based on data from 43 meteorology stations. We used a generalized additive modeling approach to construct the models, with elevation, easting, and northing co-ordinates, as the predictors. The effects of elevation were modeled mainly by linear or second-order polynomial functions, whereas the effects of location were modeled nonparametrically by thin plate regression splines. We used the leave-one-out cross-validation to validate the models, with mean errors, mean absolute errors, and root-mean-square of normalized errors as measures of generalization reliability. Spatial interpolations were then performed at a 1-km resolution. Judging by the adjusted R 2 and the proportion of deviance explained, all fitted models performed well. Monthly mean temperatures during the period of April to August decreased linearly with increasing altitude. The observed lapse rates were between 4.9 and 5.6 °C/km, with higher rates during the late spring and summer, and lower rates in the fall and early winter. For the first 3 calendar months, the lapse rates were between 3.2 and 3.6 °C/km, but rose by 0.4–0.5 °C/km for each km of increase in elevation. The three performance measures suggested that the monthly predictions were reliable. The predictions were more reliable for regions from which most of the data were collected and less reliable at the edges of the mountain regions. Model analyses and spatial interpolation results suggested that the models realistically reflected the influence of the geographic location, terrain characteristics, dominant regional climatic features, and their interactions. Predictions also confirmed the existence of the Massenerhebung effect in central Taiwan. Our models also suggest that the alpine florae of the central and the southern regions could be more vulnerable to future winter temperature increases than that of the northern region. Although this study suffered from the lack of data from remote areas, a relatively dense network with quality stations and a good representation along the altitudinal gradient enabled this study to develop reliable models over a complex terrain. Future model improvements can be achieved by placing temporary monitoring points in less represented areas and then spatially interpolating between the temporary stations and the permanent ones.

[1]  Peter Craven,et al.  Smoothing noisy data with spline functions , 1978 .

[2]  Xiaosheng Yang,et al.  Case study: the potential of bamboo resources in mountainous China , 2000 .

[3]  Robert L. Vislocky,et al.  Generalized Additive Models versus Linear Regression in Generating Probabilistic MOS Forecasts of Aviation Weather Parameters , 1995 .

[4]  R. Macarthur Mathematical Ecology and Its Place among the Sciences. (Book Reviews: Geographical Ecology. Patterns in the Distribution of Species) , 1974 .

[5]  Martin F. Price,et al.  Forests in sustainable mountain development : a state of knowledge report for 2000 , 2000 .

[6]  M. Hutchinson,et al.  A comparison of two statistical methods for spatial interpolation of Canadian monthly mean climate data , 2000 .

[7]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[8]  M. G. Ryan,et al.  Tree and forest functioning in response to global warming. , 2001, The New phytologist.

[9]  Roger G. Barry,et al.  Mountain weather and climate , 1982 .

[10]  G. Wahba Spline models for observational data , 1990 .

[11]  J. Grieser,et al.  Statistical time series decomposition into significant components and application to European temperature , 2002 .

[12]  J. Haslett Mountain Ecology: Organism Responses to Environmental Change, an Introduction , 1997 .

[13]  Dean L. Urban,et al.  Spatial estimation of air temperature differences for landscape-scale studies in montane environments , 2003 .

[14]  H. Hsu,et al.  Observed and projected climate change in Taiwan , 2002 .

[15]  O. Phillips,et al.  Extinction risk from climate change , 2004, Nature.

[16]  John H. Campbell,et al.  Biological response to climate change on a tropical mountain , 1999, Nature.

[17]  C. Rolland Spatial and Seasonal Variations of Air Temperature Lapse Rates in Alpine Regions , 2003 .

[18]  Kirsten M. de Beurs,et al.  Kriging and thin plate splines for mapping climate variables , 2001 .

[19]  N. Stuart,et al.  A comparison among strategies for interpolating maximum and minimum daily air temperatures. Part II: The interaction between number of guiding variables and the type of interpolation method , 2001 .

[20]  A. Bloom,et al.  Snowline Altitude and Climate in the Peruvian Andes (5-17° S) at Present and during the Latest Pleistocene Glacial Maximum , 1994 .

[21]  S. Wood Thin plate regression splines , 2003 .

[22]  S. Halloy,et al.  Climate-Change Effects on Alpine Plant Biodiversity: A New Zealand Perspective on Quantifying the Threat , 2003 .

[23]  R. Tibshirani,et al.  Generalized additive models for medical research , 1986, Statistical methods in medical research.

[24]  R. L. Eubank,et al.  Diagnostics for Smoothing Splines , 1985 .

[25]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[26]  S. Wood Generalized Additive Models: An Introduction with R , 2006 .

[27]  Michael F. Hutchinson,et al.  Interpolation of Rainfall Data with Thin Plate Smoothing Splines - Part II: Analysis of Topographic Dependence , 2002 .

[28]  M. Hutchinson,et al.  The development of 1901–2000 historical monthly climate models for Canada and the United States , 2006 .

[29]  Xiaogu Zheng,et al.  Thin plate smoothing spline interpolation of daily rainfall for New Zealand using a climatological rainfall surface , 2006 .

[30]  Robert Tibshirani,et al.  The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd Edition , 2001, Springer Series in Statistics.

[31]  J. Sansom,et al.  Estimation of Long-Term Climate Information at Locations with Short-Term Data Records , 2004 .

[32]  M. Stone Cross‐Validatory Choice and Assessment of Statistical Predictions , 1976 .

[33]  Michael F. Hutchinson,et al.  Interpolation of Rainfall Data with Thin Plate Smoothing Splines - Part I: Two Dimensional Smoothing of Data with Short Range Correlation , 1998 .