用叠加谱比法研究新疆地区Lg尾波Q值的区域性变化

运用叠加谱比法,对分布于新疆地区的5个台站所记录的51个地震数据进行分析,得出与各射线路径相对应的Lg尾波Q0(Lg尾波在1Hz处的Q)值。结果表明:Q0值呈现出明显的区域性变化特点,与构造活动关系紧密。在较为稳定的塔里木地台,Q0呈现出高值,在350～450之间。在南部构造活动较为激烈的昆仑褶皱系,Q0呈现出较低的值,约为200～250。在构造活动最为激烈的帕米尔地区,Q0值在170～200之间。在北部的天山褶皱系,Q0值约为220～270。准葛尔褶皱系比天山褶皱系的Q0值略高,约为260～290。     

Polymetamorphism in garnet micaschists of the Saualpe Eclogite Unit (Eastern Alps,Austria), resolved by automated SEM methods and EMP–Th–U–Pb monazite dating

Polymetamorphic garnet micaschists from the Austroalpine Saualpe Eclogite Unit (Kärnten, Austria, Eastern Alps) display complex microstructural and mineral–chemical relationships. Automated scanning electron microscopy routines with energy dispersive X‐ray (EDX) spectral mapping were applied for monazite detection and garnet mineral–chemical characterization. When the Fe, Mg, Mn and Ca element wt% compositions are used as generic labels for garnet EDX spectra, complex zonations and porphyroblast generations can be resolved in complete thin sections for selective electron‐microprobe analyses. Two garnet porphyroblast generations and diverse monazite age populations have been revealed in low‐Ca and high‐Al‐metapelites. Garnet 1 has decreasing Mn, constant Ca and significantly increasing Mg from cores to rims. Geothermobarometry of garnet 1 assemblages signals a crystallization along a M1 prograde metamorphism at ~650 °C/6–8 kbar. Sporadic monazite 1 crystallization started at c. 320 Ma. Subsequent pervasive 300–250 Ma high‐Y and high‐Gd monazite 1 formation during decompression coincided with the intrusion of Permian and Early Triassic pegmatites. Monazite 1 crystallized along the margin of garnet 1. Coronas of apatite and allanite around the large 320–250 Ma monazite signal a retrogressive stage. These microstructures suggest a Carboniferous‐to‐Early‐Permian age for the prograde M1 event with garnet 1. Such a M1 event at an intermediate‐P/T gradient has not yet been described from the Saualpe, and preceded a Permo‐Triassic low‐P stage. The M2 event with garnet 2 postdates the corona formation around Permian monazite. Garnet 2 displays first increasing X_Ca at decreasing X_Mg, then increasing X_Ca and X_Mg, and finally decreasing X_Ca with increasing X_Mg, always at high Ca and Mg, and low Mn. This records a P–T evolution which passed through eclogite facies conditions and reached maximum temperatures at ~750 °C/14 kbar during decompression‐heating. A monazite 2 population (94–86 Ma) with lower Y and Gd contents crystallized at decreasing pressure during the Cretaceous (Eo‐Alpine) metamorphism M2 at a high‐P/T gradient. The Saualpe Eclogite Unit underwent two distinct clockwise metamorphic cycles at different P–T conditions, related to continental collisions under different thermal regimes. This led to a characteristic distribution pattern of monazite ages in this unit which is different from other Austroalpine basement areas.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Sensitivity of a hydraulic model to channel erosion uncertainty during extreme flooding

Recent research into flood modelling has primarily concentrated on the simulation of inundation flow without considering the influences of channel morphology. River channels are often represented by a simplified geometry that is implicitly assumed to remain unchanged during flood simulations. However, field evidence demonstrates that significant morphological changes can occur during floods to mobilize the boundary sediments. Despite this, the effect of channel morphology on model results has been largely unexplored. To address this issue, the impact of channel cross‐section geometry and channel long‐profile variability on flood dynamics is examined using an ensemble of a 1D–2D hydraulic model (LISFLOOD‐FP) of the ~1 : 2000 year recurrence interval floods in Cockermouth, UK, within an uncertainty framework. A series of simulated scenarios of channel erosional changes were constructed on the basis of a simple velocity‐based model of critical entrainment. A Monte‐Carlo simulation framework was used to quantify the effects of this channel morphology together with variations in the channel and floodplain roughness coefficients, grain size characteristics and critical shear stress on measures of flood inundation. The results showed that the bed elevation modifications generated by the simplistic equations reflected an approximation of the observed patterns of spatial erosion that enveloped observed erosion depths. The effect of uncertainty on channel long‐profile variability only affected the local flood dynamics and did not significantly affect the friction sensitivity and flood inundation mapping. The results imply that hydraulic models generally do not need to account for within event morphodynamic changes of the type and magnitude of event modelled, as these have a negligible impact that is smaller than other uncertainties, e.g. boundary conditions. Instead, morphodynamic change needs to happen over a series of events to become large enough to change the hydrodynamics of floods in supply limited gravel‐bed rivers such as the one used in this research. Copyright © 2014 John Wiley & Sons, Ltd.     

Integrating isolated and riparian wetland modules in the PHYSITEL/HYDROTEL modelling platform: model performance and diagnosis

Mathematical modelling is a well‐accepted framework to evaluate the effects of wetlands on stream flow and watershed hydrology in general. Although the integration of wetland modules into a distributed hydrological model represents a cost‐effective way to make this assessment, the added value brought by landscape‐specific modules to a model's ability to replicate basic hydrograph characteristics remains unclear. The objectives of this paper were the following: (i) to present the adaptation of PHYSITEL (a geographic information system) to parameterize isolated and riparian wetlands; (ii) to describe the integration of specific isolated wetland and riparian wetland modules into HYDROTEL, a distributed hydrological model; and (iii) to evaluate the performance of the updated modelling platform with respect to the capacity of replicating various hydrograph characteristics. To achieve this, two sets of simulations were performed (with and without wetland modules), and the added value was assessed at three river segments of the Becancour River watershed, Quebec, Canada, using six general goodness‐of‐fit indicators and 14 water flow criteria. A sensitivity analysis of the wetland module parameters was performed to characterize their impact on stream flows of the modelled watershed. Results of this study indicate the following: (i) integration of specific wetland modules can slightly increase the capacity of HYDROTEL to replicate basic hydrograph characteristics; and (ii) the updated modelling platform allows for the explicit assessment of the impact of wetlands (e.g. typology and location) on watershed hydrology. Copyright © 2015 John Wiley & Sons, Ltd.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

基于光谱角匹配算法的水体信息提取研究

针对水体信息提取时,传统算法不能较好地解决山体阴影被误提为水体的问题,引入多用于高光谱数据处理的光谱角匹配算法,以Land-sat8 OLI(band1-band7)为数据源,以黄河小浪底水库周边区域为实验区,开展水体信息提取研究,结果显示:水体指数法和波段关系提取的结果中有大量的山体阴影信息,提取精度不超过30%,而光谱角匹配法提取结果受阴影的影响较小,提取精度在99%以上。     

高分二号卫星影像融合方法探析

图像融合技术是卫星影像处理过程中的一个重要步骤,它可以将全色影像的高空间分辨率特点与多光谱影像的高光谱分辨率特点有机地结合起来,使得融合后的影像更利于后续解译。由于融合处理方法的原理以及获取影像的传感器特性不同,不同的卫星影像适合于不同的融合方法。针对高分二号卫星影像的特点,本文采用Gram-Schmidt融合法、FIHS变换方法、HPF方法、Pansharp方法进行试验。实验结果表明:HPF方法所得融合结果在光谱特性保持方法性能最优,而Pansharp方法融合结果在信息量的保持以及清晰度方面效果更优,因此,在实际应用中,应根据具体需求进行方法选择。     

Characterisation of carbonate minerals from hyperspectral TIR scanning using features at 14 000 and 11 300 nm

Rapid characterisation of carbonate phases in hyperspectral reflectance spectra acquired from drill core material has important implications for mineral exploration and resource modelling. Major infrared active features of carbonates lie in the thermal region around 6500 nm, 11 300 nm and 14 000 nm, with the latter two features being most useful for differentiating mineral species. A scatter diagram of the wavelength of the 14 000 nm feature vs that of the 11 300 nm feature, powerfully differentiates carbonates. Although the wavelength of the 11 300 nm peak is easily measured, the 14 000 nm trough and peak are commonly weak and their wavelengths can confidently be used only after filtering the spectra, e.g. selecting only those with the trough and peak separated by 175–230 nm, typical of common carbonates. The method is demonstrated with drillhole 120R from the Rosebery polymetallic VHMS deposit in western Tasmania, which has been scanned with the HyLogger-3 system. A 14 000–11 300 plot shows a high degree of clustering of the drillhole 120R data close to the library spectra of calcite, dolomite, Fe-dolomite, ankerite, kutnohorite, rhodochrosite, Fe-rhodochrosite and siderite. The interpreted compositions of the carbonate spectral populations strongly correlate with the chemical populations of 144 analysed carbonates and provide a highly resolved spatial framework for interpreting carbonate alteration.     

顾及周期项误差和起点偏差修正的北斗卫星钟差预报

根据北斗卫星导航系统星载原子钟自身的物理特性,采用武汉大学IGS数据中心发布的2016年1月1日至2016年11月1日共306天的事后钟差产品进行谱分析。分析结果表明:北斗卫星导航系统的3类卫星钟都存在一定的周期特性;其中GEO和IGSO卫星钟的主周期相对较为明显;GEO卫星钟的主周期依次为12、24、8和6h;IGSO的主周期为24、12、8和6h;而MEO的3个主周期为12.9、6.4和24h。依据各类原子钟的周期特性,同时对各天之间钟差的起始点偏差进行修正,并利用修正模型对北斗卫星钟差进行预报和精度分析。试验结果表明,改进的预报模型能显著提升北斗卫星的钟差预报精度,北斗卫星钟差24h、12h、6h的平均预报精度分别能达到6.55ns、3.17ns和1.76ns。