An inverse modeling study in Fram Strait. Part I: dynamics and circulation

In order to reconstruct the circulation in the northern Greenland Sea, between 77°N and 81°N, and the exchanges with the Arctic Ocean through Fram Strait, a variational inverse model is applied to the density field observed in summer 1984 during the MIZEX 84 experiment. An estimate of the three-dimensional large-scale pressure field is obtained in which the solution is decomposed into a limited number of vertical modes and the mode amplitudes are described by piece-wise polynomials on a finite-element grid. The solution should be consistent with a frictional depth-integrated vorticity balance and with the density data. The global model parameters are tuned to ensure agreement between the retrieved geostrophic velocity and independent currentmeter data. In a companion paper (Schlichtholz and Houssais, 1999b), the same method, but without dynamical constraint, is applied to the same hydrographic dataset to perform a detailed water mass analysis and to estimate individual water mass transports.A comprehensive picture of the summer geostrophic circulation in Fram Strait is obtained in which northward recirculations in the East Greenland Current (EGC) and various recirculations from the West Spitsbergen Current (WSC) to the EGC are identified. It is suggested that the branch of the WSC following the upper western slope of the Yermak Plateau turns westward beyond 81°N and recirculates southward along the lower slope, then merging with a westward recirculating branch south of 79°N. At 79°N, a southward net transport of 6.5 Sv is found in the EGC which, combined with a northward net transport of only 1.5 Sv in the WSC, results in a fairly large outflow of 5 Sv from the Arctic Ocean to the Greenland Sea.The inverse solutions show that, in summer, the local induction of vorticity by the wind stress curl or by meridional advection of planetary vorticity should be small, so that, in the EGC and in the WSC, the vorticity balance is mainly achieved between the bottom pressure torque and dissipation of vorticity through bottom friction. A substantial barotropic flow associated with along-slope potential energy gradients is indeed identified on both sides of the strait.     

Nonsteady dynamics of turbulent axisymmetric jets in stratified fluid: Part 1. Experimental study

The nonsteady dynamics of turbulent axially symmetrical jets in stratified fluid with a great density gradient (pycnocline) is studied in laboratory experiments in the large test reservoir with artificial temperature stratification at the Institute of Applied Physics, Russian Academy of Sciences (IAP RAS) simulating wastewater flows from underwater collectors in coastal cities. An underwater survey of jets in the pycnocline region has been carried out along with synchronous measurements of internal waves. It is shown that a jet quasiperiodically oscillates in a vertical plane, thus effectively generating internal waves when interacting with the pycnocline. It is ascertained that the axially symmetrical mode of jet oscillations is a source of internal waves.     

A modelling study of ecosystem dynamics and nutrient cycling in the Humber plume, UK

The European Regional Seas Ecosystem Model (ERSEM) has been coupled with a two-dimensional depth-averaged transport model of the Humber plume region and run to simulate 1988–1989. Simulations of the spatial and temporal variations in chlorophyll-a, nitrate, phosphate and suspended particulate matter distributions in winter, spring and summer show how the development of the spring bloom and subsequent maintenance of primary production is controlled by the physicochemical environment of the plume zone. Results are also shown for two stations, one characterised by the high nutrient and suspended matter concentrations of the plume and the other by the relatively low nutrient and sediment concentrations of the offshore waters. The modelled net primary production at the plume site was 105 g C m⁻² a⁻¹ and 127 g C m⁻² a⁻¹ offshore. Primary production was controlled by light limitation between October and March and by the availability of nutrients during the rest of the year. The phytoplankton nutrient demand is met by in-situ recycling processes during the summer. The likely effect of increasing and decreasing anthropogenic riverine inputs of nitrate and phosphate upon ecosystem function was also investigated. Modelling experiments indicate that increasing the nitrogen to silicate ratio in freshwater inputs increased the production of non-siliceous phytoplankton in the plume. The results of this model have been used to calculate the annual and quarterly mass balances describing the usage of inorganic nitrogen, phosphate and silicate within the plume zone for the period of the NERC North Sea survey (September 1988 to October 1989). The modelled Humber plume retains 3.9% of the freshwater dissolved inorganic nitrogen, 2.2% of the freshwater phosphate and 1.3% of the freshwater silicate input over the simulated seasonal cycle. The remainder is transported into the southern North Sea in either dissolved or particulate form. The reliability of these results is discussed.     

Numerical modeling of the arctic ocean deepwater formation: Part I. Idealized study

In this study the capability of numerical models to transport anomalously dense bottom waters from the areas of their generation on the shelf to deepwater areas is considered. These numerical models have a number of mechanisms that are favorable for this transport, but unrelated directly with the real physics of that motion. A more accurate account of the resupply of oceanic deep waters requires either a very fine resolution for the numerical model or an efficient parametrization. Two such parametrizations are proposed in this study. The first is based on amplified lateral diffusion exchange if slope instability exists. The second is related to an additional procedure allowing dense waters to leak under waters of lesser density along a shelf slope with their displacement to higher layers. The first parametrization is shown to be less efficient because it leads to a significant interaction between anomaly waters and ambient waters and to a quick loss of their original features. The second procedure is more preferable because it leads to a significant resupply of dense deep waters at the expense of dense waters generated at a distance significantly far from the shelf break. With the help of the given parametrization, during a numerical experiment these waters can be transported considerable distances away from their place of generation to the shelf boundary and deepwater areas of the ocean with a slight change in their original density characteristics.     

An offshore eddy in the California current system Part I: Interior dynamics

From January 9 to 17, 1981, detailed observations of the horizontal and vertical structure beneath one of the quasi-permanent semi-stationary mesoscale offshore eddy signatures in the California Current System (CCS) discussed by Bernstein, Breaker and Whritner (1977), Burkov and Pavlova (1980), and Simpson (1982) were made. The vertical sections of temperature and density show the presence of three-layer system. A subsurface warm-core eddy, whose diameter is about 150 km at the 7°C isotherm, is the dominant feature. A warm surface layer, which extends to a depth of 75 m, lies over the eddy. Between the warm surface layer and the subsurface warm-core eddy, there is a cold-core region which extends to a depth of about 200 m. There is a high degree of symmetry about the vertical axis of rotation. Vertical sections of salinity and dissolved oxygen are entirely different from sections of temperature and density. Diagrams of water mass characteristics confirm that the core of the eddy, found between 250–600 m, consists of inshore water from the California Undercurrent (CU). Below about 700 m, local waters from the Deep Poleward Flow (DPF) have been incorporated into the eddy. The observed distributions of properties (T, S, δ_θ, O₂) are inconsistent with a single, local generation process for the eddy system. Radial distributions of angular velocity, normalized gradient velocity and relative vorticity support the use of a Gaussian radial height field as an initial condition in eddy models. Possible reasons why CCS eddies may differ dynamically from Gulf Stream rings are given in the text. At the time the observations were made, the system as a whole was in near geostrophic balance. Local geostrophic balance, however, cannot explain the observed distribution of properties and structure. The observed symmetry in the structure of the eddy system, chemical evidence (Simpson, 1984), biological distributions (Haury, 1984) and satellite images of the CC (Koblinsky, Simpson and Dickey, 1984) suggest that lateral entrainment of warm (oceanic) and cold (coastal) water into the upper two layers of the three-layer system by the subsurface eddy is a likely generation mechanism for the cold-core region. The coastal origin of the frontal structure along the northeastern quadrant and the oceanic origin of the frontal structure along the southwestern quadrant of the eddy system further support lateral entrainment as a generation mechanism for the cold core. This entrainment makes the CCS eddy system different from cold-core rings in the Gulf Stream and rather similar to some warm-core eddies found in the East Australian Current. The presence of CU water in the core of this eddy raises the question of how CU water was transported from the continental slope. Eddy generation mechanisms, other than baroclinic instability of the CC, may be required to explain the distribution, persistence, and core composition of offshore mesoscale eddies in the CCS. There is evidence that barotropic, in addition to baroclinic, processes may be important.     

Synthesis of forced oscillations,Part I: Tidal dynamics and the influence of the loading and self-attraction effect

Free oscillations in the period range from 9 to 40 h and their corresponding adjoint solutions are used to synthesize semidiurnal and diurnal tides of second degree. The free oscillations are determined with explicit consideration of frictional terms and the full loading and self-attraction effect. This allows for a detailed analysis of the loading and self-attraction effect on the dynamics of ocean tides, e.g. a physical explanation is given for the induced phase delay computed by ocean tide models.Further, the synthesis gives a spectral composition of certain well known tidal features and pairs of free oscillations are identified, diminishing their contribution either on a global or local scale.     

Macrofaunal density, biomass and composition of estuarine sediments and their relationship to the river plume of the Rhone River (NW Mediterranean)

This study investigates species community patterns and sediment relationships of benthic macrofauna off the Rhone river delta. Along a WE transect, changes in diversity, density, biomass and trophic structure were coupled with the mean extension of the river dilution plume. Species number and diversity were at a minimum off the river mouth, below the area of the minimum surface salinity. Density decreased by a factor 2 and biomass by a factor 5 from the fluvial to the marine system. These features are due to high overload of terrestrial organic matter in the river prodelta as evidenced by the carbon isotopic signature of surface sediment and by pigment content. On the basis of a non-metric MDS analysis and of Dufrêne and Legendre method (1997), groups of stations and characteristic species associated were identified. These species, mainly small capitellids, spionids, lumbrinerids and sternaspids, correspond to a successional dynamic in response to changes in sedimentation conditions, mainly in organic matter quality of the surface sediment. The succession observed in space was similar to described in macrotidal estuaries and off other deltaic systems all over the world and to that observed in time following the Rhone river severe flood events. Results suggest that organic matter quality is an important factor with regard to benthic macrofauna successions and recovery mechanisms following disturbances. The differences observed between the Rhone deltaic system and the general model of relations between shelf processes and discharge of large rivers are attributed to a more regular supply of organic material from terrestrial origin on the Rhone continental shelf.     

Design of a prototype ocean current turbine—Part I: mathematical modeling and dynamics simulation

The “C-Plane” is a submerged variable depth ocean current turbine that is tethered to the sea floor and uses sustained ocean currents to produce electricity. As part of the development of a $\frac{1}{30}$th scale physical model of the C-Plane, a mathematical model and dynamics simulation of the prototype was developed and is presented in this paper. This three-dimensional mathematical model represents the C-Plane as a rigid body with moveable control surfaces that is moored with three linear elastic cable elements. Gravitational, buoyancy, hydrodynamic, cable, gyroscopic, and inertial forces are included and a PC-based dynamics simulation is created. The simulation demonstrates that the C-Plane is stable and capable of changing depth in all expected operating conditions. The C-Plane prototype can fly level from a height of 3 to 6 m using the configuration suggested in this paper. The maximum ascent rates of the C-Plane with a water speed of 0.3 m/s are 0.015 m/s when the pitch is fixed at 0° and 0.030 m/s when the pitch is fixed at 4°. The maximum descent rates of the C-Plane are 0.018 m/s when the pitch is held at 0° and 0.031 m/s if the pitch is held at −4°.     

Running GCM physics and dynamics on different grids: algorithm and tests

An 'alternate grid' (GridAlt) technique is presented, which allows the different components of a general circulation model's governing equations to be computed on distinct grids chosen for that component or process. In the implementation presented here, the tendencies of state variables from the physical parametrizations are computed on a vertical grid with very fine resolution near the surface, whereas the remaining terms in the equations of motion are computed using an Eta coordinate with coarser vertical resolution.
Results from a suite of aquaplanet experiments show that much of the benefit of increased vertical resolution in the whole model can be realized by enhancing the vertical resolution of the 'physics grid' using GridAlt. The benefit is realized in the fields which are computed directly in the physical parametrizations, and in the vertical structure of the relative humidity and mass streamfunction. Results from a suite of realistically configured simulations demonstrated an impact of GridAlt that was similar to its impact in the simplified simulations, as well as an improved response to El Niño Southern Oscillation forcing. It is concluded that the present implementation of GridAlt offers a practical way to allow GCMs to better capture the near-surface structure of the atmosphere.     

Submarine debris flow impact on pipelines — Part I: Experimental investigation

Estimating the impact forces exerted by a submarine debris flow on a pipeline is a challenge, and there is room for considerably more work to advance the state of the art. To this end, an experimental program was performed to investigate the impact on two pipeline installation scenarios: 1) suspended pipeline and 2) laid-on-seafloor pipeline. The results and observations from the experimental investigation are discussed. The definition of Reynolds number was modified for non-Newtonian fluids and an ad hoc method was developed to estimate the drag force exerted by an impact perpendicular to the pipe axis. The method may be used in prototype situations to estimate the drag force from submarine debris flow impact on pipelines. The experimental program was complemented by Computational Fluid Dynamics (CFD) analyses, the details of which are discussed in the accompanying paper.     

Surf zone dynamics simulated by a Boussinesq type model. Part I. Model description and cross-shore motion of regular waves

This is the first of three papers on the modelling of various types of surf zone phenomena. In this first paper, part I, the model is presented and its basic features are studied for the case of regular waves. The model is based on two-dimensional equations of the Boussinesq type and it features improved linear dispersion characteristics, possibility of wave breaking, and a moving boundary at the shoreline. The moving shoreline is treated numerically by replacing the solid beach by a permeable beach characterized by an extremely small porosity. Run-up of nonbreaking waves is verified against the analytical solution for nonlinear shallow water waves. The inclusion of wave breaking is based on the surface roller concept for spilling breakers using a geometrical determination of the instantaneous roller thickness at each point and modelling the effect of wave breaking by an additional convective momentum term. This is a function of the local wave celerity, which is determined interactively. The model is applied to cross-shore motions of regular waves including various types of breaking on plane sloping beaches and over submerged bars. Model results comprise time series of surface elevations and the spatial variation of phase-averaged quantities such as the wave height, the crest and trough elevations, the mean water level, and the depth-averaged undertow. Comparisons with physical experiments are presented. The phaseaveraged balance of the individual terms in the momentum and energy equation is determined by time-integration and quantities such as the cross-sectional roller area, the radiation stress, the energy flux and the energy dissipation are studied and discussed with reference to conventional phase-averaged wave models. The companion papers present cross-shore motions of breaking irregular waves, swash oscillations and surf beats (part II) and nearshore circulations induced by breaking of unidirectional and multidirectional waves (part III).     

Winter dynamics on the northern Portuguese shelf. Part 2: bottom boundary layers and sediment dispersal

The northern Portuguese shelf, between 41°N and 42°N, is characterised by the presence of a mid-shelf mud deposit, the Douro Mud Patch (DMP). Observations conducted between July 1996 and June 1999 under the framework of project OMEX II–II, are used to examine the impact of the dynamic processes in the lower water column, particularly on the DMP. The typical wave conditions observed during the winter maintain a highly energetic environment, capable of promoting the erosion of fine sediments at mid-shelf depths. However, the bottom nepheloid layers generated beneath these waves only extend a few meters above the bottom and are contained within rocky outcrops with similar heights that fringe the outer shelf. Each year there are about ten storms, mostly associated with southerly winds that create downwelling conditions over the shelf. The waves associated with these storms produce shear velocities over 3 cm/s at mid-shelf and bottom nepheloid layers which extend a few tens of meters above the bottom and spread offshore, over the outer shelf and upper slope. A rough estimate suggests that these events account for an offshore export of about 20×10⁶ kg of fine sediments each year (equivalent to 1–2% of sediments trapped at DMP).     

Experimental study of fluid transport processes in the matrix system of the European organic-rich shales: I. Scandinavian Alum Shale

This contribution presents results from a laboratory study investigating the fluid (gas/water) transport properties in the matrix system of the Scandinavian Alum Shale. The maturity of the organic matter of the shale samples ranged between 0.5 and 2.4% vitrinite reflectance (VR_r). Gas (He, Ar, CH₄) and water flow properties were determined at effective stresses ranging between 5 and 30 MPa and a temperature of 45 °C. The effects of different controlling factors/parameters on the fluid conductivity including permeating fluid, moisture content, anisotropy, heterogeneity, effective stress, pore pressure, and load cycling were analyzed and discussed. Pore volume measurements by helium expansion were conducted under controlled “in situ” effective stress conditions on a limited number of plugs drilled parallel and perpendicular to bedding.For Alum Shale the intrinsic permeability coefficients measured parallel and perpendicular to bedding (6·10⁻²²–8·10⁻¹⁸ m²) were within the range previously reported for other shales and mudstones. Permeability coefficients were strongly dependent on permeating fluid, moisture content, anisotropy, effective stress and other sample-to-sample variations. The intrinsic/absolute permeabilities measured with helium were consistently, higher (up to five times) than those measured with argon and methane. Permeability coefficients (He, CH₄) measured on a dry sample were up to six times higher than those measured on an “as-received” sample, depending on effective stress. The effect of moisture on measured permeability coefficients became more significant as effective stress increased. Permeability coefficients (He, CH₄) measured parallel to bedding were up to more than one order of magnitude higher than those measured perpendicular to bedding. Parallel to bedding, all samples showed a nonlinear reduction in permeability with increasing effective stress (5–30 MPa). The stress dependence of permeability could be well described by an exponential relationship.     

Allochthonous salt, structure and stratigraphy of the north-eastern Gulf of Mexico. Part I: Stratigraphy

Major sequence boundaries associated with eustatic sea level changes are correlated to the general stratigraphy of the north-eastern Gulf of Mexico. The details of a Middle Cretaceous Flooding Surface, marking a major break in sedimentation, are documented. The sequence stratigraphic work provides an example of the ‘stratigraphic signature of the Neogene’. Three major episodes of sediment accumulation are represented by: (1) Late Jurassic (150.5 Ma) to Middle Cretaceous (94 Ma) aggradation and progradation of sediments with significant sediment accumulation in the present shelf and slope areas; (2) an extended period of starved sedimentation during 94-30 Ma corresponding to Middle Cretaceous flooding events (93.5 and 91.5 Ma) and the lack of sediment supply; and (3) since Late Oligocene time, unusually rapid sedimentation rates that characterize the deep water study area. These patterns of sediment accumulation directly affect the formation of allochthonous salt in the study area.     

Photo-oxidative behaviour of crude oils relative to sea pollution: Part I. Comparative study of various crude oils and model systems

The photo-oxidation of three crude oils from Algeria, the Middle East and Venezuela, respectively, is studied under conditions simulating natural ageing (λ > 300 nm, temperatures between 16 and 38°C) of oil films on the sea surface. The kinetic behaviour differs noticeably from one to another; this could be essentially due to the presence or absence of sulphur compounds. Powerful photosensitizers or photophysical stabilizers have no influence on the course of oxidation, whereas a radical trap (hindered amine) lowers the reaction rate noticeably. The study of model systems such as distillates or a synthetic mixture show the importance of polynuclear aromatics in photo-initiation and of alkyl branched aromatics in the radical propagation of oxidation chains.     

Communication over Doppler spread channels. Part I: Channel andreceiver presentation

Scattering functions from several experiments demonstrate that acoustic underwater channels are doubly spread. Receivers used on these channels to date have difficulty with large Doppler spreads. A receiver to perform coherent communication over Doppler spread channels is presented in this first paper of two. The receiver contains a channel tracker and a linear decoder. The tracker operates by means of a modified recursive least squares algorithm which makes use of frequency-domain filters called Doppler lines. The decoder makes use of the channel tracker coefficients in order to perform minimum mean square error decoding. This first paper treats theoretical aspects whereas the second part presents implementation issues and results     

Doppler radar radial winds in HIRLAM. Part I: observation modelling and validation

An observation operator for Doppler radar radial wind measurements is developed further in this article, based on the earlier work and considerations of the measurement characteristic. The elementary observation operator treats radar observations as point measurements at pre-processed observation heights. Here, modelling of the radar pulse volume broadening in vertical and the radar pulse path bending due to refraction is included to improve the realism of the observation modelling.
The operator is implemented into the High Resolution Limited Area Model (HIRLAM) limited area numerical weather prediction (NWP) system. A data set of circa 133 000 radial wind measurements is passively monitored against the HIRLAM six-hourly background values in a 1-month experiment. No data assimilation experiments are performed at this stage. A new finding is that the improved modelling reduces the mean observation minus background (OmB) vector wind difference at ranges below 55 km, and the standard deviation of the radial wind OmB difference at ranges over 25 km.
In conclusion, a more accurate and still computationally feasible observation operator is developed. The companion paper (Part II) considers optimal super-observation processing of Doppler radar radial winds for HIRLAM, with general applicability in NWP.     

The Monagas Fold-Thrust Belt of Eastern Venezuela. Part I: Structural and thermal modeling

The Eastern Venezuelan Basin (EVB) contains one of the largest hydrocarbon accumulations in the world. Main petroleum targets are buried structures of the Monagas Fold-Thrust Belt (MFTB) which forms the northeastern edge of the EVB. The objective of this study is to integrate the seismic and well data that has been acquired over the last 10 years across the MFTB and EVB, to create an updated structural model. Three regional cross sections 60-75 km long are presented across an area of 4000 km².Five structural domains are described: Amarilis, Furrial, Jusepín, Cotoperí and Pirital. They are characterized by thrusts and high-angle reverse faults. Structural style changes along strike are related to variations in depth of detachment levels and to the strike-slip component of the deformation. We have estimated a shortening between 43 and 59 km that increases eastward over a distance of 40 km.The evolution of the MFTB is divided in four episodes based on stratigraphic, structural and thermal maturity evidences: (1) Oligocene-early Miocene initial movement of Pirital thrust. (2) Early Miocene simultaneous movement on Pirital, Furrial and Cotoperí thrusts. (3) Middle Miocene increases in velocity and change in geometry of Pirital thrust, during an out of sequence period of thrusting. (4) Late Miocene to Holocene minor thrust activity. This evolution is consistent with the oblique convergence between the Caribbean and South American plates and the convergence between North and South America that affected Eastern Venezuela during the Cenozoic.By analyzing the along-strike variations in structural style, new exploratory opportunities have been identified. Under the Orocual and Santa Bárbara fields two untested duplex structures are proposed; they were developed during the middle Miocene. Other prospective hydrocarbon traps are associated to oblique transpressive faults that create anticline structures.