Geochemical characterisation and erosion stability of fine‐grained groyne field sediments of the Middle Elbe River

Because of the significant reduction of flow velocity relative to the main current there is an accumulation of sediment material in the centre of the groyne fields over time (years to decades) up to a maximum of 1.5 m. Based on the element‐ and compound‐specific depth functions of this fine‐grained material with high organic matter content the sedimentation history of the Elbe River sediments concerning the environmentally relevant nutrients and pollutants can be reconstructed. It could be shown that the fresh sediments are less polluted than the deeper older ones. During the extreme flood of the river Elbe in August 2002 the critical erosion shear stress in the groyne fields was high enough to remobilise a considerable portion of the sediment depot. The hazard to the environment is high due to the fact that a significant part of the remobilised polluted sediments was distributed over the recent floodplain which is normally used as pasture. Thus these groyne field sediments represent a temporary danger for the water quality of the river Elbe. But the hazard of the solids sedimented in the recent floodplain of the river can be viewed as long‐lasting.     

Large earthquake doublets and fault plane heterogeneity in the northern Solomon Islands subduction zone

In the Solomon Islands and New Britain subduction zones, the largest earthquakes commonly occur as pairs with small separation in time, space and magnitude. This doublet behavior has been attributed to a pattern of fault plane heterogeneity consisting of closely spaced asperities such that the failure of one asperity triggers slip in adjacent asperities. We analyzed body waves of the January 31, 1974,M _w =7.3, February 1, 1974,M _w =7.4, July 20, 1975 (1437)M _w =7.6 and July 20, 1975 (1945),M _w =7.3 doublet events using an iterative, multiple station inversion technique to determine the spatio-temporal distribution of seismic moment release associated with these events. Although the 1974 doublet has smaller body wave moments than the 1975 events, their source histories are more complicated, lasting over 40 seconds and consisting of several subevents located near the epicentral regions. The second 1975 event is well modeled by a simple point source initiating at a depth of 15 km and rupturing an approximate 20 km region about the epicenter. The source history of the first 1975 event reveals a westerly propagating rupture, extending about 50 km from its hypocenter at a depth of 25 km. The asperities of the 1975 events are of comparable size and do not overlap one another, consistent with the asperity triggering hypothesis. The relatively large source areas and small seismic moments of the 1974 doublet events indicate failure of weaker portions of the fault plane in their epicentral regions. Variations in the roughness of the bathymetry of the subducting plate, accompanying subduction of the Woodlark Rise, may be responsible for changes in the mechanical properties of the plate interface.To understand how variations in fault plane coupling and strength affect the interplate seismicity pattern, we relocated 85 underthrusting earthquakes in the northern Solomon Islands Are since 1964. Relatively few smaller magnitude underthrusting events overlap the Solomon Islands doublet asperity regions, where fault coupling and strength are inferred to be the greatest. However, these asperity regions have been the sites of several previous earthquakes withM _s 7.0. The source regions of the 1974 doublet events, which we infer to be mechanically weak, contain many smaller magnitude events but have not generated any otherM _s 7.0 earthquakes in the historic past. The central portion of the northern Solomon Islands Arc between the two largest doublet events in 1971 (studied in detail bySchwartz et al., 1989a) and 1975 contains the greatest number of smaller magnitude underthrusting earthquakes. The location of this small region sandwiched between two strongly coupled portions of the plate interface suggest that it may be the site of the next large northern Solomon Islands earthquake. However, this region has experienced no known earthquakes withM _s 7.0 and may represent a relatively aseismic portion of the subduction zone.     

Purines,pyrimidines and organic carbon in lake sediments — A comparison of sediments from lakes of varying degrees of eutrophication

Sediment cores from the Great Lakes of North America and the Bodensee (Lake Constance) were analysed for purine and pyrimidine contents and total organic carbon. Lakes subject to recent cultural eutrophication showed steep increases in both organic carbon and purine and pyrimidine concentrations in the recently deposited sediments. Analytical results for Lake Huron (oligotrophic) and for a highly eutrophic area of the Bodensee (Gnadensee) appeared remarkably similar, in that only gradual, linear increases were observed in the total organic carbon curves. The curves for total purines and pyrimidines were also less steep for these cores as compared to the others studied. A possible explanation is that both areas represent lakes with no serious, recent changes in productivity. Uracil was found to decrease in concentration faster than any other purine or pyrimidine in the first few centimeters of all cores, supporting previous suggestions of more rapid turnover of this pyrimidine in sediments.     

Crustal structure of the Levant Basin, eastern Mediterranean

A seismic refraction/wide-angle reflection experiment was undertaken in the Levant Basin, eastern Mediterranean. Two roughly east–west profiles extend from the continental shelf of Israel toward the Levant Basin. The northern profile crosses the Eratosthenes Seamount and the southern profile crosses several distinct magnetic anomalies. The marine operation used 16 ocean bottom seismometers deployed along the profiles with an air gun array and explosive charges as energy sources. The results of this study strongly suggest the existence of oceanic crust under portions of the Levant Basin and continental crust under the Eratosthenes Seamount. The seismic refraction data also indicate a large sedimentary sequence, 10–14 km thick, in the Levant Basin and below the Levant continental margin. Assuming the crust is of Cretaceous age, this gives a fairly high sedimentation rate. The sequence can be divided into several units. A prominent unit is the 4.2 km/s layer, which is probably composed of the Messinian evaporites. Overlying the evaporitic layer are layers composed of Plio–Pleistocene sediments, whose velocity is 2.0 km/s. The refraction profiles and gravity and magnetic models indicate that a transition from a two layer continental to a single-layer oceanic crust takes place along the Levant margin. The transition in the structure along the southern profile is located beyond the continental margin and it is quite gradual. The northern profile, north of the Carmel structure, presents a different structure. The continental crust is much thinner there and the transition in the crustal structure is more rapid. The crustal thinning begins under western Galilee and terminates at the continental slope. The results of the present study indicate that the Levant Basin is composed of distinct crustal units and that the Levant continental margin is divided into at least two provinces of different crustal structure.     

Automated calibration applied to watershed‐scale flow simulations

Complexity in simulating the hydrological response in large watersheds over long times has prompted a significant need for procedures for automatic calibration. Such a procedure is implemented in the basin‐scale hydrological model (BSHM), a physically based distributed parameter watershed model. BSHM simulates the most important basin‐scale hydrological processes, such as overland flow, groundwater flow and stream–aquifer interaction in watersheds. Here, the emphasis is on estimating the groundwater parameters with water levels in wells and groundwater baseflows selected as the calibration targets. The best set of parameters is selected from within plausible ranges of parameters by adjusting the values of hydraulic conductivity, storativity, groundwater recharge and stream bed permeability. The baseflow is determined from stream flow hydrographs by using an empirical scheme validated using a chemical approach to hydrograph separation. Field studies determined that the specific conductance for components of the composite hydrograph were sufficiently unique to make the chemical approach feasible. The method was applied to the Big Darby Creek Watershed, Ohio. The parameter set selected for the groundwater system provides a good fit with the estimated baseflow and observed water well data. Copyright © 1999 John Wiley & Sons, Ltd.     

Relationships between physical‐geochemical soil properties and erodibility of streambanks among different physiographic provinces of Tennessee,USA

Erosion of cohesive soils in fluvial environments is dependent on physical, geochemical and biological properties, which govern inter‐particle attraction forces and control detachment rates from stream beds and banks. Most erosion rate models are based on the excess shear stress equation where the soil erodibility coefficient (k_d) is multiplied by the difference between the boundary hydraulic shear stress (τ_b) and the soil critical shear stress (τ_c). Both k_d and τ_c are a function of soil properties and must be obtained through in situ field or laboratory testing. Many studies have generated predictive relationships for k_d and τ_c derived from various soil properties. These studies typically were conducted in watersheds within a single physiographic region with a common surficial geology and/or investigated a limited number of soil properties, particularly geochemical properties. With widely reported differences in relationships between τ_c and soil properties, this study investigated differences in predictive relationships for τ_c among different physiographic provinces in Tennessee, USA. Erodibility parameters were determined in the field using a mini‐jet test device. Among these provinces, statistically four unique clusters were identified from a dataset of 128 observations and these data clusters were used to develop predictive models for τ_c to identify dominant properties governing erosion. In these clusters, 16 significant physical and geochemical soil properties were identified for τ_c prediction. Among these soil properties, water content and passing #200 sieve (percentage soil less than 75 μm) were the dominant controlling parameters to predict τ_c in addition to clay percentage (< 2 μm), bulk density, and soil pore water chemistry. This study suggests that unique relationships exist for physiographic provinces that are likely due to soil physical‐geochemical processes associated with surficial geology that determine minerology of the cohesive soil. Copyright © 2017 John Wiley & Sons, Ltd.     

Slowness domain offset and traveltime approximations in layered vertical transversely isotropic media

Considering horizontally layered transversely isotropic media with vertical symmetry axis and all types of pure‐mode and converted waves we present a new wide‐angle series approximation for the kinematical characteristics of reflected waves: horizontal offset, intercept time, and total reflection traveltime as functions of horizontal slowness. The method is based on combining (gluing) both zero‐offset and (large) finite‐offset series coefficients. The horizontal slowness is bounded by the critical value, characterised by nearly horizontal propagation within the layer with the highest horizontal velocity. The suggested approximation uses five parameters to approximate the offset, six parameters to approximate the intercept time or the traveltime, and seven parameters to approximate any two or all three kinematical characteristics. Overall, the method is very accurate for pure‐mode compressional waves and shear waves polarised in the horizontal plane and for converted waves. The application of the method to pure‐mode shear waves polarised in the vertical plane is limited due to cusps and triplications. To demonstrate the high accuracy of the method, we consider a synthetic, multi‐layer model, and we plot the normalised errors with respect to numerical ray tracing.     

Slowness-domain kinematical characteristics for horizontally layered orthorhombic media. Part I: Critical slowness match

Kinematical characteristics of reflected waves in anisotropic elastic media play an important role in the seismic imaging workflow. Considering compressional and converted waves, we derive new, azimuthally dependent, slowness-domain approximations for the kinematical characteristics of reflected waves (radial and transverse offsets, intercept time and traveltime) for layered orthorhombic media with varying azimuth of the vertical symmetry planes. The proposed method can be considered an extension of the well-known ‘generalized moveout approximation’ in the slowness domain, from azimuthally isotropic to azimuthally anisotropic models. For each slowness azimuth, the approximations hold for a wide angle range, combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. We consider two cases for the wide-angle ray: a ‘critical slowness match’ and a ‘pre-critical slowness match’ studied in Parts I and II of this work, respectively. For the critical slowness match, the approximations are valid within the entire slowness range, up to the critical slowness. For the ‘pre-critical slowness match’, the approximations are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. The critical slowness match is particularly effective when the subsurface model includes a dominant high-velocity layer where, for nearly critical slowness values, the propagation in this layer is almost horizontal. Comparing the approximated kinematical characteristics with those computed by numerical ray tracing, we demonstrate high accuracy.     

Traveltime approximation in vertical transversely isotropic layered media

The well‐known asymptotic fractional four‐parameter traveltime approximation and the five‐parameter generalised traveltime approximation in stratified multi‐layer transversely isotropic elastic media with a vertical axis of symmetry have been widely used for pure‐mode and converted waves. The first three parameters of these traveltime expansions are zero‐offset traveltime, normal moveout velocity, and quartic coefficient, ensuring high accuracy of traveltimes at short offsets. The additional parameter within the four‐parameter approximation is an effective horizontal velocity accounting for large offsets, which is important to avoid traveltime divergence at large offsets. The two additional parameters in the above‐mentioned five‐parameter approximation ensure higher accuracy up to a given large finite offset with an exact match at this offset. In this paper, we propose two alternative five‐parameter traveltime approximations, which can be considered extensions of the four‐parameter approximation and an alternative to the five‐parameter approximation previously mentioned. The first three short‐offset parameters are the same as before, but the two additional long‐offset parameters are different and have specific physical meaning. One of them describes the propagation in the high‐velocity layer of the overburden (nearly horizontal propagation in the case of very large offsets), and the other characterises the intercept time corresponding to the critical slowness that includes contributions of the lower velocity layers only. Unlike the above‐mentioned approximations, both of the proposed traveltime approximations converge to the theoretical (asymptotic) linear traveltime at the limit case of very large (“infinite”) offsets. Their accuracy for moderate to very large offsets, for quasi‐compressional waves, converted waves, and shear waves polarised in the horizontal plane, is extremely high in cases where the overburden model contains at least one layer with a dominant higher velocity compared with the other layers. We consider the implementation of the proposed traveltime approximations in all classes of problems in which the above‐mentioned approximations are used, such as reflection and diffraction analysis and imaging.