TTG-type plutonic rocks formed in a modern arc batholith by hydrous fractionation in the lower arc crust

We present the geochemistry and intrusion pressures of granitoids from the Kohistan batholith, which represents, together with the intruded volcanic and sedimentary units, the middle and upper arc crust of the Kohistan paleo-island arc. Based on Al-in-hornblende barometry, the batholith records intrusion pressures from ~0.2 GPa in the north (where the volcano-sedimentary cover is intruded) to max. ~0.9 GPa in the southeast. The Al-in-hornblende barometry demonstrates that the Kohistan batholith represents a complete cross section across an arc batholith, reaching from the top at ~8–9 km depth (north) to its bottom at 25–35 km (south-central to southeast). Despite the complete outcropping and accessibility of the entire batholith, there is no observable compositional stratification across the batholith. The geochemical characteristics of the granitoids define three groups. Group 1 is characterized by strongly enriched incompatible elements and unfractionated middle rare earth elements (MREE)/heavy rare earth element patterns (HREE); Group 2 has enriched incompatible element concentrations similar to Group 1 but strongly fractionated MREE/HREE. Group 3 is characterized by only a limited incompatible element enrichment and unfractionated MREE/HREE. The origin of the different groups can be modeled through a relatively hydrous (Group 1 and 2) and of a less hydrous (Group 3) fractional crystallization line from a primitive basaltic parent at different pressures. Appropriate mafic/ultramafic cumulates that explain the chemical characteristics of each group are preserved at the base of the arc. The Kohistan batholith strengthens the conclusion that hydrous fractionation is the most important mechanism to form volumetrically significant amounts of granitoids in arcs. The Kohistan Group 2 granitoids have essentially identical trace element characteristics as Archean tonalite–trondhjemite–granodiorite (TTG) suites. Based on these observations, it is most likely that similar to the Group 2 rocks in the Kohistan arc, TTG gneisses were to a large part formed by hydrous high-pressure differentiation of primitive arc magmas in subduction zones.     

Damage phenomenon and petrophysical properties of sandstones at the Phnom Bakheng Temple (Angkor,Cambodia): first investigations and possible conservation measures

Sandstones, clay in the form of bricks and laterite are the building materials used by the Khmer to construct the imposing and magnificent temples in Southeast Asia. Many of these monuments suffer from fracturing, sanding, contour scaling, crust formation and salt weathering. The affinity to weathering is closely connected to the type of material. Two sandstone types classified as feldspathic arenite and quartz arenite of Angkor as well as two arkosic sandstones from Thailand are described and investigated in this study. Important petrophysical properties determined for the different sandstones consist of hydric expansion, thermal expansion, pore radii distribution and ultrasonic velocity. Different investigations such as capillary water uptake, surface hardness, hygroscopic water sorption, and salt resistance tests were undertaken in the laboratory to characterize the various rock types. Observations and quantified damage mapping were done onsite at the Phnom Bakheng Temple. Contour scaling in the form of weathering crusts is one of the main deterioration features observable at the Angkor monuments. Comparisons are made between the building stone, the crust material from the Phnom Bakheng Temple and fresh stone material used for restoration. Significant differences in hydric and especially in thermal expansion of the crust and sandstone have been determined. The results seem to indicate that extensional processes occur, which can be considered a force for detachment (i.e., contour scaling, flaking). In an experimental trial, the hydric and thermal expansion of the weathering crust and the building stone was significantly reduced by using a weak acid for the crust and a swelling inhibitor for the original building stone.     

The Namuskluft and Dreigratberg sections in southern Namibia (Kalahari Craton,Gariep Belt): a geological history of Neoproterozoic rifting and recycling of cratonic crust during the dispersal of Rodinia until the amalgamation of Gondwana

This paper presents combined U/Pb, Th/U and Hf isotope analyses on detrital and magmatic zircon grains together with whole-rock geochemical analyses of two basement and eight sedimentary rock samples from the Namuskluft and the Dreigratberg in southern Namibia (Gariep Belt). The sedimentary sections evolved during the Cryogenian on the SW part of the Kalahari Craton and where therefore deposited in an active rift setting during the break-up of Rodinia. Due to insufficient palaeomagnetic data, the position of the Kalahari Craton within Rodinia is still under discussion. There are possibilities to locate Kalahari along the western side of Australia/Mawsonland (Pisarevski et al. in Proterozoic East Gondwana: supercontinent assembly and break-up, Geological Society, London, 2003; Evans in Ancient Orogens and modern analogues. Geological Society, London, 2009; and others) or together with the Congo-Sao Francisco and Rio de la Plata Cratons (Li et al. in Prec Res 45: 203–2014, 2008; Frimmel et al. in Int J Earth Sci (Geol Rundsch) 100: 323–354, 2011; and others). It is sill unclear which craton rifted away from the Kalahari Craton during the Cryogenian. Although Middle to Upper Cryogenian magmatic activity is known for the SE Kalahari Craton (our working area) (Richtersveld Suite, Rosh Pinah Fm), all the presented samples show no U/Pb zircon ages younger than ca. 1.0 Ga and non-older than 2.06 Ga. The obtained U/Pb ages fit very well to the exposed basement of the Kalahari Craton (1.0–1.4 Ga Namaqua Province, 1.7–2.0 Ga Vioolsdrif Granite Suite and Orange River Group) and allow no correlation with a foreign craton such as the Rio de la Plata or Australia/Mawsonland. Lu–Hf isotopic signatures of detrital zircon point to the recycling of mainly Palaeoproterozoic and to a smaller amount of Archean crust in the source areas. εHf(t) signatures range between ?24 and +14.8, which relate to T_DM model ages between 1.05 and 3.1 Ga. Only few detrital zircon grains derived from magmas generated from Mesoproterozoic crustal material show more juvenile εHf(t) signatures of +14, +8 to +4 with T_DM model ages of 1.05–1.6 Ga. During Neoproterozoic deposition, only old cratonic crust with an inherited continental arc signature was available in the source area clearly demonstrated by Hf isotope composition of detrital zircon and geochemical bulk analysis of sedimentary rocks. The granodiorites of the Palaeoproterozoic basement underlying Namuskluft section are ca. 1.9 Ga old and show εHf(t) signatures of ?3 to ?5.5 with T_DM model ages of 2.4–2.7 Ga. These basement rocks demonstrate the extreme uplift and deep erosion of the underlying Kalahari Craton at its western margin before general subsidence during Cryogenian and Ediacaran time. The sedimentary sequence of the two examined sections (Namuskluft and Dreigratberg) proposes the presence of a basin and an increasing subsidence at the SW part of the Kalahari Craton during the Cryogenian. Therefore, we propose the initial formation of an intra-cratonic sag basin during the Lower Cryogenian that evolved later to a rift basin at the cratonic margin due to increasing crustal tension and rifting together with the opening of the Adamastor Ocean. As the zircons of the sedimentary rocks filling this basin show neither rift-related U/Pb ages nor an exotic craton as a possible source area, the only plausible sedimentary transport direction providing the found U/Pb ages would be from the E or the SE, directly from the heart of the Kalahari Craton. Due to subsidence and ongoing sedimentation from E/SE directions, the rift-related magmatic rocks were simply covered by the input of old intra-cratonic material that explains the absence of Neoproterozoic zircon grains in our samples. The geochemical analyses show the erosion of a continental arc and related sedimentary rocks with an overall felsic provenance. The source area was a deeply eroded and incised magmatic arc that evolved on continental crust, without any evidence for a passive margin. All of this can be explained by the erosion of rocks related to the Namaqua Belt, which represents one of the two major peaks of zircon U–Pb ages in all analysed samples. Therefore, the Namaqua Belt was well exposed during the Cryogenian, available to erosion and apart from the also well-exposed Palaeoproterozoic basement of the Kalahari Craton one potential source area for the sedimentary rocks in the investigated areas.     

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.     

On the effect of topography on surface wave propagation in the ambient noise frequency range

Due to the increasing popularity of analyzing empirical Green’s functions obtained from ambient seismic noise, more and more regional tomographical studies based on short-period surface waves are published. Results could potentially be biased in mountainous regions where topography is not small compared to the wavelength and penetration depth of the considered waves. We investigate the effect of topography on the propagation of short-period Rayleigh waves empirically by means of synthetic data using a spectral element code and a 3-D model with real topography. We show that topography along a profile through the studied area can result in an underestimation of phase velocities of up to about 0.7% at the shortest investigated period (3 s). Contrary to the expectation that this bias results from the increased surface distance along topography, we find that this error can be estimated by local topographic contrasts in the vicinity of the receiver alone. We discuss and generalize our results by considering topographic profiles through other mountain ranges and find that southern Norway is a good proxy to assess the topography effect. Nevertheless, topographic bias on phase velocity measurements is in general not large enough to significantly affect recovered velocity variations in the ambient noise frequency range.     

Bedload transport measurements at the Erlenbach stream with geophones and automated basket samplers

In the Erlenbach stream, a pre‐alpine steep channel in Switzerland, sediment transport has been monitored for more than 25 years. Near the confluence with the main valley river, stream flow is monitored and sediment is collected in a retention basin with a capacity of about 2000 m³. The basin is surveyed at regular intervals and after large flood events. In addition, sediment transport has been continuously monitored with piezoelectric bedload impact and geophone sensors since 1986. In 2008–2009, the measuring system in the Erlenbach stream was enhanced by installing an automatic system to obtain bedload samples. Movable metal baskets are mounted on a rail at the downstream wall of the large check dam above the retention basin, and they can be moved automatically into the flow to take bedload transport samples. The wire mesh of the baskets has a spacing of 10 mm to sample all sediment particles coarser than this size (which is about the limiting grain size detected by the geophones). The upgraded measuring system permits to obtain bedload samples over short sampling periods and to measure the grain size distribution of the transported material and its variation over time and with discharge. The analysis of calibration relationships for the geophone measuring system confirms findings from very similar measurements which were performed until 1999 with piezoelectric bedload impact sensors; there is a linear relationship between impulse counts and bedload mass passing over the sensors. Findings from flume experiments are used to discuss the most important factors which affect the calibration of the geophone signal. The bedload transport rates as measured by the moving baskets are among the highest measured in natural streams, with values of the order of several kilograms per meter per second. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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.     

A strategy for automated analysis of passive microseismic data to image seismic anisotropy and fracture characteristics

Monitoring of induced seismicity is gaining importance in a broad range of industrial operations from hydrocarbon reservoirs to mining to geothermal fields. Such passive seismic monitoring mainly aims at identifying fractures, which is of special interest for safety and productivity reasons. By analysing shear‐wave splitting it is possible to determine the anisotropy of the rock, which may be caused by sedimentary layering and/or aligned fractures, which in turn offers insight into the state of stress in the reservoir. We present a workflow strategy for automatic and effective processing of passive microseismic data sets, which are ever increasing in size. The automation provides an objective quality control of the shear‐wave splitting measurements and is based on characteristic differences between the two independent eigenvalue and cross‐correlation splitting techniques. These differences are summarized in a quality index for each measurement, allowing identification of an appropriate quality threshold. Measurements above this threshold are considered to be of good quality and are used in further interpretation. We suggest an automated inversion scheme using rock physics theory to test for best correlation of the data with various combinations of fracture density, its strike and the background anisotropy. This fully automatic workflow is then tested on a synthetic and a real microseismic data set.