Measurement of the signal to noise ratio in seismic profiling

The measurement of the signal to noise ratio in continuous seismic profiling systems is carried out by means of a small general purpose digital computer equipped with an analog to digital converter. The computer is programmed to compute and compare the power spectrum of a sample of noise to the power spectrum of a sample of noise plus an echo. Power spectra of some twenty to fifty shots are averaged to yield reliable estimates of the signal to noise ratio as a function of frequency. The observed spectra agree well with spectra calculated from signals measured near to an air-gun. This measurement of signal to noise ratio as a function of frequency provides unambiguous information for the adjustment of the necessary band-pass filter thus improving the operation of a continuous seismic profiler.     

Land cover change over the last three centuries due to human activities: The availability of new global data sets

Land use and land cover change is an important driver of global change (Turner et al., 1993). It is recognized that land use change has important consequences for global and regional climates, the global biogeochemical cycles such as carbon, nitrogen, and water, biodiversity, etc. Nevertheless, there have been relatively few comprehensive studies of global, long-term historical changes in land cover due to land use. In this paper, we review the development of global scale data sets of land use and land cover change. Furthermore, we assess the differences between two recently developed global data sets of historical land cover change due to land use. Based on historical statistical inventories (e.g. census data, tax records, land surveys, historical geography estimates, etc) and applying different spatial analysis techniques, changes in agricultural land cover (croplands, pastures) were reconstructed for the last 300 years. The two data sets indicate that cropland areas expanded from 3–4 million km² in 1700 to 15–18 million km² in 1990 (mostly at the expense of forests), while grazing land area expanded from 5 million km² in 1700 to 31 million km² in 1990 (mostly at the expense of natural grasslands). The data sets disagree most over Latin America and Oceania, and agree best over North America. Major differences in the two data sets can be explained by the use of a fractional versus Boolean approach, different modelling assumptions, and inventory data sets.     

Bifurcation modelling in a meandering gravel–sand bed river

The Rhine bifurcation at Pannerden forms the major distribution point for water supply in the Netherlands, distributing not only water and sediment but also flooding risks and navigability. Its morphological stability has been a concern for centuries. We present experiences from more than two decades of numerical morphological modelling of this bifurcation with a gravel–sand bed and a meandering planform. Successive computations have shown the importance of upstream approach conditions, the necessity to include physical mechanisms for grain sorting and alluvial roughness, and the need to assume a thicker active layer of the river bed than is suggested by laboratory flume experiments using a constant discharge. The active layer must be thicker in the model to account for river bed variations due to higher‐frequency discharge variations that are filtered out in morphological modelling. We discuss limitations in calibration and verification, but argue that, notwithstanding these limitations, 2D and 3D morphological models are valuable tools, not only for pragmatic applications to engineering problems, but also for revealing the limitations of established knowledge and understanding of the relevant physical processes. The application of numerical models to the Pannerden bifurcation appeared to reveal shortcomings in established model formulations that do not pose particular problems in other cases. This application is therefore particularly useful for setting the agenda for further research. Copyright © 2012 John Wiley & Sons, Ltd.     

Two‐dimensional controlled‐source electromagnetic interferometry by multidimensional deconvolution: spatial sampling aspects

We use numerically modelled data sets to investigate the sensitivity of electromagnetic interferometry by multidimensional deconvolution to spatial receiver sampling. Interferometry by multidimensional deconvolution retrieves the reflection response below the receivers after decomposition of the fields into upward and downward decaying fields and deconvolving the upward decaying field by the downward decaying field. Thereby the medium above the receiver level is replaced with a homogeneous half‐space, the sources are redatumed to the receiver level and the direct field is removed. Consequently, in a marine setting the retrieved reflection response is independent of any effect of the water layer and the air above. A drawback of interferometry by multidimensional deconvolution is a possibly unstable matrix inversion, which is necessary to retrieve the reflection response. Additionally, in order to correctly separate the upward and the downward decaying fields, the electromagnetic fields need to be sampled properly. We show that the largest possible receiver spacing depends on two parameters: the vertical distance between the source and the receivers and the length of the source. The receiver spacing should not exceed the larger of these two parameters. Besides these two parameters, the presence of inhomogeneities close to the receivers may also require a dense receiver sampling. We show that by using the synthetic aperture concept, an elongated source can be created from conventionally acquired data in order to overcome these strict sampling criteria. Finally, we show that interferometry may work under real‐world conditions with random noise and receiver orientation and positioning errors.     

Point‐spread functions for interferometric imaging

Interferometric redatuming is a data‐driven method to transform seismic responses with sources at one level and receivers at a deeper level into virtual reflection data with both sources and receivers at the deeper level. Although this method has traditionally been applied by cross‐correlation, accurate redatuming through a heterogeneous overburden requires solving a multidimensional deconvolution problem. Input data can be obtained either by direct observation (for instance in a horizontal borehole), by modelling or by a novel iterative scheme that is currently being developed. The output of interferometric redatuming can be used for imaging below the redatuming level, resulting in a so‐called interferometric image. Internal multiples from above the redatuming level are eliminated during this process. In the past, we introduced point‐spread functions for interferometric redatuming by cross‐correlation. These point‐spread functions quantify distortions in the redatumed data, caused by internal multiple reflections in the overburden. In this paper, we define point‐spread functions for interferometric imaging to quantify these distortions in the image domain. These point‐spread functions are similar to conventional resolution functions for seismic migration but they contain additional information on the internal multiples in the overburden and they are partly data‐driven. We show how these point‐spread functions can be visualized to diagnose image defocusing and artefacts. Finally, we illustrate how point‐spread functions can also be defined for interferometric imaging with passive noise sources in the subsurface or with simultaneous‐source acquisition at the surface.     

Lewis Cliff 86360: An Antarctic L-chondrite with a terrestrial age of 2.35 million years

Abstract— We measured the concentrations of the cosmogenic radionuclides ¹⁰Be (half-life = 1.51 × 10⁶ a), ²⁶Al (7.05 × 10⁵ a) and ³⁶Cl (3.01 × 10⁵ a) in Lewis Cliff (LEW) 86360, an L-chondrite from the Lewis Cliff stranding area, East Antarctica. In addition, the concentrations and isotopic compositions of He, Ne and Ar were measured. The combined results yield a terrestrial age of 2.35 ± 0.15 Ma. Only one other stony meteorite with a similar terrestrial age (～2 Ma) is known from the Allan Hills stranding area (ALH 88019), whereas all previously dated stony meteorites from Antarctica are younger than 1 Ma. We argue that LEW 86360 spent most of its terrestrial residence time deep inside the ice, near the base of the glacier, where ice flow rates are much lower than at the surface. The terrestrial ages of LEW 86360 and ALH 88019 are consistent with existing hypotheses concerning the stability and persistence of the East Antarctic ice sheet.     

Seismic reflection and transmission coefficients of a self-similar interface

The derivation of seismic reflection and transmission coefficients is generally based on the assumption that the medium parameters behave as step functions of depth, at least in a finite region around the interface. However, outliers observed in well logs generally behave quite differently from step functions. In this paper we represent an interface by a self-similar singularity, embedded between two homogeneous half-spaces, and we derive its frequency-dependent normal-incidence reflection and transmission coefficients. For ω  → 0 the expressions for the coefficients reduce to those for a discrete boundary between two homogeneous half-spaces; for ω → ∞ they become frequency-independent. These asymptotic expressions have a relatively simple form and depend on the singularity exponent α .
  The exact as well as the asymptotic expressions are used to evaluate the time-domain reflection and transmission responses of a self-similar interface. Finally, we use a numerical method to model the response of a smoothed version of a self-similar interface (note that the velocity of a smoothed singularity remains finite). It turns out that smoothing has hardly any effect on the response, provided that the smoothing does not affect the scales corresponding to the seismic frequency range.     

Understanding and managing the morphology of branches incising into sand-clay deposits in the Dutch Rhine Delta

In the Rhine-Meuse delta in the south-western part of the Netherlands,the morphology of the river branches is highly dependent on the erodibility of the subsoil.Erosion processes that were initiated after closure of the Haringvliet estuary branch by a dam(in 1970),caused a strong incision of several connecting branches.Due to the geological evolution of this area the lithology of the subsoil shows large variations in highly erodible sand and poorly erodible peat and clay layers.This study shows how the geological information can be used to create 3D maps of the erodibility of the sub-soil, and how this information can be used to schematize the sub-soil in computational models for morphological simulations.Local incisement of sand patches between areas with poorly erodible bed causes deep scour holes,hence increasing the risk on river-bank instability(flow slides) and damage to constructions such as groynes,quays,tunnels, and pipelines.Various types of mathematical models,ranging from 1D(SOBEK) to quasi-3D(Delft3D) have been applied to study the future development of the river bed and possible management options.The results of these approaches demonstrate that models require inclusion of a layer-bookkeeping approach for sub-soil schematization, non-uniform sediment fractions(sand-mud),tidal and river-discharge boundary conditions,and capacity-reduction transport modeling.For risk-reducing river management it has been shown how the development of the river bed can be addressed on a large scale and small scale.For instance,the use of sediment feeding and fixation of bed can be proposed for large-scale management,while monitoring and interventions at initiation of erosion can be proposed as response to small-scale developments that exceed predefined intervention levels.