Controls on the location,morphology and evolution of complex esker systems at decadal timescales,Breiðamerkurjökull,southeast Iceland

This paper uses detailed mapping of eskers to address three questions which are important for reconstructing meltwater behaviour beneath contemporary and ancient ice masses: ‘What controls the morphology of simple and complex esker systems?’, ‘How do esker systems evolve through time?’ and ‘Are esker patterns compatible with groundwater controlled hydraulic spacing of esker tunnels?’. Esker crestlines and widths are mapped on the Breiðamerkurjökull foreland for eight time slices between 1945 and 2007, from high resolution (~50 cm) aerial photography, permitting their long‐term morphological evolution to be analysed in a high level of detail. We find that complex eskers develop where meltwater and sediment is abundant, such that sediment clogs channels, forming distributary eskers. Isolated eskers are simpler and smaller and reflect less abundant meltwater and sediment, which is unable to clog channels. Eskers may take several decades to emerge from outwash deposits containing buried ice and can increase or decrease in size when ice surrounding and underlying them melts out. It has been suggested that groundwater–channel coupling dictates the spacing between eskers at Breiðamerkurjökull. Our results do not dispute this, but suggest that the routing of sediment and meltwater through medial moraines is an additional important control on esker location and spacing. These results may be used to better understand the processes surrounding esker formation in a variety of geographical settings, enabling a more detailed understanding of the operation of meltwater drainage systems in sub‐marginal zones beneath glaciers and ice sheets. Copyright © 2015 John Wiley & Sons, Ltd.     

A METHOD TO CORRECT YIELD SURFACE DRIFT IN SOIL PLASTICITY UNDER MIXED CONTROL AND EXPLICIT INTEGRATION

When applying an explicit integration algorithm in e.g. soil plasticity, the predicted stress point at the end of an elastoplastic increment of loading might not be situated on the updated current yield surface. This so-called yield surface drift could generally be held under control by using small integration steps. Another possibility, when circumstances might demand larger steps, is to adopt a drift correction method. In this paper, a drift correction method for mixed control in soil plasticity, under drained as well as undrained conditions, is proposed. By simulating triaxial tests in a Constitutive Driver, the capability and efficiency of this correction method, under different choices of implementation, have been analysed. It was concluded that the proposed drift correction method, for quite marginal additional computational cost, was able to correct successfully for yield surface drift giving results in close agreement to those obtained with a very large number of integration steps. © 1997 by John Wiley & Sons, Ltd.     

Pareto Joint Inversion of 2D Magnetotelluric and Gravity Data — Towards Practical Applications

In this paper, a Pareto inversion based global optimization approach, to obtain results of joint inversion of two types of geophysical data sets, is formulated. 2D magnetotelluric and gravity data were used for tests, but presented solution is flexible enough to be used for combination of any kind of two or more target functions, as long as misfits can be calculated and forward problems solved. To minimize dimensionality of the solution, space and introduce straightforward regularization Sharp Boundary Interface (SBI) method was applied. As a main optimization engine, Particle Swarm Optimization (PSO) was used. Synthetic examples based on a real geological model were used to test proposed approach and show its usefulness in practical applications.     

Seismic Wave Velocities in Deep Sediments in Poland: Borehole and Refraction Data Compilation

Sedimentary cover has significant influence on seismic wave travel times and knowing its structure is of great importance for studying deeper structures of the Earth. Seismic tomography is one of the methods that require good knowledge of seismic velocities in sediments and unfortunately by itself cannot provide detailed information about distribution of seismic velocities in sedimentary cover. This paper presents results of P-wave velocity analysis in the old Paleozoic sediments in area of Polish Lowland, Folded Area, and all sediments in complicated area of the Carpathian Mountains in Poland. Due to location on conjunction of three major tectonic units — the Precambrian East European Craton, the Paleozoic Platform of Central and Western Europe, and the Alpine orogen represented by the Carpathian Mountains the maximum depth of these sediments reaches up to 25 000 m in the Carpathian Mountains. Seismic velocities based on 492 deep boreholes with vertical seismic profiling and a total of 741 vertical seismic profiles taken from 29 seismic refraction profiles are analyzed separately for 14 geologically different units. For each unit, velocity versus depth relations are approximated by second or third order polynomials.     

Optimization routine for identification of model parameters in soil plasticity

The paper presents an optimization routine especially developed for the identification of model parameters in soil plasticity on the basis of different soil tests. Main focus is put on the mathematical aspects and the experience from application of this optimization routine. Mathematically, for the optimization, an objective function and a search strategy are needed. Some alternative expressions for the objective function are formulated. They capture the overall soil behaviour and can be used in a simultaneous optimization against several laboratory tests. Two different search strategies, Rosenbrock's method and the Simplex method, both belonging to the category of direct search methods, are utilized in the routine. Direct search methods have generally proved to be reliable and their relative simplicity make them quite easy to program into workable codes. The Rosenbrock and simplex methods are modified to make the search strategies as efficient and user‐friendly as possible for the type of optimization problem addressed here. Since these search strategies are of a heuristic nature, which makes it difficult (or even impossible) to analyse their performance in a theoretical way, representative optimization examples against both simulated experimental results as well as performed triaxial tests are presented to show the efficiency of the optimization routine. From these examples, it has been concluded that the optimization routine is able to locate a minimum with a good accuracy, fast enough to be a very useful tool for identification of model parameters in soil plasticity. Copyright © 2001 John Wiley & Sons, Ltd.