Groundwater mixing between different aquifer types in a complex structural setting discerned by elemental and stable isotope geochemistry

Karst aquifers are well known for their intricate stratigraphy and geologic structures, which make groundwater characterization challenging because flowpaths and recharge sources are complex and difficult to evaluate. Geochemical data, collected from ten closely spaced production wells constructed in two karst aquifers (Bangor Limestone (Mb) and Tuscumbia Limestone/Fort Payne Chert (Mftp)) in Trussville, north‐central Alabama, illustrate two distinctive groundwater end‐members: (1) higher major ion, dissolved inorganic carbon, conductivity, alkalinity concentrations, heavier δ¹³C ratios (max: −10.2 ± 0.2‰ Vienna Pee Dee Belemnite (PDB)) and lower residence times (mean: 19.5 ± 2 years, n = 2) of groundwater in the Mb aquifer and (2) lower constituent concentrations, lighter δ¹³C ratios (min: −13.4 ± 0.2‰ PDB) and longer residence times of groundwater (mean: 23.6 ± 2 years, n = 4) in the Mftp aquifer. Summer and fall data and the binary mixing model show aquifer inter‐flow mixing along solution fractures and confirms the distinctive groundwater geochemistry of the two aquifers. Lowering of static water levels over the summer (drawdown from 2 to 5.2 m) leads to more reducing groundwater conditions (lower Eh values) and slightly enriched δ¹⁸O and δD ratios during the fall [δ¹⁸O: −4.8 ± 0.1 to −5.4 ± 0.1‰ Vienna Standard Mean Oceanic Water (VSMOW), n = 9; δD: −25.4 ± 1 to −27.4 ± 1‰ VSMOW, n = 9] when compared with summer season samples (δ¹⁸O: −5.1 ± 0.1 to −5.7 ± 0.1‰ VSMOW, n = 11; δD: −25.0 ± 1 to −30.6 ± 1‰ VSMOW, n = 11). GIS analyses confirm the localized origin of recharge to the investigated aquifers. The combination of GIS, field parameters and geochemistry analyses can be successfully used to identify recharge sources, evaluate groundwater flow and transport pathways and to improve understanding of how groundwater withdrawals impact the sustainability and susceptibility to contamination of karst aquifers. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydraulic modeling of flood pulses in the Middle Bussento Karst System (MBSKS), UNESCO Cilento Global Geopark,southern Italy

Karst systems provide water for domestic and industrial uses and for generating hydropower, but they can also create fluvial hazards, such as upstream back‐flooding and downstream karst flash‐flood events. However, these hazards are difficult to foresee due to the complex recharge‐discharge processes as well as the lack of information on the inside of the system, which has often not been completely surveyed by speleologists or explored by boreholes. To overcome these difficulties, hydro‐chemical data from the monitoring system in the Middle Bussento Karst System (MBSKS), one of the first Experimental Karst Systems in southern Italy, were recorded and previously discussed. Based on shared background in flood karst hydraulic modeling, this paper describes the conceptual premises and rationale of a general‐purpose hydraulic model that is suitable both for the MBSKS and for other Mediterranean, multi‐recharge, mature, conduit‐dominated karst systems. To test the reliability of the model, simulations of time–space behavior and response are performed using natural and artificial flood pulses “as tracers”, considering a “pulse” as a significant variation in water quantity and/or quality. The results of the model explain the interactions between allogenic, autogenic, and anthropogenic recharges from differentiated sources and phreatic conduit systems. These results also clarify the overall response of karst springs at typical time scales of flood pulses. Table acronym name     

北斗系统周跳探测与修复方法探究

北斗卫星导航系统是中国自行研制的继美国GPS和俄罗斯GLONASS之后第3个成熟的卫星导航系统。进行GNSS相位数据处理的过程中,由于某种原因,载波相位观测值突然出现整周数的跳变的现象称为周跳。周跳现象的发生将会严重影响后续的数据处理,因此,有必要在解算前完成周跳探测与修复工作。本文就北斗系统周跳探测与修复课题进行了较为全面的探究,先后介绍了高次差法、多普勒法、伪距载波相位组合法、电离层残差法、三频数据组合探测法等一系列周跳探测与修复方法,可分别应用于不同频段接收机、不同观测环境、不同数据处理模式。     

Retrieval of effective leaf area index (LAIe) and leaf area density (LAD) profile at individual tree level using high density multi-return airborne LiDAR

As an important canopy structure indicator, leaf area index (LAI) proved to be of considerable implications for forest ecosystem and ecological studies, and efficient techniques for accurate LAI acquisitions have long been highlighted. Airborne light detection and ranging (LiDAR), often termed as airborne laser scanning (ALS), once was extensively investigated for this task but showed limited performance due to its low sampling density. Now, ALS systems exhibit more competing capacities such as high density and multi-return sampling, and hence, people began to ask the questions like—“can ALS now work better on the task of LAI prediction?” As a re-examination, this study investigated the feasibility of LAI retrievals at the individual tree level based on high density and multi-return ALS, by directly considering the vertical distributions of laser points lying within each tree crown instead of by proposing feature variables such as quantiles involving laser point distribution modes at the plot level. The examination was operated in the case of four tree species (i.e. Picea abies, Pinus sylvestris, Populus tremula and Quercus robur) in a mixed forest, with their LAI-related reference data collected by using static terrestrial laser scanning (TLS). In light of the differences between ALS- and TLS-based LAI characterizations, the methods of voxelization of 3D scattered laser points, effective LAI (LAIe) that does not distinguish branches from canopies and unified cumulative LAI (ucLAI) that is often used to characterize the vertical profiles of crown leaf area densities (LADs) was used; then, the relationships between the ALS- and TLS-derived LAIes were determined, and so did ucLAIs. Tests indicated that the tree-level LAIes for the four tree species can be estimated based on the used airborne LiDAR (R² = 0.07, 0.26, 0.43 and 0.21, respectively) and their ucLAIs can also be derived. Overall, this study has validated the usage of the contemporary high density multi-return airborne LiDARs for LAIe and LAD profile retrievals at the individual tree level, and the contribution are of high potential for advancing forest ecosystem modeling and ecological understanding.     

多光谱遥感影像亮度空间相位一致性特征点检测

提出了一种基于亮度空间和相位一致性理论的多光谱遥感影像特征点检测算法。首先利用参数自适应的灰度变换函数建立影像亮度空间;然后结合相位一致性方法在影像亮度空间进行候选特征点检测,并将候选特征点映射到原始影像上进行非极大值抑制;最后在尺度空间计算特征点的特征尺度值。本文方法有效结合了亮度空间特征检测和相位一致性特征检测的优势,对多光谱遥感影像的辐射变化具有较强的稳健性。试验结果证明,与传统特征点检测算法相比,本文方法在特征重复率和重复特征数量方面都具有明显的优势。     

双频BDS实时精密单点定位算法和精度分析

为了对BDS实时精密单点定位性能进行评估,该文提出了一种适用于BDS系统的实时精密单点定位算法。采用无电离组合模型作为双频实时精密单点定位的数学模型,采用电离层残差法和Melbourne-Wübbena组合实时探测相位周跳,进而单历元实时估计坐标、模糊度等参数,实现了BDS双频实时精密单点定位算法。基于此算法,采用轨道钟差产品和采样间隔为1s的观测数据,模拟实时BDS双频精密单点定位算法,并评估其定位精度。实验结果表明:BDS双频实时定位的平面精度和三维精度均为0.2m左右。     

基于GRNN模型的岩溶地区城市水生态足迹分析与预测

为揭示岩溶地区城市生态足迹的驱动因子与机制,评估和预测水资源可持续利用程度,根据岩溶地区水文地质特征及其与水生态足迹的相关性,分析了贵阳市2002-2014年水生态足迹的变化特征,并利用GRNN模型对其进行了时间序列预测。结果表明,贵阳市水生态足迹由2002年的12.79万hm2上升至2014年的15.74万hm2,整体呈波动上升趋势,其中,水量生态足迹变化特征主要由工业用水量决定,水质生态足迹变化特征主要由氨氮排放量决定。水生态承载力波动变化明显,且小于水生态足迹需水量,导致水资源可持续利用系数介于0.31~0.63之间。对贵阳市2015-2019年水生态足迹进行GRNN预测,表明水生态足迹呈现先升后降得趋势,水生态承载力变幅在4.315~8.038万hm2之间,但水生态足迹与承载力的缺口大,水资源可持续利用系数偏低。因此,在水生态承载力有限的贵阳市,降低水生态足迹才是实现水资源可持续利用的关键。     

长兴岛碳酸盐岩裂隙岩溶水特征分析

长兴岛碳酸盐岩区占整个岛的44%,受降雨、岩性、构造等条件的影响,岩溶较发育,地下水资源较丰富,且水质较好,适宜作为饮用水水源。在收集大量基础地质资料,并结合现场调查、钻探、抽水试验成果,分析了从上至下5层岩溶带的发育特征、赋水性规律及补、径、排条件。同时考虑供水水文地质条件,进行了天然资源和开采资源概算,及成井条件、水质的分析评价。     

A synthetic aperture radar sea surface distribution estimation by n-order Bézier curve and its application in ship detection

To dates,most ship detection approaches for single-pol synthetic aperture radar(SAR) imagery try to ensure a constant false-alarm rate(CFAR).A high performance ship detector relies on two key components:an accurate estimation to a sea surface distribution and a fine designed CFAR algorithm.First,a novel nonparametric sea surface distribution estimation method is developed based on n-order Bézier curve.To estimate the sea surface distribution using n-order Bézier curve,an explicit analytical solution is derived based on a least square optimization,and the optimal selection also is presented to two essential parameters,the order n of Bézier curve and the number m of sample points.Next,to validate the ship detection performance of the estimated sea surface distribution,the estimated sea surface distribution by n-order Bézier curve is combined with a cell averaging CFAR(CA-CFAR).To eliminate the possible interfering ship targets in background window,an improved automatic censoring method is applied.Comprehensive experiments prove that in terms of sea surface estimation performance,the proposed method is as good as a traditional nonparametric Parzen window kernel method,and in most cases,outperforms two widely used parametric methods,K and G0 models.In terms of computation speed,a major advantage of the proposed estimation method is the time consuming only depended on the number m of sample points while independent of imagery size,which makes it can achieve a significant speed improvement to the Parzen window kernel method,and in some cases,it is even faster than two parametric methods.In terms of ship detection performance,the experiments show that the ship detector which constructed by the proposed sea surface distribution model and the given CA-CFAR algorithm has wide adaptability to different SAR sensors,resolutions and sea surface homogeneities and obtains a leading performance on the test dataset.     

Seismic attenuation estimation using a complete ensemble empirical mode decomposition-based method

The frequency attenuation gradient method can provide important information for hydrocarbon detection. In this paper, a method using Complete Ensemble Empirical Mode Decomposition (CEEMD), Hilbert transform and the least-squares curve-fitting is proposed for seismic attenuation estimation as an effective frequency attenuation gradient estimation approach. We first use CEEMD to obtain the different Intrinsic Mode Functions (IMFs), which have a narrow band and can enhance the physical meaning of instantaneous attributes trace by trace. The time-frequency spectrum, which is computed using a Hilbert transform of each IMF, is represented as a spectrum with a single-peak that has narrow side lobes, which is conducive to frequency attenuation gradient estimation. Second, for each time sample, the frequency-amplitude spectrum of each IMF trace is extracted from the time-frequency spectrum to conduct the attenuation gradient computation. Then, the logarithm operation is performed for each IMF trace. Due to the very narrow bands of some IMFs in some seismic traces, a variable frequency window is adopted along the IMF trace according to the local data characteristics. Finally, the attenuation gradient for each IMF in a seismic trace can be computed using least-squares fitting. A different IMF reflects a seismic trace with a different spatiotemporal scale and can highlight different geologic and stratigraphic information. The correlation weighted average operation is used to highlight some useful details in seismic trace and obtains the attenuation gradient for each seismic trace. Field data examples demonstrate our method and its effectiveness. The proposed method can stably estimate the frequency attenuation gradient.