A giant three-stage submarine slide off Norway

One of the largest submarine slides known, The Storegga Slide, is located on the Norwegian continental margin. The slide is up to 450 m thick and has a total volume of about 5,600 km3. The headwall of the slide scar is 290 km long and the total run-out distance is about 800 km. The slide involved sediments of Quaternary to Early Tertiary age and occurred in three stages. Earthquakes combined with decomposition of gas hydrates are believed to be the main triggering agents for the slides. The first slide event is tentatively dated to be about 30,000 to 50,000 years B.P. and the two last major events are dated to be at 6,000 to 8,000 years B.P.     

The last British Ice Sheet: growth, maximum extent and deglaciation

The growth, maximum lateral extent and deglaciation of the last British Ice Sheet (BIS) has been reconstructed using sediment, faunal and stable isotope methods from a sedimentary record recovered from the Barra Fan, north-west Scotland. During a phase of ice sheet expansion postdating the early "warmth" of Marine Isotope Stage 3 (MIS 3), ice rafting events, operating with a cyclicity of approximately 1500 years, are interspersed between warm, carbonate-rich interstadials operating with a strong Dansgaard-Oeschger (D-O) cyclicity. The data suggest that the BIS expanded westwards to the outer continental shelf break shortly after 30 Ky BP (before present) and remained there until about 15 Ky BP. Within MIS 2, as the ice sheet grew to its maximum extent, the pronounced periodicities which characterize MIS 3 are lost from the record. The exact timing of the Last Glacial Maximum is difficult to define in this record; but maxima in Neogloboquadrina pachyderma (sinistral) ع⁸O are observed between 21-17 Ky BP. A massive discharge of ice-rafted detritus, coincident with Heinrich event 1, is observed at about 16 Ky BP. Deglaciation of the margin is complete by about 15 Ky BP and surface waters warm rapidly after this date.     

Knickpoint recession rate and catchment area: the case of uplifted rivers in Eastern Scotland

Knickpoint behaviour is a key to understanding both the landscape responses to a base‐level fall and the corresponding sediment fluxes from rejuvenated catchments, and must be accommodated in numerical models of large‐scale landscape evolution. Knickpoint recession in streams draining to glacio‐isostatically uplifted shorelines in eastern Scotland is used to assess whether knickpoint recession is a function of discharge (here represented by its surrogate, catchment area). Knickpoints are identified using DS plots (log slope versus log downstream distance). A statistically significant power relationship is found between distance of headward recession and catchment area. Such knickpoint recession data may be used to determine the values of m and n in the stream power law, E = KA^mSⁿ. The data have too many uncertainties, however, to judge definitively whether they are consistent with m = n = 1 (bedrock erosion is proportional to stream power and KPs should be maintained and propagate headwards) or m = 0·3, n = 0·7 (bedrock incision is proportional to shear stress and KPs do not propagate but degrade in place by rotation or replacement). Nonetheless, the E Scotland m and n values point to the dominance of catchment area (discharge) in determining knickpoint retreat rates and are therefore more consistent with the stream power law formulation in which bedrock erosion is proportional to stream power. Copyright © 2005 John Wiley & Sons, Ltd.     

Population abundance and seasonal migration patterns indicated by commercial catch-per-unit-effort of hakes (Merluccius capensis and M. paradoxus) in the northern Benguela Current Large Marine Ecosystem

We developed generalised additive models (GAMs) to estimate standardised time-series of population abundance indices for assessment purposes and to infer ecological and behavioural information on northern Benguela hakes, Merluccius capensis and M. paradoxus, using haul-by-haul commercial trawl catch-rate data as proxies for hake densities. The modelling indicated that individual ship identifiers should be used rather than general vessel characteristics, such as vessel size. The final models explained 79% and 68% of the variability in the commercial catch rates of M. capensis and M. paradoxus, respectively. The spatial density patterns were consistent and confirmed existing knowledge about these species in the northern Benguela system. Furthermore, seasonal migration patterns were described for the first time and were found to correspond to the known spawning areas and seasons for M. capensis and M. paradoxus. Spatial density patterns were validated using the geostatistical modelling results of fisheries-independent trawl survey data. Improved understanding of the relationships between fleet dynamics and fish movement can be achieved by taking into consideration the present catch-rate model and spatial and seasonal distribution maps. We conclude that the yearly standardised CPUE time-series are problematic as proxies for total stock abundance because of spatial coverage issues. Consequently, such CPUE data should not be used for stock-size assessments and fisheries advice concerning northern Benguela hakes until this is solved. We generally recommend the exclusion of standardised CPUE time-series from stock assessments when important and changing parts of the stock distribution cannot be targeted by the fishery, such as due to closed areas or seasons.     

Determination of lower and upper bounds of predicted production from history-matched models

We present a method to determine lower and upper bounds to the predicted production or any other economic objective from history-matched reservoir models. The method consists of two steps: 1) performing a traditional computer-assisted history match of a reservoir model with the objective to minimize the mismatch between predicted and observed production data through adjusting the grid block permeability values of the model. 2) performing two optimization exercises to minimize and maximize an economic objective over the remaining field life, for a fixed production strategy, by manipulating the same grid block permeabilities, however without significantly changing the mismatch obtained under step 1. This is accomplished through a hierarchical optimization procedure that limits the solution space of a secondary optimization problem to the (approximate) null space of the primary optimization problem. We applied this procedure to two different reservoir models. We performed a history match based on synthetic data, starting from a uniform prior and using a gradient-based minimization procedure. After history matching, minimization and maximization of the net present value (NPV), using a fixed control strategy, were executed as secondary optimization problems by changing the model parameters while staying close to the null space of the primary optimization problem. In other words, we optimized the secondary objective functions, while requiring that optimality of the primary objective (a good history match) was preserved. This method therefore provides a way to quantify the economic consequences of the well-known problem that history matching is a strongly ill-posed problem. We also investigated how this method can be used as a means to assess the cost-effectiveness of acquiring different data types to reduce the uncertainty in the expected NPV.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Laboratory and field investigations of rockburst phenomena using concurrent geotomographic imaging and acoustic emission/microseismic techniques

The progress made on three phases of a research project, started in 1986 to investigate mining induced seismicity/rockburst phenomena using concurrent geotomographic imaging and microseismic monitoring techniques, is described. Phase I is the geotomographic software development and laboratory calibration trials. Phase II is the enhancement of traditional microseismic monitoring instrumentation with a waveform acquisition system, so that source mechanism studies can be carried out on mining induced seismic events. Phase III is the field trials of the hybrid technique which will be used to monitor changing rock mass physical properties, in response to mining. Preliminary results from all three phases are given, together with an outline of current and future research planned.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal 1987.     

A mechanism for preferential H2O leakage from fluid inclusions in quartz,based on TEM observations

Preferential leakage of H₂O from fluid inclusions containing multiple gas components has been suspected in natural metamorphic rocks and has been demonstrated experimentally for synthetic H₂O-CO₂-rich inclusions in natural quartz. Knowledge of the physical and chemical characteristics of the leakage mechanism, which may be very complex, increases the value of natural fluid inclusions to metamorphic geology. It is proposed that crystal defects play a major role in nondecrepitative preferential H₂O leakage through quartz, and remain effective during metamorphism. Inclusions with either an internal overpressure or underpressure produce strain in the adjacent quartz crystal via the nucleation of many dislocations and planar defects (like Dauphiné twin boundaries). These defects allow preferential loss of H₂O from H₂O-CO₂-rich inclusions at supercritical conditions. The transport capacity of this leakage mechanism is enhanced by nucleation of small bubbles on defect structures. The nucleation of these bubbles seems to be a recovery process in strained crystals. Solubility gradients of quartz in water in a crystal with internally underpressurized inclusions may result in optical visible implosion halos in a three dimensional spatial arrangement, caused by the growth of small bubbles at the expense of the larger original fluid inclusion. Natural fluid inclusions from Naxos (Greece) are always associated with numerous interlinked dislocations. These dislocations may have been produced by plastic derormation or by crystal growth related processes (e.g. crack healing). The presence of small bubbles on these dislocations indicates that a similar leakage mechanism for H₂O must have occurred in these rocks.     

Multi-channel seismic reflection measurements in the Eurasian Basin, Arctic Ocean, from U.S. ice station Fram-IV

We present the first multi-channel seismic reflection data ever collected from the Eurasian Basin of the Arctic Ocean. The 200 km data set was acquired by a 20 channel sonobuoy array deployed at U.S. ice drift station FRAM-IV and operated for 34 days about 370 km north of Svalbard in April–May 1982. Cross array drift and ice floe rotation which may constitute the most serious obstacle to the advantage of multi-channel data acquisition did only occur to a minor degree during the experiment and render most of the data set suitable for processing using common mid-point binning.A 0.7–1.4 s (two-way traveltime) thick sedimentary section has been deposited over oceanic crust of mid-Oligocene age below the Barents Abyssal Plain. In the deepest part, sediments are infilling topographic lows which indicate predominantly turbidite deposition. Erosional truncations are only locally present in the central part of the section. Conformable bedforms deposited over gentle basement highs indicate a relatively stable bottom current regime since mid-Oligocene time. Thus the establishment of a deep water connection between the Arctic Ocean and lower latitude water masses appear to have had only minor effect on Eurasian Basin bottom current circulation.Extensive submarine slide scars on the north slope of Yermak Plateau show that mass waste have been a sediment source to the Barents Abyssal Plain.     

Structural evolution of the contact between two Penninic nappes (Zermatt-Saas zone and Combin zone,Western Alps) and implications for the exhumation mechanism and palaeogeography

The boundary zone between two Penninic nappes, the eclogite-facies to ultrahigh-pressure Zermatt-Saas zone in the footwall and the blueschist-facies Combin zone in the hanging wall, has been interpreted previously as a major normal fault reflecting synorogenic crustal extension. Quartz textures of mylonites from this fault were measured using neutron diffraction. Together with structural field observations, the data allow a refined reconstruction of the kinematic evolution of the Pennine nappes. The main results are: (1) the contact is not a normal fault but a major thrust towards northwest which was only later overprinted by southeast-directed normal faulting; (2) exhumation of the footwall rocks did not occur during crustal extension but during crustal shortening; (3) the Sesia-Dent Blanche nappe system originated from a continental fragment (Cervinia) in the Alpine Tethys ocean, and the Combin zone ophiolites from the ocean basin southeast of Cervinia; (4) out-of-sequence thrusting played a major role in the tectonic evolution of the Penninic nappes. An erratum to this article can be found at