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.     

SlugIn 1.0: A Free Tool for Automated Slug Test Analysis

The correct characterization of aquifer parameters is essential for water‐supply and water‐quality investigations. Slug tests are widely used for these purposes. While free software is available to interpret slug tests, some codes are not user‐friendly, or do not include a wide range of methods to interpret the results, or do not include automatic, inverse solutions to the test data. The private sector has also generated several good programs to interpret slug test data, but they are not free of charge. The computer program SlugIn 1.0 is available online for free download, and is demonstrated to aid in the analysis of slug tests to estimate hydraulic parameters. The program provides an easy‐to‐use Graphical User Interface. SlugIn 1.0 incorporates automated parameter estimation and facilitates the visualization of several interpretations of the same test. It incorporates solutions for confined and unconfined aquifers, partially penetrating wells, skin effects, shape factor, anisotropy, high hydraulic conductivity formations and the Mace test for large‐diameter wells. It is available in English and Spanish and can be downloaded from the web site of the Geological Survey of Spain. Two field examples are presented to illustrate how the software operates.     

Bivariate Gaussian mixture–based equivalent linearization method for stochastic seismic analysis of nonlinear structures

To address challenges in stochastic seismic analysis of nonlinear structures, this paper further develops a recently proposed Gaussian mixture–based equivalent linearization method (GM‐ELM). The GM‐ELM uses a Gaussian mixture distribution model to approximate the probabilistic distribution of a nonlinear system response. Using properties of the Gaussian mixture model, GM‐ELM can decompose the non‐Gaussian response of a nonlinear system into multiple Gaussian responses of linear single–degree of freedom oscillators. With the set of the equivalent linear systems identified by GM‐ELM, response statistics as crossing rate and first‐passage probability can be computed conveniently using theories of linear random vibration analysis. However, the original version of GM‐ELM may lead to an inaccurate estimate because of the heuristic parameters of the linear system introduced to supplement insufficient information. To overcome this limitation and define unique equivalent linear systems, this paper proposes a further developed version of GM‐ELM, which uses a mixture of bivariate Gaussian densities instead of univariate models. Moreover, to facilitate the use of elastic response spectra for estimating the mean peak responses of a nonlinear structure, a new response spectrum combination rule is proposed for GM‐ELM. Two numerical examples of hysteretic structural systems are presented in this paper to illustrate the application of the bivariate GM‐ELM to nonlinear stochastic seismic analysis. The analysis results obtained by the bivariate GM‐ELM are compared with those obtained by the univariate GM‐ELM, the conventional equivalent linearization method, the tail equivalent linearization method, and Monte Carlo simulation. The supporting source code and data are available for download at https://github.com/yisangri/GitHub‐bGM‐ELM‐code.git     

Coordinative similitude law considering fluid‐structure interaction for underwater shaking table tests

Generally, when a model is made of the same material as the prototype in shaking table tests, the equivalent material density of the scaled model is greater than that of the prototype because mass is added to the model to satisfy similitude criteria. When the water environment is modeled in underwater shaking table tests, however, it is difficult to change the density of water. The differences in the density similitude ratios of specimen materials and water can affect the similitude ratios of the hydrodynamic and wave forces with those of other forces. To solve this problem, a coordinative similitude law is proposed for underwater shaking table tests by adjusting the width of the upstream face of the model or the wave height in the model test to match the similitude ratios of hydrodynamic and wave forces with those of other forces. The designs of the similitude relations were investigated for earthquake excitation, wave excitation, and combined earthquake and wave excitation conditions. Series of numerical simulations and underwater shaking table tests were performed to validate the proposed coordinative similitude law through a comparison of coordinative model and conventional model designed based on the coordinative similitude law and traditional artificial mass simulation, respectively. The results show that the relative error was less than 10% for the coordinative model, whereas it reached 80% for the conventional model. The coordinative similitude law can better reproduce the dynamic responses of the prototype, and thus, this similitude law can be used in underwater shaking table tests.     

Stream channel erosion in a rapidly urbanizing region of the US–Mexico border: documenting the importance of channel hardpoints with Structure‐from‐Motion photogrammetry

Urbanization can lead to accelerated stream channel erosion, especially in areas experiencing rapid population growth, unregulated urban development on erodible soils, and variable enforcement of environmental regulations. A combination of field surveys and Structure‐from‐Motion (SfM) photogrammetry techniques was used to document spatial patterns in stream channel geometry in a rapidly urbanizing watershed, Los Laureles Canyon (LLCW), in Tijuana, Mexico. Ground‐based SfM photogrammetry was used to map channel dimensions with 1 to 2 cm vertical mean error for four stream reaches (100–300 m long) that were highly variable and difficult to survey with a differential GPS. Regional channel geometry curves for LLCW had statistically larger slopes and intercepts compared with regional curves developed for comparable, undisturbed reference channels. Cross‐sectional areas of channels downstream of hardpoints, such as concrete reaches or culverts, were up to 64 times greater than reference channels, with enlargement persisting, in some cases, up to 230 m downstream. Percentage impervious cover was not a good predictor of channel enlargement. Proximity to upstream hardpoint, and lack of riparian and bank vegetation paired with highly erodible bed and bank materials may account for the instability of the highly enlarged and unstable cross‐sections. Channel erosion due to urbanization accounts for approximately 25–40% of the total sediment budget for the watershed, and channel erosion downstream of hardpoints accounts for one‐third of all channel erosion. Channels downstream of hardpoints should be stabilized to prevent increased inputs of sediment to the Tijuana Estuary and local hazards near the structures, especially in areas with urban settlements near the stream channel. Copyright © 2017 John Wiley & Sons, Ltd.