石漠化是严峻的生态环境问题,云南省是我国石漠化最为严重的省份之一。本文分别获取云南省石漠化、非石漠化区各5个站点的1970~2017年逐日降水量、年平均和月平均气温资料,采用变差系数Cv值、Theil-Sen Median趋势分析、Mann-Kendall检验,通过对比分析法,探究局地气候对石漠化的响应规律。结果表明:1) 石漠化区降水Cv值高于非石漠化区,降水波动性较非石漠化区更明显;2) 石漠化区RX5和RX1有增加的趋势,但最长持续降水时间却是减少的,最长干旱时间趋势仍然增加,石漠化区降水量与降水时长、干旱程度的变化并不一致,较非石漠化地区更为反常;3) 石漠化区RX5、RX1、CDD发生突变的时间以及变化趋势达显著的时间均早于非石漠化区,且石漠化区变化趋势较非石漠化区显著,石漠化区和非石漠化区CWD均在70年代开始有突变发生,石漠化区除耿马站点外均在90年代减少趋势达到显著;4) 石漠化地区年均温高于相对阈值的年份频数大于非石漠化地区,极端气温出现的频率较非石漠化区更高,石漠化区气温的Cv值高于非石漠化区,气温的年际波动更大,非石漠化区年均气温总体上呈现先下降后上升的变化趋势,而石漠化区年均气温变化趋势在上升和下降之间反复变化,较非石漠化区变化更剧烈。 Rocky desertification is a severe ecological and environmental problem, Yunnan province is one of the most serious rocky desertification provinces in China. Daily precipitation, annual and monthly mean temperature data from 1970 to 2017 were obtained from 5 stations in rocky desertification and non-rocky desertification areas in Yunnan Province. Coefficient of variation Cv value, TheilSen Median trend analysis and Mann-Kendall test were used. The response of local climate to rocky desertification was studied by comparative analysis. The results show that: 1) The Cv value of precipitation in rocky desertification area is higher than that in non-rocky desertification area, and the precipitation fluctuation is more obvious than that in non-rocky desertification area. 2) RX5 and RX1 showed an increasing trend in the rocky desertification area, but the longest sustained precipitation time decreased, while the trend of the longest drought time still increased. Precipitation in the rocky desertification area was not consistent with the change of precipitation time and drought degree, which was more abnormal than that in the non-rocky desertification area. 3) The rocky desertification area RX5, RX1, CDD mutations of time and change trend of significant time is earlier than the non-rocky desertification area, and the trend of rocky desertification area change significantly from the non-rocky desertification area, rocky desertification area and the rocky desertification area CWD began to have a mutation in the 70 s, and the decreasing trend of rocky desertification area, in addition to the Gengma Station, was significant in the 90 s; 4) The annual average temperature in rocky desertification area is higher than that in non-rocky desertification area, and the frequency of extreme temperature is higher than that in non-rocky desertification area. The Cv value of temperature in rocky desertification area is higher than that in non-rocky desertification area, and the interannual fluctuation of temperature is larger. The annual average temperature in non-rocky desertification area shows a general trend of decreasing first and then increasing. However, the trend of annual average temperature in rocky desertification area changed repeatedly between rising and falling, which was more severe than that in non-rocky desertification area.
石漠化是严峻的生态环境问题,云南省是我国石漠化最为严重的省份之一。本文分别获取云南省石漠化、非石漠化区各5个站点的1970~2017年逐日降水量、年平均和月平均气温资料,采用变差系数Cv值、Theil-Sen Median趋势分析、Mann-Kendall检验,通过对比分析法,探究局地气候对石漠化的响应规律。结果表明:1) 石漠化区降水Cv值高于非石漠化区,降水波动性较非石漠化区更明显;2) 石漠化区RX5和RX1有增加的趋势,但最长持续降水时间却是减少的,最长干旱时间趋势仍然增加,石漠化区降水量与降水时长、干旱程度的变化并不一致,较非石漠化地区更为反常;3) 石漠化区RX5、RX1、CDD发生突变的时间以及变化趋势达显著的时间均早于非石漠化区,且石漠化区变化趋势较非石漠化区显著,石漠化区和非石漠化区CWD均在70年代开始有突变发生,石漠化区除耿马站点外均在90年代减少趋势达到显著;4) 石漠化地区年均温高于相对阈值的年份频数大于非石漠化地区,极端气温出现的频率较非石漠化区更高,石漠化区气温的Cv值高于非石漠化区,气温的年际波动更大,非石漠化区年均气温总体上呈现先下降后上升的变化趋势,而石漠化区年均气温变化趋势在上升和下降之间反复变化,较非石漠化区变化更剧烈。
石漠化,气候响应,降水变化,气温变化,云南
Yiqin Cao, Wenrong Ding*
Faculty of Geography, Yunnan Normal University, Kunming Yunnan
Received: Jun. 19th, 2022; accepted: Jul. 18th, 2022; published: Jul. 25th, 2022
Rocky desertification is a severe ecological and environmental problem, Yunnan province is one of the most serious rocky desertification provinces in China. Daily precipitation, annual and monthly mean temperature data from 1970 to 2017 were obtained from 5 stations in rocky desertification and non-rocky desertification areas in Yunnan Province. Coefficient of variation Cv value, Theil-Sen Median trend analysis and Mann-Kendall test were used. The response of local climate to rocky desertification was studied by comparative analysis. The results show that: 1) The Cv value of precipitation in rocky desertification area is higher than that in non-rocky desertification area, and the precipitation fluctuation is more obvious than that in non-rocky desertification area. 2) RX5 and RX1 showed an increasing trend in the rocky desertification area, but the longest sustained precipitation time decreased, while the trend of the longest drought time still increased. Precipitation in the rocky desertification area was not consistent with the change of precipitation time and drought degree, which was more abnormal than that in the non-rocky desertification area. 3) The rocky desertification area RX5, RX1, CDD mutations of time and change trend of significant time is earlier than the non-rocky desertification area, and the trend of rocky desertification area change significantly from the non-rocky desertification area, rocky desertification area and the rocky desertification area CWD began to have a mutation in the 70 s, and the decreasing trend of rocky desertification area, in addition to the Gengma Station, was significant in the 90 s; 4) The annual average temperature in rocky desertification area is higher than that in non-rocky desertification area, and the frequency of extreme temperature is higher than that in non-rocky desertification area. The Cv value of temperature in rocky desertification area is higher than that in non-rocky desertification area, and the interannual fluctuation of temperature is larger. The annual average temperature in non-rocky desertification area shows a general trend of decreasing first and then increasing. However, the trend of annual average temperature in rocky desertification area changed repeatedly between rising and falling, which was more severe than that in non-rocky desertification area.
Keywords:Rocky Desertification, Climate Response, Precipitation Change, Temperature Change, Yunnan
Copyright © 2022 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
全球气候变化对人类社会有着巨大影响,一直是众多组织和科研机构的重要研究内容。IPCC第五次评估报告指出,在全球气候变暖背景下,极端事件将在21世纪变得更加频繁 [
云南省是我国石漠化最为严重的地区之一,滇东北、滇中、滇东南、滇西北和滇西是集中分布区,其中,滇东南石漠化面积分布最广、石漠化程度严重 [
依据典型性、数据可得性和连续性,在云南岩溶区和非岩溶区甄选了10个站点(图1)。其中,泸西、广南、砚山站点选自中亚热带高原山地石漠化区,屏边、耿马站点则位于南亚热带中低山河谷盆地石漠化区 [
图1. 石漠化区和非石漠化区站点分布
表1为所选站点地理位置及周边环境情况,泸西、广南、屏边、砚山、耿马的石漠化比例分别为22.45%、50.31%、42.33%、37.65%和21.35%。
站点 Station | 经度 Longitude | 纬度 Latitude | 海拔 Altitude/m | 归一化植被指数NDVI | 石漠化比例 Ratio of rocky desertification/% |
---|---|---|---|---|---|
泸西 | E 103˚46' | N 24˚32' | 1704.3 | 0.528 | 22.45 |
广南 | E 105˚04' | N 24˚04' | 1249.6 | 0.531 | 50.31 |
屏边 | E 103˚41' | N 22˚59' | 1414.1 | 0.739 | 42.33 |
砚山 | E 104˚20' | N 23˚36' | 1542.2 | 0.351 | 37.65 |
耿马 | E 99˚24' | N 23˚32' | 1120.1 | 0.843 | 21.35 |
澜沧 | E 99˚56' | N 22˚34' | 1054.8 | 0.63 | // |
江城 | E 101˚51' | N 22˚35' | 1120.5 | 0.809 | // |
勐腊 | E 101˚34' | N 21˚29' | 631.9 | 0.806 | // |
思茅 | E 100˚59' | N 22˚47' | 1202.3 | 0.585 | // |
景东 | E 100˚50' | N 24˚27' | 2083.5 | 0.432 | // |
表1. 所选站点地理位置及周边环境
注:“//”表示无石漠化发育。
为客观揭示局地气候对石漠化的响应规律,依据资料的连续性、准确性及最长时段性等标准,分别获取符合条件的石漠化区和非石漠化区站点1970~2017年的逐日降水量、年平均气温和月平均气温资料进行分析。数据来自国家气象信息中心(http://data.cma.cn/site/index.html),均经过严格的质量控制,能满足研究的要求。
气候变化检测和指数联合专家组定义了极端降水指数。为探究石漠化背景下极端降水变化规律,本文采用其中的4个指数(表2)。降水和气温数据的数理统计基于Matlab完成。结合变差系数Cv值分析降水指标和月平均气温的离散程度,通过Sen斜率分析降水指标的变化趋势及显著性,采用Mann-Kendall检验法分析降水指标和年平均气温。得到石漠化区与非石漠化区的降水、气温特征及差异,通过对比分析,找到局地气候对石漠化的响应规律。
指数名称/单位 Index name/Unit | 缩写 Abbreviation | 定义 Definition |
---|---|---|
最长持续干旱时间/d | CDD | 日降水量 < 1 mm的最大持续日数 |
最长持续降水时间/d | CWD | 日降水量 ≥ 1 mm的最大持续日数 |
日最大降水量/mm | RX1 | 每年最大1日的降水量 |
5日最大降水量/mm | RX5 | 年内连续5日的最大降水量 |
表2. 选取的各极端降水指数及其定义
1) 变差系数
变差系数又称变异系数、离散系数,用于分析序列总体的离散程度,可以表示序列的相对变化(波动)程度,其计算公式如下 [
Cv = S x ¯ = 1 x ¯ 1 n ∑ t = 1 n ( x t − x ¯ ) 2 (1)
式中:S表示标准差, x ¯ 为平均数, x t 表示第t年的气温或降水数值。Cv值越高,序列的波动越大。
2) Theil-Sen Median趋势分析
Theil-Sen Median趋势分析又称Sen斜率估计,是一种稳健的非参数统计法。该方法计算效率高,对于测量误差和离群数据不敏感,常运用于长时间序列数据的趋势分析 [
β = median ( x j − x i j − i ) , ∀ j > i (2)
式中:xj和xi为时间序列数据。Β > 0表示降水呈现增加趋势;β < 0表示降水呈现减少趋势。
3) Mann-Kendall(M-K)检验
Mann-Kendall是一种非参数统计检验法,由Mann和Kendall提出并在20世纪中叶进一步完善。其优点是不需要测量值服从正态分布,也不要求趋势是线性的,并且不受缺失值和异常值的影响 [
对于具有n个样本量的时间序列 X 1 , X 2 , ⋯ , X n ,构造一秩序列:
S k = ∑ i = 1 k R i ( k = 2 , 3 , ⋯ , n ) (3)
式中:Ri表示Xi大于Xj( 1 ≤ j ≤ i )的积累数。在时间序列随机独立的假设下,定义统计量:
U F k = S k − E ( S k ) V a r ( S k ) ( k = 1 , 2 , ⋯ , n ) (4)
式中:当k = 1时,UF1= 0。E(Sk),Var(Sk)是累积数Sk的均值和方差,在 X 1 , X 2 , ⋯ , X n 相互独立,且有相同连续分布时,它们可由下式算出:
E ( S k ) = n ( n − 1 ) 4 (5)
V a r ( S k ) = n ( n − 1 ) ( 2 n + 5 ) 72 (6)
按时间序列X逆序 X n , X n − 1 , ⋯ , X 1 ,再重复上述过程,同时使 U B k = − U F k , k = n , n − 1 , ⋯ , 1 ,且当k = 1时,UB1= 0。
若UF > 0,则序列呈上升趋势,反之下降,当它们超过临界置信水平线时(检验置信水平a = 0.05时,置信水平线为±1.96),表明变化趋势显著。如果UF和UB曲线出现交点,且交点在临界线之间,那么交点对应的便是突变开始时间。
表3为降水指数CWD和RX1的Cv值,从表中可知,石漠化区泸西、耿马、砚山最长持续降水日数CWD的Cv值分别为0.44、0.38和0.44,而非石漠化区澜沧、景东、思茅则分别为0.35、0.35和0.37。可见,石漠化区最长持续降水日数的Cv值大于非石漠化区,石漠化区降水波动性较非石漠化区更明显。石漠化区泸西、耿马、砚山最大1日降水量RX1的Cv值在0.30~0.40之间,非石漠化区澜沧、景东、思茅则介于0.30~0.34之间。说明石漠化区每年最大1日降水量的波动较非石漠化区更大。
站点Station | Cv值Coefficient of variation | |
---|---|---|
CWD | RX1 | |
泸西 | 0.44 | 0.40 |
耿马 | 0.38 | 0.30 |
砚山 | 0.44 | 0.34 |
澜沧 | 0.35 | 0.30 |
景东 | 0.35 | 0.34 |
思茅 | 0.37 | 0.32 |
表3. 降水指数CWD和RX1的Cv值
站点Station | Cv值Coefficient of variation | |
---|---|---|
RX5 | RX1 | |
泸西 | 0.27 | 0.34 |
屏边 | 0.27 | 0.34 |
耿马 | 0.29 | 0.41 |
澜沧 | 0.24 | 0.41 |
江城 | 0.25 | 0.44 |
景东 | 0.30 | 0.37 |
表4. 降水指数RX5和CDD的Cv值
从表4中可知,石漠化区泸西、屏边、耿马连续5日最大降水量RX5的变差系数分别为0.27、0.27和0.29,非石漠化区澜沧、江城、景东则分别为0.24、0.25、0.30。可见,除景东外,非石漠化区澜沧、江城RX5的变差系数均小于石漠化区,总体上,石漠化区RX5波动性大于非石漠化区。石漠化区CDD变差系数总体上偏低,在最长持续干旱时间上石漠化区和非石漠化区的Cv值对比不明显。
表5为降水指标的Sen指数,从表中可知,石漠化区广南、屏边、耿马的RX5呈增加趋势,但均未达到0.05的置信水平,而泸西、砚山的RX5呈减少趋势,其中泸西达到显著减少;除景东外,非石漠化区RX5均呈减少趋势,其中思茅达到显著减少水平。石漠化区泸西RX1的Sen指数为−0.59,减少趋势达到0.05的置信水平,其余站点RX1均呈增加趋势;非石漠化区除景东的RX1增加外均呈减少趋势,但都未达到显著水平。从连续5日和最大1日降水量看,泸西均呈减少趋势且都达到0.05置信水平,是石漠化区降水减少最显著的站点。
站点Station | Sen指数Sen index | |||
---|---|---|---|---|
RX5 | RX1 | CWD | CDD | |
泸西 | −0.69* | −0.59* | −0.12* | 0.14 |
广南 | 0.12 | 0.18 | −0.12* | 0.14* |
屏边 | 0.37 | 0.19 | −0.22* | 0.00 |
砚山 | −0.26 | 0.07 | −0.13* | 0.09 |
耿马 | 0.42 | 0.27 | −0.21* | 0.35* |
澜沧 | −0.04 | −0.04 | −0.13 | 0.32* |
江城 | −0.60 | −0.10 | −0.33* | 0.41* |
勐腊 | −0.12 | −0.53 | −0.06 | 0.67* |
思茅 | −1.01* | −0.25 | −0.18* | 0.47* |
景东 | 0.25 | 0.48 | −0.14* | 0.37* |
表5. 降水指标的Sen指数
注:*表示达到0.05置信水平。
石漠化区和非石漠化区CWD均呈减少趋势,而石漠化区减少趋势均达到了0.05的置信水平,可见石漠化区最长持续降水日数的减少趋势较非石漠化区更显著。石漠化区除屏边的CDD无明显增减趋势外均增加,其中,广南和耿马CDD增加趋势达到0.05显著水平。非石漠化区CDD也呈上升趋势达到显著水平。
结合以上结果可以发现,从总体上看,虽然石漠化区最大5日降水量和最大1日量有增加的趋势,但最长持续降水时间却是减少的,最长干旱时间趋势仍然增加,可见石漠化区降水量与降水时长、干旱程度的变化并不一致,较非石漠化地区更为反常;非石漠化地区的最大5日和最大1日降水量呈减少趋势,最长持续降水时间的变化趋势也减少,最长干旱时间趋势增加,降水量与降水时长、干旱程度变化一致。
在用Sen斜率进行总体趋势分析后,进一步采用M-K检验方法,结果显示,非石漠化区最大5日降水量RX5发生突变的时间集中在80和90年代。从图2可知,江城1970~1989年间UF曲线在0值附近上下波动,RX5基本保持稳定,1989年发生突变,此后UF曲线位于0值以下,表明RX5从1989年后一直呈减少趋势,但未超过置信区间,减少趋势不显著;勐腊1970~1975年间RX5基本保持稳定,1975~1986年呈增加趋势,从1986年发生突变后RX5减少,但减少趋势并不显著。石漠化区则从70年代开始就有突变发生。其中,砚山RX5在1973年发生突变,由增加趋势转为减少,并于1976~1978年间减少趋势达到显著水平;泸西1982年往后RX5由增加趋势转为减少,在1985年发生突变,2008年以后RX5减少趋势显著。可见,石漠化区最大5日降水量发生突变的时间以及减少趋势达显著的时间均早于非石漠化区,且石漠化区减少趋势较非石漠化区显著。
非石漠化区最大1日降水量RX1发生突变的时间集中在1980年前后,石漠化区则集中在70和80年代。从Sen斜率总体趋势分析结果可知,非石漠化区景东站点和石漠化区泸西站点RX1变化趋势不同于各自同区的其他站点,进一步对它们的M-K检验结果进行分析,结果为图3,图中景东1983年发生突变,RX1由减少趋势转为增加,1990年开始达到显著水平,并在此后多次显著增加;泸西1970~1980年RX1呈减少趋势,之后UF曲线在0值附近波动,从1985年开始RX1一直呈减少趋势,1988年发生突变,并于2003年以后减少趋势都保持显著水平。作为石漠化区和非石漠化区RX1变化趋势与同区显著区别的站点,泸西最大1日降水量总体呈减少趋势且在2000年以后达到显著;景东则呈增加趋势并于90年代开始达到显著。
图2. 石漠化区和非石漠化区RX5的M-K检验
图3. 石漠化区和非石漠化区RX1的M-K检验
石漠化区和非石漠化区最长持续降水时间CWD均在70年代开始有突变发生,此外,石漠化区除耿马站点外均在90年代减少趋势达到显著。从图4可知,石漠化区泸西CWD在1972年往后一直呈减少趋势,1979年突变,1990年往后减少趋势一直保持在显著水平;非石漠化区勐腊1970~1980年间CWD先减少后增加,1977年发生一次突变,1980往后减少,1991~1996、2003~2011年间达到显著。可见,与勐腊相比,泸西CWD显著减少的持续时间更长。
图4. 石漠化区和非石漠化区CWD的M-K检验
非石漠化区除景东站点外CDD均在1980年前后发生突变,澜沧发生突变的时间为1979年,勐腊突变发生时间为1980年,江城和思茅1981年发生突变;石漠化区CDD从70年代开始至2000年以后均有突变发生,屏边和泸西发生突变的年份分别为1972年和2005年,砚山、广南、耿马均在80年代发生突变,时间分别为1980年、1983年、1988年。此外,非石漠化区除勐腊外均在2000年以后增加趋势达到显著,而石漠化区从80年代开始增加趋势就能够达到显著。可见,石漠化区CDD发生突变的时间较非石漠化区更分散,增加趋势达到显著的时间也早于非石漠化区。
选用云南省气象站1970~2017年石漠化区泸西、广南和非石漠化区澜沧、江城的年平均和月平均气温数据,筛选出夏、冬两个季节的月均气温数据。通过计算百分位值和Cv值,并采用Mann-Kendall检验,对比分析石漠化地区与非石漠化地区的气温变化特征。
通过计算所选站点年均温的第95个百分位值,发现石漠化地区高于相对阈值的年份频数大于非石漠化地区,由此可以推测石漠化地区极端气温出现的频率更高。进一步研究显示,石漠化区气温的Cv值均高于非石漠化区,从季节上看,所选站点气温的Cv值都是冬季高于夏季。其中,广南的Cv值最大,夏季为0.02、冬季为0.12,说明广南相对于其他站点气温的年际变化较大。由此可推断石漠化地区由于石漠化现象严重,植被NDVI指数较低,地表裸露导致地表接受太阳辐射的强度较大,吸热散热变化快,相比非石漠化地区夏冬两季温差变化较大。
由图5泸西和澜沧,广南和江城的M-K检验结果对比分析可知,石漠化区年均气温变化趋势在上升和下降之间反复变化,较非石漠化区变化更剧烈。总体上,非石漠化区年均气温变化趋势为先下降后上升,1970~1979年间下降,1980年往后上升;而石漠化区年均气温变化趋势在上升和下降之间反复变化。广南1970~1986年间气温升降变化趋势最为反复,1970~1987年除1971、1976~1978,1985~1986年外均上升,
图5. 石漠化区和非石漠化区年均气温M-K检验
1987往后一直呈上升趋势。泸西1970~1979年除1973年上升外均呈下降趋势,1979~1996年除1985和1986年下降外均上升,1997年往后一直呈上升趋势。此外,非石漠化区于90年代后期发生突变,其中,澜沧和江城发生突变的时间分别为1998年和1996年;石漠化区发生突变的时间集中在2000年以后,其中,泸西的年均气温于2008年突变,广南为2002年。
1) 石漠化区CWD的Cv值(0.38~0.44)大于非石漠化区(0.35~0.37),RX1的Cv值(0.30~0.40)波动较非石漠化区(0.30~0.34)更大,除景东外,非石漠化区RX5的变差系数(0.24~0.25)小于石漠化区(0.27~0.29)。石漠化区CDD变差系数总体上偏低,在最长持续干旱时间上石漠化区和非石漠化区的Cv值对比不明显。而石漠化区降水Cv值高于非石漠化区,降水波动性较非石漠化区更明显。
2) 从总体上看,虽然石漠化区RX5和RX1有增加的趋势,但最长持续降水时间却是减少的,最长干旱时间趋势仍然增加,可见石漠化区降水量与降水时长、干旱程度的变化并不一致,较非石漠化地区更为反常;非石漠化地区的最大5日和最大1日降水量呈减少趋势,最长持续降水时间的变化趋势也减少,最长干旱时间趋势增加,降水量与降水时长、干旱程度变化一致。
3) 非石漠化区RX5和RX1发生突变的时间分别集中在80~90年代和1980年前后,石漠化区则从70年代开始就有突变发生,发生突变的时间早于非石漠化区。石漠化区和非石漠化区CWD均在70年代开始有突变发生。非石漠化区除景东站点外CDD均在1980年前后发生突变,石漠化区CDD从70年代开始至2000年以后均有突变发生,石漠化区CDD发生突变的时间较非石漠化区更分散,增加趋势达到显著的时间也早于非石漠化区。
4) 石漠化区年均温高于相对阈值的年份频数大于非石漠化区,极端气温出现的频率较非石漠化区更高。石漠化区气温的Cv值均高于非石漠化区,气温的年际波动较非石漠化区更大。从季节上看,所选站点气温的Cv值都是冬季高于夏季。非石漠化区年均气温在90年代后期发生突变,总体上呈现先下降后上升的变化趋势,而石漠化区于2000年以后发生突变,年均气温变化趋势在上升和下降之间反复变化,较非石漠化区变化更剧烈。
国家自然科学基金项目(41761105)资助。
曹怡钦,丁文荣. 石漠化影响下的局地气候响应研究Local Climate Response to Rocky Desertification[J]. 气候变化研究快报, 2022, 11(04): 551-561. https://doi.org/10.12677/CCRL.2022.114058
https://doi.org/10.1016/j.earscirev.2014.01.005
https://doi.org/10.1002/ldr.592
https://doi.org/10.1038/s41598-017-19088-x
https://doi.org/10.1007/s12665-014-3847-8
https://doi.org/10.1016/j.gloplacha.2015.01.009
https://doi.org/10.1016/j.ecolind.2014.07.031