Data assimilation for a coastal area morphodynamic model: Morecambe Bay

This paper investigates the use of data assimilation in coastal area morphodynamic modelling using Morecambe Bay as a study site. A simple model of the bay has been enhanced with a data assimilation scheme to better predict large-scale changes in bathymetry observed in the bay over a 3-year period. The 2DH decoupled morphodynamic model developed for the work is described, as is the optimal interpolation scheme used to assimilate waterline observations into the model run. Each waterline was acquired from a SAR satellite image and is essentially a contour of the bathymetry at some level within the inter-tidal zone of the bay. For model parameters calibrated against validation observations, model performance is good, even without data assimilation. However the use of data assimilation successfully compensates for a particular failing of the model, and helps to keep the model bathymetry on track. It also improves the ability of the model to predict future bathymetry. Although the benefits of data assimilation are demonstrated using waterline observations, any observations of morphology could potentially be used. These results suggest that data assimilation should be considered for use in future coastal area morphodynamic models.     

Spatial variability in growth and reproduction of the Pacific oyster Crassostrea gigas (Thunberg, 1793) along the west European coast

The Pacific oyster Crassostrea gigas was introduced in Europe for commercial purposes in the mid 1960s. It was initially thought that low winter temperatures would restrain this species' reproduction and settlement; however, its present distribution in areas where no introduction has taken place suggests that natural invasion and expansion has occurred. Along the European coast, wild populations of Pacific oysters are already found from northern Germany to southern Portugal. Whether C. gigas will continue to further expand through northern waters will depend on its physiological performance. In this study, the performance of wild oyster populations has been studied in terms of growth and reproduction at three stations: La Rochelle (France; 46°N), Yerseke (Oosterschelde estuary, The Netherlands, 51°N), and Texel (Wadden Sea estuary, The Netherlands, 53°N). The French population had the lowest somatic-shell mass ratio and an increase in maximum shell length, somatic and gonadal mass was observed from France to the Netherlands. In addition, mean oocyte diameter decreased significantly from south to north. The combination of increasing gonadal mass and decreasing oocyte volume suggests an increasing reproductive output in terms of egg numbers from France to The Netherlands. Differences in temperature between locations will at least be partly responsible for the observed patterns; however, other environmental factors (such as food availability, predation pressure, sediment type and/or seston concentration) cannot be excluded. Since smaller eggs (oocytes) are thought to have a longer development time, the environmental conditions along the Dutch coast may result in increased larval dispersal and possibly in further population expansion.     

Wave runup on a concentric twin perforated circular cylinder

The regular wave interaction with a twin concentric porous circular cylinder system consisting of an inner impermeable cylinder and an outer perforated cylinder was studied through physical model and numerical model studies. The experiments were carried out on the twin concentric cylinder model in a wave flume to study the wave runup and rundown at the leading and trailing edges of the perforated cylinder. It was found that the maximum wave runup on the perforated cylinder is almost same as the incident wave height. The experimental results were used to develop the predictive formulae for the wave runup and rundown on the perforated cylinder, which can be easily used for design applications. The wave runup profiles around the perforated cylinder for different values of ka and porosities were studied numerically using Green's Identity Method. The results of the numerical study are presented and compared with the experimental measurements.     

Recruitment patterns of Serpula vermicularis L. (Polychaeta, Serpulidae) in Loch Creran, Scotland

This study aimed to contribute to conservation management of reefs of Serpula vermicularis by increasing understanding of the factors influencing larval settlement. The study was carried out in Loch Creran, which supports the most extensive known development of S. vermicularis reefs in the world. Settlement plates were deployed to examine the influence of season, depth, reef density, substrate type and orientation. Monthly deployment of plates revealed settlement of S. vermicularis to occur predominantly from mid-June to mid-October, peaking in late August to early September. Settlement of Pomatoceros spp. peaked much earlier, in late May to early June. Deployment of plates at different depths revealed a marked reduction in S. vermicularis settlement intensity between 6 and 12 m. As this corresponds with the deeper limit of the peripheral fringe of serpulid reefs in the loch, it is suggested that this limit is imposed by a depth-correlated settlement response, rather than reduction in available substrata. Comparisons of various substrata showed a preference by S. vermicularis larvae for a slate over a scallop substrate and no evidence of enhanced recruitment to occupied or unoccupied tubes of S. vermicularis, suggesting that gregarious attraction is unlikely to be a factor causing reef formation. Settlement onto the upper side of a horizontal scallop substrate was found to be insignificant in comparison with the underside or a vertically orientated scallop. Evidence for the role of light in controlling the depth and substrate-orientation preferences of S. vermicularis larvae is discussed. Based on the results of this study, recommendations are made regarding remediation of areas suffering reef damage.     

Bathymetric segregation of sea urchins on reefs of the Canarian Archipelago: Role of flow-induced forces

We examined whether adults of three species of sea urchins species (Diadema antillarum, Arbacia lixula, and Paracentrotus lividus) exhibit a consistent depth-dependent partitioning pattern on rocky reefs of the Canarian Archipelago (eastern Atlantic). Hydrodynamic experiments were carried out to quantify the resistance to flow-induced dislodgement in these three species. We tested the model that different morphology can result in habitat partitioning among these sea urchins. Abundances of D. antillarum increased with depth. In contrast, A. lixula and P. lividus showed the opposite zonation pattern, coexisting in high abundances in the shallowest depths (<5 m), and occurring in low densities in the deepest part of reefs (>7 m). Both A. lixula and P. lividus had greater adhesion-surface to body-height ratios than D. antillarum. Similarly, A. lixula and P. lividus showed a greater ability to resist flow-induced dislodgement compared with D. antillarum. The mean “velocity of dislodgement” was 300% and 50% greater for A. lixula and P. lividus, respectively, relative to D. antillarum, for any particular size. As a result, A. lixula and P. lividus are better fitted to life in high-flow environments than D. antillarum. We conclude that the risk of dislodgement by water motion likely play a relevant role in the vertical distribution patterns of these sea urchins in the eastern Atlantic.     

Geological evolution of the West Luzon Basin (South China Sea,Philippines)

Interpretation of deep 2-D multi-channel seismic data sheds insights into the geological evolution of the West Luzon Basin, Philippines. This basin is a sediment-filled trough that is located between the island of Luzon and the outer arc high of the west Luzon subduction zone. High-amplitude, low-frequency reflection bands mark the acoustic basement. The basement, at about 6 s (TWT), is dissected by normal faults with some of them being inverted in a later phase of deformation. The sedimentary successions, overlying the basement are stratified with partly chaotic structures and discontinuous reflectors. Five regional unconformities separate major stratigraphic units. Grid calculations of our seismic data reveal variations in the sedimentation pattern of the basin with a shift of the deposition centre from east to west and backwards during formation. A distinct bottom-simulating reflector is commonly observed. Because the northern boundary of the continental fragments to the South of the West Luzon Basin is unclear we speculate that the basin may be (partly) underlain by continental crust. The continental crust was affected by rifting prior to and during the opening of the South China Sea and the basin was overprinted at a later stage by a forearc structural setting when subduction was initiated.     

Depositional architecture and sequence stratigraphy of the Karoo basin floor to shelf edge succession, Laingsburg depocentre, South Africa

The Laingsburg depocentre of the SW Karoo Basin, South Africa preserves a well-exposed 1200 m thick succession of upper Permian strata that record the early filling of a basin during an icehouse climate. Uniformly fine-grained sandstones were derived from far-field granitic sources, possibly in Patagonia, although the coeval staging and delivery systems are not preserved. Early condensed shallow marine deposits are overlain by distal basin plain siltstone-prone turbidites and volcanic ashes. An order of magnitude increase in siliciclastic input to the basin plain is represented by up to 270 m of siltstone with thin sandstone turbidites (Vischkuil Formation). The upper Vischkuil Formation comprises three depositional sequences, each bounded by a regionally developed zone of soft sediment deformation and associated 20-45 m thick debrite that represent the initiation of a major sand delivery system. The overlying 300 m thick sandy basin-floor fan system (Unit A) is divisible into three composite sequences arranged in a progradational-aggradational-retrogradational stacking pattern, followed by up to 40 m of basin-wide hemipelagic claystone. This claystone contains Interfan A/B, a distributive lobe system that lies 10 m beneath Unit B, a sandstone-dominated succession that averages 150 m thickness and is interpreted to represent a toe of slope channelized lobe system. Unit B and the A/B interfan together comprise 4 depositional sequences in a composite sequence with an overall basinward-stepping stacking pattern, overlain by 30 m of hemipelagic claystone. The overlying 400 m thick submarine slope succession (Fort Brown Formation) is characterized by 10-120 m thick sand-prone to heterolithic packages separated by 30-70 m thick claystone units. On the largest scale the slope stratigraphy is defined by two major cycles interpreted as composite sequence sets. The lower cycle comprises lithostratigraphic Units B/C, C and D while the upper cycle includes lithostratigraphic Units D/E, E and F. In each case a sandy basal composite sequence is represented by an intraslope lobe (Units B/C and D/E respectively). The second composite sequence in each cycle (Units C and E respectively) is characterized by slope channel-levee systems with distributive lobes 20-30 km down dip. The uppermost composite sequence in each cycle (Units D and F respectively) are characterised by deeply entrenched slope valley systems. Most composite sequences comprise three sequences separated by thin (<5 m thick) claystones. Architectural style is similar at individual sequence scale for comparable positions within each composite sequence set and each composite sequence. The main control on stratigraphic development is interpreted as late icehouse glacio-eustasy but along-strike changes associated with changing shelf edge delivery systems and variable bathymetry due to differential substrate compaction complicate the resultant stratigraphy.     

Defining the concept of stratigraphic grade and applying it to stratal (reservoir) architecture and evolution of the slope-to-basin profile: An outcrop perspective

Stratigraphic grade is the similarity of the morphology of successive slope-to-basin profiles in a genetically related depositional system. In this article we use data collected from regional cross-sections of six depositional systems, stratal architecture derived from outcrops of the Lewis Shale (Wyoming, USA), and the Ross Sandstone (Ireland), and supplementary outcrop and subsurface data from other depositional systems to determine how stratigraphic grade relates to stratal (reservoir) architecture in deepwater systems.Four methods are developed that collectively define stratigraphic grade: (1) regional stacking patterns of fourth-order stratigraphic surfaces, (2) the relationship between the trajectory of the shelf edge (T_se) and the trajectory of the depocenter (T_d) for fourth-order stratigraphic units, (3) morphology of the slope-to-basin profiles of fourth-order stratigraphic surfaces, and (4) the similarity of the morphologies of slope-to-basin profiles of fourth-order surfaces in a system (σ_s, σ_r). Several characteristics of stratigraphic (reservoir) architecture of fourth-order stratigraphic cycles are related to stratigraphic grade: (1) longitudinal distribution of sandstone in fourth-order cycles, (2) location of maximum sandstone relative to the depocenter of fourth-order cycles, (3) lengths of fourth-order submarine fans, and (4) longitudinal and vertical distribution of architectural elements. Stratigraphic grade is thus a predictor of reservoir architecture and can thereby be used to reduce the uncertainty in the interpretation of subsurface data.The concept of stratigraphic grade is useful in understanding the stratigraphic evolution of deepwater systems. Most deepwater systems analyzed in this study initiated as out-of-grade and temporally evolved to graded systems over a time span of millions of years. Systems rarely evolve from graded to out-of-grade. First-order controls on stratigraphic grade are determined to be angle of slope, tectonically forced changes in angle of slope during deposition, and sediment supply.     

Composition, porosity, and reservoir potential of the Middle Jurassic Kashafrud Formation, northeast Iran

In the Kopet-Dagh Basin of Iran, deep-sea sandstones and shales of the Middle Jurassic Kashafrud Formation are disconformably overlain by hydrocarbon-bearing carbonates of Upper Jurassic and Cretaceous age. To explore the reservoir potential of the sandstones, we studied their burial history using more than 500 thin sections, supplemented by heavy mineral analysis, microprobe analysis, porosity and permeability determination, and vitrinite reflectance.The sandstones are arkosic and lithic arenites, rich in sedimentary and volcanic rock fragments. Quartz overgrowths and pore-filling carbonate cements (calcite, dolomite, siderite and ankerite) occluded most of the porosity during early to deep burial, assisted by early compaction that improved packing and fractured quartz grains. Iron oxides are prominent as alteration products of framework grains, probably reflecting source-area weathering prior to deposition, and locally as pore fills. Minor cements include pore-filling clays, pyrite, authigenic albite and K-feldspar, and barite. Existing porosity is secondary, resulting largely from dissolution of feldspars, micas, and rock fragments, with some fracture porosity. Porosity and permeability of six samples averages 3.2% and 0.0023 mD, respectively, and 150 thin-section point counts averaged 2.7% porosity. Reflectance of vitrinite in eight sandstone samples yielded values of 0.64-0.83%, in the early mature to mature stage of hydrocarbon generation, within the oil window.Kashafrud Formation petrographic trends were compared with trends from first-cycle basins elsewhere in the world. Inferred burial conditions accord with the maturation data, suggesting only a moderate thermal regime during burial. Some fractures, iron oxide cements, and dissolution may reflect Cenozoic tectonism and uplift that created the Kopet-Dagh Mountains. The low porosity and permeability levels of Kashafrud Formation sandstones suggest only a modest reservoir potential. For such tight sandstones, fractures may enhance the reservoir potential.     

The continent-ocean transition at the southeastern margin of the South China Sea

The South China Sea formed by magma-poor, or intermediate volcanic rifting in the Paleogene. We investigate the structure of the continent-ocean transition (COT) at its southern margin, off NW Palawan between the continental blocks of Reed Bank and the islands of Palawan and Calamian. Several surveys, recorded by the BGR from 1979 to 2008, established a comprehensive database of regional seismic lines, accompanied with magnetic and gravity profiles.We interpret two major rifted basins, extending in the NE direction across the shelf and slope, separated by a structural high of non volcanic origin.The continent-ocean transition is interpreted at the seaward limit of the continental crust, when magnetic spreading anomalies terminate some 80-100 km farther north. The area in between displays extensive volcanism - as manifest by extrusions that occasionally reach and cut the seafloor, by dykes, and by presumed basaltic lava flows - occurring after break-up.The COT is highly variable along the NW Palawan slope: One type shows a distinct outer ridge at the COT with a steep modern seafloor relief. The other type is characterised by rotated fault blocks, bounded by listric normal faults ramping down to a common detachment surface. Half-grabens developed above a strongly eroded pre-rift basement. The seafloor relief is smooth across this other type of COT.We suggest the pre-rift lithospheric configuration had major influence on the formation of the COT, besides transfer zones. Volcanic domains, confined to the north of competent crustal blocks correlate with the style of the COT.Gravity modelling revealed an extremely thinned crust across the shelf. We propose a depth-dependent extension model with crust being decoupled from mantle lithosphere, explaining the discrepancy of subsidence observed across the South China Sea region.