Modeling equilibrium bed profiles of short tidal embayments

In many tidal embayments, bottom patterns, such as the channel-shoal systems of the Wadden Sea, are observed. To gain understanding of the mechanisms that result in these bottom patterns, an idealized model is developed and analyzed for short tidal embayments. In this model, the water motion is described by the depth- and width-averaged shallow water equations and forced by a prescribed sea surface elevation at the entrance of the embayment. The bed evolves due to the divergence and convergence of suspended sediment fluxes. To model this suspended-load sediment transport, the three-dimensional advection–diffusion equation is integrated over depth and averaged over the width. One of the sediment fluxes in the resulting one-dimensional advection–diffusion equation is proportional to the gradient of the local water depth. In most models, this topographically induced flux is not present. Using standard continuation techniques, morphodynamic equilibria are obtained for different parameter values and forcing conditions. The bathymetry of the resulting equilibrium bed profiles and their dependency on parameters, such as the phase difference between the externally prescribed M₂ and M₄ tide and the sediment fall velocity, are explained physically. With this model, it is then shown that for embayments that are dominated by a net import of sediment, morphodynamic equilibria only exist up to a maximum embayment length. Furthermore, the sensitivity of the model to different morphological boundary conditions at the entrance of the embayment is investigated and it is demonstrated how this strongly influences the shape and number of possible equilibrium bottom profiles. This paper ends with a comparison between the developed model and field data for the Wadden Sea’s Ameland and Frisian inlets. When the model is forced with the observed M₂ and M₄ tidal constituents, morphodynamic equilibria can be found with embayment lengths similar to those observed in these inlets. However, this is only possible when the topographically induced suspended sediment flux is included. Without this flux, the maximum embayment length for which morphodynamic equilibria can be found is approximately a third of the observed length. The sensitivity of the model to the topographically induced sediment flux is discussed in detail.     

Seasonal variability in somatic and reproductive investment of the bivalve Scrobicularia plana (da Costa, 1778) along a latitudinal gradient

Monthly investment in soma and gonads in the bivalve Scrobicularia plana is described for three populations along its distributional range: Minho estuary, Portugal; Westerschelde estuary, The Netherlands and Buvika estuary, Norway. Seasonal cycles in body mass (BMI), somatic mass (SMI) and gonadal mass (GMI) indices were observed for all populations. In Portugal, BMI and SMI peaked in mid-autumn, while in The Netherlands both indices were at their highest in mid-spring. Norway showed a different pattern with two distinct peaks: one in mid-autumn and a second peak in spring. GMI reached maximum values in July in Portugal and Netherlands and in June in Norway. Overall, mean BMI and SMI were lower in Portugal while mean GMI was lower in Norway. The spawning period lasted the whole summer in Portugal, but was shorter (only two months) in The Netherlands and Norway. The reproductive investment in The Netherlands was significantly higher than in Portugal and Norway, with the lowest values being observed in Norway. Differences in annual cycles between populations were attributed to environmental factors, namely temperature and food availability. Temperature seems important in shaping the reproductive pattern with more northern populations showing shorter reproductive periods starting later in the year, and a lower reproductive output. In addition, winter water temperatures can explain the lower mean body and somatic mass values observed in Portugal. Food availability influenced the physiological performance of the species with peaks in somatic mass coinciding with phytoplankton blooms. This relation between physiological performance and environmental factors influences S. plana distribution, densities and even survival, with natural consequences on its commercial importance.     

Reconstructing recent sedimentation in two urbanised coastal lagoons (NSW,Australia) using radioisotopes and geochemistry

In this study, we combined grain size and geochemical analyses with radioisotope analysis of lead-210 (²¹⁰Pb), caesium-137 (¹³⁷Cs) and radiocarbon (¹⁴C) ages to reconstruct the sedimentation history of two urbanised coastal lagoons in south-east Australia. Towradgi and Fairy Lagoons were both found to exhibit slow initial sedimentation of less than 1 mm year⁻¹ prior to anthropogenic influences. Land clearing in the catchments increased runoff and erosion in the creeks feeding into the estuaries, and has resulted in progradation of fluvial material into the estuarine systems with a marked increase in sedimentation to between 2 and 7 mm year⁻¹. The upper 20–50 cm of the sediment column in both lagoons contained elevated concentrations of heavy metals such as Pb, Cu, Ni and Zn. This pollution trend was found to be consistent with the history of industrialisation and urbanisation in the region, which expanded rapidly post World War 2. The total metal concentrations were consistent with other urbanised/industrialised estuaries around the world. Despite the fairly disturbed nature of these coastal lagoons, the use of ²¹⁰Pb, ¹³⁷Cs and ¹⁴C dating in combination with bulk geochemical analyses allowed detailed reconstruction of sedimentation history.     

Geomechanical modelling to assess fault integrity at the CO2CRC Otway Project,Australia

This paper presents the first published 3D geomechanical modelling study of the CO2CRC Otway Project, located in the state of Victoria, Australia. The results of this work contribute to one of the main objectives of the CO2CRC, which is to demonstrate the feasibility of CO₂ storage in a depleted gas reservoir. With this aim in mind, a one-way coupled flow and geomechanics model is presented, with the capability of predicting changes to the in situ stress field caused by changes in reservoir pressure owing to CO₂ production and injection. A parametric study investigating the pore pressures required to reactivate key, reservoir-bounding faults has been conducted, and the results from the numerical simulation and analytical analysis are compared. The numerical simulation indicates that the critical pore fluid pressure to cause fault reactivation is 1.15 times the original pressure as opposed to 1.5 times for the comparable analytical model. Possible reasons for the differences between the numerical and analytical models can be ascribed to the higher degree of complexity incorporated in the numerical model. Heterogeneity in terms of lateral variations of hydrological and mechanical parameters, effect of topography, presence of faults and interaction between cells are considered to be the main sources for the different estimation of critical pore pressure. The numerical model, which incorporates this greater complexity, is able then to better describe the state of stress that acts in the subsurface compared with a simple 1D analytical model. Moreover, the reactivation pressures depend mainly on the state of stress described; therefore we suggest that numerical models be performed when possible.     

On the mathematical stability of stratified flow models with local turbulence closure schemes

Occasionally, numerical simulations using local turbulence closure schemes to estimate vertical turbulent fluxes exhibit small-scale oscillations in space, causing the eddy coefficients to vary over several orders of magnitude on short distances. Theoretical developments suggest that these spurious oscillations are essentially due to the way the eddy coefficients depend on the vertical gradient of the model’s variables. An instability criterion is derived based on the assumptions that the artefacts under study are due to the development of small-amplitude, small time- and space-scale perturbations of a smooth solution. The relevance of this criterion is demonstrated by applying it to a series a closure schemes, ranging from the Pacanowski–Philander formulas to the Mellor–Yamada level 2.5 model.     

Tectonics and global change - inferences from Late Cenozoic subsidence and uplift patterns in the Atlantic/ Mediterranean region

Quantitative subsidence analysis for a number of rifted basins in the northern Atlantic/Mediterranean region provides evidence for rapid phases of PlioQuaternary subsidence. The observed acceleration in tectonic subsidence occurs after a phase of general quiescence in subsidence in these basins and deviates from predictions of stretching models. The latter indicate a decay of subsidence with time after Mesozoic-Tertiary basin formatioh and a slow tectonic subsidence in Plio-Quaternary times. A possible explanation for the observed patterns of anomalous subsidence could be an increase in the level of intraplate compression in the northern Atlantic region. Intraplate stress changes in the Plio-Quaternary are related to the dynamics of African/ Eurasian collision processes and a reorganization of spreading directions in the AtlAntic, possibly reflecting a plate reorganization of global nature.
It seems that the Plio-Quaternary record reflects a period of increased levels of neotectonic activity, interplaying with periods of (de)glaciation. Stress-induced topography in the onshore parts of continental margins, coupled with the stress-induced subsidence in the offshore deeper parts of the basins, could have contributed to recent phases of uplift in Fennoscandia, augmenting the uplift induced by glacial unloading. Estimates of ice thicknesses are directly inferred from the observed uplift ignoring other driving mechanisms whereas topography plays a crucial role in the dynamics of glaciation. It is, therefore, important to quantdy the interplay of rapid tectonic uplift and subsidence phases with climatic effects during the Plio-Quaternary.     

Beggiatoa in microbial mats at hydrocarbon vents in the Gulf of Mexico and Warm Mineral Springs,Florida

Microbial mats were collected from a variety of sites near hydrocarbon vents along the slope in the northern Gulf of Mexico and, for comparison, from Warm Mineral Springs, Florida, USA. A predominant microorganism in each of the mats was the giant bacterium,Beggiatoa. Diameters of the bacterial filaments ranged from about 6 µm to approximately 200 µm. The latter organisms are the largest prokaryotic organisms yet found. All filaments over about 10 µm in diameter contained a large central vacuole, producing a cell with the cytoplasm as a cylindrical tube underlying the cytoplasmic membrane. Sulfur globules were confined to this peripheral layer. Push cores often contained pyrite tubules whose appearance is suggestive of aBeggiatoa origin. Determinations of ¹³C inBeggiatoa mats from vents along the Louisiana slope yielded values in the range of –26.6 to –27.9 (PDB), suggesting an unusually high degree of isotope fractionation (–24.9) relative to the carbon source in the ambient seawater, which is typical of sulfur-oxidizing chemoautotrophs. The presence of S^O (elemental sulfur) within cells ofBeggiatoa resulting from oxidation of H₂S supports the importance of bacterial sulfate reduction processes in the underlying vents for the sustenance of theBeggiatoa mats.