Non-linear channel–shoal dynamics in long tidal embayments

The dynamics of finite-amplitude bed forms in a tidal channel is studied with the use of an idealized morphodynamic model. The latter is based on depth-averaged equations for the tidal flow over a sandy bottom. The model considers phenomena on spatial scales of the order of the tidal excursion length. Transport of sediment mainly takes place as suspended load. The reference state of this model is characterized by a spatially uniform M₂ tidal current over a fixed horizontal bed. The temporal evolution of deviations from this reference state is governed by amplitude equations: these are a set of non-linear equations that describe the temporal evolution of bed forms. These equations are used to obtain new morphodynamic equilibria which may be either static or time-periodic. Several of these bottom profiles show strong similarity with the tidal bars that are observed in natural estuaries. The dependence of the equilibrium solutions on the value of bottom friction and channel width is investigated systematically. For narrow channels (width small compared to the tidal excursion length) stable static equilibria exist if bottom friction is slightly larger than r_cr. For channel widths more comparable to the tidal excursion length, multiple stable steady states may exist for bottom friction parameter values below r_cr. Regardless of channel width, stable time-periodic equilibria seem to emerge as the bottom friction is increased.Responsible Editor: Jens Kappenberg     

Anthropogenic iodine-129 in seawater along a transect from the Norwegian coastal current to the North Pole

Variation in the concentrations of iodine-129 (129I, T1/2=15.7 Myr), a low-level radioactive component of nuclear fuel waste, is documented in surface waters and depth profiles collected during 2001 along a transect from the Norwegian Coastal Current to the North Pole. The surface waters near the Norwegian coast are found to have 20 times higher 129I concentration than the surface waters of the Arctic Ocean. The depth profiles of 129I taken in the Arctic Ocean reveal a sharp decline in the concentration to a depth of about 300-500 m followed by a weaker gradient extending down to the bottom. A twofold increase in the 129I concentration is observed in the upper 1000 m since 1996. Based on known estimates of marine transient time from the release sources (the nuclear reprocessing facilities at La Hague, France, and Sellafield, UK), a doubling in the 129I inventory of the top 1000 m of the Arctic Ocean is expected to occur between the years 2001 and 2006. As 129I of polar mixed layer and Atlantic layer of the Arctic Ocean is ventilated by the East Greenland Current into the Nordic Seas and North Atlantic Ocean, further dispersal and increase of the isotope concentration in these regions will be encountered in the near future.     

Liquefaction mitigation in silty soils using composite stone columns and dynamic compaction

The objective of this study is to develop an analytical methodology to evaluate the effectiveness ofvibro stone column (S. C.) and dynamic compaction (D.C.) techniques supplemented with wick drains to densify and mitigate liquethctionin saturated sands and non-plastic silty soils. It includes the following: (i) develop numerical models to simulate and analyze soil densification during S.C. installation and D.C. process, and (ii) identify parameters controlling post-improvement soil density in both cases, and (iii) develop design guidelines for densification of silty soils using the above techniques. An analytical procedure was developed and used to simulate soil response during S.C. and D.C. installations, and the results were compared with available case history data. Important construction design parameters and soil properties that affect the effectiveness of these techniques, and construction design choices suitable for sands and non-plastic silty soils were identified. The methodology is expected to advance the use of S.C. and DC. in silty soils reducing the reliance on expensive field trials as a design tool. The ultimate outcome of this research will be design charts and design guidelines for using composite stone columns and composite dynamic compaction techniques in liquefaction mitigation of saturated silty soils.     

GOCE Gravity Field Processing Strategy

A processing strategy and the corresponding software architecture for the processing of GOCE (Gravity field and steady-state Ocean Circulation Explorer) observables is presented and described, with the major objective to compute a high-accuracy, high-resolution spherical harmonic model of the Earth's gravity field. The combination of two numerical solution strategies, i.e. the rigorous solution of the corresponding large normal equation systems applying parallel processing (on a PC cluster) as the core solver, and the fast semianalytic approach as a quick-look gravity field analysis (QL-GFA) tool, is proposed. Such a method fusion benefits from the advantages of the individual components: the rigorous inversion of the system providing also the full variance-covariance information, and the quickness enabling the consecutive production of intermediate gravity field solutions, for the purpose to analyse partial and incomplete data sets and to derive a diagnosis of the performance of the GOCE measurement system. The functionality and operability of the individual components are demonstrated in the framework of a closed loop simulation, which is based on a realistic mission scenario both in terms of the orbit configuration and the coloured measuring noise. Special concern is given to the accuracy of the recovered coefficients, the numerical behaviour, the required computing time, and the particular role of the individual modules within the processing chain. In the case of the core solver, it is demonstrated that the assembling and rigorous solution of large normal equation systems can be handled by using Beowulf clusters within a reasonable computing time. The application of the quick-look tool to partial data sets with short-term data gaps is demonstrated on the basis of several case studies. Additionally, the spectral analysis of the residuals of the adjustment is presented as a valuable tool for the verification of the noise characteristics of the GOCE gradiometer.     

On the Incorporation of the Effect of Crustal Structure into Empirical Strong Ground Motion Estimation

This article has two purposes. Firstly, a validation exercise of the modal summation technique for the computation of synthetic strong-motion records is performed for two regions of Europe (Umbria-Marche and south Iceland), using a variety of region specific crustal structure models, by comparing the predicted ground motion amplitudes with observed motions. It is found that the rate of decay of ground motions is well predicted by the theoretical decay curves but that the absolute size of the ground motions is underpredicted by the synthetic time-histories. This is thought to be due to the presence of low-velocity surface layers that amplify the ground motions but are not included in the crustal structure models used to compute the synthetic time-histories. Secondly, a new distance metric based on the computed theoretical decay curves is introduced which should have the ability to model the complex decay of strong ground motions. The ability of this new distance metric to reduce the associated scatter in empirically derived equations for the estimation of strong ground motions is tested. It is found that it does not lead to a reduction in the scatter but this is thought to be due to the use of crustal structure models that are not accurate or detailed enough for the regions studied. This revised version was published online in July 2006 with corrections to the Cover Date.     

River runoff evaluation for ungauged watersheds by SWAP model. 2. Application of methods of physiographic similarity and spatial geostatistics

The efficiency of the methods of spatial proximity and geostatistics, as well as physico-geographic similarity, is studied as applied to the evaluation of the key model parameters of ungauged watersheds to be used in river runoff calculation by SWAP model. The target geographic objects were 323 experimental watersheds of MOPEX project. The quality of model parameter estimates and reproduction of river runoff hydrographs was analyzed in the case of the use of different similarity methods, and the order of decisions to be made was developed for the problem of river runoff calculation from an ungauged watershed for the entire area under study.     

Reproduction of secondary data in projection pursuit transformation

A challenge when working with multivariate data in a geostatistical context is that the data are rarely Gaussian. Multivariate distributions may include nonlinear features, clustering, long tails, functional boundaries, spikes, and heteroskedasticity. Multivariate transformations account for such features so that they are reproduced in geostatistical models. Projection pursuit as developed for high dimensional data exploration can also be used to transform a multivariate distribution into a multivariate Gaussian distribution with an identity covariance matrix. Its application within a geostatistical modeling context is called the projection pursuit multivariate transform (PPMT). An approach to incorporate exhaustive secondary variables in the PPMT is introduced. With this approach the PPMT can incorporate any number of secondary variables with any number of primary variables. A necessary alteration to the approach to make this numerically practical was the implementation of a continuous probability estimator that relies on Bernstein polynomials for the transformation that takes place in the projections. Stopping criteria were updated to incorporate a bootstrap t test that compares data sampled from a multivariate Gaussian distribution with the data undergoing transformation.