基于Google Earth Engine云平台的植被覆盖度变化长时间序列遥感监测

基于Google Earth Engine遥感大数据云计算平台,以云南省南洞地下河流域为例,利用近2 000景30 m分辨率Landsat-NDVI长时间序列数据,采用像元二分模型对研究区1988-2016年的年最大植被覆盖度进行定量估算,并分别从流域整体和像元尺度分析近29 a间植被覆盖度的时空变化特征。研究结果表明:（1）南洞地下河流域大部分区域处于中等覆盖度和中高覆盖度,覆盖度随高程和坡度的增加而增大,其中年最大植被覆盖度 > 60%的区域占流域总面积的45.75%;（2）近29 a来,流域年最大植被覆盖度整体呈现不断增加的趋势,年均增长速率为0.56%,其中植被覆盖度轻微改善或是明显改善的面积占38.84%;（3）相比1988年,2016年高植被覆盖区和中高植被覆盖区面积分别增长50.51%、18.40%;而中等植被覆盖区、中低植被覆盖区和低植被覆盖区面积分别减少24.05%、47.95%和37.72%。封山育林等石漠化治理工程以及气候变化对于流域植被恢复和生态环境重建具有重要影响,其研究成果可为后续石漠化监测提供重要的基础研究数据。      

具有球形空腔饱和土分数导数 黏弹性土的动力特性

将土骨架视为具有分数阶导数本构关系的黏弹性体,基于Biot两相饱和介质模型,建立具有球形空腔饱和分数导数黏弹性土体稳态振动的控制方程。通过引入势函数,得到球对称情形下具有球形空腔饱和分数导数黏弹性土体的位移、应力和孔隙流体压力等解析表达式。考察分数导数模型参数和饱和土参数等对土体振动特性的影响,结果表明,流体压缩性对饱和土体的动力特性有显著影响,而土骨架压缩性和流-固耦合系数的影响相对较小;分数导数阶数对土体动力特性的影响与材料参数比的取值有关。同时,边界不排水条件下饱和土体的动力响应大于排水条件下饱和土体的动力响应。     

四合屯玄武岩斑晶中单个熔体包裹体元素组成及其对岩浆演化的指示

玄武岩斑晶中熔体包裹体成分特征可以推断玄武岩源区物质组成,反映岩浆形成演化过程。利用LA—ICPMS对四合屯义县组玄武岩橄榄石、单斜辉石斑晶中单个熔体包裹体的元素组成进行了分析测试。研究结果表明,橄榄石、单斜辉石斑晶中的熔体包裹体在主、微量元素含量上表现出了比全岩更大的变化范围,但微量元素分配特征总体和全岩一致。单斜辉石斑晶中包裹体的CaO含量、CaO／Al2O3比值和Cr2O3含量随着单斜辉石Mg#值的降低而降低,反映了单斜辉石结晶分离的影响,Al2O3与Sr之间的显著相关关系则记录了斜长石结晶分离作用的影响,MgO—Ni和MgO—CaO／Al2O3的变化则反映了橄榄石的分离结晶作用。包裹体元素组成变化总体受橄榄石、单斜辉石和斜长石的结晶分离作用控制。结合前人研究成果,认为四合屯玄武岩在微量元素和同位素组成上的壳源组分特征可能部分地继承自原岩（即橄榄岩＋榴辉岩部分熔融体反应形成的（橄榄）辉石岩）,而不是岩浆上升过程中受地壳岩石混染的结果。高Mg#值单斜辉石斑晶中少量高Mg馆、高Si含量,低CaO、TiO2、Al2O3和微量元素含量的熔体包裹体反映玄武岩浆上升过程中受到了S1质岩石的混染,这与义县组玄武岩下伏地层为长城系大红裕组石英岩、石英砂岩的地质特征一致。因此,高Fo橄榄石斑晶中的熔体包裹体比采用向全岩中简单添加橄榄石方式计算出的原始熔体可能更能真实反映原始熔体组成。     

西昆仑山普鲁新生代火山岩中包体的发现及其地质意义

报道了在西昆仑用于田县南普鲁火山岩首次发现的包体的产出特下与矿物组成特征，提出其包体为岩浆在高位岩浆房中分离结晶作用形成的堆积相，并估算了其原始岩浆的成分特征。结果表明：其原始岩浆和处于同一火山带的康西瓦火山岩相似，暗示了它们来自于同一源区。但在高位岩浆房中，随着分离结晶作用的进行。岩浆的温度和氧逸度逐渐降低。     

湖南金船塘矿区矽卡岩矿物及Sn元素在绿帘石中的分布特征

金船塘矿区位于千里山花岗岩岩体西侧, 与东侧柿竹园矿区均是湖南柿竹园特大型W-Sn-Mo-Bi矿床目前正在开采的重要矿区.该矿区块状矽卡岩的矿物特点是: 石榴子石以钙铝石榴子石为主, 辉石存在透辉石-钙铁辉石连续系列, 长石以钙长石大量产出为特点, 硅灰石出现了不常见的锰铁硅灰石; 绿帘石中普遍含Sn, SnO₂含量0.2%~2%.根据绿帘石中Sn元素的详细测定和计算, 查明Sn元素在绿帘石中含量与在矿物中位置具幂律函数关系, 计算的分形维为1.93, 从矿物微观的尺度确认Sn等成矿元素在运移和分布上具分形时-空结构, 从而证实自组织的临界状态是柿竹园特大型多金属矿床动力学的基本特点.      

非差模糊度整数固定解PPP新方法及实验

分析了标准模型中非差模糊度无法固定的原因,提出了非差模糊度整数解的小数偏差分离(fractional bias isolating,FBI)模型与方法。实验结果表明,使用30min的静态观测数据,即可达到水平方向mm级、高程cm级的定位精度,相比浮点解水平方向精度提高了一个数量级,高程方向也改善了50%～60%,同时,对流层天顶延迟ZPD的精度也有20%～30%的改善。     

地震数据分数阶域空间预测滤波方法研究

将分数阶傅氏变换引入地震数据空间预测滤波中,研究滤波效果与信号保真的最佳折衷.利用分数阶傅氏变换把信号旋转变换到任意1个角度时"频"平面,运用奇异值分解(SVD)方法求取空间预测滤波算子,实现预测滤波;通过数据相关系数和相对误差能量曲线评价滤波效果.研究结果表明,可以找到1个合适的分数阶域,使得分数阶域空间预测滤波能够在压制噪声和保留信号特征两方面达到最佳折衷,较一般f-x域预测滤波具有更大灵活性.     

不同分辨率影像反演植被覆盖度的参数确定与尺度效应分析

选择新疆古尔通班古特沙漠南缘的石河子地区作为研究区,针对统计回归模型需要大量实测数据,而像元二分模型(dimidiate pixel model,DPM)参数难以确定的缺点,提出了利用统计回归的结果反推像元二分模型参数来建立研究区遥感反演植被覆盖度的方法。运用该方法研究了利用TM和MODIS遥感影像数据分别反演植被覆盖度的参数化方案;同时,还研究了采用不同遥感数据源和不同反演模型时研究区植被覆盖度信息提取中存在的尺度效应。     

一维均质与非均质土柱中溶质迁移的分数微分对流-弥散模拟

分数微分对流-弥散方程(FADE)是模拟溶质迁移问题的新理论,但应用FADE来模拟溶质迁移时能否克服弥散的尺度效应尚待验证。利用长土柱实验资料结合FADE的解析解拟合推求FADE的弥散系数,并分析其与尺度之间的相关关系。研究结果表明,FADE的弥散系数具有随尺度增大而增大的现象,且均质土柱中FADE的弥散系数尺度效应小于非均质土柱中弥散系数尺度效应。在均质土柱中,弥散系数与尺度之间成指数相关关系,在非均质土柱中,弥散系数与尺度之间成幂相关关系。考虑了弥散系数分别与迁移时间和迁移距离呈线性递增两种相关关系,进而分别构建了3种考虑弥散尺度效应的FADE模型,并提出了求解的差分方法。利用上述3种考虑弥散尺度效应的FADE来模拟和预测不同空间位置处的溶质迁移过程。结果表明,对均质土柱中的溶质迁移可得到较好的模拟结果;对于非均质土柱,其模拟结果与实测结果仍然存在一定的差异。     

Using a shifted Grünwald-Letnikov scheme for the Caputo derivative to study anomalous solute transport in porous medium

We propose an extension of the shifted Grünwald-Letnikov method to solve fractional partial differential equations in the Caputo sense with arbitrary fractional order derivative α and with an advective term. The method uses the relation between Caputo and Riemann-Liouville definitions, the shifted Grünwald-Letnikov, and the traditional backward and forward finite difference method. The stability of the method is investigated for the implicit and explicit scheme with homogeneous boundary conditions, and a stability criterion is found for the advective-dispersive equation. An application of the method is used to solve contaminant diffusion and advective-dispersive problems. The numerical solution for the fractional diffusion and fractional advection-dispersion is compared with their respective analytical solutions for different time and space grid refinements. The diffusion simulation exhibited a good fit between the analytical and numerical solutions, with the explicit scheme going from stable to unstable as the time and space refinement changes. The fractional advection-dispersion application produced small deviations from the analytical solution. These deviations, however, are analogous to the numerical dispersions encountered in conventional finite difference solutions of the advection-dispersion equation. The new method is also compared with the traditional L2 method. Notably, an example that involves asymmetrical fractional conditions, a fractional diffusivity that depends on time, and a source term show how the methods compare. Overall, this study assesses the quality and easiness of use of the numerical method.