云南气候变化高分辨率模拟与RCP4.5情景预估

利用高分辨率区域气候模式CCLM（COSMO Model in Climate Mode）对云南省气温和降水的模拟资料,采用多种评价指标,对比分析了试验期（1961~2005年）的模拟结果与同期25个气象站的观测值,并对RCP45情景下的近期（2030~2040年）气温可能变化趋势进行了预估。结果表明：（1）CCLM区域气候模式能够较好地模拟云南省气温的演变趋势,而对降水的模拟能力相对较弱;（2）RCP45情景下,年平均气温（T）、最高气温（Tmax）和最低气温（Tmin）将呈现一致上升趋势,2030~2040年比基准期1986~2005年上升幅度均为12℃;（3）2030~2040年,云南省暖事件发生的可能性将增加,冷事件可能有所减少     

21世纪三峡库区极端气温指数的情景预估

利用气候模式诊断与比较计划(PCMDI）提供的9个新一代气候系统模式的模拟结果，通过多模式集合方法预估分析了3种排放情景(高排放A2、中等排放A1B和低排放B1）下21世纪三峡库区3种极端气温指数的可能变化。结果表明：21世纪三峡库区气温年较差呈震荡的趋势，主要将以增大为主。暖夜指数和热浪指数都将显著增加。整个21世纪，库区气温年较差将增加04~08℃，暖夜指数将增加133%~174%，热浪指数将增加85~133 d。分阶段来看，21世纪前期，气温年较差将增加02~06℃；暖夜指数将增加51%~73%，热浪指数将增加31~41 d；21世纪中期，气温年较差将增加04~10℃，暖夜指数将增加136%~189%，热浪指数将增加77~121 d；21世纪后期，气温年较差将增加06~14℃，暖夜指数将增加191%~289%，热浪指数将增加143~237 d     

长江流域逐月气温空间插值方法的探讨

长江流域地形地貌特征复杂。为了探索适合长江流域的逐月气温数据空间插值方法，在考虑海拔、经纬度、坡度、坡向对气温影响和没有考虑这些地形因子对气温影响的两种情况下，利用流域151个气象站2007年逐月气温与对应的站点地形因子进行回归分析，并与回归直线截距相加产生栅格化的回归气温值，同时对回归多项式的残差分别运用反距离权重法（IDW）、普通克立格法（OK）和样条函数法（SPLINE）气温进行了空间插值，然后将栅格化的回归气温值与残差的插值结果相加得到空间化的逐月气温数据，并利用交叉检验方法对插值精度进行了评估。结果表明：考虑了地形因子影响的3种插值方法的精度都有比较明显的提高，对于普通克立格法，平均绝对误差（MAE）从103℃降到060℃，均方根误差（RMSE）从238℃降到123℃；对于反距离权重法，MAE从110℃降到065℃，RMSE从248℃降到138℃；对于样条插值法，MAE从124℃降到074℃，RMSE从266℃降到151℃。考虑了地形因子影响的空间插值方法整体上要优于没有考虑地形因子的空间插值方法，其中，基于地形因子的普通克里格插值方法结果相对较好     

南京江北新区城市规划对区域热环境影响的WRF模拟研究

区域热环境问题是影响城市生态环境的重要因素之一,如何通过合理的城市规划从而最大限度降低热环境问题带来的负面影响,是城市规划需要解决的关键问题之一。以南京江北新区为研究对象,基于城市规划下垫面类型数据,利用WRF模式(Weather Research and Forecasting model,WRF)来模拟近地表气温变化,进而分析研究区域热环境对城市规划的响应。研究表明:(1)高时空分辨率土地利用数据的输入能够显著提高WRF模式模拟的近地表气温精度,在此基础上耦合城市冠层模式,可以更准确地反映城市内部气温变化情况;(2)区域热环境受城市规划中土地利用类型变化的影响显著,城镇用地增加区域的气温变化幅度在不同时段有差异,其中气温升高幅度最大为1. 8℃,出现在20∶00～22∶00时段;相比于自然地表,城镇用地呈现出升温快和降温慢的特点;(3)与2013年相比,2030年研究区热岛范围扩大,热岛强度增强。因此,初步建议在城市规划时,除有效控制城镇用地增长外,需考虑该用地的分散布局,将其对热环境造成的影响降到最低,更有利于江北新区生态规划建设。     

城、郊能量及辐射平衡特征观测分析

城市化发展导致城、郊能量平衡及辐射平衡存在明显差异。利用城、郊夏季涡度相关系统观测数据,分析了南京城、郊在晴天条件下的辐射平衡、能量平衡、反照率以及储热项变化特征。结果表明：（1）夏季晴天条件下,城、郊辐射收支分配有显著不同。城市辐射陷阱作用导致城区向上短波辐射日均值小于郊区,城、郊分别为246 W/m2和395 W/m2,城市辐射陷阱效应导致截留的短波辐射日平均值可达237 W/m2。同时城市地表向上长波辐射日均值大于郊区自然地表,城、郊分别为5051 W/m2和4883 W/m2,这与城、郊地区不同的地表温度直接相关;（2）城、郊能量平衡分配的方式有显著不同,城市地表储热及感热为能量分配的主要方式,感热全天为正,始终处于不稳定层结。郊区潜热为能量分配的主要方式,夜间存在稳定层结条件;（3）城区总体日平均反照率为011,郊区草地日平均反照率为017,城区水泥楼顶日平均反照率为023。城区铁塔上日均反照率代表整体城市下垫面反照率的特征,城市辐射陷阱作用导致其小于郊区反照率。而水泥楼顶的反照率代表水泥板单一下垫面城市特征,因此大于前2者的反照率日均值。另外,郊区日平均反照率的日变化比城区表现出更明显的不对称性;（4）白天城区储热项占净辐射的56%,而郊区仅占净辐射的7%,并与已有试验中夏季白天储热项数据进行了比较,城市储热项特征与观测所处站点的水泥（硬地）面积所占比重以及是否考虑人为热的影响密切相关     

近40年气候变化对江西自然植被净第一性生产力的影响

根据全球气候变化的趋势,采用植被净第一性生产力模型,对江西省南昌、吉安、赣州3地近40年气候变化对自然植被净第一性生产力（NPP）的影响进行研究,并模拟了3地自然植被NPP在未来气候3种水热条件下的变化趋势。此外还以1980年江西全省自然植被NPP为例分析了自然植被NPP的区域分布特征,结果表明：3地近40年自然植被NPP平均值分别为1319 、1311和1320 t/hm2〖DK1〗·a,总体上都呈上升的趋势。 当年均气温增加2℃且降水量增加20%时,NPP值增加了149%~1585%;随着年均气温增加2℃且降水量减少20%,NPP减少了477%~516%;当年均气温增加2℃且降水量不变时,NPP增加了530%~569%。江西自然植被NPP区域分布特征由东、南、西3个方向向北呈放射状分布,随着地形由高山向丘陵、平原的方向变化而减小。     

江汉平原马铃薯秋播密度对农艺性状及产量的影响

马铃薯（Solanum tuberosum）是一种重要的经济作物,它可以用作粮食、蔬菜、饲料及工业原料作物,江汉平原是马铃薯发展的新区。费乌瑞它是一个适宜于低山平原的早熟、优质马铃薯品种。探索马铃薯高产、优质的栽培技术规程是获得马铃薯最大经济效益、促进马铃薯生产持续发展的重要前提,为了探讨马铃薯秋播的最佳密度,设置了4 000穴/667 m2、7 000穴/667 m2、10 000穴/667 m2、13 000穴/667 m2、16 000穴/667 m2、19 000穴/667 m2等6个播种密度,研究密度对马铃薯费乌瑞它的主要农艺性状及产量的影响。试验表明：播种密度与生育期呈极显著正相关（y=0.000 8x+80514,r1=0990 7*[KG-*2]*）,密度每增加300穴/667 m2,生育期延长24 d;密度与株高（r2）呈显著负相关（r2=-0895 1*）,与出苗率（r3）、主茎数（r4）、单株块茎数（r5）和商品薯率（r6）呈极显著负相关（r3=-0941 1*[KG-*2]*,r4=-0992 4*[KG-*2]*,r5=-0964 4*[KG-*2]*,r6=-0937 9*[KG-*2]*）;密度与单株块茎重（r7）呈二次曲线关系（y=-1E-06x2+0016 8x+13978,r7=-0837 0*）,当密度为8 498穴/667 m2时,单株块茎重达到最大值;密度与单产（r8）为二次曲线关系（y= -8E-06x2+0208x+15761,r8=-0726）,当密度为12 838穴/667 m2时,单产达到最大值;密度与商品薯产量为二次曲线关系（y=-6E-06x2+0141 9x+19395,r8=-0767*）,当密度为1 1087穴/667 m2时,商品薯产量达最大值;10 000穴/667 m2播种密度的马铃薯单产和商品薯产量居首,分别为1 6309 kg/667 m2和1 109 kg/667 m2,综合性状优良;13 000穴/667 m2的效果次之;单产随后的各处理依次为19 000穴/667 m2（1 3046 kg/667 m2）、16 000穴/667 m2（1 2113 kg/667 m2）、7 000穴/667 m2（1 2105 kg/667 m2）、4 000穴/667 m2（8091 kg/667 m2）;商品薯随后的各处理依次为13 000穴/667 m2（8473 kg/667 m2）、16 000穴/667 m2（7679 kg/667 m2）、19 000穴/667 m2（634 kg/667 m2）、4 000穴/667 m2（6181 kg/667 m2）;4 000穴/667 m2的商品薯率最高,19 000穴/667 m2的商品薯率最底,各处理的商品薯率依次为4 000穴/667 m2（7640%）、7 000穴/667 m2（740%）、10 000穴/667 m2（679%）、13 000穴/667 m2（639%）、16 000穴/667 m2（585%）、19 000穴/667 m2（488%     

上海城市草坪CO2排放通量特征及其影响因子

草坪是城市绿地系统的重要组成部分，为量化其在城市碳平衡中的作用，采用静态箱 气相色谱法对上海结缕草草坪CO2排放通量进行了研究。结果表明：CO2排放通量具有明显的日变化和季节变化特征，日变化表现为昼高夜低的特点，最大排放通量一般出现在正午前后，最小在清晨，最大日较差为1 1877 mg/（m2·h）（7月1日），最小892 mg/（m²·h）（1月3日）。CO₂排放通量的季节变化呈单峰型，夏季7月最大，春秋季次之，冬季1月最小，与温度的季节变化一致。观测期间年平均CO2排放通量为6194 mg/（m2·h）。对上海城市结缕草草坪观测期间，温度是CO2排放通量变化的主要影响因素。研究成果不仅可为城市草坪碳排放提供基础数据，也可为城市碳循环研究和管理提供参考.     

三峡水库上空甲烷浓度时空变化及与水库甲烷通量的关系

2010年，对三峡水库水-气界面上空05m处和岸边甲烷的浓度进行全年的观测，并采用静态浮箱-气相色谱法测定水库水-气界面上的甲烷通量，研究了三峡水库上空的甲烷浓度的背景值及其与甲烷排放强度之间的关系。结果表明：三峡水库上空甲烷浓度的年平均值为2216±0224 mL/m3，岸边甲烷浓度的年平均值为2211±0206 mL/m3，两者差异不显著；除夏季时云阳上空的甲烷浓度较高外（2850 mL/m3），其余地点上空的甲烷浓度都接近年平均值；三峡水库上空的甲烷浓度存在明显的空间差异，上游地区（云阳）上空的甲烷浓度（231±033 mL/m3）显著高于下游地区（秭归：214±013 mL/m3，巫山：220±018 mL/m3），坝后河流（三斗坪：221±016 mL/m3）上空的甲烷浓度高于坝前秭归处，但差异不显著。这一格局与水库中甲烷通量的空间变异基本相同，即从上游到下游（云阳、巫山、秭归）江面上甲烷通量依次为0454，0260和0115 mg/（m2·h），呈下降趋势，坝后三斗坪处的甲烷通量（0280 mg/（m2·h））显著高于坝前秭归处通量，这表明：水库的水-气界面上的甲烷排放强度是影响水库上空甲烷浓度变化的主要因素     

基于AVHRR和DEM的重庆城市热岛效应分析

为研究复杂地形上城市热岛(UHI)的特点和变化规律，给城市生态环境分析提供依据，使用重庆市2006年14个时相的NOAA/AVHRR（美国国家海洋大气局/改进的甚高分辨率辐射计）卫星遥感数据和数字高程模型（DEM），得到研究区域的亮温场，并参照热力学概念和模型，根据当地高程，以海拔 300 m 为基准高度将亮温场计算成位温场。结果表明：（1）计算后山脉、河谷处的温度场普遍改变；山脉植被对下垫面的降温作用明显；白天UHI范围较大，长江和嘉陵江河谷对UHI的切隔作用显著；夜间嘉陵江两侧的UHI趋为一体，范围缩小；（2）7月 25～26日存在多个UHI中心，下午和凌晨最大强度分别为64℃和22℃，平均值分别为605℃和205℃，昼夜平均为405℃；（3）夏、秋季白天和夜间平均最大强度分别为523℃和329℃；（4）重庆UHI强度白天大于夜间，与成都的同类研究结果相同；但夏、秋季的平均最大强度为426℃，小于成都的650℃。     

重庆缙云山地区森林土壤酸化特征

通过在2010～2012年8月份对重庆缙云山4种典型林分土壤进行野外取样和室内分析,研究了酸沉降对研究地森林土壤酸化的影响,并从多角度深入分析了4种林分土壤在时间和空间上的酸化机理。结果表明：（1）在常年酸沉降的影响下,缙云山森林土壤酸化现象非常严重,pH值较低,表层土壤pH值为387~447,4种林分土壤pH值都随土层加深递增;（2）4种林分有机质含量相差较大,灌木林最高（377 g/kg）,毛竹林最低（1144 g/kg）,且随土层增加而降低,这主要是受枯落物类型和来源的影响;（3）由于长期酸雨淋溶作用,4种林分土壤CEC含量均较低为1276~1973 cmol/kg,表层含量最高,随土层加深递减;（4）4种林分Al3+含量均较高,阔叶林（4097 mmol/kg）>灌木林（3132 mmol/kg）>针阔混交林（3006 mmol/kg）>毛竹林（2615 mmol/kg）。Al3+含量与土壤pH值呈显著负相关关系     

基于水动力模型及遥感水边线方法的潮滩高程反演

由于潮滩范围广,变化频繁,缺乏常规的地形测量资料,给潮滩开发建设带来极大不便。以崇明东滩为实验区,依据长江口1999~2004年多时相遥感影像的光谱特征,对不同潮情的影像采用不同的波段提取水边线信息。同时利用水动力模型软件（Delft3D）模拟卫星过境时刻潮位并赋值给相应的水边线,对具有高程信息的水边线利用地理信息系统空间分析功能进行空间插值处理,生成潮滩的数字高程模型。将初始高程作为模型输入重新运行水动力模型,并将水边线模拟结果与卫星影像提取的水边线对比进一步修正潮滩地形,直到提取水边线与模拟水边线满足精度要求为止。最后利用崇明东滩2005年4条实测剖面与对应年份高程模拟结果进行对比检验。结果表明,该方法能够较高精度地反演潮滩的地形信息,其平均误差分别达到050 m（σ2=031）、019 m（σ2=004）、015 m（σ2=003）及013 m（σ2=003）。     

长株潭城市群人为热排放对城市热环境影响研究

采用"自上而下"的能源清单法,研究了长株潭城市群2017年人为热排放量及其时空分布特征;同时,利用中尺度天气预报模式(Weather Research and Forecasting,WRF)及其耦合的UCM对夏冬两季分别进行数值模拟试验,定量分析了不同类型建成区人为热排放增温效应的季节性差异.结果表明:(1)长株潭2017年人为热排放总量为3.49x1017J/a,平均人为热排放强度为14.22 W/m2,工业、建筑、交通和人体新城代谢对人为热的贡献率分别为48.15％、40％、11.3％和0.55％.(2)人为热的引入使得城市群主城区夏季和冬季的平均温度分别升高了约0.7℃和1.5℃;冬季增温效果是夏季的2倍.(3)不同密度分区人为热排放导致的增温效果不同,总的来说,工业/商业区＞高密度住宅区＞低密度住宅区.     

三峡库区近50年来的气温变化趋势

利用三峡库区及周边32个气象站点1960~2006年的气温资料，运用线性趋势分析、累积距平分析和〖WTBX〗t〖WTBZ〗检验等统计学方法研究了三峡库区近50 a来的气温变化趋势。研究结果表明：（1）近50 a来三峡库区气温变化总体上呈显著上升趋势，增温率为013℃/10 a；其中1960’s~1980’s末存在一个缓慢降温过程，1980’s末后快速增温。（2）三峡库区各季节平均气温变化过程与年气温变化过程相似，总体上也呈上升趋势，春、夏、秋、冬四季平均气温增温率分别为010、0005、019和021℃/10 a，其中冬季气温上升对库区年平均气温上升的贡献率最大。（3）年均气温跃变出现在1996年，春、夏、秋、冬四季气温跃变点分别出现在1996、1993、1997、1996年，季节气温的跃变与年均气温跃变具有较好的同步性。（4）三峡库区偏暖和显著偏暖年份都发生在1996年以后，其中1998〖JP2〗和2006年为异常偏暖年份；偏冷年份基本出现在1990’s以前，尤其集中在1980’s，但无显著偏冷和异常偏冷年份。〖     

利用微气象测定仪MINCER分析热害条件下的水稻冠层微气象特征

采用模型模拟水稻热害过程及其可能的产量损失是水稻耐高温研究的重要一环,但较精确地测定水稻冠层的微气象条件还存在诸多困难。采用最新研制的水稻冠层微气象测定仪Mincer实地测定法研究了水稻花期热害受害下的冠层微气象特征,并将相关数据与邻近国家基本气象观测站的数据比较,从而形成了水稻冠层、水稻田田面上部(距田面2倍株高,约1.9m)和空旷地(距地面1.5m)等的对比。结果表明,供试品种受到花期热害后,结实率降低4.0%,达到显著受害水平。在此条件下,日平均气温值由水稻田冠层、水稻田上部与空旷地相比依次升高,水稻冠层与水稻田上部的温度值一般分别比空旷地低2.05、0.92℃,且在连续高温的部分时段出现冠层温度的异常升高;相对湿度值反之,由冠层内部、水稻田上部与空旷地依次降低,水稻冠层与水稻田上部的相对湿度值分别比空旷地低高14.55%、5.97%(绝对值)。此外,高温天气下开花前后水稻田上部与冠层温湿度差的变化幅度都小于常温天气,且出现的时间段早于常温。这些结果,对设定水稻高温受害鉴定条件、模拟水稻高温受害过程均具有重要指导意义。     

湘中丘陵区不同经营目标毛竹林年龄结构、生物量分配与碳贮量格局

经营模式对毛竹林生物量、碳贮量具有重要影响。研究了湘中丘陵区毛竹笋用林(Ⅰ)、笋材兼用林 (Ⅱ)和材用林(Ⅲ) 3种不同经营目标下的竹林年龄结构、生物量分配及碳贮量格局。结果表明：应减少1～2 a、增加5～6 a生竹的留养比例，控制达到1～2 a、3～4 a、5～6 a各占1/3左右的立竹年龄结构。不同层次生物量表现为乔木层＞凋落物层＞林下植被层，毛竹笋用林经营有利于增加乔木层生物量。乔木层生物量及所占比例分别为5183～5566 t/hm2、8895%～9293%，林下植被层生物量及所占比例分别为154～258 t/hm2、258%～443%，凋落物层生物量及所占比例分别为269～386 t/hm2、449%～662%。毛竹林总碳贮量排队顺序为Ⅱ（14263 t/hm2）＞Ⅰ（13389 t/hm2）＞Ⅲ（13004 t/hm2），笋材兼用林有利于提高竹林碳贮能力。不同层次碳贮量排列顺序总体均表现为土壤层＞乔木层＞凋落物层＞林下植被层。湘中丘陵区毛竹林生物量、碳贮量较低，应提高集约经营水平     

环形热岛格局演变过程的遥感分析

城市热岛效应的空间格局在一定程度上可以反映城市规划和管理的成败得失，也是城市大气污染的驱动因素之一。〖JP2〗使用2000、2003、2004、2006、2008和2010年2～3月下午成都的6次NOAA/AVHRR（National Oceanic and Atmospheric Administration/Advanced Very High Resolution Radiometer）卫星遥感数据反演亮温，采用均值 标准差法将热岛强度分为7个等级，用面积加权平均法将研究区强度分为高温、中温和低温3类。结果表明：（1）热岛空间格局呈现巨大变异，由前期的中心型演变为中期的环型，后期是热岛环的形态和高温中心的调整阶段；（2）研究区的强度为前期强，2004年迅速减弱，高温类的强度2004年比2003年减弱259℃，2004年之后的平均减弱率仅为024℃/a；（3）城市东北部和西南部高温区的温度梯度出现逆向升降趋势。由于地表温度对城市气温有重要影响，热岛形态的变化是人们从更复杂角度分析大气污染机制、设计热岛数值模型和进行城市规划的客观依据     

基于LUR模型探究城市景观格局对PM2.5浓度的影响——以长株潭城市群为例

随着城市化进程的加快，空气污染问题已成为中国最主要城市问题之一，严重影响公众健康。当前微观尺度下空气监测点周围景观格局对PM25浓度影响的研究较少，以长株潭城市群为例，选取地形、污染源、人口、道路交通、土地利用与城市景观格局6大类预测变量，其中城市景观格局选取边缘密度、连续度、形状指数、斑块平均面积、蔓延度、均匀度指数7个景观指数，运用逐步线性回归模型，探究城市景观格局对PM25浓度的影响。研究结果显示：（1）所选取的土地景观格局指数可以解释研究区PM25浓度的732%的变异，模型拟合较好；（2）影响PM25浓度的土地利用类型包括建设用地、林地、草地与水体。微观尺度下城市各类型景观格局中连续度和形状指数对PM25影响显著，建设用地连续度越高，分布越集聚，PM25浓度越高；水体形状指数越小，形状越简单规则，越易降低PM25浓度；（3）城市整体景观格局中，景观聚集程度与景观多样性等因素对PM25浓度产生重要影响。减少景观内各类型斑块的离散分布，使各景观类型均匀分布于整体景观内，有助于降低PM25浓度。研究结果可为未来大气防治与城市规划提供参考依据。〖HJ1〗〖HJ〗〖JP+1〗     

基于开源GIS的长江流域晴日辐射反演

以开源GIS软件GRASS为平台，首次通过辐射过程模型r.sun和编程语言Shell实现长江流域500 m分辨率的晴日辐射反演（包括直接辐射、散射辐射和反射辐射），为了解长江流域辐射时空分布规律和相关宏观气候、生态建模奠定基础。结论如下：（1）精度验证表明，晴日辐射反演结果较理想；相关晴日指数的误差指标MPE、MAPE和RMSE分别为37%、75%和64%，r.sun模型可进一步用于中国其它区域；（2）长江流域晴日辐射年均日值为2424±245 MJ/(m2·d)，其各组分差异显著，直接、散射、反射辐射分别占总辐射的842%、156%和02%；（3）晴日辐射空间分布梯度明显，沿三级地形阶梯向东辐射值及变异逐渐减小，海拔变化引起的大气衰减因素在决定晴日辐射分布时起重要作用；（4）晴日辐射季节分布不对称，月际变化总体呈倒“U”形分布；直接辐射和散射辐射变化趋势基本一致，但二者比例随大气浑浊度的变化波动剧烈。     

基于CA Markov的土壤有机碳储量空间格局重建研究——以泛长三角地区为例

选取了泛长三角地区的高程、坡度、坡向、GDP、人口、温度、降雨量、河流距离、城市距离和海岸线距离10个驱动因子作为土地利用变化驱动力，基于这些驱动力因子，利用Logistic分析制作土地适宜性图集；运用CA Markov模型模拟得到2000年、1995年、1990年、1985年、1980年、1975年6期的土地利用格局；采用全数检测的方法对模拟结果进行精度检验；在此基础上，得到了泛长三角地区土壤有机碳储量，进而重建历史土壤有机碳储量空间格局。研究结果表明：（1）模型模拟精度为8873%，说明CA Markov模型对历史土地利用模拟具有一定可行性；（2）泛长三角地区历史土壤有机碳储量从今往前逐期递增，且明显高于2010年土壤有机碳储量水平，1995、1990、1985、1980、1975年土壤有机碳储量分别为3235 1、3241 9、3247 6、3251 8、3255 3；（3）各时期泛长三角地区土壤有机碳储量空间格局呈现出南高北低的趋势，土地利用类型与土壤有机碳储量密切相关，城市建成区有机碳储量低于其他区域，土壤有机碳储量增加区域主要出现在其高值区周边；（4）耕地与林地面积的增加，建设用地面积减少，导致了1995~1975年间有机碳储量的增加。基于CA Markov模型可以有效重建历史土壤有机碳储量的空间格局