Understory and small trees contribute importantly to stemflow of a lower montane cloud forest

Stemflow (Sf) measurements in tropical rain and montane forests dominated by large trees rarely include the understory and small trees. In this study, contributions of lower (1‐ to 2‐m height) and upper (>2‐m height and <5‐cm diameter at breast height [DBH]) woody understory, small trees (5 < DBH < 10 cm), and canopy trees (>10‐cm DBH) to Sf per unit ground area (Sf_a) of a Mexican lower montane cloud forest were quantified for 32 days with rainfall (P) during the 2014 wet season. Rainfall, stemflow yield (Sf_y), vegetation height, density, and basal area were measured. Subsequently, stemflow funneling ratios (SFRs) were calculated, and three common methods to scale up Sf_y from individual trees to the stand level (tree‐Sf_y correlation, P‐Sf_y correlation, and mean‐Sf_y extrapolation) were used to calculate Sf_a. Understory woody plants, small trees, and upper canopy trees represented 96%, 2%, and 2%, respectively, of the total density. Upper canopy trees had the lowest SFRs (1.6 ± 0.5 Standard Error (SE) on average), although the lower understory had the highest (36.1 ± 6.4). Small trees and upper understory presented similar SFRs (22.9 ± 5.4 and 20.2 ± 3.9, respectively). Different Sf scaling methods generally yielded similar results. Overall Sf_a during the study period was 22.7 mm (4.5% of rainfall), to which the understory contributed 70.1% (15.9 mm), small trees 10.6% (2.4 mm), and upper canopy trees 19.3% (4.4 mm). Our results strongly suggest that for humid tropical forests with dense understory of woody plants and small trees, Sf of these groups should be measured to avoid an underestimation of overall Sf at the stand level.     

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

A new approach to estimate canopy evaporation and canopy interception capacity from evapotranspiration and sap flow measurements during and following wetting

Canopy interception and its evaporation into the atmosphere during irrigation or a rainfall event are important in irrigation scheduling, but are challenging to estimate using conventional methods. This study introduces a new approach to estimate the canopy interception from measurements of actual total evapotranspiration (ET) using eddy covariance and estimation of the transpiration from measurements of sap flow. The measurements were conducted over a small‐scale sprinkler‐irrigated cotton field before, during and after sprinkler irrigation. Evaporation and sap flow dynamics during irrigation show that the total ET during irrigation increased significantly because of the evaporation of free intercepted water while transpiration was suppressed almost completely. The difference between actual ET and transpiration (sap flow) during and immediately following irrigation (post irrigation) represents the total canopy evaporation while the canopy interception capacity was calculated as the difference between actual ET and transpiration (sap flow) during drying (post irrigation) following cessation of the irrigation. The canopy evaporation of cotton canopy was calculated as 0.8 mm, and the interception capacity was estimated to be 0.31 mm of water. The measurement uncertainty in both the non‐dimensional ET and non‐dimensional sap flow was shown to be very low. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

运用ENVI实现SPOT-5卫星影像的RVI提取

运用ENVI软件处理长沙地区的SPOT-5影像,将全色影像和多波段影像进行融合,分析融合后的影像,再选择感兴趣区域进行影像裁剪;采用ENVI的波段运算(Band3/Band2)提取该地区比值植被指数并进行分析。结果表明:比值植被指数对植被覆盖率较高的区域有非常高的敏感度;在植被覆盖率小于50%的时候,其敏感度明显降低,在绿色植被覆盖区域的比值植被指数远大于1,在裸露地表、建筑物、水体等没有植被覆盖区域的比值植被指数则在1附近,而比值植被指数大于2的则出现在高覆盖健康绿色植被上。     

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

针对水体信息提取时,传统算法不能较好地解决山体阴影被误提为水体的问题,引入多用于高光谱数据处理的光谱角匹配算法,以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。     

多源卫星重力场模型精度与频谱分析

针对新一代卫星重力探测技术对地球重力场的频谱贡献问题,该文提出了一种基于GPS/水准数据获取多源卫星重力场模型频谱变化特征的方法。采用GPS/水准外符合检验,有效分析评估了多源卫星重力场模型在中国东、西部地区的精度水平。研究结果表明,以CHAMP、GRACE和GOCE卫星为代表的高-低卫星跟踪卫星、低-低卫星跟踪卫星和卫星重力梯度技术,对地球重力场的频谱贡献分别集中在600km以上的长波和中长波、300km以上的中波、200~350km之间的中短波部分。