利用伴随法反演风应力系数

基于无结构网格有限体积法的正模型和伴随数据同化模型,建立了风应力系数反演表达式,将风应力系数的确定变成依赖于区域内观测值的数值迭代过程,实现了风应力系数确定的自动化,并成功进行了概化河口一系列风应力系数反演的伴随同化实验,反演得到的结果精度较高,实验结果表明：1）伴随法反演风应力系数是合理可行的;2）充分的观测数据能够提高反演的精度;3）风速在10 m/s或15 m/s的情况下,对于风应力系数反演的结果影响不大。     

稻虾轮作模式面源污染排放特征及其能值评估

为探明稻虾轮作模式面源污染排放特征并合理评价该模式的环境可持续性,通过对江汉平原稻虾轮作模式小龙虾养殖排放尾水中总氮(TN)、总磷(TP)、COD和氨氮(NH₃-N)浓度进行监测,对稻虾轮作模式稻田养殖小龙虾的排污系数进行了估算,并采用等标污染负荷法进行了主要污染物解析。运用能值分析方法对稻虾轮作模式进行了包含面源污染的能值评估,对比单季稻模式,对其可持续发展能力进行了定量评价。结果表明：江汉平原稻虾轮作模式小龙虾养殖排放尾水中TN、TP、COD和NH₃-N的浓度范围分别为0.53～5.36、0.12～0.70、6.60～78.39和0.34～1.75 mg/L,TN、TP和COD平均排放浓度高于《地表水环境质量标准》(GB 3838-2002)Ⅲ类水质标准。等标污染负荷法分析结果表明TN的等标负荷比最高,是稻虾轮作模式面源污染控制的关键污染物。稻虾轮作模式小龙虾养殖排放尾水中TN、TP、COD和NH₃-N的排污系数分别为2.994、0.458、35.132和1.405 kg/t,表明稻虾轮作模式面源污染排放系数较低,对...     

新疆阿尔泰可可托海伟晶岩型稀有金属矿床中磷灰石地球化学特征及意义

磷灰石可以记录和保存岩浆和热液活动的信息。可可托海伟晶岩型稀有金属矿床磷灰石发育,为研究该矿床伟晶岩成岩成矿过程提供了优良的条件。已有对可可托海伟晶岩型稀有金属矿床磷灰石的研究集中在其稀土元素特征,较少讨论其对伟晶岩成岩成矿过程的制约。本文选取可可托海伟晶岩型稀有金属矿床富矿伟晶岩脉（3号脉）和相对贫矿伟晶岩脉（1号、2b号和3a号脉）中的磷灰石作为研究对象,进行磷灰石岩相学和地球化学研究。岩相学分析表明,磷灰石主要与钠长石、石英、白云母、锰铝榴石等伴生。EPMA分析显示,磷灰石F含量为3.67%~4.41%,Cl含量小于0.67%,较低的Cl含量表明伟晶岩熔体出溶的流体Cl含量较低;大部分磷灰石MnO含量为4.67%~8.71%,但2b号脉磷灰石MnO含量变化较大（1.23%~14.28%）,这是由于2b号脉磷灰石具有分带结构,暗示其遭受后期热液作用,促使磷灰石溶解-再沉淀,导致MnO含量发生较大变化。LA-ICP-MS分析显示,贫矿伟晶岩脉磷灰石的稀土元素含量较低（< 180×10^-6）;相反,富矿伟晶岩脉磷灰石的稀土元素含量较高（>700×10^-6）,并具有明显的四分组效应（TE_1-3平均值为1.7）。1号脉和3a号脉磷灰石均显示轻稀土元素富集,反映其形成过程中有含Cl热液的参与。3号脉磷灰石显示强烈Eu负异常和Ce正异常,而2b号脉磷灰石显示强烈Ce负异常和中等Eu负异常,这种Eu、Ce异常的差异可能与岩浆-热液阶段大量流体出溶密切相关。磷灰石的沉淀将导致热液中HF含量的降低,促使磷灰石周围铌钽矿结晶和Nb、Ta进入磷灰石中。可见,在伟晶岩形成过程中,磷灰石并非保持稳定,其分带结构和主微量成分变化记录了后期热液活动,暗示后期热液活动对伟晶岩的成矿具有重要作用。

     

Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada

Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (S_ke and S_kv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

Effects of surface flux parameterization on the numerically simulated intensity and structure of Typhoon Morakot (2009)

The effects of surface flux parameterizations on tropical cyclone(TC) intensity and structure are investigated using the Advanced Research Weather Research and Forecasting(WRF-ARW) modeling system with high-resolution simulations of Typhoon Morakot(2009).Numerical experiments are designed to simulate Typhoon Morakot(2009) with different formulations of surface exchange coefficients for enthalpy(C_K) and momentum(C_D) transfers,including those from recent observational studies based on in situ aircraft data collected in Atlantic hurricanes.The results show that the simulated intensity and structure are sensitive to C_K and C_D,but the simulated track is not.Consistent with previous studies,the simulated storm intensity is found to be more sensitive to the ratio of C_K/C_D than to C_K or C_D alone.The pressure-wind relationship is also found to be influenced by the exchange coefficients,consistent with recent numerical studies.This paper emphasizes the importance of C_D and C_K on TC structure simulations.The results suggest that C_D and C_K have a large impact on surface wind and flux distributions,boundary layer heights,the warm core,and precipitation.Compared to available observations,the experiment with observed C_D and C_K generally simulated better intensity and structure than the other experiments,especially over the ocean.The reasons for the structural differences among the experiments with different C_D and C_K setups are discussed in the context of TC dynamics and thermodynamics.     

M-Quadtree索引:一种基于改进四叉树编码方法的云存储环境下空间索引方法

为了解决基于"键-值"模型的云存储环境仅支持简单的关键字查询,不支持多维空间查询的问题,提出了一种新的分布式空间索引方法——M-Quadtree索引。在索引构建过程中,设计了一种基于改进四叉树的空间数据划分方法,该方法规定了叶节点区域的最小数据量,通过四叉树叶节点的再合并,解决了划分后各子区域间存储量不平衡的问题,并且满足了MapReduce并行化要求。给出了MapReduce框架下M-Quadtree索引的快速构建、查询与更新算法,并在搭建的Hadoop平台进行了关键参数对索引效率的影响以及不同规模数据下索引的创建、查询和更新试验。与现有分布式空间索引的对比试验及分析结果表明,M-Quadtree索引在数据存储量负载均衡、算法并行化和空间查询效率等方面表现得更好。     

Effect of sensor modelling methods on computation of 3-D coordinates from Cartosat-1 stereo data

The orbital and the rational polynomial coefficients (RPC) models are the two most commonly used models to compute a three-dimensional coordinates from an image stereo-pair. But it is still confusing that with the identical user provided inputs, which one of these two models provides more accurate digital elevation model (DEM), especially for mountainous terrain. This study aimed to find out the answer by evaluating the impact of used models on the vertical accuracy of DEM extracted from Cartosat-1 stereo data. We used high-accuracy photogrammetric DEM as the reference DEM. Apart from general variations in statistics, surprisingly in a few instances, both the DEMs provided contrasting results, thus proving the significance of this study. The computed root mean square errors and linear error at 90% (LE90) were lower in case of RPC DEM for various classes of slope, aspect and land cover, thus suggesting its better relative accuracy.     

基于最优分割分级法的月球撞击坑分级及其演化分析

撞击坑是月球表面广泛分布的重要构造形态,占据了月球表面的大部分面积。撞击坑的直径差别很大,从几微米到数百千米,其退化程度与形成年代具有密切关系。为了研究不同地质年代形成的撞击坑直径大小及其演化规律,需采用量化分级方法对大小不同的撞击坑进行定量分级和统计分析。本文在月表撞击坑数据库LU60645GT和Lunar_Impact_Crater_Database(2011)的基础上,结合数据库中撞击坑的直径、深度和年代信息,利用最优分割分级法对撞击坑直径进行定量化分级,并根据分级结果,综合分析撞击坑几何形态特征及其演化规律。研究结果表明,撞击坑形态特征的演化与年代有密切的关系。在相同级别、相同地体下,撞击坑形成的年代越早,其形态特征的精细结构退化程度越明显,只保留了大体的几何形状;而在不同级别、相同地体、相同年代下的撞击坑形态特征则由简单逐渐变为复杂,坑物质也逐渐变得复杂。     

西天山巴音布鲁克地区早古生代成矿地质环境:岩浆岩及其时代和元素同位素约束

西天山位于中亚造山带(CAOB)西南部,是其重要组成部分。CAOB晚古生代金属成矿环境和过程尤为典型,但早古生代成矿地质环境还不甚清楚。新疆巴音布鲁克地区出露(原定时代晚志留世)巴音布鲁克组火山岩夹浅海相碎屑岩和灰岩,是认识西天山早古生代成矿地质环境的难得对象。巴音布鲁克组出露于Nikolaev-那拉提山北缘断裂与Atbash-Inylchek-那拉提山南缘断裂之间的中天山,在巴音布鲁克地区典型发育,火山岩包括玄武岩、玄武安山岩、英安岩、流纹岩及相应的火山碎屑岩,其中侵入岩有正长斑岩和花岗闪长岩。LA-ICP-MS测得玄武安山岩、英安岩、正长斑岩、花岗闪长岩、流纹岩锆石U-Pb年龄分别为455.6±8.1Ma、444.5±1.9Ma、441.4±1.6Ma、455.4±5.3Ma、424±1.9Ma,岩浆活动于晚奥陶-早志留世,喷出和侵入时代接近,原定巴音布鲁克组地层时代晚志留世应改为晚奥陶-早志留世。这些岩浆岩具有相似的稀土元素地球化学特征,微量元素相比原始地幔均亏损Nb、Ta、P、Ti。玄武岩正的ε(t)=+1.6~+6.7,低的(~(87)~Sr/~(86)NdSr)i=0.70377~0.70489,指示岩浆源区具有亏损地幔特征,弱的Zr-Hf负异常,低的Th/Nb比值,较窄的同位素变化范围暗示地壳混染并不显著,微量元素及铅同位素特征(~(206)Pb/~(204)Pb=18.26~18.77,~(207)Pb/~(204)Pb=15.63~15.69,~(208)Pb/~(204)Pb=38.21~38.34)表明岩浆源区可能是俯冲流体及洋底沉积物交代的地幔楔橄榄岩部分熔融成因。西天山巴音布鲁克地区早古生代岩浆岩应是南天山洋晚奥陶-早志留世向北向中天山陆块之下俯冲在中天山-伊犁板块南缘活动大陆边缘的岩浆产物,指示了陆缘岩浆弧环境。这种陆缘弧环境有利于斑岩铜金成矿系统发育,值得高度关注相关铜金矿的地质找矿。