Application of a geographically-weighted regression analysis to estimate net primary production of Chinese forest ecosystems

Aim The objective of this paper is to obtain a net primary production (NPP) regression model based on the geographically weighted regression (GWR) method, which includes spatial non‐stationarity in the parameters estimated for forest ecosystems in China. Location We used data across China. Methods We examine the relationships between NPP of Chinese forest ecosystems and environmental variables, specifically altitude, temperature, precipitation and time‐integrated normalized difference vegetation index (TINDVI) based on the ordinary least squares (OLS) regression, the spatial lag model and GWR methods. Results The GWR method made significantly better predictions of NPP in simulations than did OLS, as indicated both by corrected Akaike Information Criterion (AIC_c) and R². GWR provided a value of 4891 for AIC_c and 0.66 for R², compared with 5036 and 0.58, respectively, by OLS. GWR has the potential to reveal local patterns in the spatial distribution of a parameter, which would be ignored by the OLS approach. Furthermore, OLS may provide a false general relationship between spatially non‐stationary variables. Spatial autocorrelation violates a basic assumption of the OLS method. The spatial lag model with the consideration of spatial autocorrelation had improved performance in the NPP simulation as compared with OLS (5001 for AIC_c and 0.60 for R²), but it was still not as good as that via the GWR method. Moreover, statistically significant positive spatial autocorrelation remained in the NPP residuals with the spatial lag model at small spatial scales, while no positive spatial autocorrelation across spatial scales can be found in the GWR residuals. Conclusions We conclude that the regression analysis for Chinese forest NPP with respect to environmental factors and based alternatively on OLS, the spatial lag model, and GWR methods indicated that there was a significant improvement in model performance of GWR over OLS and the spatial lag model.     

基于不同空间模型的马尾松林生态系统碳密度研究

以江西省马尾松林生态系统为研究对象,基于样地调查及样品碳含量测定结果计算其碳密度,并选取立地、植被及气象等方面的15个因子,采用多元线性逐步回归方法筛选出对生态系统碳密度影响显著的因子,然后分别利用最小二乘模型(OLS)、空间误差模型(SEM)、空间滞后模型(SLM)和地理加权回归模型(GWR)构建生态系统碳密度与其影响因子之间的关系模型,筛选出最优的拟合模型。结果表明:对马尾松林生态系统碳密度影响显著的因子分别为海拔、坡度、土层厚度、胸径、年均温度和年均降水量。4种模型拟合结果均显示碳密度与坡度呈负相关,与海拔、土层厚度、胸径呈正相关。模型的决定系数(R²)由大到小分别为GWR(0.8043)>SEM(0.6371)>SLM(0.6364)>OLS(0.6321),模型均方误差(MSE)与赤池信息准则(AIC)最大的均为OLS模型,最小的均为GWR模型;残差检验表明GWR模型能有效降低模型残差的空间自相关性。综合分析得出GWR模型的拟合效果最优,更适用于江西省马尾松林生态系统碳密度的估测。     

Model selection and information theory in geographical ecology

Aim   Although parameter estimates are not as affected by spatial autocorrelation as Type I errors, the change from classical null hypothesis significance testing to model selection under an information theoretic approach does not completely avoid problems caused by spatial autocorrelation. Here we briefly review the model selection approach based on the Akaike information criterion (AIC) and present a new routine for Spatial Analysis in Macroecology (SAM) software that helps establishing minimum adequate models in the presence of spatial autocorrelation.
Innovation    We illustrate how a model selection approach based on the AIC can be used in geographical data by modelling patterns of mammal species in South America represented in a grid system ( n  = 383) with 2° of resolution, as a function of five environmental explanatory variables, performing an exhaustive search of minimum adequate models considering three regression methods: non-spatial ordinary least squares (OLS), spatial eigenvector mapping and the autoregressive (lagged-response) model. The models selected by spatial methods included a smaller number of explanatory variables than the one selected by OLS, and minimum adequate models contain different explanatory variables, although model averaging revealed a similar rank of explanatory variables.
Main conclusions    We stress that the AIC is sensitive to the presence of spatial autocorrelation, generating unstable and overfitted minimum adequate models to describe macroecological data based on non-spatial OLS regression. Alternative regression techniques provided different minimum adequate models and have different uncertainty levels. Despite this, the averaged model based on Akaike weights generates consistent and robust results across different methods and may be the best approach for understanding of macroecological patterns.     

A comparison of simultaneous autoregressive and generalized least squares models for dealing with spatial autocorrelation

Aim  In their recent paper, Kissling & Carl (2008 ) recommended the spatial error simultaneous autoregressive model (SAR_err) over ordinary least squares (OLS) for modelling species distribution. We compared these models with the generalized least squares model (GLS) and a variant of SAR (SAR_vario). GLS and SAR_vario are superior to standard implementations of SAR because the spatial covariance structure is described by a semivariogram model.
Innovation  We used the complete datasets employed by Kissling & Carl (2008 ), with strong spatial autocorrelation, and two datasets in which the spatial structure was degraded by sample reduction and grid coarsening. GLS performed consistently better than OLS, SAR_err and SAR_vario in all datasets, especially in terms of goodness of fit. SAR_vario was marginally better than SAR_err in the degraded datasets.
Main conclusions  GLS was more reliable than SAR-based models, so its use is recommended when dealing with spatially autocorrelated data.     

The biogeography of invasive alien plants in California: an application of GIS and spatial regression analysis

The spatial distribution of invasive alien plants has been poorly documented in California. However, with the increased availability of GIS software and spatially explicit data, the distribution of invasive alien plants can be explored. Using bioregions as defined in Hickman (1993 ), I compared the distribution of invasive alien plants (n = 78) and noninvasive alien plants (n = 1097). The distribution of both categories of alien plants was similar with the exception of a higher concentration of invasive alien plants in the North Coast bioregion. Spatial autocorrelation analysis using Moran's I indicated significant spatial dependence for both invasive and noninvasive alien plant species. I used both ordinary least squares (OLS) and spatial autoregressive (SAR) models to assess the relationship between alien plant species distribution and native plant species richness, road density, population density, elevation, area of sample unit, and precipitation. The OLS model for invasive alien plants included two significant effects; native plant species richness and elevation. The SAR model for invasive alien plants included three significant effects; elevation, road density, and native plant species richness. The SAR model for noninvasive alien plants resulted in the same significant effects as invasive alien plants. Both invasive and noninvasive alien plants are found in regions with low elevation, high road density, and high native‐plant species richness. This is in congruity with previous spatial pattern studies of alien plant species. However, the similarity in effects for both categories of alien plants alludes to the importance of autecological attributes, such as pollination system, dispersal system and differing responses to disturbance in the distribution of invasive plant species. In addition, this study emphasizes the critical importance of testing for spatial autocorrelation in spatial pattern studies and using SAR models when appropriate.     

Analysis of determinants of mammalian species richness in South America using spatial autoregressive models

Classically, hypotheses concerning the distribution of species have been explored by evaluating the relationship between species richness and environmental variables using ordinary least squares (OLS) regression. However, environmental and ecological data generally show spatial autocorrelation, thus violating the assumption of independently distributed errors. When spatial autocorrelation exists, an alternative is to use autoregressive models that assume spatially autocorrelated errors. We examined the relationship between mammalian species richness in South America and environmental variables, thereby evaluating the relative importance of four competing hypotheses to explain mammalian species richness. Additionally, we compared the results of ordinary least squares (OLS) regression and spatial autoregressive models using Conditional and Simultaneous Autoregressive (CAR and SAR, respectively) models. Variables associated with productivity were the most important at determining mammalian species richness at the scale analyzed. Whereas OLS residuals between species richness and environmental variables were strongly autocorrelated, those from autoregressive models showed less spatial autocorrelation, particularly the SAR model, indicating its suitability for these data. Autoregressive models also fit the data better than the OLS model (increasing R² by 5–14%), and the relative importance of the explanatory variables shifted under CAR and SAR models. These analyses underscore the importance of controlling for spatial autocorrelation in biogeographical studies.     

基于地理加权回归模型探究环境异质性对秦岭大熊猫空间利用的影响

环境异质性对野生动物分布的影响具有明显的空间不均匀性。传统分析中多采用经典线性回归模型来量化野生动物分布与环境变量之间的关系,难以准确反映物种-环境关系的空间异质特征。地理加权回归(GWR)是近年来提出的一种新的空间分析方法,通过将空间结构嵌入线性回归模型中,以此来探测空间关系的非均匀性。以秦岭大熊猫为例,应用GWR模型分析大熊猫空间分布与环境异质性特征之间的潜在关系,并同经典的全局最小二乘回归法(OLS)进行比较。结果表明,GWR模型的AIC、R²和校正R²均显著优于OLS模型,GWR模型的局部回归系数估计能够更加深刻地揭示大熊猫空间分布与环境变量间的复杂空间关系,且GWR模型能够为物种的科学保护提供更加有效的理论支撑。因此,GWR模型可为探究物种-环境关系的空间异质特征提供一种新的方法,在物种栖息地选择与利用研究中具有一定的应用前景。     

Red herrings revisited: spatial autocorrelation and parameter estimation in geographical ecology

There have been numerous claims in the ecological literature that spatial autocorrelation in the residuals of ordinary least squares (OLS) regression models results in shifts in the partial coefficients, which bias the interpretation of factors influencing geographical patterns. We evaluate the validity of these claims using gridded species richness data for the birds of North America, South America, Europe, Africa, the ex‐USSR, and Australia. We used richness in 110×110 km cells and environmental predictor variables to generate OLS and simultaneous autoregressive (SAR) multiple regression models for each region. Spatial correlograms of the residuals from each OLS model were then used to identify the minimum distance between cells necessary to avoid short‐distance residual spatial autocorrelation in each data set. This distance was used to subsample cells to generate spatially independent data. The partial OLS coefficients estimated with the full dataset were then compared to the distributions of coefficients created with the subsamples. We found that OLS coefficients generated from data containing residual spatial autocorrelation were statistically indistinguishable from coefficients generated from the same data sets in which short‐distance spatial autocorrelation was not present in all 22 coefficients tested. Consistent with the statistical literature on this subject, we conclude that coefficients estimated from OLS regression are not seriously affected by the presence of spatial autocorrelation in gridded geographical data. Further, shifts in coefficients that occurred when using SAR tended to be correlated with levels of uncertainty in the OLS coefficients. Thus, shifts in the relative importance of the predictors between OLS and SAR models are expected when small‐scale patterns for these predictors create weaker and more unstable broad‐scale coefficients. Our results indicate both that OLS regression is unbiased and that differences between spatial and nonspatial regression models should be interpreted with an explicit awareness of spatial scale.     

福州市土壤铬含量高光谱预测的GWR模型研究

通过系统分析不同光谱分辨率和光谱变换对土壤铬高光谱预测模型的不确定性影响,筛选出最优的光谱分辨率及光谱变量进行土壤铬含量预测的地理权重回归(GWR)模型构建,利用该模型进行福州市土壤铬含量预测,并将预测结果与普通最小二乘法回归(OLS)结果进行比较分析,探讨GWR模型在土壤铬高光谱预测中的适用性及局限性。结果表明:(1)在10 nm分辨率尺度下,以土壤全铬含量为因变量,反射率的二阶微分和反射率倒数的二阶微分为自变量构建的GWR模型对土壤铬预测的效果最好。GWR模型的R~2和调节R~2分别为0.821和0.716,较OLS模型分别提高了0.529和0.450,而AIC值为720.703,较OLS模型减少了22个单位,残差平方和仅为OLS模型的1/4,说明GWR模型的预测效果较OLS模型有了显著提高。(2)土壤铬预测模型的精度受光谱分辨率影响。对于OLS预测模型来说,3 nm分辨率的模型预测效果最好,而对于GWR预测模型来说,10nm分辨率的模型不仅预测效果最好,其相较于OLS模型的改善作用显著,为土壤铬含量GWR预测的最佳光谱分辨率。(3)光谱的一阶微分变换可以有效增强土壤铬的光谱特征,而其余的光谱变换对土壤铬的光谱特征则未起到增强作用,但可以很好地提高模型的预测效果。(4)研究得出土壤铬GWR模型预测的最佳光谱分辨率为10 nm,为EO-1 Hyperion影像的光谱分辨率,而且随着采样点的增加,GWR模型的预测效果趋于稳定,适合空间异质性大的区域尺度土壤铬预测。故该模型与高光谱影像结合,实现模型从实验室尺度向区域尺度的推广,为格网尺度土壤铬的空间预测提供可能。     

Spatial autocorrelation and red herrings in geographical ecology

Aim Spatial autocorrelation in ecological data can inflate Type I errors in statistical analyses. There has also been a recent claim that spatial autocorrelation generates ‘red herrings’, such that virtually all past analyses are flawed. We consider the origins of this phenomenon, the implications of spatial autocorrelation for macro‐scale patterns of species diversity and set out a clarification of the statistical problems generated by its presence. Location To illustrate the issues involved, we analyse the species richness of the birds of western/central Europe, north Africa and the Middle East. Methods Spatial correlograms for richness and five environmental variables were generated using Moran's I coefficients. Multiple regression, using both ordinary least‐squares (OLS) and generalized least squares (GLS) assuming a spatial structure in the residuals, were used to identify the strongest predictors of richness. Autocorrelation analyses of the residuals obtained after stepwise OLS regression were undertaken, and the ranks of variables in the full OLS and GLS models were compared. Results Bird richness is characterized by a quadratic north–south gradient. Spatial correlograms usually had positive autocorrelation up to c. 1600 km. Including the environmental variables successively in the OLS model reduced spatial autocorrelation in the residuals to non‐detectable levels, indicating that the variables explained all spatial structure in the data. In principle, if residuals are not autocorrelated then OLS is a special case of GLS. However, our comparison between OLS and GLS models including all environmental variables revealed that GLS de‐emphasized predictors with strong autocorrelation and long‐distance clinal structures, giving more importance to variables acting at smaller geographical scales. Conclusion Although spatial autocorrelation should always be investigated, it does not necessarily generate bias. Rather, it can be a useful tool to investigate mechanisms operating on richness at different spatial scales. Claims that analyses that do not take into account spatial autocorrelation are flawed are without foundation.     

Information-Theoretic Model Selection and Model Averaging for Closed-Population Capture-Recapture Studies

Specification of an appropriate model is critical to valid statistical inference. Given the “true model” for the data is unknown, the goal of model selection is to select a plausible approximating model that balances model bias and sampling variance. Model selection based on information criteria such as AIC or its variant AIC_c, or criteria like CAIC, has proven useful in a variety of contexts including the analysis of open-population capture-recapture data. These criteria have not been intensively evaluated for closed-population capture-recapture models, which are integer parameter models used to estimate population size (N), and there is concern that they will not perform well. To address this concern, we evaluated AIC, AIC_c, and CAIC model selection for closed-population capture-recapture models by empirically assessing the quality of inference for the population size parameter N. We found that AIC-, AIC_c-, and CAIC-selected models had smaller relative mean squared errors than randomly selected models, but that confidence interval coverage on N was poor unless unconditional variance estimates (which incorporate model uncertainty) were used to compute confidence intervals. Overall, AIC and AIC_c outperformed CAIC, and are preferred to CAIC for selection among the closed-population capture-recapture models we investigated. A model averaging approach to estimation, using AIC, AIC_c, or CAIC to estimate weights, was also investigated and proved superior to estimation using AIC-, AIC_c-, or CAIC-selected models. Our results suggested that, for model averaging, AIC or AIC_c should be favored over CAIC for estimating weights.     

Incorporating spatial autocorrelation may invert observed patterns

Though still often neglected, spatial autocorrelation can be a serious issue in ecology because the presence of spatial autocorrelation may alter the parameter estimates and error probabilities of linear models. Here I re-analysed data from a previous study on the relationship between plant species richness and environmental correlates in Germany. While there was a positive relationship between native plant species richness and an altitudinal gradient when ignoring the presence of spatial autocorrelation, the use of a spatial simultaneous liner error model revealed a negative relationship. This most dramatic effect where the observed pattern was inverted may be explained by the environmental situation in Germany. There the highest altitudes are in the south and the lowlands in the north that result in some locally or regionally inverted patterns of the large-scale environmental gradients from the equator to the north. This study therefore shows the necessity to consider spatial autocorrelation in spatial analyses.     

Effects of spatial autocorrelation,natal philopatry and phenotypic plasticity on the heritability of laying date

We investigated the effect of spatial autocorrelation on heritability (h²) estimates of laying date and clutch size in a population of great tits Parus major. We found that h² of laying date, but not clutch size, declined significantly with increasing distance between the nestbox of mothers and daughters. This decline was caused by a decreasing effect of spatial autocorrelation in laying date, rather than by the existence of genotype–environment interactions (GEI). After correcting for the effect of spatial autocorrelation, h² of laying date was low (0.16 ± 0.07), but significant, and surprisingly consistent with increasing distance between parental and offspring environments. The h² of clutch size was not much affected by spatial autocorrelation. Most previously published estimates of the heritability of laying date include various degrees of common environment effects, which can bias estimates both upwards and downwards. We suggest that using techniques that take spatial autocorrelation into account might be a fruitful approach to estimate h² of traits that show a high degree of plasticity.     

Accounting for multi‐scale spatial autocorrelation improves performance of invasive species distribution modelling (iSDM)

Aim Analyses of species distributions are complicated by various origins of spatial autocorrelation (SAC) in biogeographical data. SAC may be particularly important for invasive species distribution models (iSDMs) because biological invasions are strongly influenced by dispersal and colonization processes that typically create highly structured distribution patterns. We examined the efficacy of using a multi‐scale framework to account for different origins of SAC, and compared non‐spatial models with models that accounted for SAC at multiple levels. Location We modelled the spatial distribution of an invasive forest pathogen, Phytophthora ramorum, in western USA. Methods We applied one conventional statistical method (generalized linear model, GLM) and one nonparametric technique (maximum entropy, Maxent) to a large dataset on P. ramorum occurrence (n = 3787) to develop four types of model that included environmental variables and that either ignored spatial context or incorporated it at a broad scale using trend surface analysis, a local scale using autocovariates, or multiple scales using spatial eigenvector mapping. We evaluated model accuracies and amounts of explained spatial structure, and examined the changes in predictive power of the environmental and spatial variables. Results Accounting for different scales of SAC significantly enhanced the predictive capability of iSDMs. Dramatic improvements were observed when fine‐scale SAC was included, suggesting that local range‐confining processes are important in P. ramorum spread. The importance of environmental variables was relatively consistent across all models, but the explanatory power decreased in spatial models for factors with strong spatial structure. While accounting for SAC reduced the amount of residual autocorrelation for GLM but not for Maxent, it still improved the performance of both approaches, supporting our hypothesis that dispersal and colonization processes are important factors to consider in distribution models of biological invasions. Main conclusions Spatial autocorrelation has become a paradigm in biogeography and ecological modelling. In addition to avoiding the violation of statistical assumptions, accounting for spatial patterns at multiple scales can enhance our understanding of dynamic processes that explain ecological mechanisms of invasion and improve the predictive performance of static iSDMs.     

基于带宽优选地理加权回归模型的深圳市植被碳储量反演

植被碳储量估测是自然资源监测的重要内容,遥感技术结合地面样地进行反演可以获得区域范围内植被碳储量的空间连续分布,弥补了传统人工抽样调查估测的不足。然而,现有的参数和非参数遥感估测模型大多忽略了样地数据的变异与空间自相关关系。研究以Landsat 8 OLI影像为数据源提取遥感变量,结合植被碳储量实测调查数据,利用最小信息准则(AICc)、最大空间自相关距离(MSAD)和交叉验证(CV)分别确定最优带宽,组合Gaussian、Bi-square和Exponential核函数构建地理加权回归(GWR)模型估算深圳市植被碳储量,并与多元线性回归(MLR)进行比较,选择最优模型绘制深圳市植被碳储量空间分布图。研究结果表明,GWR模型整体精度优于MLR模型,GWR模型的决定系数(R~2)均高于MLR模型,且均方根误差(RMSE)和平均绝对误差(MAE)显著降低。带宽和核函数的选择对GWR模型估测结果产生了显著影响。以CV确定带宽、Exponential为核函数组合构建的GWR模型效果最佳,其R~2为0.697,RMSE为10.437 Mg C/hm~2,相比其它模型精度上升了13.87%—32....     

Biological and molecular comparison between localized and systemic acquired resistance induced in tobacco by a Phytophthora megasperma glycoprotein elicitin

We have compared localized (LAR) and systemic (SAR) acquired resistance induced in tobacco by a hypersensitive response (HR) inducing Phytophthora megasperma glycoprotein elicitin. Three different zones were taken into account: LAR, SAR_T and SAR_S. The LAR zone was 5–10 mm wide and surrounded the HR lesion. SAR_T was the tissue of the elicitor-treated leaf immediately beyond the LAR zone. The systemic leaf was called SAR_S. Glycoprotein-treated plants showed enhanced resistance to challenge infection by tobacco mosaic virus (TMV). Disease resistance was similar in SAR_T and SAR_S, and higher in LAR. The expression pattern, in glycoprotein-treated plants, of acidic and basic PR1, PR2, PR3 and PR5 proteins and of O-methyltransferases (OMT), enzymes of the phenylpropanoid pathway, was similar to that in TMV-infected plants. OMT was stimulated in LAR but not in SAR_T and SAR_S. The four classes of acidic and basic PR proteins accumulated strongly in LAR. Reduced amounts of acidic PR1, PR2, PR3 and only minute amounts of basic PR2 and PR3 accumulated in SAR_T and SAR_S. In glycoprotein-treated plants, expression of the acidic and basic PR proteins in LAR and SAR of transgenic NahG and ETR tobacco plants and in LAR of plants treated with inhibitors of salicylic acid accumulation and of ethylene biosynthesis indicated a salicylic acid-dependent signalling pathway for acidic isoform activation and an ethylene-dependent signalling pathway for basic isoform activation.     

Local–regional boundary shifts in oribatid mite (Acari: Oribatida) communities: species–area relationships in arboreal habitat islands of a coastal temperate rain forest, Vancouver Island, Canada

Aim This study investigates the species–area relationship (SAR) for oribatid mite communities of isolated suspended soil habitats, and compares the shape and slope of the SAR with a nested data set collected over three spatial scales (core, patch and tree level). We investigate whether scale dependence is exhibited in the nested sampling design, use multivariate regression models to elucidate factors affecting richness and abundance patterns, and ask whether the community composition of oribatid mites changes in suspended soil patches of different sizes. Location Walbran Valley, Vancouver Island, Canada. Methods A total of 216 core samples were collected from 72 small, medium and large isolated suspended soil habitats in six western redcedar trees in June 2005. The relationship between oribatid species richness and habitat volume was modelled for suspended soil habitat isolates (type 3) and a nested sampling design (type 1) over multiple spatial scales. Nonlinear estimation parameterized linear, power and Weibull function regression models for both SAR designs, and these were assessed for best fit using R² and Akaike's information criteria (ΔAIC) values. Factors affecting oribatid mite species richness and standardized abundance (number per g dry weight) were analysed by anova and linear regression models. Results Sixty‐seven species of oribatid mites were identified from 9064 adult specimens. Surface area and moisture content of suspended soils contributed to the variation in species richness, while overall oribatid mite abundance was explained by moisture and depth. A power‐law function best described the isolate SAR (S = 3.97 × A^0.12, R² = 0.247, F_1,70 = 22.450, P < 0.001), although linear and Weibull functions were also valid models. Oribatid mite species richness in nested samples closely fitted a power‐law model (S = 1.96 × A^0.39, R² = 0.854, F_1,18 = 2693.6, P < 0.001). The nested SAR constructed over spatial scales of core, patch and tree levels proved to be scale‐independent. Main conclusions Unique microhabitats provided by well developed suspended soil accumulations are a habitat template responsible for the diversity of canopy oribatid mites. Species–area relationships of isolate vs. nested species richness data differed in the rate of accumulation of species with increased area. We suggest that colonization history, stability of suspended soil environments, and structural habitat complexity at local and regional scales are major determinants of arboreal oribatid mite species richness.     

基于地理加权回归拓展模型的天然次生林碳储量空间分布

为精准获取区域尺度天然次生林的碳储量及其空间分布格局,以吉林省汪清林业局浪溪林场的天然次生林为研究对象,基于165块局级固定样地,以林分因子、地形因子和土壤因子为影响因子,将普通地理加权回归模型(GWR)作为基础,从空间维度、参数异质性特征和残差空间自相关性3个方面进行改进,构建7类拓展模型,即地理海拔加权回归模型(G...     

基于GWR模型的东江水质空间分异与水生态功能分区验证

流域水生态功能分区研究是我国正在开展的一项重要工作,如何验证分区结果的合理性,是当前亟待解决的问题。采用地理加权回归(GWR)模型评估流域特征对东江水质的影响,验证水质及流域影响空间差异是否与一二级水生态功能分区结果吻合,并对比了GWR模型与普通最小二乘(OLS)模型性能,讨论了GWR在分区验证方面的应用价值及不足。结果显示:1)水质指标以及GWR模型局部解释率(Local R~2)均在一二级水生态功能分区间存在显著差异;2)相比OLS模型,GWR模型校正R~2更高,残差空间自相关指数Moran′s I更低。研究表明东江水生态功能分区结果能合理反映水陆耦合关系,有效解释水质空间差异。此外建议选择总氮(TN)而非溶解氧(DO)和总磷(TP)作为分区验证指标。GWR模型在分区结果验证中具有广泛应用前景。降低数据空间自相关影响及改善距离测度方法是未来GWR模型研究的难点问题。