Short‐ and long‐term evapotranspiration rates at ecological restoration sites along a large river receiving rare flow events

Many large rivers around the world no longer flow to their deltas, due to ever greater water withdrawals and diversions for human needs. However, the importance of riparian ecosystems is drawing increasing recognition, leading to the allocation of environmental flows to restore river processes. Accurate estimates of riparian plant evapotranspiration (ET) are needed to understand how the riverine system responds to these rare events and achieve the goals of environmental flows. In 2014, historic environmental flows were released into the Lower Colorado River at Morelos Dam (Mexico); this once perennial but now dry reach is the final stretch to the mighty Colorado River Delta. One of the primary goals was to supply native vegetation restoration sites along the reach with water to help seedlings establish and boost groundwater levels to foster the planted saplings. Patterns in ET before, during, and after the flows are useful for evaluating whether this goal was met and understanding the role that ET plays in this now ephemeral river system. Here, diurnal fluctuations in groundwater levels and Moderate Resolution Imaging Spectroradiometer (MODIS) data were used to compare estimates of ET specifically at 3 native vegetation restoration sites during 2014 planned flow events, and MODIS data were used to evaluate long‐term (2002–2016) ET responses to restoration efforts at these sites. Overall, ET was generally 0–10 mm d^?1 across sites, and although daily ET values from groundwater data were highly variable, weekly averaged estimates were highly correlated with MODIS‐derived estimates at most sites. The influence of the 2014 flow events was not immediately apparent in the results, although the process of clearing vegetation and planting native vegetation at the restoration sites was clearly visible in the results.     

Effects of fracture intersections on seismic dispersion: theoretical predictions versus numerical simulations

The detection and characterisation of domains of intersecting fractures are important goals in several disciplines of current interest, including exploration and production of unconventional reservoirs, nuclear waste storage, CO₂ sequestration, and groundwater hydrology, among others. The objective of this study is to propose a theoretical framework for quantifying the effects of fracture intersections on the frequency‐dependent elastic properties of fluid‐saturated porous and fractured rocks. Three characteristic frequency regimes for fluid pressure communication are identified. In the low‐frequency limit, fractures are in full pressure communication with the embedding porous matrix and with other fractures. Conversely, in the high‐frequency limit, fractures are hydraulically isolated from the matrix and from other fractures. At intermediate frequencies, fractures are hydraulically isolated from the matrix porosity but can be in hydraulic communication with each other, depending on whether fracture sets are intersecting. For each frequency regime, the effective stiffness coefficients are derived using the linear‐slip theory and anisotropic Gassmann equations. Explicit mathematical expressions for the two characteristic frequencies that separate the three frequency regimes are also determined. Theoretical predictions are then applied to two synthetic 2D samples, each containing two orthogonal fracture sets: one with and another without intersections. The resulting stiffness coefficients, Thomsen‐style anisotropy parameters, and the transition frequencies show good agreement with corresponding numerical simulations. The theoretical results are applicable not only to 2D but also to 3D fracture systems and are amenable to being employed in inversion schemes designed to characterise fracture systems.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

Water quality and presence of pesticides in a tropical coastal wetland in southern Mexico

The water quality of the Pozuelos-Murillo lagoon system in southern Mexico was evaluated during three periods between March and October 2002, with particular emphasis on the detection of organochlorine and organophosphorous pesticide residues in water and sediments. Physicochemical and microbiological parameters of water were also measured and integrated in a water quality index (WQI). Multivariate analysis was used to regionalise the lagoon system as a function of the behaviour of all measured parameters. Solid phase micro-extraction followed by gas chromatography (SPME-GC) was used for pesticide analysis. The concentration of phosphorous was found to be higher than that of nitrogenous compounds. This, besides a deficiency in dissolved oxygen and a high organic matter concentration (as COD), reflects eutrophication processes in some areas of the system. Measured levels of faecal coliforms and oils and greases were above the limits established by Mexican law and comparable to the concentrations reported for other highly polluted systems in Mexico. Residues of DDD (2.0 microg L(-1)) in water and DDE (247 ng g(-1)) and endosulfan I (814 ng g(-1)) sediments were detected by SPME-GC. The spatial distribution of these contaminants implies major potential risks because the most polluted sites were found to be those with the highest fishing activity. Although in general the WQI is on acceptable levels (65-80%), some contamination problems are evident.     

Extracting water-related features using reflectance data and principal component analysis of Landsat images

This study aimed to map water features using a Landsat image rather than traditional land cover. We involved the original bands, spectral indices and principal components (PCs) of a principal component analysis (PCA) as input data, and performed random forest (RF) and support vector machine (SVM) classification with water, saturated soil and non-water categories. The aim was to compare the efficiency of the results based on various input data. Original bands provided 93% overall accuracy (OA) and bands 4–5–7 were the most informative in this analysis. Except for MNDWI (modified normalized differenced water index, with 98% OA), the performance of all water indices was between 60 and 70% (OA). The PCA-based approach conducted on the original bands resulted in the most accurate identification of all classes (with only 1% error in the case of water bodies). We therefore show that both water bodies and saturated soils can be identified successfully using this approach.     

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.     

Identifying slope processes over time and their imprint in soils of medium‐high mountains of Central Europe (the Karkonosze Mountains,Poland)

Soils in mountainous areas are often polygenetic, developed in slope covers that relate to glacial and periglacial activities of the Pleistocene and Holocene and reflect climatic variations. Landscape development during the Holocene may have been influenced by erosion/solifluction that often started after the Holocene climatic optimum. To trace back soil evolution and its timing, we applied a multi‐methodological approach. This approach helped us to outline scenario of soil transformation. According to our results, some aeolian input must have occurred in the late Pleistocene. During that time and the early Holocene, the soils most likely had features of Cryosols or Leptosols. Physico‐chemical and mineralogical analyses have indicated that the material was denudated (between late Boreal to the Atlantic) from the ridge and upper‐slope positions forming a colluvium at mid‐slope positions. Later, during the Sub‐Boreal, mass wasting of the remains of silt material deposited at the end of the Pleistocene age on the ridge top seems to have occurred. In addition, the cool and moist conditions caused the deposition of a colluvium at the lower‐slope positions. The next phase was characterized by the transformation of Leptosols/Cambisols into Podzols at upper‐slope or shoulder positions and to Albic Cambisols at mid‐slope positions. During the Sub‐Boreal period, Stagnosols started to form at the lower part of the slope catena. Overall, the highest erosion rates were calculated at the upper‐slope position and the lowest rates at mid‐slope sites. Berylium‐10 (¹⁰Be) data showed that the Bs, BC/C were covered during the Holocene by a colluvium with a different geological composition which complicated the calculation of erosion or accumulation rates. The interpretation of erosion and accumulation rates in such multi‐layered materials may, therefore, be hampered. However, the multi‐methodological reconstruction we applied shed light on the soil and landscape evolution of the eastern Karkonosze Mountains. Copyright © 2017 John Wiley & Sons, Ltd.     

Transport of moderately sorted gravel at low bed shear stresses: The role of fine sediment infiltration

A reliable estimation of sediment transport in gravel‐bed streams is important for various practical engineering and biological studies (e.g., channel stability design, bed degradation/aggradation, restoration of spawning habitat). In the present work, we report original laboratory experiments investigating the transport of gravel particles at low bed shear stresses. The laboratory tests were conducted under unsteady flow conditions inducing low bed shear stresses, with detailed monitoring of the bed topography using a laser scanner. Effects of bed surface arrangements were documented by testing loose and packed bed configurations. Effects of fine sediments were examined by testing beds with sand, artificial fine sand or cohesive silt infiltrated in the gravel matrix. Analysis of the experimental data revealed that the transport of gravel particles depends upon the bed arrangement, the bed material properties (e.g., size and shape, consolidation index, permeability) and the concentration of fine sediments within the surface layer of moving grains. This concentration is directly related to the distribution of fine particles within the gravel matrix (i.e., bottom‐up infiltration or bridging) and their transport mode (i.e., bedload or suspended load). Compared to loose beds, the mobility of gravel is reduced for packed beds and for beds clogged from the bottom up with cohesive fine sediments; in both cases, the bed shear stress for gravel entrainment increases by about 12%. On the other hand, the mobility of gravel increases significantly (bed shear stress for particle motion decreasing up to 40%) for beds clogged at the surface by non‐cohesive sand particles. Copyright © 2017 John Wiley & Sons, Ltd.     

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.     

Estimating external magnetic field differences at high geomagnetic latitudes from a single station

Providing an accurate estimate of the magnetic field on the Earth's surface at a location distant from an observatory has useful scientific and commercial applications, such as in repeat station data reduction, space weather nowcasting or aeromagnetic surveying. While the correlation of measurements between nearby magnetic observatories at low and mid‐latitudes is good, at high geomagnetic latitudes () the external field differences between observatories increase rapidly with distance, even during relatively low magnetic activity. Thus, it is of interest to describe how the differences (or errors) in external magnetic field extrapolation from a single observatory grow with distance from its location. These differences are modulated by local time, seasonal and solar cycle variations, as well as geomagnetic activity, giving a complex temporal and spatial relationship. A straightforward way to describe the differences are via confidence intervals for the extrapolated values with respect to distance. To compute the confidence intervals associated with extrapolation of the external field at varying distances from an observatory, we used 695 station‐years of overlapping minute‐mean data from 37 observatories and variometers at high latitudes from which we removed the main and crustal fields to isolate unmodelled signals. From this data set, the pairwise differences were analysed to quantify the variation during a range of time epochs and separation distances. We estimate the 68.3%, 95.4% and 99.7% confidence levels (equivalent to the 1σ, 2σ and 3σ Gaussian error bounds) from these differences for all components. We find that there is always a small non‐zero bias that we ascribe to instrumentation and local crustal field induction effects. The computed confidence intervals are typically twice as large in the north–south direction compared to the east‐west direction and smaller during the solstice months compared to the equinoxes.