Onshore and offshore seismic and lithostratigraphic analysis of a deeply incised Quaternary buried valley system in the Northern Netherlands

High-resolution seismic data (onshore and offshore), geophysical borehole data as well as detailed lithofacies from airlift boreholes were acquired in northern Netherlands on and around the island of Ameland. Marine and land seismic data combined with information from land boreholes have been explored with the objective of providing a sedimentary model. Qualitative seismic facies analysis of the valley fill commonly shows a thin unit with high amplitude reflectors at the base. Thick units of variable seismic facies (transparent to high amplitude) occur higher up in the sequence. Onlap is common at mid–upper levels within the sandy valley fill (with clay in mm layering), and a transparent seismic facies, corresponding to firm clays, is common at the top. Almost all lithological unit boundaries recognised within core parameters correspond with seismic unconformities within error margins. Subunits contain multiple cyclical trends in gamma ray and grain size. Cyclical trends show lower order fluctuations in gamma radiation on a scale of less than 1 m. Gamma-ray pattern variability between units, e.g. in general coarsening-up or fining-up units, suggests migration of subaqueous outwash fans or ice margin fluctuations. Seismic results could support a headward excavation and backfilling process suggested by Praeg [Morphology, stratigraphy and genesis of buried Elsterian tunnel valleys in the southern North Sea basin [PhD thesis]: University of Edinburgh, 207 pp.; Journal of Applied Geophysics, (this volume)] as being responsible for the formation of buried valleys. On a lithological scale, a more complicated, detailed and cyclical pattern arises. Catastrophic processes are considered unlikely as being responsible for the infill sequence because of the observed small-scale facies variability and because of the presence of diamicton layers. Diamicton layers at the base of basal unconformities as well as higher in sequence could suggest subglacial deformation by grounded ice before and during the valley-fill process.     

Towards a 4D topographic view of the Norwegian sea margin

The present-day topography/bathymetry of the Norwegian mainland and passive margin is a product of complex interactions between large-scale tectonomagmatic and climatic processes that can be traced back in time to the Late Silurian Caledonian Orogeny. The isostatic balance of the crust and lithosphere was clearly influenced by orogenic thickening during the Caledonian Orogeny, but was soon affected by post-orogenic collapse including overprinting of the mountain root, and was subsequently affected by a number of discrete extensional events eventually leading to continental break-up in Early Eocene time. In the mid-Jurassic the land areas experienced deep erosion in the warm and humid climate, forming a regional paleic surface. Rift episodes in the Late Jurassic and Early Cretaceous, with differential uplift along major fault zones, led to more pronounced topographic contrasts during the Cretaceous, and thick sequences of clastic sediments accumulated in the subsiding basins on the shelf. Following renewed extension in the Late Cretaceous, a new paleic surface developed in the Paleocene. Following break-up the margin has largely subsided thermally, but several Cenozoic shortening events have generated positive contraction structures. On the western side of the on-shore drainage divide, deeper erosion took place along pre-existing weakness zones, creating the template of the present day valleys and fjords. In the Neogene the mainland and large portions of the Barents Sea were uplifted. It appears that this uplift permitted ice caps to nucleate and accumulate during the Late Pliocene northern hemisphere climatic deterioration. The Late Pliocene to Pleistocene glacial erosion caused huge sediment aprons to be shed on to the Norwegian Sea and Barents Sea margins. Upon removal of the ice load the landmass adjusted isostatically, and this still continues today.     

The use of potential field data in revealing the basement structure in sub-basaltic settings: an example from the Møre margin, offshore Norway

This study investigates the utility of the potential fields (gravity and magnetics) in volcanic settings as observed on the Møre margin. Synthetic models are used to investigate the effect of volcanics on the gravity and magnetic fields. The focus is on detecting sub-basaltic basement structures. The methods applied to the models are Euler deconvolution on magnetic data, gravity gradients and integrated 3D gravity and magnetic forward modelling. The same methods are used on the Møre margin and the results compared to the synthetic models. The Euler deconvolution on the magnetic signal does provide limited depth solutions in the volcanic environment and the use of different observation levels does not enhance the results. Forward gravity and magnetic models provide a valuable tool to estimate both the basalt and sub-basaltic sedimentary thickness but are limited by the ambiguity inherent in potential field methods. The use of gravity gradients significantly decreases the available model solutions and provides boundary detection even in sub-basaltic settings.     

Gravity Spectra from the Density Distribution of Earth’s Uppermost 435 km

The Earth masses reside in a near-hydrostatic equilibrium, while the deviations are, for example, manifested in the geoid, which is nowadays well determined by satellite gravimetry. Recent progress in estimating the density distribution of the Earth allows us to examine individual Earth layers and to directly see how the sum approaches the observed anomalous gravitational field. This study evaluates contributions from the crust and the upper mantle taken from the LITHO1.0 model and quantifies the gravitational spectra of the density structure to the depth of 435 km. This is done without isostatic adjustments to see what can be revealed with models like LITHO1.0 alone. At the resolution of 290 km (spherical harmonic degree 70), the crustal contribution starts to dominate over the upper mantle and at about 150 km (degree 130) the upper mantle contribution is nearly negligible. At the spatial resolution \(<150\,\hbox {km},\) the spectra behavior is driven by the crust, the mantle lid and the asthenosphere. The LITHO1.0 model was furthermore referenced by adding deeper Earth layers from ak135, and the gravity signal of the merged model was then compared with the observed satellite-only model GOCO05s. The largest differences are found over the tectonothermal cold and old (such as cratonic), and over warm and young areas (such as oceanic ridges). The misfit encountered comes from the mantle lid where a velocity–density relation helped to reduce the RMS error by 40%. Global residuals are also provided in terms of the gravitational gradients as they provide better spatial localization than gravity, and there is strong observational support from ESA’s satellite gradiometry mission GOCE down to the spatial resolution of 80–90 km.     

Integrated 3D density modelling and segmentation of the Dead Sea Transform

A 3D interpretation of the newly compiled Bouguer anomaly in the area of the “Dead Sea Rift” is presented. A high-resolution 3D model constrained with the seismic results reveals the crustal thickness and density distribution beneath the Arava/Araba Valley (AV), the region between the Dead Sea and the Gulf of Aqaba/Elat. The Bouguer anomalies along the axial portion of the AV, as deduced from the modelling results, are mainly caused by deep-seated sedimentary basins (D > 10 km). An inferred zone of intrusion coincides with the maximum gravity anomaly on the eastern flank of the AV. The intrusion is displaced at different sectors along the NNW–SSE direction. The zone of maximum crustal thinning (depth 30 km) is attained in the western sector at the Mediterranean. The southeastern plateau, on the other hand, shows by far the largest crustal thickness of the region (38–42 km). Linked to the left lateral movement of approx. 105 km at the boundary between the African and Arabian plate, and constrained with recent seismic data, a small asymmetric topography of the Moho beneath the Dead Sea Transform (DST) was modelled. The thickness and density of the crust suggest that the AV is underlain by continental crust. The deep basins, the relatively large intrusion and the asymmetric topography of the Moho lead to the conclusion that a small-scale asthenospheric upwelling could be responsible for the thinning of the crust and subsequent creation of the Dead Sea basin during the left lateral movement. A clear segmentation along the strike of the DST was obtained by curvature analysis: the northern part in the neighbourhood of the Dead Sea is characterised by high curvature of the residual gravity field. Flexural rigidity calculations result in very low values of effective elastic lithospheric thickness (t _e < 5 km). This points to decoupling of crust in the Dead Sea area. In the central, AV the curvature is less pronounced and t _e increases to approximately 10 km. Curvature is high again in the southernmost part near the Aqaba region. Solutions of Euler deconvolution were visualised together with modelled density bodies and fit very well into the density model structures. An erratum to this article can be found at     

Mismatch of geophysical datasets in distal rifted margin studies

The interpretation protocols for defining offshore rifted margin architecture normally include seismic‐reflection analysis supplemented by refraction and/or potential field modelling to help constrain sedimentary, basement and Moho geometries at depth and/or the presence of magmatic material. Interpretation of modern high‐resolution long‐offset reflection profiles shows that significant mismatches may arise between the structural observations made from these data and the common translation of density, magnetic or velocity values into specific rock types made by geophysical models. We illustrate this problem with three examples from the Mid‐Norwegian rifted system, and discuss the implications with respect to the geological interpretation.     

The crustal structure of the Eastern Alps from a combination of 3D gravity modelling and isostatic investigations

The transition between European and Adriatic crust is an important feature related to the plate collision that formed the European Alps. The diversity of seismic and geological information allows the construction of two alternative 3D density models, which both match the observed gravity field. Different seismic experiments suggest a thickness for the Adriatic crust between 30 and 40 km. The thick crust model requires an unusually dense lower crust (>3050 kg/m³) to reproduce the observed Bouguer anomaly. To evaluate the two alternative models, the isostatic implications of the geometry and density distribution within both 3D models are investigated, using local (Airy) and regional (Vening Meinesz) isostasy.Airy isostatic investigations show that the Eastern Alps are not isostatically compensated and the residuals correlate strongly with exposed geological formations. Subsequently, subsurface loading is an important factor controlling isostatic processes. The different geometry and densities in the two 3D models imply different loading at the crust–mantle boundary. The subsurface loads calculated from the 3D density models were used to estimate regional isostasy by a convolution method. In general, small rigidity values (D<10×10²¹ Nm) are determined for the Eastern Alpine lithosphere. In the model with a 40-km-thick Adriatic crust, high flexural rigidities are inferred for the Adriatic plate (>100×10²¹ Nm), but these values are unusual for an active orogenic region. The results point to the interfingering of European and Adriatic crust that results in the squeezing of European crust between Adriatic crust and mantle with additional contamination by mantle material.     

A discussion of structural and thermal control of magnetic anomalies on the mid-Norwegian margin

We discuss the correlation between the depth extent of magnetic sources, the Curie temperature depth and crustal structures on the mid-Norwegian margin. Spectral methods can be used to estimate the depth extent of magnetic sources, even if the bottom is located in the lower crust, however, only with limited resolution. The bottom of the magnetic surfaces is often regarded to represent the depth to the Curie isotherm. However, comparison with a 3D model based on the interpretation of potential field and seismic reflection data and thermal modelling shows that the depth extent of the magnetic sources is merely controlled by the overall geometry of the crystalline crust and not the temperature distribution. The observed changes in the magnetic field between the inner and outer part of the mid-Norwegian margin appears not to reflect, as previously assumed, the depth to the Curie temperature but the geometry of the basement and lower crust. Our 3D model of the mid-Norwegian margin reveals a basement configuration that involves a basement with different petrophysical properties, which can be connected with lithological basement units of onshore Norway.     

New gravity maps of the Eastern Alps and significance for the crustal structures

The deep seismic profile Transalp crosses, from north to south, Germany, Austria and Italy. The gravity measurements for each country were made by national agencies with different reference systems and data reduction methods. Within the frame of the Transalp-project a comprehensive database of the Eastern Alps was compiled covering an area of 3.5° by 4° in longitude and latitude (275 by 445 km), respectively. To increase the data coverage in the south Alpine area two gravity surveys were carried out, resulting in 469 areally distributed new stations, of which 215 have been measured with the intent to improve the geoid in the area of the planned Brenner Basistunnel (BBT). The resulting gravity database is the best in terms of resolution and data quality presently available for the Eastern Alps. Here the free air, Bouguer and isostatic gravity fields are critically discussed. The spatial density of existing gravity stations in the three countries is discussed. On the Italian side of the Alps the spatial density is rather sparse compared to the Austrian side. The Bouguer-gravity field varies between − 190 * 10^− 5 m/s² and + 25 * 10^− 5 m/s², with the minimum located along the Alpine high topographic chain, but with a small offset (a few tens of km) to the greatest topographic elevation, showing that the Airy-type local isostatic equilibrium does not fully apply here. The maximum of the Bouguer anomaly has an elongated shape of 100 by 50 km located between the towns of Verona and Vicenza and covers the Venetian Tertiary Volcanic Province (VTVP), a feature not directly related to the plate collision in the Eastern Alps. The gravity high is only partly explainable by high-density magmatic rocks and requires also a deeper source, like a shallowing of the Moho. The isostatic residual anomalies (Airy model) are in the range ± 50 * 10^− 5 m/s², with the greatest positive anomaly corresponding to the location of the VTVP, indicating here under-compensation of masses. At last a discussion of a 2D density model based on reflection seismic data and receiver functions is made.