The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

Parameter sensitivity of three Kalman filter schemes for assimilation of water levels in shelf sea models

In applications of data assimilation algorithms, a number of poorly known assimilation parameters usually need to be specified. Hence, the documented success of data assimilation methodologies must rely on a moderate sensitivity to these parameters. This contribution presents a parameter sensitivity study of three well known Kalman filter approaches for the assimilation of water levels in a three dimensional hydrodynamic modelling system. The filters considered are the ensemble Kalman filter (EnKF), the reduced rank square root Kalman filter (RRSQRT) and the steady Kalman filter. A sensitivity analysis of key parameters in the schemes is undertaken for a setup in an idealised bay. The sensitivity of the resulting root mean square error (RMSE) is shown to be low to moderate. Hence the schemes are robust within an acceptable range and their application even with misspecified parameters is to be encouraged in this perspective. However, the predicted uncertainty of the assimilation results are sensitive to the parameters and hence must be applied with care. The sensitivity study further demonstrates the effectiveness of the steady Kalman filter in the given system as well as the great impact of assimilating even very few measurements.     

Trends in the solubility of iron,aluminium, manganese and phosphorus in aerosol collected over the Atlantic Ocean

The solubility of iron, aluminium, manganese and phosphorus has been determined in aerosol samples collected between 49°N and 52°S during three cruises conducted in the Atlantic Ocean as part of the European Union funded IRONAGES programme. Solubilities (defined at pH 4.7) determined for Fe and Al in samples of Saharan dust were significantly lower (medians 1.7% and 3.0%, respectively) than the solubilities of these metals in aerosols from other source regions (whole dataset medians 5.2% and 9.0%, respectively). Mn solubility also varied with aerosol source, but the median solubility of Mn in Saharan dust was very similar to the median for the dataset as a whole (55% and 56%, respectively). The observed solubility of aerosol P was ∼ 32%, with P solubility in Saharan aerosol perhaps as low as 10%. Laboratory studies have indicated that aerosol Fe solubility is enhanced by acid processing. No relationship could be found between Fe solubility and the concentrations of acid species (non-seasalt SO₄²⁻, NO₃⁻) nor the net acidity of the aerosol, so we are unable to confirm that this process is significant in the atmosphere. In terms of the supply of soluble Fe to oceanic ecosystems on a global scale, the observed higher solubility for Fe in non-Saharan aerosols is probably not significant because the Sahara is easily the dominant source of Fe to the Atlantic. On a smaller scale however, higher solubility for aerosol Fe may alter our understanding of Fe cycling in regions such as the remote Southern Ocean.     

Stress Sharing Behavior and Its Mechanism during Consolidation Process in Composite Ground Improved by Sand Compaction Piles with Low Replacement Area Ratio

In order to accurately design a sand compaction pile (SCP) with low replacement area ratio, it is important to understand the mechanical interaction between the sand pile and clay ground and its mechanism during consolidation process in composite ground. In this article, therefore, a series of numerical analyses on composite ground improved by SCP with low replacement area ratio were carried out. The applicability of numerical analyses, in which an elasto-viscoplastic consolidation finite element method was applied, were confirmed by comparing the results obtained from a series of laboratory model tests with the composite ground improved by SCP. Through the results of the numerical analyses, mechanical behavior of the sand pile and clay in composite ground during consolidation is elucidated, together with a stress sharing mechanism between sand pile and clay.     

Large aperture seismic imaging at a convergent margin: Techniques and results from the Costa Rica seismogenic zone

In March and April 1995 a cooperative German, Costa Rican, and United States research team recorded onshore-offshore seismic data sets along the Pacific margin of Costa Rica using the R/V Ewing. Off the Nicoya Peninsula we used a linear array of ocean bottom seismometers and hydrophones (OBS/H) with onshore seismometers extending across much of the isthmus. In the central area we deployed an OBS/H areal array consisting of 30 instruments over a 9 km by 35-km area and had land stations on the Nicoya Peninsula adjacent to this marine array and also extending northeast on the main Costa Rican landmass. Our goal in these experiments was to determine the crustal velocity structure along different portions of this convergent margin and to use the dense instrument deployments to create migrated reflection images of the plate boundary zone and the subducting Cocos Plate. Our specific goal in the central area was to determine whether a subducted seamount is present at the location of the 1990, M 7 earthquake off the Nicoya Peninsula and can thus be linked to its nucleation. Subsequently we have processed the data to improve reflection signals, used the data to calculate crustal velocity models, and developed several wide-aperture migration techniques, based on a Kirchhoff algorithm, to produce reflection images. Along the northern transect we used the ocean bottom data to construct a detailed crustal velocity model, but reflections from the plate boundary and top and bottom of the subducting Cocos plate are difficult to identify and have so far produced poor images. In contrast, the land stations along this same transect recorded clear reflections from the top of the subducting plate or plate boundary, within the seismogenic zone, and we have constructed a clear image from this reflector beneath the Nicoya shelf. Data from the 3-D seismic experiment suffer from high-amplitude, coherent noise (arrivals other than reflections), and we have tried many techniques to enhance the signal to noise ratio of reflected arrivals. Due to the noise, an apparent lack of strong reflections from the plate boundary zone, and probable structural complexity, the resulting 3-D images only poorly resolve the top of the subducting Cocos Plate. The images are not able to provide compelling evidence of whether there is a subducting seamount at the 1990 earthquake hypocenter. Our results do show that OBS surveys are capable of creating images of the plate boundary zone and the subducting plate well into the seismogenic zone if coherent reflections are recorded at 1.8 km instrument spacing (2-D) and 5 km inline by 1 km crossline spacing for 3-D acquisition. However, due to typical high amplitude coherent noise, imaging results may be poorer than expected, especially in unfavorable geologic settings such as our 3-D survey area. More effective noise reduction in acquisition, possibly with the use of vertical hydrophone arrays, and in processing, with advanced multiple removal and possibly depth filtering, is required to achieve the desired detailed images of the seismogenic plate boundary zone.     

A hydrodynamic model of a two-part underwater towed system

A three-dimensional model of a two-part underwater towed system is studied. In the model, the governing equations of cables are established based on the Ablow and Schechter method. The boundary conditions for the two-part underwater towed system are derived. The six-degrees-of-freedom equations of motion for submarine simulations are adopted to predict the hydrodynamic performance of a towed vehicle. The established governing equations for the system are then solved using a central finite difference method. In this paper several algorithms are used to solve this special form of finite difference equations. The results in this paper indicate that the two-part underwater towed system improves the dynamic behavior of the towed vehicle and is an easy way to decouple the towing ship motion from the towed vehicle. Because the model uses an implicit time integration, it is stable for large time steps and is an effective algorithm for simulation of a large-scale underwater towed system.