Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Exploring the interaction of surface roughness and slope gradient in controlling rates of soil loss from sloping farmland on the Loess Plateau of China

Surface roughness and slope gradient are two important factors influencing soil erosion. The objective of this study was to investigate the interaction of surface roughness and slope gradient in controlling soil loss from sloping farmland due to water erosion on the Loess Plateau, China. Following the surface features of sloping farmland in the plateau region, we manually prepared rough surfaces using four tillage practices (contour drilling, artificial digging, manual hoeing, and contour plowing), with a smooth surface as the control measure. Five slope gradients (3°, 5°, 10°, 15°, and 20°) and two rainfall intensities (60 and 90 mm/hr) were considered in the artificial rainfall simulation experiment. The results showed that the runoff volume and sediment yield increased with increasing slope gradient under the same tillage treatment. At gentle slope gradients (e.g., 3° and 5°), the increase in surface roughness prevented the runoff and sediment production, that is, the surface roughness reduced the positive effect of slope gradient on the runoff volume and sediment yield to a certain extent. At steep slope gradients, however, the enhancing effect of slope gradient on soil erosion gradually increased and surpassed the reduction effect of surface roughness. This study reveals the existence of a critical slope gradient that influences the interaction of surface roughness and slope gradient in controlling soil erosion on sloping farmland. If the slope gradient is equal to or less than the critical value, an increase in surface roughness would decrease soil erosion. Otherwise, the increase in surface roughness would be ineffective for preventing soil erosion. The critical slope gradient would be smaller under higher rainfall intensity. These findings are helpful for us to understand the process of soil erosion and relevant for supporting soil and water conservation in the Loess Plateau region of China.     

Geochemical and isotopic characterization of trace fossil infillings: New insights on tracemaker activity after the K/Pg impact event

Geochemical and isotopic analyses of the Cretaceous–Paleogene (K/Pg) boundary deposits were conducted at the Caravaca section (External Subbetic, southeast of Spain) in order to evaluate the recovery of the macrobenthic tracemaker community and the bioturbational disturbance. Samples from the infilling material of several lower Danian dark-colored trace fossils (Chondrites, Planolites, Thalassinoides and Zoophycos) located in the uppermost 8-cm of the light upper Maastrichtian strata, as well as samples from the host sedimentary rock of these trace fossils, were analyzed and compared with data from the lower Danian deposits. The values of element ratios indicative of extraterrestrial contamination (Cr/Al, Co/Al and Ni/Al) are higher in the infilling trace fossil material than in the upper Maastrichtian and lower Danian deposits, which suggests a contribution of the ejecta layer. Regarding the isotope composition, the δ¹³C values are lower in the infilling material than in the Maastrichtian host sedimentary rocks surrounding the traces, while the δ¹⁸O are higher in the infilling material. The geochemical and isotopic compositions of the infilling material evidence the unconsolidated character of the sediment, including the red boundary layer. Softground conditions confirm a relatively rapid recovery by the macrobenthic tracemaker community, starting a few millimeters above the K/Pg boundary layer. The mixture of the infilling material of the trace fossils moreover reveals a significant macrobenthic tracemaker activity affecting K–Pg boundary transition sediments that may have significantly altered original signatures.     

DEM simulations of methane hydrate exploitation by thermal recovery and depressurization methods

Methane hydrate (MH, also called fiery ice) exists in forms of pore filling, cementing and load-bearing skeleton in the methane hydrate bearing sediment (MHBS) and affects its mechanical behavior greatly. To study the changes of macro-scale and micro-scale mechanical behaviors of MHBS during exploitation by thermal recovery and depressurization methods, a novel 2D thermo-hydro-mechanical bonded contact model was proposed and implemented into a platform of distinct element method (DEM), PFC2D. MHBS samples were first biaxially compressed to different deviator stress levels to model different in-situ stress conditions. With the deviator stress maintained at constant, the temperature was then raised to simulate the thermal recovery process or the pore water pressure (i.e. confining pressure for MH bond) was decreased to simulate the depressurization process. DEM simulation results showed that: during exploitation, the axial strain increased with the increase of temperature (in the thermal recovery method) or decrease of pore water pressure (in the depressurization method); sample collapsed during MH dissociation if the deviator stress applied was larger than the compression strength of a pure host sand sample; sample experienced volume contraction but its void ratio was slightly larger than the pure host sand sample at the same axial strain throughout the test. By comparison with the laboratory test results, the new model was validated to be capable of reproducing the exploitation process by thermal recovery and depressurization methods. In addition, some micro-scale parameters, such as contact distribution, bond distribution, and averaged pure rotation rate, were also analyzed to investigate their relationships with the macroscopic responses.     

利用引力梯度数据计算径向梯度的优化方法

根据重力梯度观测各分量的方差及协方差信息,提出了利用GOCE梯度数据计算径向重力梯度的优化方法。首先给出了径向重力梯度的计算方法,并深入分析了误差传播规律,通过建立相应的条件极值问题,给出了计算径向重力梯度最优组合因子的方法;通过模拟数据验证了本文所提出的优化因子的优越性。实际数据计算表明:相对于传统方法,采用优化组合因子可使反演所得引力位模型的累积大地水准面精度在250阶时提高约2 cm。由于径向重力梯度不仅可以用于地球引力场模型的求解,也可直接应用于地球物理问题的讨论,因此本文所提出的优化方法也可对部分地球动力学问题的讨论提供方便。     

云蒙湖流域不同坡度的土壤侵蚀特征

针对定量分析土壤侵蚀在各坡度等级上的空间分布研究较少的现状,该文选用通用的土壤流失预报方程,对云蒙湖流域1986—2010年间的土壤水力侵蚀状况进行了定量的估算,以探讨不同坡度上的土壤侵蚀特征,并进一步分析了土壤侵蚀变化与人类活动的关系。分析得出:土壤侵蚀强度发生在人类活动比较频繁的区域上(8~25°坡度)更为严重;2010年比1986年强度以上所占比例在15°坡度等级上相对更低,在15°坡度等级上有所增加;云蒙湖流域主要土壤侵蚀量发生在25°坡度上;2010年比1986年耕地面积减少、林地和居民用地面积增加是土壤侵蚀降低的主要因素。     

Reactive transport impacts on recovered freshwater quality during multiple partially penetrating wells (MPPW-)ASR in a brackish heterogeneous aquifer

The use of multiple partially penetrating wells (MPPW) during aquifer storage and recovery (ASR) in brackish aquifers can significantly improve the recovery efficiency (RE) of unmixed injected water. The water quality changes by reactive transport processes in a field MPPW-ASR system and their impact on RE were analyzed. The oxic freshwater injected in the deepest of four wells was continuously enriched with sodium (Na⁺) and other dominant cations from the brackish groundwater due to cation exchange by repeating cycles of ‘freshening’. During recovery periods, the breakthrough of Na⁺ was retarded in the deeper and central parts of the aquifer by ‘salinization’. Cation exchange can therefore either increase or decrease the RE of MPPW-ASR compared to the RE based on conservative Cl⁻, depending on the maximum limits set for Na⁺, the aquifer’s cation exchange capacity, and the native groundwater and injected water composition. Dissolution of Fe and Mn-containing carbonates was stimulated by acidifying oxidation reactions, involving adsorbed Fe²⁺ and Mn²⁺ and pyrite in the pyrite-rich deeper aquifer sections. Fe²⁺ and Mn²⁺ remained mobile in anoxic water upon approaching the recovery proximal zone, where Fe²⁺ precipitated via MnO₂ reduction, resulting in a dominating Mn²⁺ contamination. Recovery of Mn²⁺ and Fe²⁺ was counteracted by frequent injections of oxygen-rich water via the recovering well to form Fe and Mn-precipitates and increase sorption. The MPPW-ASR strategy exposes a much larger part of the injected water to the deeper geochemical units first, which may therefore control the mobilization of undesired elements during MPPW-ASR, rather than the average geochemical composition of the target aquifer.     

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.     

反渗透海水淡化技术现状和展望

反渗透膜技术于20世纪60年代取得突破性进展,促使反渗透海水淡化在近50年间高速发展,淡化产能自20世纪90年代起激增.海水反渗透(SWRO)淡化已成为目前投资最省、成本最低的利用海水制备饮用水的过程.文中主要对海水反渗透淡化的发展状况进行了介绍,如膜和组件的改进,关键设备高压泵和能量回收装置效率的提高,多种工艺过程的不断发展,包括预处理和后处理的新工艺,以及对环境的影响和相应对策等.反渗透技术的发展也推动了其他膜分离技术的进步,并扩展其应用领域.预计在不久的将来,膜技术在海水淡化和水再利用、扩大和保护水资源、循环经济、清洁生产、改造传统产业、节能减排及提高人民生活水平等方面发挥的作用会越来越显著.     

Possibilities of land use change analysis in a mountainous rural area: a methodological approach

Decreasing population density is a current trend in the European Union, and causes a lower environmental impact on the landscape. However, besides the desirable effect on the regeneration processes of semi-natural forest ecosystems, the lack of traditional management techniques can also lead to detrimental ecological processes. In this study we investigated the land use pattern changes in a micro-region (in North-Eastern Hungary) between 1952 and 2005, based on vectorised land use data from archive aerial photos. We also evaluated the methodology of comparisons using GIS methods, fuzzy sets and landscape metrics. We found that both GIS methods and statistical analysis of landscape metrics resulted in more or less the same findings. Differences were not as relevant as was expected considering the general tendencies of the past 60 years in Hungary. The change in the annual rate of forest recovery was 0.12%; settlements extended their area by an annual rate of 3.04%, while grasslands and arable lands had a net loss in their area within the studied period (0.60% and 0.89%, respectively). The kappa index showed a smaller similarity (~60%) between these dates but the fuzzy kappa and the aggregation index, taking into account both spatial and thematic errors, gave a more reliable result (~70–80% similarity). Landscape metrics on patch and class level ensured the possibility of a detailed analysis. We arrived at a similar outcome but were able to verify all the calculations through statistical tests. With this approach we were able to reveal significant (p < 0.05) changes; however, effect sizes did not show large magnitudes. Comparing the methods of revealing landscape change, the approach of landscape metrics was the most effective approach, as it was independent of spatial errors and ensuring a multiple way of interpretation.