Catastrophic detachment and high-velocity long-runout flow of Kolka Glacier, Caucasus Mountains, Russia in 2002

In September 2002, a catastrophic geomorphic event occurred in the Caucasus Mountains, southern Russia, in which almost the entire mass of Kolka Glacier detached from its bed, accelerated to a very high velocity (max. 65–80 m/s), and traveled a total distance of 19 km downstream as a glacier-debris flow. Based on the interpretation of satellite imagery obtained only 8.5 h before the event occurred, the analysis of seismograms from nearby seismic stations, and subsequent detailed field observations and measurements, we suggest that this remarkable event was not a response to impulse loading from a rock avalanche in the mountainside above the glacier, or to glacier surging, but due entirely to the static and delayed catastrophic response of the Kolka glacier to ice and debris loading over a period of months prior to the September 20 detachment. We reconstruct the glacier-debris flow using field observations in conjunction with the interpretation of seismographs from nearby seismic stations and successfully simulate the behaviour (runout, velocity, and deposition) of the post-detachment glacier-debris flow using a three-dimensional analytical model. Our demonstration of a standing-start hypothesis in the 2002 Kolka Glacier detachment has substantial implications for glacier hazard assessment and risk management strategies in valleys downstream from unstable debris-covered glaciers in the mountain regions of the world.     

GNSS Software Receivers

Software Global Navigation Satellite Systems (GNSS) receivers are those that implement signal correlation processing not in hardware, but in their software. The main problem for the development of real-time software (SW) multichannel GNSS receivers is the tremendous amount of calculations to perform signal correlation. The article reviews recent developments of SW GNSS receivers. The emphasis is made on the computationally effective correlation processing algorithms and the optimization of processing allocation to the receiver's hardware (HW) and SW. An architecture is suggested that implements the PRN signals despreading in a special HW preprocessor while all the other correlation processing functions are still kept in SW. The combination of the most time-consuming processing in HW, and all signal structure-dependent processing in SW, enables unique flexibility of sophisticated GNSS receiver design based on inexpensive digital signal processors. ? 2000 John Wiley & Sons, Inc.     

Shifted hyperbola moveout approximation revisited

The shifted hyperbola approximation is widely used in seismic applications. Mostly, this approximation is applied to reflection moveout in multilayered media. The traditional domain for this application is the t ‐ x domain. In this paper, we discuss the use of this approximation in the τ ‐ p and t ‐ p domains. The accuracy of the shifted hyperbola approximation defined in different domains is illustrated by analytical and numerical examples.     

Source–receiver two‐way wave extrapolation for prestack exploding‐reflector modelling and migration

Most modern seismic imaging methods separate input data into parts (shot gathers). We develop a formulation that is able to incorporate all available data at once while numerically propagating the recorded multidimensional wavefield forward or backward in time. This approach has the potential for generating accurate images free of artiefacts associated with conventional approaches. We derive novel high‐order partial differential equations in the source–receiver time domain. The fourth‐order nature of the extrapolation in time leads to four solutions, two of which correspond to the incoming and outgoing P‐waves and reduce to the zero‐offset exploding‐reflector solutions when the source coincides with the receiver. A challenge for implementing two‐way time extrapolation is an essential singularity for horizontally travelling waves. This singularity can be avoided by limiting the range of wavenumbers treated in a spectral‐based extrapolation. Using spectral methods based on the low‐rank approximation of the propagation symbol, we extrapolate only the desired solutions in an accurate and efficient manner with reduced dispersion artiefacts. Applications to synthetic data demonstrate the accuracy of the new prestack modelling and migration approach.     

A fast algorithm for 3D azimuthally anisotropic velocity scan

The conventional velocity scan can be computationally expensive for large‐scale seismic data sets, particularly when the presence of anisotropy requires multiparameter scanning. We introduce a fast algorithm for 3D azimuthally anisotropic velocity scan by generalizing the previously proposed 2D butterfly algorithm for hyperbolic Radon transforms. To compute semblance in a two‐parameter residual moveout domain, the numerical complexity of our algorithm is roughly as opposed to of the straightforward velocity scan, with N being the representative of the number of points in a particular dimension of either data space or parameter space. Synthetic and field data examples demonstrate the superior efficiency of the proposed algorithm.     

Stratigraphic coordinates: a coordinate system tailored to seismic interpretation

In certain seismic data processing and interpretation tasks such as spiking deconvolution, tuning analysis, impedance inversion, and spectral decomposition, it is commonly assumed that the vertical direction is normal to reflectors. This assumption is false in the case of dipping layers and may therefore lead to inaccurate results. To overcome this limitation, we propose a coordinate system in which geometry follows the shape of each reflector and the vertical direction corresponds to normal reflectivity. We call this coordinate system stratigraphic coordinates. We develop a constructive algorithm that transfers seismic images into the stratigraphic coordinate system. The algorithm consists of two steps. First, local slopes of seismic events are estimated by plane‐wave destruction; then structural information is spread along the estimated local slopes, and horizons are picked everywhere in the seismic volume by the predictive‐painting algorithm. These picked horizons represent level sets of the first axis of the stratigraphic coordinate system. Next, an upwind finite‐difference scheme is used to find the two other axes, which are perpendicular to the first axis, by solving the appropriate gradient equations. After seismic data are transformed into stratigraphic coordinates, seismic horizons should appear flat, and seismic traces should represent the direction normal to the reflectors. Immediate applications of the stratigraphic coordinate system are in seismic image flattening and spectral decomposition. Synthetic and real data examples demonstrate the effectiveness of stratigraphic coordinates.