Approximating constant‐Q seismic propagation in the time domain

In this study, we investigate the accuracy of approximating constant‐Q wave propagation by series of Zener or standard linear solid (SLS) mechanisms. Modelling in viscoacoustic and viscoelastic media is implemented in the time domain using the finite‐difference (FD) method. The accuracy of numerical solutions is evaluated by comparison with the analytical solution in homogeneous media. We found that the FD solutions using three SLS relaxation mechanisms as well as a single SLS mechanism, with properly chosen relaxation times, are quite accurate for both weak and strong attenuation. Although the RMS errors of FD simulations using a single relaxation mechanism increase with increasing offset, especially for strong attenuation (Q = 20), the results are still acceptable for practical applications. The synthetic data of the Marmousi‐II model further illustrate that the single SLS mechanism, to model constant Q, is efficient and sufficiently accurate. Moreover, it benefits from less computational costs in computer time and memory.     

Angular and Frequency-Dependent Wave Velocity and Attenuation in Fractured Porous Media

Wave-induced fluid flow generates a dominant attenuation mechanism in porous media. It consists of energy loss due to P-wave conversion to Biot (diffusive) modes at mesoscopic-scale inhomogeneities. Fractured poroelastic media show significant attenuation and velocity dispersion due to this mechanism. The theory has first been developed for the symmetry axis of the equivalent transversely isotropic (TI) medium corresponding to a poroelastic medium containing planar fractures. In this work, we consider the theory for all propagation angles by obtaining the five complex and frequency-dependent stiffnesses of the equivalent TI medium as a function of frequency. We assume that the flow direction is perpendicular to the layering plane and is independent of the loading direction. As a consequence, the behaviour of the medium can be described by a single relaxation function. We first consider the limiting case of an open (highly permeable) fracture of negligible thickness. We then compute the associated wave velocities and quality factors as a function of the propagation direction (phase and ray angles) and frequency. The location of the relaxation peak depends on the distance between fractures (the mesoscopic distance), viscosity, permeability and fractures compliances. The flow induced by wave propagation affects the quasi-shear (qS) wave with levels of attenuation similar to those of the quasi-compressional (qP) wave. On the other hand, a general fracture can be modeled as a sequence of poroelastic layers, where one of the layers is very thin. Modeling fractures of different thickness filled with CO₂ embedded in a background medium saturated with a stiffer fluid also shows considerable attenuation and velocity dispersion. If the fracture and background frames are the same, the equivalent medium is isotropic, but strong wave anisotropy occurs in the case of a frameless and highly permeable fracture material, for instance a suspension of solid particles in the fluid.     

Temporal dynamics of soil water balance components in a karst range in southeastern Spain: estimation of potential recharge

Abstract

This paper analyses the temporal dynamics of soil water balance components in a representative recharge area of the Sierra de Gádor (Almeria, southeastern Spain) in two hydrological years. Two approaches are used to estimate daily potential recharge (PR): Approach 1 based on deriving PR from the water balance as the difference between measurements of rainfall (P) and actual evapotranspiration (E) obtained by eddy covariance; and Approach 2 with PR obtained from the dynamic pattern of the soil moisture (θ) recorded at two depths in the site's thin soil (average 0.35 m thickess). For the hydrological year 2003/04, which was slightly drier than the 30-year average, E accounted for 64% of rainfall and occurred mainly in late spring and early summer. The PR estimated by Approach 1 was 181 ± 18 mm year^-1 (36% of rainfall), suggesting an effective groundwater recharge in the study area. In the unusually dry hydrological year 2004/05, E was about 215 mm year^-1, close to the annual rainfall input, and allowing very little (8 ± 12 mm year^-1) PR according to Approach 1. Estimation of PR based on Approach 2 resulted in PR rates lower than those found by Approach 1, because Approach 2 does not take into account the recharge that occurs through preferential flow pathways (cracks, joints and fissures) which were not monitored with the θ probes. Moreover, using Approach 2, the PR estimates differed widely depending on the time scale considered: with daily mean θ data, PR estimation was lower, especially in late spring, while θ data at 30 min resolution yielded a more reliable prediction of the fraction of total PR resulting from the downward movement of soil water by gravity.

Citation Cantón, Y., Villagarcía, L., Moro, M. J., Serrano-Ortíz, P., Were, A., Alcalá, F. J., Kowalski, A. S., Solé-Benet, A., Lázaro, R. & Domingo, F. (2010) Temporal dynamics of soil water balance components in a karst range in southeastern Spain: estimation of potential recharge. Hydrol. Sci. J. 55(5), 737–753.     

Runoff and sediment yield from land under various uses in a Mediterranean mountain area: long‐term results from an experimental station

In this study we analyzed runoff and sediment yield from land under various traditional and current land uses in Mediterranean mountain areas, using long‐term data from an experimental station in the Aísa Valley, Central Spanish Pyrenees. Monitoring at this station has provided 20 years of data that can help explain the hydrological and geomorphological changes that have been observed at larger spatial scales, and also the changes that have occurred to some of the most characteristic landscapes of the Mediterranean middle mountains. In spite of the problems associated with the use of small experimental plots, the results obtained are consistent with other studies in the Mediterranean region, and confirm the strong influence of land use changes on runoff generation and sediment yield. The results indicate that: (i) cereal cultivation on steep slopes (both alternating cereal cultivation and fallow on sloping fields and shifting agriculture on the steepest slopes) represents a major problem for soil conservation. This explains the occurrence throughout the Mediterranean mountains of many degraded hillslopes, which show evidence of sheet wash erosion, rilling, gullying and shallow landsliding; (ii) farmland abandonment has led to a marked reduction in runoff and sediment yield as a consequence of rapid plant recolonization, particularly by dense shrubs; (iii) the natural transformation of abandoned fields into grazing meadows has reduced runoff and sediment yield. Land use trends in the Mediterranean mountains are mainly characterized by generalized farmland abandonment and a decrease in livestock pressure. From a hydrological and geomorphological point of view the main consequences have been a reduction in overland flow from the hillslopes, and a reduction in sediment sources, with differences up to one order of magnitude in sediment yield from dense shrub cover and grazing meadow areas compared with areas under shifting agriculture. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaporation from shallow groundwater in closed basins in the Chilean Altiplano

Abstract

The hydrological cycle in arid and semi-arid climates is highly controlled by evaporation. The correct quantification of this process is essential for improving the accuracy of water balance estimates, especially in closed basins. The objective of this paper is to characterize evaporation rates from shallow groundwater using the chamber approach in six closed basins in the Altiplano of northern Chile. Measurements were made at 49 locations with water-table depths ranging from 0.09 m to 3.3 m. Estimated daily evaporation rates appeared to be strongly related to groundwater depth and soil texture. In particular, the highest rates were recorded in areas with high groundwater tables and coarse-grained soils. Evaporation curves were derived by fitting exponential and power relationships as functions of the groundwater depths that we proposed to use in the study area. An application of these curves for the Salar de Pedernales basin produced an estimated evaporation flow of 530 L s^-1, using the average curve.

Citation Johnson, E., Yáñez, J., Ortiz, C. & Muñoz, J. (2010) Evaporation from shallow groundwater in closed basins in the Chilean Altiplano. Hydrol. Sci. J. 55(4), 624–635.     

Geophysical-Archaeological Survey in Lake Tequesquitengo,Morelos, Mexico

In August 2009, a marine geophysical survey was conducted in Lake Tequesquitengo (located in the state of Morelos, Mexico) to delineate the extent of the remains of a small town that has been submerged since the mid 19th century. The survey consists of the acquisition and mapping of magnetic, single beam bathymetric and side-scan sonar data. A dual receiver marine GPS navigation system was used to position the boat during the survey. Except for the larger structural remains that are visible on the side scan sonar images, the magnetic anomaly map proved to be most useful in delineating the extent of the town. These anomalies exhibit short wavelength components in the area surrounding a submerged church, with the shortest wavelength components being confined to the area immediately east of the church. These short wavelength components are only observed near the church; therefore, we propose that they delineate the buried remnants of the submerged town.     

Hydro-climatic variability and forest fires in Mexico’s northern temperate forests

Global warming is likely modifying the hydrological cycle of forested watersheds. This report set as objectives to: a) assess the hydrological variables interception loss, I, potential and actual evapo-transpiration, E, Et, runoff, Q, and soil moisture content, θ; b) evaluate whether these variables are presenting consistent trends or oscillations that can be associated to global warming or climate variability; and c) relate θ to the number of wildfires and the burned area in Durango, Mexico. A mass balance approach estimated daily variables of the water cycle using sub-models for I and Et to calculate Q and θ for a time series from 1945 to 2007. Regression and auto-regressive and moving averaging (ARIMA) techniques evaluated the statistical significance of trends. The cumulative standardized z value magnified and ARIMA models projected statistically similar monthly and annual time series data of all variables of the water cycle. Regression analysis and ARIMA models showed monthly and annual P, I, E, and Et, Q, and θ do not follow consistent up or downward linear tendencies over time with statistical significance; they rather follow oscillations that could be adequately predicted by ARIMA models (r² ≥ 0.70). There was a consistent statistical association (p ≤ 0.05) of θ with the number of wildfires and the area burned regardless of the different spatial scales used in evaluating these variables. The analysis shows seasonal variability is increasing over time as magnifying pulses of dryness and wetness, which may be the response of the hydrological cycle to climate change. Further research must center on using longer time series data, testing seasonal variability with additional statistical analysis, and incorporating new variables in the analysis.     

Effect of galvanic distortions on the series and parallel magnetotelluric impedances and comparison with other responses

The series and parallel impedances of the magnetotelluric tensor are appraised in relation to their relative immunity to galvanic electric distortions. The distorted responses are modeled using the Groom-Bailey decomposition of the tensor in terms of twist, shear, statics and strike direction. These four parameters and the undistorted responses are normally considered as unknowns, and are obtained from field data through the solution of an inverse problem. In the present work we use the decomposition as a forward model to simulate distorted sounding curves. Starting with undistorted 2-d te and tm responses, the tensor is distorted by assuming arbitrary values of twist, shear, static and strike direction. By default, both series and parallel responses are immune to the strike direction because they are invariants under rotation. In addition, series responses are immune to twist and shear and parallel responses only to twist. The dependence of the latter on shear is in the form of a real factor that shifts downwards the amplitude curves. On the other hand, the effect of statics on both series and parallel responses is more complicated than that on the impedance tensor because it cannot be accounted for by a simple shift of the curves. On the whole, there is a positive balance on the part of the series and parallel impedances over the te and tm responses because some of the distortions are filtered out by the invariants. It is shown that invariance is not sufficient to be immune to any of the distortions. The example chosen is Eggers’ eigenvalues, which are immune only to the by-the-fault strike direction. Invariance is not necessary either, as evidenced by the phase tensor, whose elements depend on strike but are immune to all distortions. The derivations are illustrated using soundings from the synthetic coprod2s1 and field-recorded coprod2 and bc87 data sets.     

Announcements

Abstract

This paper investigates for a 25-year period the sediment distribution in a semi-arid Brazilian basin (2 × 10⁴ km²) with a network containing more than 4000 surface reservoirs. The methodology is based on rating curves and fitted parameters derived from field data. The results showed that suspended load corresponded to 70% of the total sediment yield (148 t km^-2 year^-1). The relatively low contribution of the suspended load (compared with other semi-arid regions) was attributed to the impact of the numerous upstream reservoirs, which retained 235 t km^-2 year^-1. The micro (<1 hm³), small (1–10 hm³), medium-sized (10–50 hm³), and large or strategic (>50 hm³) reservoirs responded to, respectively, 5, 17, 30 and 48% of the total sediment retention by the reservoir network. This indicates that retention in the non-strategic reservoirs has a positive impact on water availability, since siltation of the strategic reservoirs would be expected to more than double if only such reservoirs existed.

Citation Lima Neto, I. E., Wiegand, M. C. &; de Araújo, J. C. (2011) Sediment redistribution due to a dense reservoir network in a large semi-arid Brazilian basin. Hydrol. Sci. J. 56(2), 319–333.     

Comparison of Different Techniques of Tilt Testing and Basic Friction Angle Variability Assessment

A relevant parameter for estimating discontinuity shear strength is the basic friction angle, usually derived from different types of tilt tests. However, the tilt tests described in the literature produce varying basic friction angle values. From a large number of different types of tilt tests on different kinds of rocks, it was possible to conclude that the mechanisms of sliding along cylinder generatrixes and planar surfaces are quite different, and that tests based on sliding along generatrixes are not appropriate for determining reliable basic friction angle values for discontinuity planes. Tests on small specimens are also not recommended, for geometry reasons and because ensuring reliable stress conditions is difficult. To quantify the natural variability in tilt testing, large specimens of the same granite were tested. The results revealed coefficients of variation for the basic friction angle in the range of 5–10?%, a variability which is no greater than that found for other rock mechanics parameters. This observation enables to forward some recommendations concerning the appropriate number of tests needed to obtain reliable results.