Cenozoic alongslope processes and sedimentation on the NW European Atlantic margin

Based on studies of sediment accumulations deposited from-and erode by-alongslope flowing ocean currents on the European continental margin from Porcupine (Ireland) to Lofoten (Norway), the evolution of the Cenozoic paleocirculation was reconstructed as part of the STRATAGEM project. There is evidence of ocean current-controlled erosion and deposition in the Rockall Trough, in the Faeroe-Shetland Channel and on the Vøring Plateau since the late Eocene, although the circulation pattern remains ambiguous. The late Palaeogene flow in the Rockall Trough was almost probably driven by southerly-derived Tethyan Outflow Water. The extent and strength of any northerly-derived flow is uncertain. From the early Neogene (early-mid-Miocene), there was a massive regional expansion of contourite drift development both in the North Atlantic and in the Norwegian-Greenland Sea. This was most probably related to the development of the Faroe Conduit, the opening of the Fram Strait and the general subsidence of the Greenland-Scotland Ridge. These may have combined to cause a considerable acceleration in the exchange and overflow of deep waters between the Arctic and Atlantic Oceans. An early late Neogene (late early Pliocene) regional erosional event has been ascribed to a vigorous pulse of bottom-current activity, most probably the result of a global reorganisation of ocean currents associated with the closure of the Central American Seaway. During the late Neogene, contourites and sediment drifts developed in deep-water basins, between units of glacigenic sediments as well as infill of several paleo-slide scars. These sediments were derived from areas of bottom-current erosion as well as from the development of Plio-Pleistocene prograding sediment wedges, incorporating the extensive sediment supply derived from shelf-wide ice sheets. Presently a profound winnowing prevails along the shelf and upper slope due to the inflowing currents of Atlantic water. Depocentres of sediments derived from the winnowing are located (locally) in lower slope embayments and in slide scars.     

Oblique sand ridges in confined tidal channels due to Coriolis and frictional torques

Ocean Dynamics - The role of the Coriolis effect in the initial formation of bottom patterns in a tidal channel is studied by means of a linear stability analysis. The key finding is that the...     

Modelling grain-size distributions. A comparison of two models and their numerical solution

Three main diffusion-based models are currently used to study grain-size distributions. In this paper, two of these approaches — perfect sorting and imperfect sorting — are compared in a parameter study. First, the numerical solution of the imperfect-sorting model is extensively discussed, and numerical tests are performed. Then, the two sedimentation models are compared for a basin under varying conditions. For some of the imposed variations, predictions of both models differ markedly due to the different approach. The position of the gravel front in the perfect sorting model depends on gravel input and proximal accommodation space. The position of grain-size boundaries in the imperfect-sorting model is strongly controlled by gravel input, the position of the basin axis and the difference in diffusivities. As a result, those two models may predict gravel progradation for different situations. Both models suggest that gravel progradation should always be coupled with sedimentation rates in order to suggest an explanation of gravel progradation observed in the geological record. Simulations with the imperfect sorting model show that this criterion may also fail, showing that a unique interpretation of gravel progradation may be impossible.     

Influence of topography on tide propagation and amplification in semi-enclosed basins

An idealized model for tide propagation and amplification in semi-enclosed rectangular basins is presented, accounting for depth differences by a combination of longitudinal and lateral topographic steps. The basin geometry is formed by several adjacent compartments of identical width, each having either a uniform depth or two depths separated by a transverse topographic step. The problem is forced by an incoming Kelvin wave at the open end, while allowing waves to radiate outward. The solution in each compartment is written as the superposition of (semi)-analytical wave solutions in an infinite channel, individually satisfying the depth-averaged linear shallow water equations on the f plane, including bottom friction. A collocation technique is employed to satisfy continuity of elevation and flux across the longitudinal topographic steps between the compartments. The model results show that the tidal wave in shallow parts displays slower propagation, enhanced dissipation and amplified amplitudes. This reveals a resonance mechanism, occurring when the length of the shallow end is roughly an odd multiple of the quarter Kelvin wavelength. Alternatively, for sufficiently wide basins, also Poincaré waves may become resonant. A transverse step implies different wavelengths of the incoming and reflected Kelvin wave, leading to increased amplitudes in shallow regions and a shift of amphidromic points in the direction of the deeper part. Including the shallow parts near the basin’s closed end (thus capturing the Kelvin resonance mechanism) is essential to reproduce semi-diurnal and diurnal tide observations in the Gulf of California, the Adriatic Sea and the Persian Gulf.     

Extreme storm surge modelling in the North Sea

This study shows that storm surge model performance in the North Sea is mostly unaffected by the application of temporal variations of surface drag due to changes in sea state provided the choice of a suitable constant Charnock parameter in the sea-state-independent case. Including essential meteorological features on smaller scales and minimising interpolation errors by increasing forcing data resolution are shown to be more important for the improvement of model performance particularly at the high tail of the probability distribution. This is found in a modelling study using WAQUA/DCSMv5 by evaluating the influence of a realistic air-sea momentum transfer parameterization and comparing it to the influence of changes in the spatial and temporal resolution of the applied forcing fields in an effort to support the improvement of impact and climate analysis studies. Particular attention is given to the representation of extreme water levels over the past decades based on the example of the Netherlands. For this, WAQUA/DCSMv5 is forced with ERA-Interim reanalysis data. Model results are obtained from a set of different forcing fields, which either (i) include a wave-state-dependent Charnock parameter or (ii) apply a constant Charnock parameter (α _{C h} =?0.032) tuned for young sea states in the North Sea, but differ in their spatial and/or temporal resolution. Increasing forcing field resolution from roughly 79 to 12 km through dynamically downscaling can reduce the modelled low bias, depending on coastal station, by up to 0.25 m for the modelled extreme water levels with a 1-year return period and between 0.1 m and 0.5 m for extreme surge heights.     

A comparison of management planning and implementation in three Indonesian marine protected areas

A comparison of three marine protected areas in Indonesia illustrates many of the problems in management planning and implementation. Major factors influencing management are proximity to urban areas, the protected area's economic potential, current uses and the resolution of jurisdictional disputes. Isolation does not prevent resource depletion. Improving the economic conditions of local residents and gaining their commitment to environmentally sustainable commercial and subsistence resource uses within and near protected areas are vital to implementing area management.     

Meridional overturning circulation: stability and ocean feedbacks in a box model

A box model of the inter-hemispheric Atlantic meridional overturning circulation is developed, including a variable pycnocline depth for the tropical and subtropical regions. The circulation is forced by winds over a periodic channel in the south and by freshwater forcing at the surface. The model is aimed at investigating the ocean feedbacks related to perturbations in freshwater forcing from the atmosphere, and to changes in freshwater transport in the ocean. These feedbacks are closely connected with the stability properties of the meridional overturning circulation, in particular in response to freshwater perturbations. A separate box is used for representing the region north of the Antarctic circumpolar current in the Atlantic sector. The density difference between this region and the north of the basin is then used for scaling the downwelling in the north. These choices are essential for reproducing the sensitivity of the meridional overturning circulation observed in general circulation models, and therefore suggest that the southernmost part of the Atlantic Ocean north of the Drake Passage is of fundamental importance for the stability of the meridional overturning circulation. With this configuration, the magnitude of the freshwater transport by the southern subtropical gyre strongly affects the response of the meridional overturning circulation to external forcing. The role of the freshwater transport by the overturning circulation (M _ov ) as a stability indicator is discussed. It is investigated under which conditions its sign at the latitude of the southern tip of Africa can provide information on the existence of a second, permanently shut down, state of the overturning circulation in the box model. M _ov will be an adequate indicator of the existence of multiple equilibria only if salt-advection feedback dominates over other processes in determining the response of the circulation to freshwater anomalies. M _ov is a perfect indicator if feedbacks other than salt-advection are negligible.     

Simulating ground water-lake interactions: approaches and insights

Approaches for modeling lake-ground water interactions have evolved significantly from early simulations that used fixed lake stages specified as constant head to sophisticated LAK packages for MODFLOW. Although model input can be complex, the LAK package capabilities and output are superior to methods that rely on a fixed lake stage and compare well to other simple methods where lake stage can be calculated. Regardless of the approach, guidelines presented here for model grid size, location of three-dimensional flow, and extent of vertical capture can facilitate the construction of appropriately detailed models that simulate important lake-ground water interactions without adding unnecessary complexity. In addition to MODFLOW approaches, lake simulation has been formulated in terms of analytic elements. The analytic element lake package had acceptable agreement with a published LAKI problem, even though there were differences in the total lake conductance and number of layers used in the two models. The grid size used in the original LAKI problem, however, violated a grid size guideline presented in this paper. Grid sensitivity analyses demonstrated that an appreciable discrepancy in the distribution of stream and lake flux was related to the large grid size used in the original LAKI problem. This artifact is expected regardless of MODFLOW LAK package used. When the grid size was reduced, a finite-difference formulation approached the analytic element results. These insights and guidelines can help ensure that the proper lake simulation tool is being selected and applied.     

Approximate solutions for radial travel time and capture zone in unconfined aquifers

Radial time-of-travel (TOT) capture zones have been evaluated for unconfined aquifers with and without recharge. The solutions of travel time for unconfined aquifers are rather complex and have been replaced with much simpler approximate solutions without significant loss of accuracy in most practical cases. The current "volumetric method" for calculating the radius of a TOT capture zone assumes no recharge and a constant aquifer thickness. It was found that for unconfined aquifers without recharge, the volumetric method leads to a smaller and less protective wellhead protection zone when ignoring drawdowns. However, if the saturated thickness near the well is used in the volumetric method a larger more protective TOT capture zone is obtained. The same is true when the volumetric method is used in the presence of recharge. However, for that case it leads to unreasonableness over the prediction of a TOT capture zone of 5 years or more.