Identification of potential subsidence related to pumping in the Almería basin (SE Spain)

Ground subsidence of detrital deposits in the Almería basin (SE Spain) was studied using the remote sensing technique of Differential Interferometry SAR (DInSAR). This basin is one of the most arid in Europe, receiving an average rainfall of 250 mm per year. Over the last 60 years the region has experienced an enormous agricultural and urban expansion, whose water demand has been largely supplied from groundwater, leading to the current situation of overexploitation of water resources. This paper outlines the likely relationship between groundwater abstraction and subsidence. To this end, 34 ERS and Envisat images, taken between 2003 and 2009, were analysed to estimate ground surface deformations, and hence, compared with water table variations measured in a number of piezometers in the basin. The analysis shows a clear parallelism between the variations in piezometric level and deformation of the ground surface. In addition, the zones of greatest subsidence coincide with those areas where groundwater abstractions are concentrated. Subsidence over the examined period varies from 10 to 30 mm, with extreme values as high as 50 mm, which translates to a rate of between 1·7 and 5 mm/year, reaching maximum rates of 8 mm/year at some points. Given such subsidence rates, damage to urban infrastructures are, for the moment, incipient. Copyright © 2011 John Wiley & Sons, Ltd.     

Scientific management of Mediterranean coastal zone: a hybrid ocean forecasting system for oil spill and search and rescue operations

The oil spill from Prestige tanker showed the importance of scientifically based protocols to minimize the impacts on the environment. In this work, we describe a new forecasting system to predict oil spill trajectories and their potential impacts on the coastal zone. The system is formed of three main interconnected modules that address different capabilities: (1) an operational circulation sub-system that includes nested models at different scales, data collection with near-real time assimilation, new tools for initialization or assimilation based on genetic algorithms and feature-oriented strategic sampling; (2) an oil spill coastal sub-system that allows simulation of the trajectories and fate of spilled oil together with evaluation of coastal zone vulnerability using environmental sensitivity indexes; (3) a risk management sub-system for decision support based on GIS technology. The system is applied to the Mediterranean Sea where surface currents are highly variable in space and time, and interactions between local, sub-basin and basin scale increase the non-linear interactions effects which need to be adequately resolved at each one of the intervening scales. Besides the Mediterranean Sea is a complex reduced scale ocean representing a real scientific and technological challenge for operational oceanography and particularly for oil spill response and search and rescue operations.     

Acoustic emission processes occurring during high-pressure sand compaction

Granular materials submitted to uniaxial compression undergo pore space reduction due to mechanisms such as particle rearrangement and grain crushing. These changes in the internal structure of the material release energy in the form of elastic waves that can be detected by sensors sensitive to acoustic emission. In this study, Acoustic emission monitoring with a wavelet-based signal processing technique is used to identify the various mechanisms occurring during high-pressure sand compaction. Particle movement, grain failure, friction between grains and the surface of the compression cell and intergranular friction are studied. Acoustic emission data recorded during these simplified tests are used to characterize each phenomenon. Wavelet transform analyses allow the identification of useful features, making possible frequency discrimination among sliding, rolling, friction and grain fragmentation processes. For instance, we observe that at low stress, grain flow is characterized by the lowest centroid and peak frequencies, while at greater stresses, intergranular friction and grain fragmentation have the higher values. In the tests performed, the stress–strain evolution and final condition of the tested sand are broadly consistent, irrespective of the condition employed: continuous, stepwise or even variable loading rate or temperature. However, Acoustic emission data manifest much more complex behaviour (including thermal, load-rate dependency and delayed fragmentation phenomena) than that suggested by stress–strain relationships. At low stress levels, grain flow (sliding/rolling) is a relevant strain-accommodation mechanism, but so is crushing due to the effect of concentrated forces at the grain contact level. At high stresses, when crushing is generalized, intergranular friction is also a relevant phenomenon due to the increase in the coordination number produced by previous fragmentation.     

Comparison of stress-dependent geophysical,hydraulic and mechanical properties of synthetic and natural sandstones for reservoir characterization and monitoring studies

Synthetic rock samples can offer advantages over natural rock samples when used for laboratory rock physical properties studies, provided their success as natural analogues is well understood. The ability of synthetic rocks to mimic the natural stress dependency of elastic wave, electrical and fluid transport properties is of primary interest. Hence, we compare a consistent set of laboratory multi-physics measurements obtained on four quartz sandstone samples (porosity range 20–25%) comprising two synthetic and two natural (Berea and Corvio) samples, the latter used extensively as standards in rock physics research. We measured simultaneously ultrasonic (P- and S-wave) velocity and attenuation, electrical resistivity, permeability and axial and radial strains over a wide range of differential pressure (confining stress 15–50 MPa; pore pressure 5–10 MPa) on the four brine saturated samples. Despite some obvious physical discrepancies caused by the synthetic manufacturing process, such as silica cementation and anisotropy, the results show only small differences in stress dependency between the synthetic and natural sandstones for all measured parameters. Stress dependency analysis of the dry samples using an isotropic effective medium model of spheroidal pores and penny-shaped cracks, together with a granular cohesion model, provide evidence of crack closure mechanisms in the natural sandstones, seen to a much lesser extent in the synthetic sandstones. The smaller grain size, greater cement content, and cementation under oedometric conditions particularly affect the fluid transport properties of the synthetic sandstones, resulting in lower permeability and higher electrical resistivity for a similar porosity. The effective stress coefficients, determined for each parameter, are in agreement with data reported in the literature. Our results for the particular synthetic materials that were tested suggest that synthetic sandstones can serve as good proxies for natural sandstones for studies of elastic and mechanical properties, but should be used with care for transport properties studies.     

Weather types,runoff and sediment yield in a Mediterranean mountain landscape

The relationship between the synoptic weather types (WTs), runoff and sediments in a Mediterranean mountain landscape was analysed. The study was performed between 2005 and 2012 using one of the most complete and extensive daily databases of rainfall, runoff and sediment recorded in the Spanish Pyrenees, coupled with WTs defined from the National Meteorological Center/National Center for Atmospheric Research (NMC/NCAR) 40‐Year Reanalysis Project ? ? Correction added on 18 December 2013, after online publication: EMULATE project was replaced with NMC/NCAR 40‐Year Reanalysis Project.
. The results show that the three wettest WTs accounted for 30% of rainy days and 46% of rainfall, but comprised only 13% of total daily records. To obtain a much more robust association between WTs and rainfall an analysis was carried out using a longer rainfall record (1989–2011). The analyses confirmed that the results obtained from the Araguás catchment are representative of a longer time span. The cyclonic, north‐westerly and westerly WTs play an important role in runoff generation, coinciding with the wettest WTs. Extreme floods are commonly associated also with south‐westerly and westerly airflows, whereas less flooding was generated under cyclonic circulations. Sediment transport was concentrated in 2.9% of total time mainly related to westerly WTs. Seasonal differences exist in WT frequency. In winter and spring north‐west and West are the most prominent WTs related to rainfall, water and sediment yield, although in spring cyclonic frequency was higher. During autumn north‐west and south‐west were the most frequent, but sediment yield was produced nearby under south‐west flow. In summer the WTs that produce sediment are north and west. A magnitude–frequency analysis shows the different behaviour of WTs according to their efficiency in producing runoff and sediment. A study with different monitored areas around the Iberian Peninsula and the Mediterranean basin would be very valuable in providing information for hydrological and sediment behaviour under the current conditions of global climate change. Copyright © 2013 John Wiley & Sons, Ltd.     

On matrix diffusion: formulations, solution methods and qualitative effects

 Matrix diffusion has become widely recognized as an important transport mechanism. Unfortunately, accounting for matrix diffusion complicates solute-transport simulations. This problem has led to simplified formulations, partly motivated by the solution method. As a result, some confusion has been generated about how to properly pose the problem. One of the objectives of this work is to find some unity among existing formulations and solution methods. In doing so, some asymptotic properties of matrix diffusion are derived. Specifically, early-time behavior (short tests) depends only on φ _m ² R _m D _m  / L _m ², whereas late-time behavior (long tracer tests) depends only on φ _m R _m , and not on matrix diffusion coefficient or block size and shape. The latter is always true for mean arrival time. These properties help in: (a) analyzing the qualitative behavior of matrix diffusion; (b) explaining one paradox of solute transport through fractured rocks (the apparent dependence of porosity on travel time); (c) discriminating between matrix diffusion and other problems (such as kinetic sorption or heterogeneity); and (d) describing identifiability problems and ways to overcome them. Received, October 1997 · Revised, November 1997 · Accepted, December 1997     

Study of the interaction between U(VI) and the anoxic corrosion products of carbon steel

The objective of this study was to investigate the removal mechanism of U(VI) from groundwater by magnetite as the main product of anoxic steel corrosion. For this purpose, a systematic sequence of batch experiments was conducted to focus the active role of magnetite in the reduction of U under different conditions. Results indicated that under anoxic conditions U(VI) was sorbed at the magnetite surface, whereas under reducing conditions at different H₂(g) pressures, U was present in tetravalent form as amorphous UO₂.     

Vulnerability of simple reinforced concrete buildings to damage by rockfalls

A procedure is presented for investigating the response of reinforced concrete buildings to rockfall impact. The method considers a single rock hit on the basement columns, and it includes four steps: (a) calculation of the probability of a rock impact on a member of the load-bearing system, taking into account the block size and the design of the structure; (b) evaluation of the response of one or more structural elements to the hit based on element capacity; (c) in the case of structural element failure, assessment of the robustness of the whole structural system, calculating the potential for progressive collapse; and (d) calculation of a damage index (DI), which is the ratio of structural elements that fail to the total number of structural elements. The proposed method is applied to a reinforced concrete building for a range of rockfall paths and intensities. The analysis has been carried out for a 2-m-diameter block and velocities < 3.5 m/s. The possible damage range is found to be highly variable, with DI values ranging from 0.01 to 1 depending on the impact location and block velocity.