青藏高原净生态系统碳交换量的遥感评估模型研究

青藏高原作为世界海拔最高的区域,是全球气候变化的敏感区之一。定量估算这一区域的净生态系统碳交换量（NEE）有利于理解陆地生态系统碳平衡对未来气候变化的响应。本文构建了一个模拟该地区NEE动态变化的净碳收支模型（NCBM）。该模型由来源于MODIS影像的增强型植被指数（EVI）、陆地表面水分指数（LSWI）以及来源于地面观测的空气温度和短波辐射共同驱动,并利用青藏高原地区的3种植被类型（包括高寒灌丛、高寒湿地和高寒草甸）的碳通量长期观测数据对模型进行了校准和验证。结果表明,在模型校准站点年,NCBM模型可以模拟NEE观测值81%的变化,均方根误差（RMSE）为0.03molC/m2/d,模型效率（EF）为0.81。在模型验证站点年,NCBM模型可以预测NEE观测值84%的变化,RMSE为0.03molC/m2/d,EF为0.81。在大多数情况下,NCBM模型可以清晰地模拟各植被类型的NEE季节和年际变化。此外,NCBM模型因为结构简单,模型驱动变量易于获取等优势,具有在区域尺度上模拟NEE时空变化的潜力。但是该模型还需要进一步的改进和发展,特别需要提高对植被非常稀疏地区NEE变化的模拟能力。     

Is the soil moisture precipitation feedback enhanced by heterogeneity and dry soils? A comparative study

The interaction between the land surface and the atmosphere is a crucial driver of atmospheric processes. Soil moisture and precipitation are key components in this feedback. Both variables are intertwined in a cycle, that is, the soil moisture – precipitation feedback for which involved processes and interactions are still discussed. In this study the soil moisture – precipitation feedback is compared for the sempiternal humid Ammer catchment in Southern Germany and for the semiarid to subhumid Sissili catchment in West Africa during the warm season, using precipitation datasets from the Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), from the German Weather Service (REGNIE) and simulation datasets from the Weather Research and Forecasting (WRF) model and the hydrologically enhanced WRF-Hydro model. WRF and WRF-Hydro differ by their representation of terrestrial water flow. With this setup we want to investigate the strength, sign and variables involved in the soil moisture – precipitation feedback for these two regions. The normalized model spread between the two simulation results shows linkages between precipitation variability and diagnostic variables surface fluxes, moisture flux convergence above the surface and convective available potential energy in both study regions. The soil moisture – precipitation feedback is evaluated with a classification of soil moisture spatial heterogeneity based on the strength of the soil moisture gradients. This allows us to assess the impact of soil moisture anomalies on surface fluxes, moisture flux convergence, convective available potential energy and precipitation. In both regions the amount of precipitation generally increases with soil moisture spatial heterogeneity. For the Ammer region the soil moisture – precipitation feedback has a weak negative sign with more rain near drier patches while it has a positive signal for the Sissili region with more rain over wetter patches. At least for the observed moderate soil moisture values and the spatial scale of the Ammer region, the spatial variability of soil moisture is more important for surface-atmosphere interactions than the actual soil moisture content. Overall, we found that soil moisture heterogeneity can greatly affect the soil moisture – precipitation feedback.     

Enhanced identification of a hydrologic model using streamflow and satellite water storage data: A multicriteria sensitivity analysis and optimization approach

Hydrologic model development and calibration have continued in most cases to focus only on accurately reproducing streamflows. However, complex models, for example, the so‐called physically based models, possess large degrees of freedom that, if not constrained properly, may lead to poor model performance when used for prediction. We argue that constraining a model to represent streamflow, which is an integrated resultant of many factors across the watershed, is necessary but by no means sufficient to develop a high‐fidelity model. To address this problem, we develop a framework to utilize the Gravity Recovery and Climate Experiment's (GRACE) total water storage anomaly data as a supplement to streamflows for model calibration, in a multiobjective setting. The VARS method (Variogram Analysis of Response Surfaces) for global sensitivity analysis is used to understand the model behaviour with respect to streamflow and GRACE data, and the BORG multiobjective optimization method is applied for model calibration. Two subbasins of the Saskatchewan River Basin in Western Canada are used as a case study. Results show that the developed framework is superior to the conventional approach of calibration only to streamflows, even when multiple streamflow‐based error functions are simultaneously minimized. It is shown that a range of (possibly false) system trajectories in state variable space can lead to similar (acceptable) model responses. This observation has significant implications for land‐surface and hydrologic model development and, if not addressed properly, may undermine the credibility of the model in prediction. The framework effectively constrains the model behaviour (by constraining posterior parameter space) and results in more credible representation of hydrology across the watershed.     

基于投影变换的线阵旋转扫描相机标定方法

针对线阵旋转扫描相机标定数据难以获取及模型解算不稳定的问题,提出了一种基于投影变换的线阵旋转扫描相机标定方法。首先构建了一种新的线阵旋转扫描相机成像模型,该模型可以根据相机旋转平台参数、线阵旋转扫描相机成像平面与其切平面的位置关系将旋转扫描图像投影为框幅式图像;然后根据构建出的成像模型,采用直接线性变换法和非线性优化方法进行相机标定和参数优化。试验结果表明,此方法简单易行,具有较高的精度。     

河北遵化汤泉地热田成因模式

河北遵化汤泉地热资源丰富,阐明其成因模式对于该地热田的进一步开发和热水资源的可持续利用具有一定的指导意义。基于地温测量和水文地球化学分析等方法对其进行了系统研究,结果表明该地热田属中低温对流型地热系统。地热田在其以北的山区接受大气降水补给,补给高程下限约为935 m。地热流体在深循环过程中于正常的大地热流背景下被围岩逐渐加热,热储温度约为130 ℃,循环深度约为4 800 m。在汤泉福泉宫至疗养院一带,构造切割花岗岩体,使得深部地热流体沿破碎带上涌,混合并加热赋存于蓄水构造中的片麻岩裂隙水,形成汤泉地区的地热异常。     

含水层中海水入侵的数学模型

提出了考虑过渡带的适合潜水含水层和承压含水层的海水入侵数学模型,导出了相应的方程和定解条件.经实例验证,模型是合理和可靠的.     

模拟地震的弹簧滑块模型的混沌运动

我们构造了一种自由度为２的双弹簧滑块系统，讨论模型参数的不均匀性和滑块间的相互作用与复杂现象的关系。假定滑块在所受力小于其静摩擦力时保持静止，因此，只有发生了滑动，两滑块所受力之差才改变。如果我们以两滑块受力之差为变量，就可以构造一种一维的映射。这种映射直观地显示了系统的演化，反映了方程终态解的形态。结果表明，用混沌运动来描述地震现象可能是较为合理的。     

Analysis of snowpack dynamics during the spring melt season for a sub‐alpine site using point measurements and numerical modeling

Snowpack dynamics through October 2014–June 2017 were described for a forested, sub‐alpine field site in southeastern Wyoming. Point measurements of wetness and density were combined with numerical modeling and continuous time series of snow depth, snow temperature, and snowpack outflow to identify 5 major classes of distinct snowpack conditions. Class (i) is characterized by no snowpack outflow and variable average snowpack temperature and density. Class (ii) is characterized by short durations of liquid water in the upper snowpack, snowpack outflow values of 0.0008–0.005 cm hr^?1, an increase in snowpack temperature, and average snow density between 0.25–0.35 g cm^?3. Class (iii) is characterized by a partially saturated wetness profile, snowpack outflow values of 0.005–0.25 cm hr^?1, snowpack temperature near 0 °C, and average snow density between 0.25–0.40 g cm^?3. Class (iv) is characterized by strong diurnal snowpack outflow pattern with values as high as 0.75 cm hr^?1, stable snowpack temperature near 0 °C, and stable average snow density between 0.35–0.45 g cm^?3. Class (v) occurs intermittently between Classes (ii)–(iv) and displays low snowpack outflow values between 0.0008–0.04 cm hr^?1, a slight decrease in temperature relative to the preceding class, and similar densities to the preceding class. Numerical modeling of snowpack properties with SNOWPACK using both the Storage Threshold scheme and Richards' equation was used to quantify the effect of snowpack capillarity on predictions of snowpack outflow and other snowpack properties. Results indicate that both simulations are able to predict snow depth, snow temperature, and snow density reasonably well with little difference between the 2 water transport schemes. Richards' equation more accurately simulates the timing of snowpack outflow over the Storage Threshold scheme, especially early in the melt season and at diurnal timescales.     

中川花岗岩构造岩浆活动特征与成矿作用

中川花岗岩在空间上呈同心圆环状分布，由三个单元组成，三个单元之间，岩石矿物学、地球化学成分作有规律性地变化，这种变化与圈层结构变化相适应。构造岩浆演化具有同源性，连续性、脉动性；岩体上升侵位方式为具有多次脉动性的复式气球膨胀体；在空间上与金、铀关系密切，但为成金仅提供热源，而为成铀既提供物源，又提供热源。     

GPS卫星运动状态算法及其软件实现

确定卫星的位置和速度是GPS定位导航的基础,以二体问题为基础,介绍了卫星运动状态的计算方法,建立了计算卫星实时位置、速度的数学模型,并运用C 语言实现了算法,证实了算法的正确性。