A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

The thermal structure of the Dvurechenskii mud volcano and its implications for gas hydrate stability and eruption dynamics

The sediment temperature distribution at mud volcanoes provides insights into their activity and into the occurrence of gas hydrates. If ambient pressure and temperature conditions are close to the limits of the gas hydrate stability field, the sediment temperature distribution not only limits the occurrence of gas hydrates, but is itself influenced by heat production and consumption related to the formation and dissociation of gas hydrates. Located in the Sorokin Trough in the northern Black Sea, the Dvurechenskii mud volcano (DMV) was in the focus of detailed investigations during the M72/2 and M73/3a cruises of the German R/V Meteor and the ROV Quest 4000 m in February and March 2007. A large number of in-situ sediment temperature measurements were conducted from the ROV and with a sensor-equipped gravity corer. Gas hydrates were sampled in pressurized cores using a dynamic autoclave piston corer (DAPC). The thermal structure of the DMV suggests a regime of fluid flow at rates decreasing from the summit towards the edges of the mud volcano, accompanied by intermittent mud expulsion at the summit. Modeled gas hydrate dissociation temperatures reveal that the gas hydrates at the DMV are very close to the stability limits. Changes in heat flow due to variable seepage rates probably do not result in changes in sediment temperature but are compensated by gas hydrate dissociation and formation.     

Reproduction in the sea pen Funiculina quadrangularis (Anthozoa: Pennatulacea) from the west coast of Scotland

The sea pen Funiculina quadrangularis (Pallas, 1766) is a species of conservation concern in Scottish coastal waters, due to its restricted geographical distribution and high sensitivity to demersal fishing activities. Reproduction in F. quadrangularis was investigated in a population located in southern Loch Linnhe, west Scotland. This was accomplished through the analysis of trends in oocyte size-frequency distribution and relative fecundity over a 12-month period. Funiculina quadrangularis is dioecious and the study population exhibited a sex ratio of 1:1. Oogenesis in female F. quadrangularis is characterised by the maintenance of a large pool of asynchronously developing oocytes throughout the year, of which a small proportion (<10%) mature with increasing sychronicity and are spawned in midwinter. The reasons for this distinct pattern of oogenesis and winter spawning remain unclear, although the potential influence of environmental cues and the role of endogenous factors in relation to this sea pen's deep-sea habit are discussed. Whilst the duration of oogenesis is prolonged (>12 months), it is proposed that spawning is a brief and synchronous annual event. Relative fecundity is high and is independent of colony size, varying between approximately 500–2000 oocytes per 1 cm rachial midsection. This measure of fecundity exhibited pronounced seasonality and was significantly lower during the post-spawning winter months. Total fecundity in F. quadrangularis is considered to be high; although a small proportion of the total number of oocytes is spawned annually, this is compensated for by large colony size. Funiculina quadrangularis produces large oocytes (>800 μm), indicative of the production of lecithotrophic larvae.     

Impact of dynamic feedbacks between sedimentation,sea-level rise,and biomass production on near-surface marsh stratigraphy and carbon accumulation

Salt marshes accrete both organic and inorganic sediments. Here we present analytical and numerical models of salt marsh sedimentation that, in addition to capturing inorganic processes, explicitly account for above- and belowground organic processes including root growth and decay of organic carbon. The analytical model is used to examine the bias introduced by organic processes into proxy records of sedimentation, namely ¹³⁷Cs and ²¹⁰Pb. We find that accretion rates estimated using ²¹⁰Pb will be less than accretion rates estimated using the ¹³⁷Cs peak in steadily accreting marshes if (1) carbon decay is significant and (2) data for ²¹⁰Pb extend below the ¹³⁷Cs peak. The numerical model expands upon the analytical model by including belowground processes such as compaction and root growth, and by explicitly tracking the evolution of aboveground biomass and its effect on sedimentation rates. Using the numerical model we explore how marsh stratigraphy responds to sediment supply and the rate of sea-level rise. It is calibrated and tested using an extensive data set of both marsh stratigraphy and measurements of vegetation dynamics in a Spartina alterniflora marsh in South Carolina, USA. We find that carbon accumulation in marshes is nonlinearly related to both the supply of inorganic sediment and the rate of sea-level rise; carbon accumulation increases with sea-level rise until sea-level rise reaches a critical rate that drowns the marsh vegetation and halts carbon accumulation. The model predicts that changes in carbon storage resulting from changing sediment supply or sea-level rise are strongly dependent on the background sediment supply: if inorganic sediment supply is reduced in an already sediment poor marsh the storage of organic carbon will increase to a far greater extent than in a sediment-rich marsh, provided that the rate of sea-level rise does not exceed a threshold. These results imply that altering sediment supply to estuaries (e.g., by damming upstream rivers or altering littoral sediment transport) could lead to significant changes in the carbon budgets of coastal salt marshes.     

Evaluation of policy options for the live reef food fish trade in the province of Palawan,Western Philippines

The live reef food fish (LRFF) trade has an estimated annual retail value of US$1 billion in the Asia-Pacific region and at least US$30 million in the Philippines. Palawan Province is the trade's center accounting for as much as 55% of the commodity's total export. Concern about the LRFF trade has caused the Palawan Provincial government to take action to regulate the trade. Four key policy options to address the short-term policy objective of developing a sustainable LRFF industry through the reduction of threats associated with destructive fishing and overfishing were proposed. These four policy options and associated analyses were presented to stakeholders at three separate policy consultations. The stakeholders reached a consensus that a regulated LRFF trade was the most sensible option. As a result, the Provincial Government of Palawan took action to enact “The Palawan Live Reef Fish Ordinance of 2006”, which provides for the sustainable and integrated regulation of the LRFF industry. This provincial ordinance will assist in crafting a national policy for the LRFF trade and may also be useful in revising the 1998 Philippine Fisheries Code. The key challenge is how to translate the ordinance provisions into effective on-the-ground actions.     

The effect of a topographic high on the morphological stability of a two-inlet bay system

The cross-sectional stability of two tidal inlets connecting the same back-barrier lagoon to the ocean is investigated. The condition for equilibrium is the cross-sectional area tidal prism relationship. In an earlier study [Van de Kreeke, J., 1990. Can multiple inlets be stable? Estuarine, Coastal and Shelf Science 30: 261–273.], using the same equilibrium condition, it was concluded that where two inlets connect the same basin to the ocean ultimately one inlet will close. One of the major assumptions in that study was that the water level in the basin fluctuated uniformly. In hindsight this assumption might be too restrictive. For example, in the Wadden Sea the back barrier lagoon consists of a series of basins, rather than one single basin, separated by topographic highs. These topographic highs limit but do not exclude the exchange of water between the sub-basins. For this reason in the present study, a topographic high in the form of a weir was added, separating the back-barrier lagoon in two sub-basins. The water level in the sub-basins, rather than in the back-barrier as a whole, is assumed to fluctuate uniformly. For this schematization the hydrodynamic equations are solved using a finite difference method. The results, together with the equilibrium condition, yield the equilibrium flow curve for each of the inlets. The intersections of the two equilibrium flow curves represent combinations of cross-sectional areas for which both inlets are in equilibrium. The stability of the equilibriums was investigated using a non-linear stability analysis resulting in a flow diagram. Calculations were carried out for different inlet and weir characteristics. Sinussoidal tides were the same for both inlets. The results show that for relatively large wetted cross-sectional areas over the topographic high, approaching the situation of a single basin, there are no combinations of inlet cross-sectional areas for which both inlets are in a stable equilibrium. This supports the conclusion in the earlier study. For relatively small wetted cross-sectional areas over the topographic high there is one set of stable equilibriums. In that case the two-inlet bay system approaches that of two single-inlet bay systems.     

Analysis of threshold and incipient conditions for sediment movement

Prediction of threshold conditions and incipient motion is the essential issue for the study of sediment transport. This work compares existing empirical threshold curves proposed for Shields diagram, a method based on the concept of probability of sediment movement, and an empirical method based on movability number. These methods are used to predict the incipient motion conditions for experimental runs taken from various studies. Most of the experimental data, used in this work, have not been used before in derivation of alternative formulations for Shields diagram and other methods. The empirical threshold curves based on the Shields entrainment function was the least successful at predicting the measured incipient motion conditions, while the use of the movability number gives good predictions of critical shear velocity compared with experimental data.     

Lagrangian measurements and modelling of fluid advection in the inner surf and swash zones

New laboratory and field data are presented on fluid advection into the swash zone. The data illustrate the region of the inner surf zone from which sediment can be directly advected into the swash zone during a single uprush, which is termed the advection length. Experiments were conducted by particle tracking in a Lagrangian reference frame, and were performed for monochromatic breaking waves, solitary bores, non-breaking solitary waves and field conditions. The advection length is normalised by the run-up length to give an advection ratio, A, and different advection ratios are identified on the basis of the experimental data. The data show that fluid enters the swash zone from a region of the inner surf zone that can extend a distance seaward of the bore collapse location that is approximately equal to half of the run-up length. This region is about eight times wider than the region predicted by the classical swash solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. Journal of Fluid Mechanics 16, 113–125], as illustrated by Pritchard and Hogg [Pritchard, D., Hogg, A.J., 2005. On the transport of suspended sediment by a swash event on a plane beach. Coastal Engineering 52, 1–23]. Measured advection ratios for periodic waves show no significant trend with Iribarren number, consistent with self-similarity in typical swash flows. The data are compared to recent characteristic solutions of the non-linear shallow water wave (NLSW) equations and both finite difference and finite volume solutions of the NLSW equations.     

Morphodynamic response of a two-barred beach to a shoreface nourishment

The response to a shoreface nourishment of the two-bar system at Noordwijk (the Netherlands) is analyzed based on a daily data set of time-exposure video images collected during about 6 years, complemented with topographic and bathymetric surveys. The 1.7 Mm³ nourishment, implemented as a 3 km alongshore bump seaward of the outer bar, migrated more than 300 m onshore in 4 years before losing its integrity. Furthermore, the nourishment interrupted the autonomous seaward migration of both bars for the entire duration of the study period and, allaying earlier fears, did not intensify the three-dimensional patterns in the bars, such as the crescentic plan-shape and rip channels. The nourishment did result in clear head effects on both flanks, with the bar becoming discontinuous and the flank section decaying or becoming attached to an offshore-located bar, while the section of bar landward of the nourishment became attached to a landward-located bar. This sequence of morphologies is known as bar switching. Each switching episode took almost one year to complete and can therefore not be ascribed to individual wave events. We suspect that shoreface nourishments enhance the possibility of bar switching by creating alongshore variability in the position and depth of the outer bar and in its cross-shore migration rate and direction. The Noordwijk nourishment did not influence the shoreline position as its trend did not undergo distinctive variations after 1998. Differences in the response of the Noordwijk sandbars to the shoreface nourishment compared with other Dutch nourishments are attributed to the location and size (volume per unit length) of the nourishment with respect to the sandbars and to the median grain size of the nourished material.