Glacial dispersal studies using indicator minerals and till geochemistry around two eastern Finland kimberlites

Diamondiferous kimberlites occur in eastern Finland, in the areas of Kaavi–Kuopio and Kuhmo. Active diamond exploration has been ongoing in the country for over two decades, but the Karelian craton still remains under explored given its size and potential. In order to develop techniques that can be applied to diamond exploration in glaciated terrains, the Geological Survey of Finland (GTK) carried out a detailed heavy mineral and geochemical survey of Quaternary till in 2001–2003 around two of the known kimberlitic bodies in Finland, Pipe 7 in Kaavi and Dyke 16 in Kuhmo. The mineralogical and geochemical signatures of these two kimberlites were studied in the basal till deposited down-ice from the targets. The kimberlites were selected to represent two different types in terms of shape, size, age and petrology, as well as showing contrasting country rocks and Quaternary deposits. Till samples up to 60 kg in weight were taken by excavator and by drill rig. Kimberlitic indicator mineral grains (0.25–1.0 mm) were concentrated using a GTK modified 3″Knelson Concentrator. Fine fractions (< 0.063 mm) of selected samples were analyzed by XRF and ICP-MS. The indicator grains down-ice from Pipe 7 form a well-defined fan in the basal till that can be followed for at least 2 km with a maximum concentration at 1.2 km distance from the pipe. Another kimberlitic body discovered during the study 300 m down-ice from Pipe 7 demonstrates that there are in fact at least two superimposed indicator fans. The results do not rule out the possibility of even more undiscovered kimberlitic sources in the area. In contrast, the indicator dispersal trail from Dyke 16 is shorter (1 km) and less well-defined than that at Kaavi, mainly due to the lower indicator content in the kimberlite itself and subsequently in till, as well as a large population of background chromites in till. The latter population is likely having been derived from the Archean Näätäniemi serpentinite massif and the associated ultramafic metavolcanics of the Kuhmo greenstone belt, located ca. 30 km up-ice from the sampling area. The indicator maximum at Seitaperä dyke swarm occurs immediately down-ice from the kimberlite, after which the concentration drops rapidly. Results of this study contribute to the overall understanding of the Quaternary history of the Kaavi and Kuhmo areas, and more importantly, provide key information to diamond exploration in these particular regions and also elsewhere in glaciated terrains.     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Likely effects of climate change on groundwater availability in a Mediterranean region of Southeastern Spain

Groundwater resources are typically the main fresh water source in arid and semi‐arid regions. Natural recharge of aquifers is mainly based on precipitation; however, only heavy precipitation events (HPEs) are expected to produce appreciable aquifer recharge in these environments. In this work, we used daily precipitation and monthly water level time series from different locations over a Mediterranean region of Southeastern Spain to identify the critical threshold value to define HPEs that lead to appreciable aquifer recharge in this region. Wavelet and trend analyses were used to study the changes in the temporal distribution of the chosen HPEs (≥20 mm day^?1) over the observed period 1953–2012 and its projected evolution by using 18 downscaled climate projections over the projected period 2040–2099. The used precipitation time series were grouped in 10 clusters according to similarities between them assessed by using Pearson correlations. Results showed that the critical HPE threshold for the study area is 20 mm day^?1. Wavelet analysis showed that observed significant seasonal and annual peaks in global wavelet spectrum in the first sub‐period (1953–1982) are no longer significant in the second sub‐period (1983–2012) in the major part of the ten clusters. This change is because of the reduction of the mean HPEs number, which showed a negative trend over the observed period in nine clusters and was significant in five of them. However, the mean size of HPEs showed a positive trend in six clusters. A similar tendency of change is expected over the projected period. The expected reduction of the mean HPEs number is two times higher under the high climate scenario (RCP8.5) than under the moderate scenario (RCP4.5). The mean size of these events is expected to increase under the two scenarios. The groundwater availability will be affected by the reduction of HPE number which will increase the length of no aquifer recharge periods (NARP) accentuating the groundwater drought in the region. Copyright © 2016 John Wiley & Sons, Ltd.     

Measurement of stream cross section using ground penetration radar with Hilbert–Huang transform

This study presents a new method to measure stream cross section without having contact with water. Compared with conventional measurement methods which apply instruments such as sounding weight, ground penetration radar (GPR), used in this study, is a non‐contact measurement method. This non‐contact measurement method can reduce the risk to hydrologists when they are conducting measurements, particularly in high flow period. However, the original signals obtained by using GPR are very complex, different from studies in the past where the measured data were mostly interpreted by experts with special skill or knowledge of GPR so that the results obtained were less objective. This study employs Hilbert–Huang transform (HHT) to process GPR signals which are difficult to interpret by hydrologists. HHT is a newly developed signal processing method that can not only process the nonlinear and non‐stationary complex signals, but also maintain the physical significance of the signal itself. Using GPR with HHT, this study establishes a non‐contact stream cross‐section measurement method with the ability to measure stream cross‐sectional areas precisely and quickly. Also, in comparison with the conventional method, no significant difference in results is found to exist between the two methods, but the new method can considerably reduce risk, measurement time, and manpower. It is proven that the non‐contact method combining GPR with HHT is applicable to quickly and accurately measure stream cross section. Copyright © 2013 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

基于遥感的GEF海河项目实施效果评价

GEF海河项目于2011年完工,该项目在华北地区16个县市实施,惠及2 000多万人。本文从ET、地下水超采量削减、污染物排放量控制、地表水水资源量与利用量以及出入境水质水量5个方面对GEF海河项目实施效果进行了评价。结果表明,全流域ET值减少了68.4 mm,COD和NH3-N入河总量分别减少了36.8%和8.5%、入海总量分别减少了46.1%和36.4%,地表水资源量增加了8.94%,总用水量减少了1.99×108m3,生态环境用水量增加了7.77×108m3,河湖环境明显改善。     

EMD-Wavelet-ICA耦合模型及其在GPS坐标序列降噪中的应用

在分析比较经验模态分解(EMD)、小波变换(Wavelet)和独立分量分析(ICA)优缺点的基础上,提出一种新的EMD-Wavelet-ICA耦合模型。该模型充分利用了EMD的自适应性,对原始信号进行分解获得不同频率的模态函数(IMF),采用标准化模量的累计均值对IMF进行尺度划分;进而分别采用Wavelet和ICA对高频和低频IMF进行降噪,将降噪后的IMF进行多尺度重构,获得降噪后的信号;采用信噪比、标准差、偏差和相关系数等指标对降噪效果进行评价。仿真数据和GPS坐标序列的处理结果表明:与EMD模型和EMD-ICA模型相比,新模型的标准差、偏差均有不同程度的减小;信噪比和相关系数有一定程度的增大,可以获得更好的降噪效果。     

热红外与可见光数据的云边缘高度匹配

高分辨率数据中云高度的差异性突显,特别是边缘处高度在云阴影识别和地表辐射估算等方面成为需要考虑的重要因素。热红外数据获取云高度分辨率较低、缺乏细部差异性特征,为解决这一问题,首先将对应的热红外和可见光数据进行特征点配准,再将基于热红外数据计算的云高度重采样至高分辨率,然后以基于欧式距离变换的围线搜索方法及距离加权将热红外云边缘高度匹配至对应的可见光图像,最后根据云阴影的相似度匹配方法确定真实云高度。结果表明,算法在遵循热红外云高信息分布变化规律的同时,可以得到较准确的高分辨率云边缘高度,一定程度上解决了热红外技术获取云高在分辨率上的局限,扩展了其在云高反演方面的作用。