Optimal spherical spline filters for the analysis and comparison of regional-scale tomographic models

Advances in seismic tomography lead to increasingly detailed models of the Earth that are often represented on irregular and resolution-adaptive grids. To take full advantage of such models, their assessment must progress beyond a purely visual analysis, and tools must become available for their quantitative comparison.We present a method for the spectral analysis and comparison of multi-scale tomographic models. The method is applicable to irregular grids on the sphere, and is more efficient that filters based on spherical-harmonic expansions or convolution integrals. The combination of a spherical spline representation of tomographic information with Abel-Poisson scaling enables the construction of targetted spatial filters by solving a nonlinear inverse problem for appropriate weighting coefficients. This can be readily achieved with a simulated annealing approach for the limited number of weights. Once suitable filters have been generated they can be employed to address issues such as the patterns of small-scale heterogeneity, transitional structures and comparison of independent models from a region.We illustrate our method in a series of applications where we use different bandpass filters to detect differences in the distribution of small-scale heterogeneity beneath central and eastern Europe, and to compare several recent tomographic models of the Australian region.     

Magnetic signature of hydrocarbon-contaminated soils and sediments at the former oil field H?nigsen, Germany

Magnetic properties of hydrocarbon (HC) containing soils and sediments from two sites (Site A and B) of the former oil-field H?nigsen were analyzed in order to determine whether magnetic methods can be employed to delineate HC contamination of soils and sediments. Magnetic parameters such as magnetic susceptibility and induced isothermal remanent magnetizations, as well as soil and sediment properties such as pH, iron content and water content, HC content and most probable number counts of iron-metabolizing microorganisms were determined. The magnetic concentration-dependent parameters for HC contaminated samples were 25 times higher in soils from Site A than in sediment samples from Site B. However, at Site B the magnetic susceptibility was still four times higher in comparison to lithologically similar non-contaminated sediment samples from a third Site C. Newly formed magnetite containing mainly single domain particles was responsible for the magnetic enhancement, whereas superparamagnetic grains represented only a minor component. Site A had an acidic pH compared to neutral pH at Site B, and a higher crystalline and bioavailable total iron content. Nevertheless, Site B samples contained significant numbers of both iron(II)-oxidizing and iron(III)-reducing microorganisms indicating that microbial iron cycling might have taken place at this site and potentially played a role for iron mineral transformation, including magnetite (trans)formation. The content of total non-polar hydrocarbons (TNPH) at Site A was one order of magnitude higher than at Site B. Only at Site A magnetic susceptibility correlated well with TNPH. Our results demonstrate that HC contaminated samples had an enhanced magnetite content compared to non-contaminated soils and sediments. Therefore, magnetic methods may provide a quick and cost-effective way to assess HC contamination in soils and sediments. However, more field sites and laboratory investigations are needed to reveal the complex nature of the processes involved.     

Larval dispersal and vertical migration behaviour – a simulation study for short dispersal times

Current speed often varies with depth, so vertical movements of larvae are expected to have profound effects on dispersal velocity and therefore dispersal potential. Systematic behaviours are expected to have strong effects on dispersal. However, reliable information on the presence of vertical migrations in larvae is scarce, but the few well investigated empirical examples justify a detailed simulation study and an analysis of potential effects. We present a spatially explicit 3D hydrodynamic model that incorporates biological information in the form of active particles advected in a Lagrangian fashion. The set‐up is designed to analyze the sensitivity of dispersal distances to variation in vertical behaviour of larvae. We simulated short (4 days) pelagic larval durations (PLDs) to determine whether behaviour might be important over short dispersal periods. We found that sinusoidal behaviours (slow vertical migration) in or out of phase with tides did not significantly change the dispersal patterns compared to those of larvae that remained at the surface. By contrast, a quadratic pattern of behaviour resulting in rapid vertical migration, in or out of phase with tides, had dramatic effects on both distance and direction of dispersal. The resulting dispersal kernels were found to be multimodal due to the interaction between tidal and meteorological components in flow. Incorporating biological information on larval migrations in Lagrangian simulation of dispersal will be important in estimates of connectivity and forecasting marine reserve networks.     

An advanced signal processing technique for deriving grain size information of bedload transport from impact plate vibration measurements

A reliable characterization of bedload transport is required to gauge the engineering and theoretical issues related to the dynamics of sediments transport in rivers. However, while significant advances have been made in the development of monitoring techniques, robust quantitative predictive relationships have proven difficult to derive. In this article, we develop a dedicated signal processing technique aimed at improving the usage of impact plate measurements for material transport characterization. Our set‐up consists of a piezoelectric hydrophone mounted on the bottom side of a stainless steel plate, thus acting as a ‘sediment vibration sensor’. While the classical analysis with such systems is usually limited to rather simple procedures, such as impact counting, a large amount of useful information is contained in the actual waveform of the impact signal, which conveys the force and the contact time that the bedload imposes on the plate. An advanced signal processing technique called ‘first arrival atomic decomposition’ is used to improve the characterization of bedload transport by analysing the amplitude and frequency attributes of each single impact. This new processing approach proves to be well suited for bedload transport monitoring using plate systems and allows us to establish a relationship between the median grain size (D₅₀) and the impact signal properties. This link is first observed and validated with controlled flume experiments and then applied to continuous impact records in a small gravel‐bed river during a flood event. The estimated D₅₀ offers a novel possibility to observe the time‐varying grain size distribution of bedload transport. Copyright © 2014 John Wiley & Sons, Ltd.     

Complexity and non‐linearity in earth surface processes – concepts,methods and applications

Complexity has long been recognized and is increasingly becoming mainstream in geomorphology. However, the relative novelty of various concepts and techniques associated to it means that ambiguity continues to surround complexity. In this commentary, we present and discuss a variety of recent contributions that have the potential to help clarify issues and advance the use of complexity in geomorphology. Copyright © 2015 John Wiley & Sons, Ltd.     

Incorporating structural uncertainty of hydrological models in likelihood functions via an ensemble range approach

ABSTRACT

Model ensembles are possibly the most powerful tool to assess uncertainties in runoff predictions stemming from inadequacies in model structure. But in many applications little knowledge is gained about the specific weaknesses of the individual models. Here we introduce the ensemble range approach (ERA). Compared to other ensemble techniques, ERA is primarily intended to facilitate hydrological reasoning about model structural uncertainty. This is attempted by separate modelling of data uncertainty and structural uncertainty with two different error density functions that are combined in one likelihood function. The width of the structural error density is in accordance with the range of runoff predictions calculated by a small model ensemble at each individual time step. Albeit not the only choice, this study is restricted on the use of a modified beta density to represent structural uncertainty. The performance of ERA is assessed in some synthetic and real data case studies. Ensembles of two structurally identical models are applied, made possible by estimating the parameters of ERA and both models simultaneously.

Editor D. Koutsoyiannis; GUEST editor S. Weijs     

Hybrid testing and modeling of a suspended zipper steel frame

An analytical and experimental study has been conducted to evaluate the seismic performance of a three‐story suspended zipper steel frame. The frame was concentrically braced and had zipper struts to transfer the unbalanced forces induced on the beams due to the buckling of the lower‐story braces. The experimental study was conducted with the hybrid test technique, in which only the bottom‐story braces of the three‐story frame were physically tested, while the behavior of the rest of the frame was modeled using a general structural analysis software. The paper discusses issues pertinent to the calibration of the computer model for the analytical substructure as well as for the entire frame, including the selection of an appropriate damping matrix, and the modeling of the buckling behavior of the braces and bracing connections. The analytical model of the entire frame was validated with the hybrid tests and was able to accurately capture the material and geometric nonlinearities that developed when the braces yielded and buckled. This study has demonstrated the usefulness of hybrid testing in improving analytical models and modeling assumptions and providing information that cannot be obtained from an analytical study alone. The results have shown that the suspended zipper frame can distribute the brace nonlinearity over the first two stories as intended in the design and will not have catastrophic failure under the design‐level earthquakes considered in this study, despite the significant inelastic deformations. Copyright © 2009 John Wiley & Sons, Ltd.