Colonization and evolution of lakes on the central Norwegian coast following deglaciation and land uplift 9500 to 7800 years B.P.

Invertebrate colonization of lakes following the uplift of land from the sea was studied in four lakes, currently situated between 39 and 24 m a.s.l., on the central Norwegian coast. The lakes were isolated from the sea between 9500 and 7700 years B.P. Animal and algal remains picked from core samples showed that the first colonizers preserved as fossils were usually members of the Chironomidae, Daphnidae/Chydoridae, Acarina, Porifera (Ephydatia mülleri and Spongilla lacustris), Bryozoa (Cristatella mucedo and Plumatella spp.) and Charophyta (Chara sp.). Of the chironomids, the genus Chironomus was present in the oldest lacustrine layers of all four lakes, but other genera recorded at the marine/lacustrine boundary were Dicrotendipes, Procladius (?), Einfeldia, Microtendipes, and Glyptotendipes. Remains of the caddis fly family Limnephilidae were also present in the earliest lacustrine sediments in Kvennavatnet and Kvernavatnet. The oldest invertebrate fauna is typical for mesotrophic lakes. However, chironomids and mites have been present in this area from at least about 10?500 years B.P. A diverse chironomid community was established between 300 and 800 years after isolation from the sea at Kvernavatnet on the island of Hitra, while only between 80 and 120 years passed before a comparably diverse community developed at Kvennavatnet on the mainland coast. A similar development of the invertebrate fauna occurred in Kvennavatnet, Kvernavatnet and Storkuvatnet. However, Litjvatnet deviates greatly from the ‘normal’ pattern because a tsunami disturbed the bottom sediments and fauna. The tsunami, a gigantic sea wave, was caused by a submarine slide from the Norwegian continental slope. It reached Litjvatnet, today located 24 m a.s.l., but was not traced in Storkuvatnet at 30 m a.s.l. This event happened about 7200 years B.P.     

Relative Height Estimation by Cross-Correlating Ground-Range Synthetic Aperture Sonar Images

The relative height of the seafloor can be estimated by using two vertically displaced receivers. In this paper, we propose techniques to improve the accuracy of the estimated height. Our results are based on the use of synthetic aperture sonar (SAS) imaging, which implies coherent addition of complex images acquired from a moving platform. The SAS processing improves the along-track (or azimuth) resolution, as well as the signal-to-noise ratio (SNR), which in turn improves the estimated height accuracy. We show that the shift of the effective center frequency induced by coherent, frequency-dependent scattering affect the time-delay estimates from complex cross correlations, and we propose a correction technique for broadband signals with uneven magnitude spectra. To reduce the effect of geometrical decorrelation and increase the coherence between the images, we beamform the sonar images onto an a priori estimate of the seafloor height before correlating. We develop a mathematical model for the imaging geometry. Finally, we demonstrate our proposed estimators by providing relative seafloor height estimates from real aperture and SAS images, obtained during the InSAS-2000 experiment at Elba Island in Italy. In particular, we demonstrate that the SAS image quality is significantly improved by inclusion of the height estimates as a priori information.     

Late glacial fans in the eastern Skagerrak; depositional environment interpreted from swath bathymetry and seismostratigraphy

The origin of acoustically transparent fan deposits overlying glacial till and ice-proximal sediments on the southern margin of the Norwegian Channel has been studied using high-resolution seismic-reflection profiles and multibeam bathymetry. The first deposits overlying glacigenic sediments are a series of stacked, acoustically transparent submarine fans. The lack of glaciomarine sediments below and between individual fans indicates that deposition was rapid and immediately followed the break up of the Late Weichselian ice cover. The fans are overlain by stratified glaciomarine sediments and Holocene mud. Because of the uniformity of this drape, the upper surface of the fan deposits is mimicked at the present seafloor, and the bathymetric images clearly show the spatial relationship of the fans to bedrock ridges and the presence of braided channel-levee systems on the surface of the youngest fans. The acoustically transparent character of the fan deposits indicates that they comprise silt and clay, and their lobate form and lack of internal stratification indicates that they were deposited by debris flows. The channel-levee morphology indicates deposition from more watery hyperconcentrated fluid flows. The fan sediments were either derived from 1) erosion of Mid Weichselian lake deposits in southern Skagerrak or 2) from Late glacial ice-margin lake deposits, ponded against the Norwegian Channel ice stream, which collapsed catastrophically when the lateral support was removed as the ice disintegrated. Fans composed almost exclusively of fine-grained sediment need not, therefore, rule out an origin in a deglacial setting relatively close to the former margins of glaciers and ice sheets.     

Formation of De Geer moraines and implications for deglaciation dynamics

De Geer moraines are very common in the Møre area, western Norway. These moraines occur below the marine limit and outside the Younger Dryas ice limit and occupy tributaries that connect the main fjords through the mountain passes. During deglaciation, ice in these tributaries flowed to the major ice streams. Sections across three De Geer moraines show that the ridges are composed of diamictons and fine-grained sediment, partly in stacked sequences. The diamicton units are interpreted as being composed of water-lain tills, lodgements tills and subaqueous flow deposits. The fine-grained sediment is though to have formed in a proglacial marine environment. Clast fabric of diamictons and deformation structures in underlying sands show that depositional directions for diamicton units and the direction of deformation for the sands is perpendicular to the ridge crests. Mainly based on this evidence, the ridges are thought to have formed by push at the glacier grounding line. The formation of transverse ridges (relative to ice flow) do occur in basal crevasses on modern glaciers, as do swarms of ridges along the front of retreating glaciers. The first mechanism of deposition does not seem to explain the ridges studied in the present paper and hence the importance of this process in the formation of De Geer moraines is questioned. The De Geer moraines were deposited by ice lobes advancing from one main fjord into another; therefore by studying the drainage pattern of the tributary lobes and their sequence of deglaciation, many features of the style of deglaciation of the ice sheet across the area can be determined. The northwestern part of the area was deglaciated earliest. After that, deglaciation proceeded to the southwest parallel to the coast. Subsequently the outer and the central part of Romsdalsfjorden were deglaciated causing ice to drain towards this fjord from both the north and south. The last fjord to be deglaciated was Storfjorden in the south.     

Reliable estimation of orbit errors in spaceborne SAR interferometry

An approach to improve orbital state vectors by orbit error estimates derived from residual phase patterns in synthetic aperture radar interferograms is presented. For individual interferograms, an error representation by two parameters is motivated: the baseline error in cross-range and the rate of change of the baseline error in range. For their estimation, two alternatives are proposed: a least squares approach that requires prior unwrapping and a less reliable gridsearch method handling the wrapped phase. In both cases, reliability is enhanced by mutual control of error estimates in an overdetermined network of linearly dependent interferometric combinations of images. Thus, systematic biases, e.g., due to unwrapping errors, can be detected and iteratively eliminated. Regularising the solution by a minimum-norm condition results in quasi-absolute orbit errors that refer to particular images. For the 31 images of a sample ENVISAT dataset, orbit corrections with a mutual consistency on the millimetre level have been inferred from 163 interferograms. The method itself qualifies by reliability and rigorous geometric modelling of the orbital error signal but does not consider interfering large scale deformation effects. However, a separation may be feasible in a combined processing with persistent scatterer approaches or by temporal filtering of the estimates.     

Fault displacement modelling using 3D vector fields

In history matching and sensitivity analysis, flexibility in the structural modelling is of great importance. The ability to easily generate multiple realizations of the model will have impact both on the updating workflow in history matching and uncertainty studies based on Monte Carlo simulations. The main contribution to fault modelling by the work presented in this paper is a new algorithm for calculating a 3D displacement field applicable to a wide range of faults due to a flexible representation. This gives the possibility to apply this field to change the displacement and thereby moving horizons and fault lines. The fault is modelled by a parametric format where the fault has a reference plane defined by a centre point, dip and strike angles. The fault itself is represented as a surface defined by a function z = f(x,y), where x, y and z are coordinates local to the reference plane, with the z-axis being normal to the plane. The displacement associated with the fault outside the fault surface is described by a 3D vector field. The displacement on the fault surface can be found by identifying the intersection lines between horizons and the fault surface (fault lines), and using kriging techniques to fill in values at all points on the surface. Away from the fault surface the displacement field is defined by a monotonic decreasing function which ensures zero displacement at a specified distance from the fault. An algorithm is developed where the displacement can be increased or decreased according to user-defined parameters. This means that the whole displacement field is changed and points on horizons around the fault can be moved accordingly by applying the modified displacement field on them. The interaction between several faults influencing the same points is taken care of by truncation rules and the ordering of the faults. The method is demonstrated on a realistic synthetic case based on a real reservoir.     

A quantitative study on the use of converted waves for sub-basalt imaging

The idea of imaging beneath a high-velocity layer using converted waves has been popular since 1990. Because these wave types have their maximum amplitudes at mid- to far-offsets, the search for pure P-waves at the highly multiple-contaminated near-offsets can be avoided. For the Atlantic Margin, with buried thin-layered basalts, our quantitative study shows that the initial single-layered approach is not viable. Even in an unrealistic ideal geological setting, the amplitude of the symmetrical PSP-mode is far too weak to be recognized on towed streamer data. Furthermore, in the far-offset window, where locally converted waves have their strongest amplitudes, there is a multitude of other reflections, refractions and interbedded multiples, which have similar moveouts and, often, higher amplitudes. Without the removal of these events, a reliable image of the subsurface cannot be produced. We show that even if this problem were solved, it would be far easier to use the P-wave reflection from beneath the basalt at near-offsets. Our study shows that this wave type is by far the strongest response. A borehole-derived model using a thin-layered basalt sequence reveals that the strongest locally converted wave has an asymmetrical path and is 10 times weaker. All our results indicate that the pure P-modes provide the best chance of imaging sub-basalt sedimentary interfaces.