Star Cluster Simulations: the State of the Art

This paper concentrates on four key tools for performing star cluster simulations developed during the last decade which are sufficient to handle all the relevant dynamical aspects. First we discuss briefly the Hermite integration scheme which is simple to use and highly efficient for advancing the single particles. The main numerical challenge is in dealing with weakly and strongly perturbed hard binaries. A new treatment of the classical Kustaanheimo-Stiefel two-body regularization has proved to be more accurate for studying binaries than previous algorithms based on divided differences or Hermite integration. This formulation employs a Taylor series expansion combined with the Stumpff functions, still with one force evaluation per step, which gives exact solutions for unperturbed motion and is at least comparable to the polynomial methods for large perturbations. Strong interactions between hard binaries and single stars or other binaries are studied by chain regularization which ensures a non-biased outcome for chaotic motions. A new semi-analytical stability criterion for hierarchical systems has been adopted and the long-term effects on the inner binary are now treated by averaging techniques for cases of interest. These modifications describe consistent changes of the orbital variables due to large Kozai cycles and tidal dissipation. The range of astrophysical processes which can now be considered by N-body simulations include tidal capture, circularization, mass transfer by Roche-lobe overflow as well as physical collisions, where the masses and radii of individual stars are modelled by synthetic stellar evolution. This revised version was published online in July 2006 with corrections to the Cover Date.     

Emplacement mechanisms of sill complexes: Information from the geochemical architecture of the Golden Valley Sill Complex,South Africa

The well-exposed Golden Valley Sill Complex, Karoo basin, South Africa, consists of four large sills (ca. 100 m thick; long axes: 13–24 km), one small sill (55–80 m thick; long axis: 4 km; forming an appendix to one of the large sills), and two large dykes (15–20 m thick; 25 and 70 km long), plus some minor intrusions. Field mapping shows physical connections between the small sill and one of the large sills, but no other connection between the large sills, or connections between the sills and the large dykes.     

The 3D facies architecture and petrophysical properties of hyaloclastite delta deposits: An integrated photogrammetry and petrophysical study from southern Iceland

Hyaloclastites develop where lava interacts with water resulting in deposits that have a unique and often complex range of petrophysical properties. A combination of eruptive style and emplacement environment dictates the size, geometry and distribution of different hyaloclastite facies and their associated primary physical properties such as porosity, permeability and velocity. To date, links between the 3D facies variability within these systems and their petrophysical properties remain poorly understood. Hjörleifshöfði in southern Iceland presents an exceptional outcrop exposure of an emergent hyaloclastite sequence >1 km wide by >200 m high and enables an investigation of the distribution of the hyaloclastite deposits at seismic scale. Within this study we present a photogrammetry-based 3D model from part of this recent hyaloclastite delta and incorporate previous work by Watton et al. (Journal of Volcanology and Geothermal Research, 2013, 250, 19) to undertake detailed facies interpretation and quantification. Laboratory petrophysical analyses were performed on 34 core plugs cut from key field facies samples, including P- and S-wave velocity, density, porosity and permeability at both ambient and confining pressure. Integration of the 3D model with the petrophysical data has enabled the production of pseudo-wireline logs and property distribution maps which demonstrate the variability of physical properties within hyaloclastite sequences at outcrop to seismic scale. Through comparison of our data with examples of older buried hyaloclastite sequences we demonstrate that the wide-ranging properties of young hyaloclastites become highly uniform in older sequences making their identification by remote geophysical methods for similar facies variations more challenging. Our study provides an improved understanding of the petrophysical property distribution within hyaloclastite sequences and forms a valuable step towards improving the understanding of similar subsurface sequences and their implications for imaging and fluid flow.     

Nested intrashelf platform clinoforms—Evidence of shelf platform growth exemplified by Lower Cretaceous strata in the Barents Sea

Two nested clinoform set types of different scales and steepness are mapped and analysed from high-resolution seismic data. Restoration of post-depositional faulting reveals a persistent pattern of small-scale, high-angle clinoforms contained within platform-scale, low-angle clinothems, showing a combined overall progradational depositional system. The large clinoforms lack a well-defined platform edge, and show a gradual increase in dip from topset to foreset. A consistent recurring stratal pattern is evident from the architecture, and is considered a result of interplay between relative sea-level change and autocyclic switching of sediment delivery focal points that brought sediment to the platform edge. This un-interrupted succession records how intra-shelf platforms prograde. Quantitative clinoform analysis may assist in determining the most influential depositional factors. Post-depositional uplift and erosion requires restoration with re-burial to maximum burial depth. Backstripping, decompaction and isostatic correction was performed assuming a range of lithologic compositions, as no wells test the lithology. Nearby wells penetrate strata basinward of the clinoforms, proving mudstone content above 50%, which in turn guide restoration values. Typical restored platform heights are 250–300 m, with correspondingly sized platform-scale clinoform heights. Typical large-scale clinoform foreset dip values are 1.3°–2.4°. Small-scale clinothems are typically 100 m thick, with restored foreset dip angles at 4.4° - > 10°. The results suggest that intrashelf platform growth occurs in pulses interrupted by draping of strata over its clinoform profile. The resultant architecture comprises small-scale clinoforms nested within platform-scale clinothems.     

Rift–shear architecture and tectonic development of the Ghana margin deduced from multichannel seismic reflection and potential field data

The Ghana margin displays one of the best-known transform margins. Studies of the margin have provided the framework for a number of conceptual models aimed at understanding transform margin development worldwide. However, the deep structure of the margin is poorly known as knowledge is based only on wide-angle refraction measurements obtained from two separate localities on the margin. Consequently, complexities in the rift–shear margin architecture have been overlooked by current interpretations of margin development. Based on combined analysis of a detailed grid of ∼2710 km multichannel (MCS) lines and potential field data, we provide new insights into the structural architecture and tectonic development of the Ghana margin. In particular, we outline the deep structure of the entire margin using a series of 2D gravity modelled transects constrained by MCS and published wide-angle data. Our study reveals more complex rift–shear margin architecture than previously envisaged. We demonstrate that the main transform boundary representing the continental extension of the Romanche Fracture Zone, is actually composed of two distinct margin segments, i.e., the ENE–WSW trending sheared margin segment of the Cote d’Ivoire-Ghana Ridge and the NE–SW trending rift-influenced sheared margin segment of the Ghana Platform. These segments evolved under varying stress regimes, and during different time intervals. West of the transform margin, divergent rifting during the Early Cretaceous initiated the development of the Deep Ivorian Basin, essentially, as a single major pull-apart structure. However, east of the shear zone, oblique rifting resulted in the development of the Eastern Ghana Slope Basin as a composite of at least two coalescing pull-apart basins displaced along strike-slip faults. Our structural interpretation of the transform boundary geometry shows that the ridge and platform margin segments were each subjected to separate thermal influences from two different migrating spreading centres. Tectonic uplift of the ridge began through transpression during mid-Albian time following a change in relative direction of plate motion from NE–SW to ENE–WSW. However, the ridge uplift was amplified by thermal heating from a previously undocumented spreading centre whose progressive westward migration along the ridge followed closely after the Albian transpressional phase. The structural architecture of the Ghana margin resulted from a combination of factors, notably, pre-existing basement structure, plate boundary geometry, the relative direction of plate motion and thermal heating.     

On the joint distribution of heights and periods of sea waves

In the design of offshore structures, a consistent method of choosing the height and period for the adopted design wave may be crucial for a reliable design. With this objective, the joint probability density function for the height and period is considered in this paper. At first the adequacy of a simple theoretical model for the joint distribution under stationary conditions is investigated, using measurements achieved during storms in Northern North Sea. In a slightly modified form the model is found to be of a reasonable accuracy as far as the highest waves are of interest. Design curves regarding simultaneous values of heights and periods are estimated by means of this model and their possible impact for design is discussed.     

Formation of glauconie in foraminiferal shells on the continental shelf off Norway

Glauconie inside shells of foraminifers, assumed to be formed during the timespan from the Last Glaciation to the present, is reported. The process is thought to take place at water temperatures of 6–7°C. The nature of the infilling process is suggested to be neoformation of clay minerals rather than transformation of pre-infilled clay.     

Chronology and extent of the Lofoten–Vesterålen sector of the Scandinavian Ice Sheet from 26 to 16 cal. ka BP

The interplay between the onshore and offshore areas during the Last Glacial Maximum and the deglaciation of the Scandinavian Ice Sheet is poorly known. In this paper we present new results on the glacial morphology, stratigraphy and chronology of Andøya, and the glacial morphology of the nearby continental shelf off Lofoten–Vesterålen. The results were used to develop a new model for the timing and extent of the Scandinavian Ice Sheet in the study area during the local last glacial maximum (LLGM) (26 to 16 cal. ka BP). We subdivided the LLGM in this area into five glacial events: before 24, c. 23 to 22.2, 22.2 to c. 18.6, 18 to 17.5, and 16.9–16.3 cal. ka BP. The extent of the Scandinavian Ice Sheet during these various events was reconstructed for the shelf areas off Lofoten, Vesterålen and Troms. Icecaps survived in coastal areas of Vesterålen–Lofoten after the shelf was deglaciated and off Andøya ice flowed landwards from the shelf. During the LLGM the relative sea level was stable until 18.5 cal. ka BP, and thereafter there was a sea‐level drop on Andøya. Thus, relative sea level (i.e. a sea level rise) does not seem to be a driving mechanism for ice‐margin retreat in this area but the fall in sea level may have had some importance for the grounding episodes on the banks during deglaciation. The positions of the grounding zone wedges (GZWs) in the troughs are related to the morphology as they are often located where the troughs narrow.     

Deglaciation of the western margin of the Barents Sea Ice Sheet — A swath bathymetric and sub-bottom seismic study from the Kveithola Trough

Kveithola Trough, an E–W trending glacial trough in the NW Barents Sea, was surveyed for the first time during the EGLACOM cruise of R/V OGS-Explora in summer 2008. Swath bathymetry shows that the seafloor is characterized by E–W trending mega-scale glacial lineations (MSGL) that record a fast flowing ice stream draining the Svalbard/Barents Sea Ice Sheet (SBIS) during the Last Glacial Maximum (LGM). MSGL are overprinted by transverse sediment ridges about 15 km apart which give rise to a staircase axial profile of the trough. Such transverse ridges are interpreted to be grounding-zone wedges (GZWs) formed by deposition of subglacial till during episodic ice stream retreat. Sub-bottom (CHIRP) and multi-channel reflection seismic data show that the present-day morphology is largely inherited from the palaeo-seafloor topography at the time of deposition of the transverse ridges, overlain by a draping glaciomarine unit which in places is over 15 m thick. Our data allow the reconstruction of depositional processes which accompanied deglaciation of the Spitsbergen Bank area. The sedimentary drape deposited on top of the GZWs is suggested to have accumulated at a very high rate, (on average in the order of 1–1.5 m ka^?1) and therefore may potentially preserve a high-resolution palaeoclimatic record of deglaciation and post-glacial conditions in this sector of the Barents Sea.     

Relationships between shelf‐edge trajectories and sediment dispersal along depositional dip and strike: a different approach to sequence stratigraphy

Analysis of shelf‐edge trajectories in prograding successions from offshore Norway, Brazil, Venezuela and West Africa reveals systematic changes in facies associations along the depositional dip. These changes occur in conjunction with the relative sea‐level change, sediment supply, inclination of the substratum and the relief of the margin. Flat and ascending trajectories generally result in an accumulation of fluvial and shallow marine sediments in the topset segment. Descending trajectories will generally result in erosion and bypass of the topset segment and deposition of basin floor fans. An investigation of incised valley fills reveals multiple stages of filling that can be linked to distinct phases of deepwater fan deposition and to the overall evolution of the margin. In the case of high sediment supply, like the Neogene Niger and Orinoco deltas, basin floor fans may develop systematically even under ascending trajectory styles. In traditional sequence stratigraphic thinking, this would imply the deposition of basin floor fans during a period of relative sea‐level highstand. Facies associations and sequence development also vary along the depositional strike. The width and gradient of the shelf and slope show considerable variations from south to north along the Brazilian continental margin during the Cenozoic. During the same time interval, the continental shelf may display high or low accommodation conditions, and the resulting stacking patterns and facies associations may be utilized to reconstruct palaeogeography and for prediction of lithology. Application of the trajectory concept thus reveals nuances in the rock record that would be lost by the application of traditional sequence stratigraphic work procedures. At the same time, the methodology simplifies the interpretation in that less importance is placed on interpretation and labelling of surface boundaries and systems tracts.