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.     

Model for multiblock columns subjected to base excitation

A new model is presented for multiblock columns subjected to earthquakes, which contains an impact and an opening model. Both in the impact and in the opening model, all the possible opening configurations are investigated because it was found that in many practical cases, unexpected patterns may occur. The model is purely mechanical: assuming rigid blocks and classical (inelastic) impact. The effect of energy dissipation during impact was investigated. Using our model in accordance with the literature, it was found that monolithic blocks are more vulnerable to overturning than multiblock systems.     

Depth of closure: New calculation method based on sediment data

Obtaining depth of closure (DoC) in an accurate manner is a fundamental issue for coastal engineering, since good results for coastal structures and beach nourishment depend mainly on DoC. Currently, there are two methods for obtaining the DoC, mathematical formulations and profile surveys. However, these methods can incur important errors if one does not take into account the characteristics and morphology of the area, or if one does not have a sufficiently long time series. In this work the DoC is obtained from the break in the trend of the sediment with the depth, that is, in general with the increase of the depth a decrease in the size of the sediment takes place. However, at one point this tendency changes and the size increases, and then decreases again. When comparing the point where the minimum sediment size occurs before the increase, it is observed that the error incurred is small compared to other methods. If the Standard Deviation of Depth Change (SDDC) method is considered as the most accurate method, the error incurred by the proposed method is less than 7%. In addition, it can be seen that the dispersion of the sediment method always occurs outside the zone of bar movement. Whereas in the methods of profiles survey (using 2 cm precision profiles), sometimes the DoC is obtained within the active zone of bar movement. In addition, where the relative minimum of the median sediment size is found, and the sizes of 0.063 and 0.125 mm predominate in the composition of the sample. Therefore, this new method allows the precise location of the DoC to be obtained in a fast and simple way. Furthermore, this method has the advantage that it is not affected by the modifications that may be experienced by both the study area and the cross-shore beach profile.     

Simulating rainfall time-series: how to account for statistical variability at multiple scales?

Daily rainfall is a complex signal exhibiting alternation of dry and wet states, seasonal fluctuations and an irregular behavior at multiple scales that cannot be preserved by stationary stochastic simulation models. In this paper, we try to investigate some of the strategies devoted to preserve these features by comparing two recent algorithms for stochastic rainfall simulation: the first one is the modified Markov model, belonging to the family of Markov-chain based techniques, which introduces non-stationarity in the chain parameters to preserve the long-term behavior of rainfall. The second technique is direct sampling, based on multiple-point statistics, which aims at simulating a complex statistical structure by reproducing the same data patterns found in a training data set. The two techniques are compared by first simulating a synthetic daily rainfall time-series showing a highly irregular alternation of two regimes and then a real rainfall data set. This comparison allows analyzing the efficiency of different elements characterizing the two techniques, such as the application of a variable time dependence, the adaptive kernel smoothing or the use of low-frequency rainfall covariates. The results suggest, under different data availability scenarios, which of these elements are more appropriate to represent the rainfall amount probability distribution at different scales, the annual seasonality, the dry-wet temporal pattern, and the persistence of the rainfall events.     

Hybrid denoising-jittering data pre-processing approach to enhance multi-step-ahead rainfall–runoff modeling

Successful modeling of stochastic hydro-environmental processes widely relies on quantity and quality of accessible data and noisy data might effect on the functioning of the modeling. On the other hand in training phase of any Artificial Intelligence based model, each training data set is usually a limited sample of possible patterns of the process and hence, might not show the behavior of whole population. Accordingly in the present article first, wavelet-based denoising method was used in order to smooth hydrological time series and then small normally distributed noises with the mean of zero and various standard deviations were generated and added to the smoothed time series to form different denoised-jittered training data sets, for Artificial Neural Network (ANN) and Adaptive Neuro-Fuzzy Inference System (ANFIS) modeling of daily and multi-step-ahead rainfall–runoff process of the Milledgeville station of the Oconee River and the Pole Saheb station of the Jighatu River watersheds, respectively located in USA and Iran. The proposed hybrid data pre-processing approach in the present study is used for the first time in modeling of time series and especially in modeling of hydrological processes. Furthermore, the impacts of denoising (smoothing) and noise injection (jittering) have been simultaneously investigated neither in hydrology nor in any other engineering fields. To evaluate the modeling performance, the outcomes were compared with the results of multi linear regression and Auto Regressive Integrated Moving Average models. Comparing the achieved results via the trained ANN and ANFIS models using denoised-jittered data showed that the proposed data pre-processing approach which serves both denoising and jittering techniques could improve performance of the ANN and ANFIS based single-step-ahead rainfall–runoff modeling of the Milledgeville station up to 14 and 12% and of the Pole Saheb station up to 22 and 16% in the verification phase. Also the results of multi-step-ahead modeling using the proposed data pre-processing approach showed improvement of modeling for both watersheds.     

Estimation of carbon stock for greenhouse gas emissions from hydropower reservoirs

The identification and accurate quantification of sources or sinks of greenhouse gas (GHG) have become a key challenge for scientists and policymakers working on climate change. The creation of a hydropower reservoir, while damming a river for power generation, converts the terrestrial ecosystems into aquatic and subsequently aerobic and anaerobic decomposition of flooded terrestrial soil organic matter resulting in the emission of significant quantity of GHG to the atmosphere. Tropical/subtropical hydropower reservoirs are more significant sources of GHG compared to boreal or temperate one. This paper aims to estimate the emission factor (gCO₂eq./kWh) and net GHG emission from Koteshwar hydropower reservoir in Uttarakhand, India. Further, estimated GHG are compared with those from global reservoirs located in the same eco-region so that its impact could be timely minimized/mitigated. Results have shown that emission factor and net GHG emission of Koteshwar reservoir are, respectively, estimated as 13.87 gCO₂eq./kWh and 167.70 Gg C year^?1 which are less than other global reservoirs located in the same eco-region. This information could be helpful for the hydropower industries to construct reservoirs in tropical eco-regions.     

Declustering experimental variograms by global estimation with fourth order moments

The variogram is a key parameter for geostatistical estimation and simulation. Preferential sampling may bias the spatial structure and often leads to noisy and unreliable variograms. A novel technique is proposed to weight variogram pairs in order to compensate for preferential or clustered sampling . Weighting the variogram pairs by global kriging of the quadratic differences between the tail and head values gives each pair the appropriate weight, removes noise and minimizes artifacts in the experimental variogram. Moreover, variogram uncertainty could be computed by this technique. The required covariance between the pairs going into variogram calculation, is a fourth order covariance that must be calculated by second order moments. This introduces some circularity in the calculation whereby an initial variogram must be assumed before calculating how the pairs should be weighted for the experimental variogram. The methodology is assessed by synthetic and realistic examples. For synthetic example, a comparison between the traditional and declustered variograms shows that the declustered variograms are better estimates of the true underlying variograms. The realistic example also shows that the declustered sample variogram is closer to the true variogram.     

Distributions-oriented wind forecast verification by a hidden Markov model for multivariate circular–linear data

Winds from the North–West quadrant and lack of precipitation are known to lead to an increase of PM10 concentrations over a residential neighborhood in the city of Taranto (Italy). In 2012 the local government prescribed a reduction of industrial emissions by 10% every time such meteorological conditions are forecasted 72 h in advance. Wind forecasting is addressed using the Weather Research and Forecasting (WRF) atmospheric simulation system by the Regional Environmental Protection Agency. In the context of distributions-oriented forecast verification, we propose a comprehensive model-based inferential approach to investigate the ability of the WRF system to forecast the local wind speed and direction allowing different performances for unknown weather regimes. Ground-observed and WRF-forecasted wind speed and direction at a relevant location are jointly modeled as a 4-dimensional time series with an unknown finite number of states characterized by homogeneous distributional behavior. The proposed model relies on a mixture of joint projected and skew normal distributions with time-dependent states, where the temporal evolution of the state membership follows a first order Markov process. Parameter estimates, including the number of states, are obtained by a Bayesian MCMC-based method. Results provide useful insights on the performance of WRF forecasts in relation to different combinations of wind speed and direction.