Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Impedance of flexible suction caissons

The dynamic response of offshore wind turbines is affected by the properties of the foundation and the subsoil. The aim of this paper is to evaluate the dynamic soil–structure interaction of suction caissons for offshore wind turbines. The investigations include evaluation of the vertical and coupled sliding–rocking vibrations, influence of the foundation geometry and examination on the properties of the surrounding soil. The soil is simplified as a homogenous linear viscoelastic material and the dynamic stiffness of the suction caisson is expressed in terms of dimensionless frequency‐dependent coefficients corresponding to different degrees of freedom. The dynamic stiffness coefficients for the skirted foundation are evaluated using a three‐dimensional coupled boundary element/finite element model. Comparisons with known analytical and numerical solutions indicate that the static and dynamic behaviours of the foundation are predicted accurately using the applied model. The analysis has been carried out for different combinations of the skirt length, Poisson's ratio of the subsoil and the ratio of the soil stiffness to the skirt stiffness. Copyright © 2007 John Wiley & Sons, Ltd.     

Determination of gravity anomaly at sea level from inversion of satellite gravity gradiometric data

Gravity gradients can be used to determine the local gravity field of the Earth. This paper investigates downward continuation of all elements of the disturbing gravitational tensor at satellite level using the second-order partial derivatives of the extended Stokes formula in the local-north oriented frame to determine the gravity anomaly at sea level. It considers the inversion of each gradient separately as well as their joint inversion. Numerical studies show that the gradients T_zz, T_xx, T_yy and T_xz have similar capability of being continued downward to sea level in the presence of white noise, while the gradient T_yz is considerably worse than the others. The bias-corrected joint inversion process shows the possibility of recovering the gravity anomaly with 1 mGal accuracy. Variance component estimation is also tested to update the observation weights in the joint inversion.     

Use of biliary PAH metabolites as a biomarker of pollution in fish from the Baltic Sea

During field campaigns of the BEEP project (Biological Effects of Environmental Pollution in Marine Coastal Ecosystems) in 2001-2002, metabolites of polycyclic aromatic hydrocarbons (PAHs) were determined in bile samples from three fish species, flounder (Platichthys flesus), perch (Perca fluviatilis) and eelpout (Zoarces viviparus), from four separate areas in the Baltic Sea. Two determination methods were applied: fixed wavelength fluorescence (FF) for pyrene-type metabolites and high-performance liquid chromatography with fluorescence detection (HPLC). There was a good correlation between the FF method and 1-OH pyrene determined by HPLC. Normalisation of the FF data for absorbance at 380 nm or bile protein concentrations greatly increased variance in one third and decreased it in two thirds of the cases and resulted in a loss of significant differences (protein normalisation) between the sampling stations, but normalisation of the HPLC data had little effect on the results. The biliary PAH metabolite content was usually higher in males than in females. In perch and eelpout the biliary PAH contents were at similar levels, whereas in flounder the levels were lower. The sampling areas arranged in decreasing order of biliary PAH contents were: Wismar Bay > Gulf of Gdansk > Lithuanian coast > Kvadofj?rden (reference area). It is concluded that FF with un-normalised data is a reliable and simple method for monitoring purposes and only one sex of a selected species should be used.     

The termination of the last major phase of aeolian sand movement,coastal dunefields,Denmark

Optically stimulated luminescence (OSL) dating of sand samples from stabilized (or inactive) coastal dunes in Denmark provides information on the age of the termination phase of the last major aeolian activity period. A total of 26 sand samples were taken from four different coastal dunefields around the North Sea, Skagerrak and Kattegat coasts of Denmark. The OSL dates indicate that the last major phase of aeolian activity terminated between ad 1860 and 1905. Most of the dunes examined in this study were active around 1820, during a period documented to have been very stormy. A dune management scheme started around 1792, and this no doubt was a major cause of dunefield stabilization, but an overall decline of storminess, particularly spring and summer storminess, around the end of the 19th century must also have contributed to the increasing inactivity of coastal dunes. The new OSL dates on aeolian sand movement agree well with historical data and data from topographic maps on dune movement. This agreement supports the observation from earlier work that OSL dating of recent aeolian sand movement is accurate over the last few decades to centuries. Copyright © 2005 John Wiley & Sons, Ltd.     

Evaluation of Combined Direct-Push Methods Used for Aquifer Model Generation

Most established methods to characterize aquifer structure and hydraulic conductivities of hydrostratigraphical units are not capable of delivering sufficient information in the spatial resolution that is desired for sophisticated numerical contaminant transport modeling and adapted remediation design. With hydraulic investigation methods based on the direct-push (DP) technology such as DP slug tests, DP injection logging, and the hydraulic profiling tool, it is possible to rapidly delineate hydrogeological structures and estimate their hydraulic conductivity in shallow unconsolidated aquifers without the need for wells. A combined application of these tools was used for the investigation of a contaminated German refinery site and for the setup of hydraulic aquifer models. The quality of DP investigation and the models was evaluated by comparisons of tracer transport simulations using these models and measured breakthroughs of two natural gradient tracer tests. Model scenarios considering the information of all tools together showed good reproduction of the measured breakthroughs, indicating the suitability of the approach and a minor impact of potential technical limitations. Using the DP slug tests alone yielded significantly higher deviations for the determined hydraulic conductivities compared to considering two or three of the tools. Realistic aquifer models developed on basis of such combined DP investigation approaches can help optimize remediation concepts or identify flow regimes for aquifers with a complex structure.     

Wave-induced steady streaming,mass transport and net sediment transport in rough turbulent ocean bottom boundary layers

The interaction between two important mechanisms which causes streaming has been investigated by numerical simulations of the seabed boundary layer beneath both sinusoidal waves and Stokes second order waves, as well as horizontally uniform bottom boundary layers with asymmetric forcing. These two mechanisms are streaming caused by turbulence asymmetry in successive wave half-cycles (beneath asymmetric forcing), and streaming caused by the presence of a vertical wave velocity within the seabed boundary layer as earlier explained by Longuet-Higgins. The effect of wave asymmetry, wave length to water depth ratio, and bottom roughness have been investigated for realistic physical situations. The streaming induced sediment dynamics near the ocean bottom has been investigated; both the resulting suspended load and bedload are presented. Finally, the mass transport (wave-averaged Lagrangian velocity) has been studied for a range of wave conditions. The streaming velocities beneath sinusoidal waves (Longuet-Higgins streaming) is always in the direction of wave propagation, while the streaming velocities in horizontally uniform boundary layers with asymmetric forcing are always negative. Thus the effect of asymmetry in second order Stokes waves is either to reduce the streaming velocity in the direction of wave propagation, or, for long waves relative to the water depth, to induce a streaming velocity against the direction of wave propagation. It appears that the Longuet-Higgins streaming decreases as the wave length increases for a given water depth, and the effect of wave asymmetry can dominate, leading to a steady streaming against the wave propagation. Furthermore, the asymmetry of second order Stokes waves reduces the mass transport (wave-averaged Lagrangian velocity) as compared with sinusoidal waves. The boundary layer streaming leads to a wave-averaged transport of suspended sediments and bedload in the direction of wave propagation.