Distribution characteristics and controlling factors of typical heavy metals in Huanghe River estuary,China

Journal of Oceanology and Limnology - The geochemical characteristics and potential controlling factors of colloidal Zn, Cd, and Pb in Huanghe (Yellow) River estuary (HRE), China were investigated....     

Direct Conservative Domain in the Continuous Galerkin Method for Groundwater Models

The continuous Galerkin finite element method is commonly considered locally nonconservative because a single element with fluxes computed directly from its potential distribution is unable to conserve its mass and fluxes across edges that are discontinuous. Some literature sources have demonstrated that the continuous Galerkin method can be locally conservative with postprocessed fluxes. This paper proposes the concept of a direct conservative domain (DCD), which could conserve mass when fluxes are computed directly from the potential distribution. Also presented here is a method for modifying the advection fluxes to obtain different conservative domains from the DCDs. Furthermore, DCDs are used to analyze the local conservation of several postprocessing algorithms, for which DCDs provide the theoretical basis. The local conservation of DCDs and the proposed method are illustrated and verified by using a hypothetical 2‐D model.     

A regime shift in sediment export from a coastal watershed during a record wet winter,California: Implications for landscape response to hydroclimatic extremes

Small, steep watersheds are prolific sediment sources from which sediment flux is highly sensitive to climatic changes. Storm intensity and frequency are widely expected to increase during the 21st century, and so assessing the response of small, steep watersheds to extreme rainfall is essential to understanding landscape response to climate change. During record winter rainfall in 2016–2017, the San Lorenzo River, coastal California, had nine flow peaks representing 2–10‐year flood magnitudes. By the third flood, fluvial suspended sediment showed a regime shift to greater and coarser sediment supply, coincident with numerous landslides in the watershed. Even with no singular catastrophic flood, these flows exported more than half as much sediment as had a 100‐year flood 35 years earlier, substantially enlarging the nearshore delta. Annual sediment load in 2017 was an order of magnitude greater than during an average‐rainfall year, and 500‐fold greater than in a recent drought. These anomalous sediment inputs are critical to the coastal littoral system, delivering enough sediment, sometimes over only a few days, to maintain beaches for several years. Future projections of megadroughts punctuated by major atmospheric‐river storm activity suggest that interannual sediment‐yield variations will become more extreme than today in the western USA, with potential consequences for coastal management, ecosystems, and water‐storage capacity. The occurrence of two years with major sediment export over the past 35 years that were not associated with extremes of the El Niño Southern Oscillation or Pacific Decadal Oscillation suggests caution in interpreting climatic signals from marine sedimentary deposits derived from small, steep, coastal watersheds, to avoid misinterpreting the frequencies of those cycles. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.     

Deterioration processes affecting prehistoric rock art engravings in granite in NW Spain

Prehistoric rock art sites are endangered despite conservation efforts. The lack of scientific documentation regarding weathering agents affecting rock art and the absence of specific diagnostic protocols hinder the development of conservation strategies. The aim of this research was to investigate active deterioration processes in a granite petroglyph site located in Mougás (Galicia, NW Spain) by characterizing the granite, conducting a geotechnical study of the outcrop and describing and analysing the main weathering processes. Two main deterioration factors were identified. First, water favours block disjunction at the massif scale and causes pitting and surface erosion at the millimetre scale that affects the readability of the engravings. Second, high temperatures associated with wildfires cause mineral transformations that increase the susceptibility of the rock to weathering. Identifying deterioration factors is a first step in developing appropriate preventive conservation measures, which should aim to reduce rock contact time with water (technically affordable in the short term) and to reduce the probability of wildfire occurrence (technically more complex and possibly with longer‐term results). Copyright © 2018 John Wiley & Sons, Ltd.     

Assessment of soil factors controlling ephemeral gully erosion on agricultural fields

The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion. Copyright © 2018 John Wiley & Sons, Ltd.     

Mesocosm experiments identifying hotspots of groundwater upwelling in a water column by fibre optic distributed temperature sensing

Lacustrine groundwater discharge (LGD) can substantially impact ecosystem characteristics and functions. Fibre optic distributed temperature sensing (FO‐DTS) has been successfully used to locate groundwater discharge into lakes and rivers at the sediment–water interface, but locating groundwater discharge would be easier if it could be detected from the more accessible water surface. So far, it is not clear if how and under which conditions the LGD signal propagates through the water column to the water surface–atmosphere interface, and what perturbations and signal losses occur along this pathway. In the present study, LGD was simulated in a mesocosm experiment. Under winter conditions, water with temperatures of 14 to 16 °C was discharged at the bottom of a 10 × 2.8‐m mesocosm. Water within this mesocosm ranged from 4.0 to 7.4 °C. Four layers (20, 40, 60, and 80 cm above the sediment) of the 82 cm deep mesocosm were equipped with FO‐DTS for tracing thermal patterns in the mesocosm. Aims are (a) to test whether the positive buoyancy of relatively warm groundwater imported by LGD into shallow water bodies allows detection of LGD at the lake's water surface–atmosphere interface by FO‐DTS, (b) to analyse the propagation of the temperature signal from the sediment‐water interface through the water column, and (c) to learn more about detectability of the signal under different discharge rates and weather conditions. The experiments supported the benchmarking of scale dependencies and robustness of FO‐DTS applications for measuring upwelling into aquatic environments and revealed that weather conditions can have important impacts on the detection of upwelling at water surface–atmosphere interfaces at larger scales.     

Operation of two major reservoirs of Iran under IPCC scenarios during the XXI century

We assess the effects of prospective climate change until 2100 on water management of two major reservoirs of Iran, namely, Dez (3.34 × 10⁹ m³) and Alavian (6 × 10⁷ m³). We tune the Poly‐Hydro model suited for simulation of hydrological cycle in high altitude snow‐fed catchments. We assess optimal operation rules (ORs) for the reservoirs using three algorithms under dynamic and static operation and linear and non‐linear decision rules during control run (1990–2010 for Dez and 2000–2010 for Alavian). We use projected climate scenarios (plus statistical downscaling) from three general circulation models, EC‐Earth, CCSM4, and ECHAM6, and three emission scenarios, or representative concentration pathways (RCPs), RCP2.6, RCP4.5, and RCP8.5, for a grand total of nine scenarios, to mimic evolution of the hydrological cycle under future climate until 2100. We subsequently test the ORs under the future hydrological scenarios (at half century and end of century) and the need for reoptimization. Poly‐Hydro model when benchmarked against historical data well mimics the hydrological budget of both catchments, including the main processes of evapotranspiration and streamflows. Teaching–learning‐based optimization delivers the best performance in both reservoirs according to objective scores and is used for future operation. Our projections in Dez catchment depict decreased precipitation along the XXI century, with ?1% on average (of the nine scenarios) at half century and ?6% at the end of century, with changes in streamflows on average ?7% yearly and ?13% yearly, respectively. In Alavian, precipitation would decrease by ?10% on average at half century and ?13% at the end of century, with streamflows ?14% yearly and ?18% yearly, respectively. Under the projected future hydrology, reservoirs' operation would provide lower performance (i.e., larger lack of water) than now, especially for Alavian dam. Our results provide evidence of potentially decreasing water availability and less effective water management in water stressed areas like Northern Iran here during this century.     

A multicatchment analysis of headwater and downstream temperature effects from contemporary forest harvesting

Stream temperature is a key physical water‐quality parameter, controlling many biological, chemical, and physical processes in aquatic ecosystems. Maintenance of cool stream temperatures during summer is critical for high‐quality aquatic habitat. As such, transmission of warm water from small, nonfish‐bearing headwater streams after forest harvesting could cause warming in downstream fish‐bearing stream reaches with negative consequences. In this study, we evaluate (a) the effects of contemporary forest management practices on stream temperature in small, headwater streams, (b) the transmission of thermal signals from headwater reaches after harvesting to downstream fish‐bearing reaches, and (c) the relative role of lithology and forest management practices in influencing differential thermal responses in both the headwater and downstream reaches. We measured summer stream temperatures both preharvest and postharvest at 29 sites—12 upstream sites (4 reference, 8 harvested) and 17 downstream sites (5 reference, 12 harvested)—across 3 paired watershed studies in western Oregon. The 7‐day moving average of daily maximum stream temperature (T_7DAYMAX) was greater during the postharvest period relative to the preharvest period at 7 of the 8 harvested upstream sites. Although the T_7DAYMAX was generally warmer in the downstream direction at most of the stream reaches during both the preharvest and postharvest period, there was no evidence for additional downstream warming related to the harvesting activity. Rather, the T_7DAYMAX cooled rapidly as stream water flowed into forested reaches ~370–1,420 m downstream of harvested areas. Finally, the magnitude of effects of contemporary forest management practices on stream temperature increased with the proportion of catchment underlain by more resistant lithology at both the headwater and downstream sites, reducing the potential for the cooling influence of groundwater.     

The effects of initial soil moisture conditions on swale flow hydrographs

The effects of soil water content (SWC) on the formation of run‐off in grass swales draining into a storm sewer system were studied in two 30‐m test swales with trapezoidal cross sections. Swale 1 was built in a loamy fine‐sand soil, on a slope of 1.5%, and Swale 2 was built in a sandy loam soil, on a slope of 0.7%. In experimental runs, the swales were irrigated with 2 flow rates reproducing run‐off from block rainfalls with intensities approximately corresponding to 2‐month and 3‐year events. Run‐off experiments were conducted for initial SWC (SWC_ini) ranging from 0.18 to 0.43 m³/m³. For low SWC_ini, the run‐off volume was greatly reduced by up to 82%, but at high SWC_ini, the volume reduction was as low as 15%. The relative swale flow volume reductions decreased with increasing SWC_ini and, for the conditions studied, indicated a transition of the dominating swale functions from run‐off dissipation to conveyance. Run‐off flow peaks were reduced proportionally to the flow volume reductions, in the range from 4% to 55%. The swale outflow hydrograph lag times varied from 5 to 15 min, with the high values corresponding to low SWC_ini. Analysis of swale inflow/outflow hydrographs for high SWC_ini allowed estimations of the saturated hydraulic conductivities as 3.27 and 4.84 cm/hr in Swales 1 and 2, respectively. Such estimates differed from averages (N = 9) of double‐ring infiltrometer measurements (9.41 and 1.78 cm/hr). Irregularities in swale bottom slopes created bottom surface depression storage of 0.35 and 0.61 m³ for Swales 1 and 2, respectively, and functioned similarly as check berms contributing to run‐off attenuation. The experimental findings offer implications for drainage swale planning and design: (a) SWC_ini strongly affect swale functioning in run‐off dissipation and conveyance during the early phase of run‐off, which is particularly important for design storms and their antecedent moisture conditions, and (b) concerning the longevity of swale operation, Swale 1 remains fully functional even after almost 60 years of operation, as judged from its attractive appearance, good infiltration rates (3.27 cm/hr), and high flow capacity.     

MoisturEC: A New R Program for Moisture Content Estimation from Electrical Conductivity Data

Noninvasive geophysical estimation of soil moisture has potential to improve understanding of flow in the unsaturated zone for problems involving agricultural management, aquifer recharge, and optimization of landfill design and operations. In principle, several geophysical techniques (e.g., electrical resistivity, electromagnetic induction, and nuclear magnetic resonance) offer insight into soil moisture, but data‐analysis tools are needed to “translate” geophysical results into estimates of soil moisture, consistent with (1) the uncertainty of this translation and (2) direct measurements of moisture. Although geostatistical frameworks exist for this purpose, straightforward and user‐friendly tools are required to fully capitalize on the potential of geophysical information for soil‐moisture estimation. Here, we present MoisturEC, a simple R program with a graphical user interface to convert measurements or images of electrical conductivity (EC) to soil moisture. Input includes EC values, point moisture estimates, and definition of either Archie parameters (based on experimental or literature values) or empirical data of moisture vs. EC. The program produces two‐ and three‐dimensional images of moisture based on available EC and direct measurements of moisture, interpolating between measurement locations using a Tikhonov regularization approach.