Influence of shallow infiltration on time-lapse ERT: Experience of advanced interpretation

Previous time-lapse Electrical Resistivity Tomography (ERT) studies have experienced difficulties in reconstructing reliable calculated resistivity changes in the subsurface. Increases or decreases of resistivity appear in the calculated ERT image where no changes were noted in the subsurface, leading to erroneous hydrological interpretations of the geophysical results. In this article, we investigate how a variation of actual resistivity with time and at shallow depth can influence time-lapse ERT results and produce resistivity artefacts at depth. We use 1 and 2-D numerical modelling to simulate infiltration scenarios. Using a standard time-lapse inversion, we demonstrate the resistivity artefact production according to the electrode spacing parameter. We used an advanced inversion methodology with a decoupling line at shallow depth to attenuate or remove resistivity artefacts. We also applied this methodology to a field data set obtained in a semi-arid environment in Burkina Faso, West Africa. Here, time-lapse ERT shows several resistivity artefacts of calculated resistivity if a standard inversion is used. We demonstrate the importance of a dense sampling of shallow resistivity variations at shallow depth. Advanced interpretation allows us to significantly attenuate or remove the resistivity artefact production at intermediate depth and produce reliable interpretation of hydrological processes.     

Contribution of geophysical surveys to groundwater modelling of a porous aquifer in semiarid Niger: An overview

Subsurface geophysical surveys were carried out using a large range of methods in an unconfined sandstone aquifer in semiarid south-western Niger for improving both the conceptual model of water flow through the unsaturated zone and the parameterization of numerical a groundwater model of the aquifer. Methods included: electromagnetic mapping, electrical resistivity tomography (ERT), resistivity logging, time domain electromagnetic sounding (TDEM), and magnetic resonance sounding (MRS). Analyses of electrical conductivities, complemented by geochemical measurements, allowed us to identify preferential pathways for infiltration and drainage beneath gullies and alluvial fans. The mean water content estimated by MRS (13%) was used for computing the regional groundwater recharge from long-term change in the water table. The ranges in permeability and water content obtained with MRS allowed a reduction of the degree of freedom of aquifer parameters used in groundwater modelling.     

Variability of aerosol properties during the 2007–2010 spring seasons over central Europe

Aerosol optical properties have been studied for spring seasons when increased values of PM10 are registered. Measurements of aerosol optical properties were taken by collocated lidar and sun-photometers at Belsk, Poland, and Minsk, Belarus. A significant increase of registered aerosol optical thickness (AOT) was found during episodes with elevated PM10 concentrations. An increase of AOT at 1020 nm amounted to 50% in the case of Minsk and 18% in the case of Belsk, while an increase of AOT at 400 nm was 66% and 33%, respectively. We noted an increase of Ångström exponent by 6% at both stations and no significant increase of single scattering albedo. The LIDAR measurements together with NAAPS model results and backtrajectory analysis suggest that both the biomass burning products and the Saharan dust are responsible for increased PM10 concentrations and large AOT values during spring time. The smoke aerosol is transported over Central Europe mainly in the boundary layer, increasing both PM10 concentration and AOT. The dust aerosol transported in the free troposphere slightly affects the AOT values only. Statistically significant correlation between PM10 concentration and AOT was found during reporting period.     

Isotropic and Nonisotropic Components of Earthquakes and Nuclear Explosions on the Lop Nor Test Site, China

—?We test the hypothesis that the existence of an observable non-zero isotropic component of seismic moment can be used as a discriminant to distinguish nuclear explosions from shallow earthquakes. We do this by applying the method described herein to a small set of data recorded between 1990 and 1996 following events (seven nuclear explosions, three earthquakes) that occurred on the Lop Nor test site in Western China. We represent each source as a sum of an isotropic component at the surface and a nonisotropic, double-couple component at an estimated depth. The explosions all possess a significant non-zero isotropic component and the estimated depth of the double-couple component of the moment tensor, presumably the result of tectonic release, lies between about 0 and 3?km. For the earthquakes studied, the isotropic component is indistinguishable from zero and the depths of the sources are estimated at 3, 17 and 31?km. The data set we have studied, although still very small, suggests that certain source characteristics (namely, double-couple depth and the ratio of the isotropic to nonisotropic components of seismic moment) may prove useful in discriminating explosions from shallow earthquakes. Further work is needed to determine whether these observations hold for explosions at other test sites, to investigate a much larger set of shallow earthquakes located in regions of interest, and to study the robustness of the estimated source parameters as source magnitude and the number of observing stations decrease.     

Caveats in Multi-modal Inversion of Seismic Surface Wavefields

We consider several examples demonstrating that the formal modal representation of surface wavefields often does not describe adequately observable wave parameters, such as the phase and group velocity dispersion of higher modes. The main reason for this is the existence in the medium of several waveguides or weakly coupled wavefields in the same waveguide. In such cases the separation of neighboring higher modes may be impossible, and observed dispersion curves may significantly differ from the ones predicted by the theory. From the example related to the studies of the crustal and upper mantle structure we found that the difficulty in the separation of first and second crustal higher modes can be overcome by applying a special inversion procedure. This procedure ignores the existence of a low velocity layer in the upper mantle when fitting the observable higher-mode dispersion curve to the one predicted by the model.     

Rayleigh Wave Group Velocity Tomography of Siberia, China and the Vicinity

—Rayleigh waves are used in a tomographic inversion to obtain group velocity maps of East Asia (40° E–160° E and 20° N–70° N). The period range studied is 30 to 70 seconds. Seismograms used for this study were recorded at CDSN stations, at a temporary broadband seismic array in Tibet, at several SRO stations, and Kirnos-equipped stations established in Asia by the former Soviet Union, in Siberia, in the Sakhalin and in Mongolia. Altogether more than 1200 paths were available in the tomographic inversion. The study area includes the Angara craton, the geologically ancient core of Asia, and the subsequently accreted units, the Altaids (a Paleozoic collision complex), the Sino-Korean platform (a chain of Archaen terranes separated by belts of active structures), the south China platform (a collage of Precambrian, Paleozoic and Mesozoic metamorphic and igneous terranes), as well as the Tibetan plateau (an active tectonic feature created in late Cenozoic through collision of the Indian subcontinent and the Asian continent). Many of these main units are recognizable in the tomographic images as distinctive units; Tibet appears as a prominent low velocity (about ?15% from the average) structure, with western and central Tibet often appearing as the areas with the lowest velocities, the Central Asian fold-belt, and the Angara craton are consistently high group velocity areas. Some lesser tectonic features are also recognizable. For example, Lake Baikal is seen as a high velocity feature at periods greater than 40 seconds. However, the high group velocity feature does not stop near the southern end of Lake Baikal; it extends south-southwestward across Mongolia. The North China Plain, a part of the platform where extensional tectonics dominate, is an area of high velocities as a result of relatively thin crust. The south China block, the least tectonically active region of China, is generally an area of high velocity. For periods longer than 40 seconds, a NNE trending high group velocity gradient clearly exists in eastern China; the velocities are noticeably higher in the east. From the group velocity maps, average dispersion curves at twelve locations were determined and inverted to obtain velocity structures. Main results of group velocity inversion include: (1) a Tibetan crust of around 60?km thick, with low crustal and upper mantle shear velocities, at 3.3?km/s and 4.2?km/s, respectively; (2) with the Moho constrained at 40–43?km, the Angara craton and the Central Asian foldbelt have a V _S in excess of 4.6?km/s; (3) relatively low shear velocities are obtained for tectonically active areas. In many parts of the study area, where Precambrian basement is exposed, the process in the crust and upper mantle due to recent tectonic activities have modified the crust and upper mantle velocity structures under the Precambrian terranes, they are no longer underlain by high velocity crust and mantle.