Refine
Has Fulltext
- no (124) (remove)
Year of publication
- 2011 (124) (remove)
Document Type
- Article (111)
- Doctoral Thesis (10)
- Review (2)
- Preprint (1)
Is part of the Bibliography
- yes (124)
Keywords
- Site effects (4)
- Tibetan Plateau (4)
- Holocene (3)
- Wave propagation (3)
- Benzene (2)
- Climate change (2)
- Crustal structure (2)
- Garnet (2)
- Geochronology (2)
- Land use (2)
Institute
- Institut für Geowissenschaften (124) (remove)
The Karaburun Peninsula, which is considered part of the Anatolide-Tauride Block of Turkey, contains clastic and carbonate sequences deposited on the northern margin of Gondwana. The Palaeozoic clastic sequence, which is intruded by the Early Triassic granitoid and tectonically overlies a Mesozoic melange sequence, can be divided into three subunits: a lower clastic subunit consisting of a sandstone-shale alternation, an upper clastic subunit consisting of black chert-bearing shales, sandstone and conglomerate, and a Permo-Carboniferous carbonate subunit. The lower Triassic Karaburun I-type granitoid has a high initial Sr-87/Sr-86 ratio (0.709021-0.709168), and low Nd-143/Nd-144 ratio (0.512004-0.512023) and epsilon Nd (-5.34 to -5.70) isotopic values. Geochronological data indicate a crystallization (intrusion) age of 247.1 +/- 2.0 Ma (Scythian). Geochemically, the acidic magmatism reflects a subduction-related continental-arc basin tectonic setting, which can be linked to the opening of the northern branch of Neo-Tethys as a continental back-arc rifting basin on the northern margin of Gondwana. This can be related to the closure through southward subduction of the Palaeotethys Ocean beneath Gondwana.
A parameterization for the motion of ice-shelf fronts on a Cartesian grid in finite-difference land-ice models is presented. The scheme prevents artificial thinning of the ice shelf at its edge, which occurs due to the finite resolution of the model. The intuitive numerical implementation diminishes numerical dispersion at the ice front and enables the application of physical boundary conditions to improve the calculation of stress and velocity fields throughout the ice-sheet-shelf system. Numerical properties of this subgrid modification are assessed in the Potsdam Parallel Ice Sheet Model (PISM-PIK) for different geometries in one and two horizontal dimensions and are verified against an analytical solution in a flow-line setup.
Aiming at the stimulation of intrinsic microbial activity, pulses of pure oxygen or pressurized air were recurrently injected into groundwater polluted with chlorobenzene. To achieve well-controlled conditions and intensive sampling, a large, vertical underground tank was filled with the local unconfined sandy aquifer material. In the course of two individual gas injections, one using pure oxygen and one using pressurized air, the mass transfer of individual gas species between trapped gas phase and groundwater was studied. Field data on the dissolved gas composition in the groundwater were combined with a kinetic model on gas dissolution and transport in porous media. Phase mass transfer of individual gas components caused a temporary enrichment of nitrogen, and to a lower degree of methane, in trapped gas leading to the formation of excess dissolved nitrogen levels downgradient from the dissolving gas phase. By applying a novel gas sampling method for dissolved gases in groundwater it was shown that dissolved nitrogen can be used as a partitioning tracer to indicate complete gas dissolution in porous media.
A poorly understood lag time of 15-20 m.y. exists between the initial Arabia-Eurasia continental collision in late Eocene to early Oligocene time and the acceleration of tectonic and sedimentary processes across the collision zone in the early to late Miocene. The late Eocene to Miocene-Pliocene clastic and shallow-marine sedimentary rocks of the Kond, Eyvanekey, and Semnan Basins in the Alborz Mountains (northern Iran) offer the possibility to track the evolution of this orogen in the framework of collision processes. A transition from volcaniclastic submarine deposits to shallow-marine evaporites and terrestrial sediments occurred shortly after 36 Ma in association with reversals in sediment provenance, strata tilting, and erosional unroofing. These events followed the termination of subduction arc magmatism and marked a changeover from an extensional to a contractional regime in response to initiation of continental collision with the subduction of stretched Arabian lithosphere. This early stage of collision produced topographic relief associated with shallow foreland basins, suggesting that shortening and tectonic loading occurred at low rates. Starting from the early Miocene (17.5 Ma), flexural subsidence in response to foreland basin initiation occurred. Fast sediment accumulation rates and erosional unroofing trends point to acceleration of shortening by the early Miocene. We suggest that the lag time between the initiation of continental collision (36 Ma) and the acceleration of regional deformation (20-17.5 Ma) reflects a two-stage collision process, involving the "soft" collision of stretched lithosphere at first and "hard" collision following the arrival of unstretched Arabian continental litho sphere in the subduction zone.
End member mixing analysis (EMMA) is a commonly applied method to identify and quantify the dominant runoff producing sources of water. It employs tracers to determine the dimensionality of the hydrologic system. Many EMMA studies have been conducted using two to six tracers, with some of the main tracers being Ca, Na, Cl(-), water isotopes, and alkalinity. Few studies use larger tracer sets including minor trace elements such as Li, Rb, Sr, and Ba. None of the studies has addressed the question of the tracer set size and composition, despite the fact that these determine which and how many end members (EM) will be identified. We examine how tracer set size and composition affects the conceptual model that results from an EMMA. We developed an automatic procedure that conducts EMMA while iteratively changing tracer set size and composition. We used a set of 14 tracers and 9 EMs. The validity of the resulting conceptual models was investigated under the aspects of dimensionality, EM combinations, and contributions to stream water. From the 16,369 possibilities, 23 delivered plausible results. The resulting conceptual models are highly sensitive to the tracer set size and composition. The moderate reproducibility of EM contributions indicates a still missing EM. It also emphasizes that the major elements are not always the most useful tracers and that larger tracer sets have an enhanced capacity to avoid false conclusions about catchment functioning. The presented approach produces results that may not be apparent from the traditional approach and it is a first step to add the idea of statistical significance to the EMMA approach.
The seismicity pattern along the San Andreas fault near Parkfield and Cholame, California, varies distinctly over a length of only fifty kilometres. Within the brittle crust, the presence of frictionally weak minerals, fault-weakening high fluid pressures and chemical weakening are considered possible causes of an anomalously weak fault northwest of Parkfield(1-4). Non-volcanic tremor from lower-crustal and upper-mantle depths(5-7) is most pronounced about thirty kilometres southeast of Parkfield and is thought to be associated with high pore-fluid pressures at depth(8). Here we present geophysical evidence of fluids migrating into the creeping section of the San Andreas fault that seem to originate in the region of the uppermost mantle that also stimulates tremor, and evidence that along-strike variations in tremor activity and amplitude are related to strength variations in the lower crust and upper mantle. Interconnected fluids can explain a deep zone of anomalously low electrical resistivity that has been imaged by magnetotelluric data southwest of the Parkfield-Cholame segment. Near Cholame, where fluids seem to be trapped below a high-resistivity cap, tremor concentrates adjacent to the inferred fluids within a mechanically strong zone of high resistivity. By contrast, sub-vertical zones of low resistivity breach the entire crust near the drill hole of the San Andreas Fault Observatory at Depth, northwest of Parkfield, and imply pathways for deep fluids into the eastern fault block, coincident with a mechanically weak crust and the lower tremor amplitudes in the lower crust. Fluid influx to the fault system is consistent with hypotheses of fault-weakening high fluid pressures in the brittle crust.
The aim of this paper is to estimate the Hurst parameter of Fractional Gaussian Noise (FGN) using Bayesian inference. We propose an estimation technique that takes into account the full correlation structure of this process. Instead of using the integrated time series and then applying an estimator for its Hurst exponent, we propose to use the noise signal directly. As an application we analyze the time series of the Nile River, where we find a posterior distribution which is compatible with previous findings. In addition, our technique provides natural error bars for the Hurst exponent.
Tectonic and metamorphic data for the Central Alps (Switzerland/Italy) are used to discuss this classic example of a Barrovian metamorphic terrain, notably the evolution of its thermal structure in space and time. Available P-T-t data indicate variable contributions of advective and conductive heat transport during collision and subsequent cooling and exhumation. Some areas experienced a prolonged period of partial melting while other areas, at the same time, show but moderate heating. The Barrow-type metamorphic field gradient observed in the final orogen is the result of two distinct tectonic processes, with their related advective and conductive heat transport processes. The two tectonic processes are (1) accretion of material within a subduction channel related to decompression and emplacement of high-pressure units in the middle crust and (2) wedging and related nappe formation in the continental lower plate. The second process postdates the first one. Wedging and underthrusting of continental lower plate material produces heat input into lower crustal levels, and this process is responsible for predominantly conductive heat transport in the overlying units. The interacting processes lead to different maximum temperatures at different times, producing the final Barrovian metamorphic field gradient. The south experienced rapid cooling, whereas the north shows moderate cooling rates. This discrepancy principally reflects differences in the temperature distribution in the deeper crust prior to cooling. Differences in the local thermal gradient that prevailed before the cooling also determined the relationships between cooling rate and exhumation rate in the different areas. Citation: Berger, A., S. M. Schmid, M. Engi, R. Bousquet, and M. Wiederkehr (2011), Mechanisms of mass and heat transport during Barrovian metamorphism: A discussion based on field evidence from the Central Alps (Switzerland/northern Italy), Tectonics, 30, TC1007, doi:10.1029/2009TC002622.
The Upper Cretaceous La Cova limestones (southern Pyrenees, Spain) host a rich and diverse larger foraminiferal fauna, which represents the first diversification of K-strategists after the mass extinction at the Cenomanian-Turonian boundary.
The stratigraphic distribution of the main taxa of larger foraminifera defines two assemblages. The first assemblage is characterised by the first appearance of lacazinids (Pseudolacazina loeblichi) and mean-dropsinids (Eofallotia simplex), by the large agglutinated Montsechiana montsechiensis, and by several species of complex rotalids (Rotorbinella campaniola, Iberorotalia reicheli, Orbitokhatina wondersmitti and Calcarinella schaubi). The second assemblage is defined by the appearance of Lacazina pyrenaica, Palandrosina taxyae and Martiguesia cyclamminiformis.
A late Coniacian-early Santonian age was so far accepted for the La Cova limestones, based on indirect correlation with deep-water fades bearing planktic foraminifers of the Dicarinella concavata zone. Strontium isotope stratigraphy, based on many samples of pristine biotic calcite of rudists and ostreids, indicates that the La Cova limestones span from the early Coniacian to the early-middle Santonian boundary. The first assemblage of larger foraminifera appears very close to the early-middle Coniacian boundary and reaches its full diversity by the middle Coniacian. The originations defining the second assemblage are dated as earliest Santonian: they represent important bioevents to define the Coniacian-Santonian boundary in the shallow-water facies of the South Pyrenean province.
By means of the calibration of strontium isotope stratigraphy to the Geological Time Scale, the larger foraminiferal assemblages of the La Cova limestones can be correlated to the standard biozonal scheme of ammonites, planktonic foraminifers and calcareous nannoplankton. This correlation is a first step toward a larger foraminifera standard biozonation for Upper Cretaceous carbonate platform facies.
Well-preserved primary contact relationships between a Late Proterozoic metasedimentary and the metagranitic core and Palaeozoic cover series of the Menderes Massif have been recognized in the eastern part of the Cine submassif on a regional-scale. Metaconglomerates occur as laterally discontinuous channel-fill bodies close the base of the metaquartzarenite directly above the basement. The pebbles in the metaconglomerates consist mainly of different types of tourmaline-rich leucocratic granitoids, tourmalinite and schist in a sandy matrix. Petrographic features, geochemical compositions and zircon radiometric ages (549.6 +/- 3.7-552.3 +/- 3.1 Ma) of the diagnostic clasts of the metaconglomerates (e.g. leucocratic granitoids and tourmalinites) show excellent agreement with their in situ equivalents (549.0 +/- 5.4 Ma) occurring in the Pan-African basement as stocks and veins.
The correlation between clasts in the metaconglomerates and granitoids of the basement suggests that the primary contact between the basement and cover series is a regional unconformity (supra-Pan-African Unconformity) representing deep erosion of the Pan-African basement followed by the deposition of the cover series. Hence the usage of 'core-cover' terminology in the Menderes Massif is valid. Consequently, these new data preclude the views that the granitic precursors of the leucocratic orthogneisses are Tertiary intrusions.
Recent studies have shown that rhizosphere hydraulic properties may differ from those of the bulk soil. Specifically, mucilage at the root-soil interface may increase the rhizosphere water holding capacity and hydraulic conductivity during drying. The goal of this study was to point out the implications of such altered rhizosphere hydraulic properties for soil-plant water relations. We addressed this problem through modeling based on a steady-rate approach. We calculated the water flow toward a single root assuming that the rhizosphere and bulk soil were two concentric cylinders having different hydraulic properties. Based on our previous experimental results, we assumed that the rhizosphere had higher water holding capacity and unsaturated conductivity than the bulk soil. The results showed that the water potential gradients in the rhizosphere were much smaller than in the bulk soil. The consequence is that the rhizosphere attenuated and delayed the drop in water potential in the vicinity of the root surface when the soil dried. This led to increased water availability to plants, as well as to higher effective conductivity under unsaturated conditions. The reasons were two: (i) thanks to the high unsaturated conductivity of the rhizosphere, the radius of water uptake was extended from the root to the rhizosphere surface; and (ii) thanks to the high soil water capacity of the rhizosphere, the water depletion in the bulk soil was compensated by water depletion in the rhizosphere. We conclude that under the assumed conditions, the rhizosphere works as an optimal hydraulic conductor and as a reservoir of water that can be taken up when water in the bulk soil becomes limiting.
The role of bedrock fractures and rock mass strength is often considered a primary influence on the efficiency of surface processes and the morphology of landscapes. Quantifying bedrock characteristics at hillslope scales, however, has proven difficult. Here, we present a new field-based method for quantifying the depth and apparent density of bedrock fractures within the shallow subsurface based on seismic refraction surveys. We examine variations in subsurface fracture patterns in both Fiordland and the Southern Alps of New Zealand to better constrain the influence of bedrock properties in governing rates and patterns of landslides, as well as the morphology of threshold landscapes. We argue that intense tectonic deformation produces uniform bedrock fracturing with depth, whereas geomorphic processes produce strong fracture gradients focused within the shallow subsurface. Additionally, we argue that hillslope strength and stability are functions of both the intact rock strength and the density of bedrock fractures, such that for a given intact rock strength, a threshold fracture-density exists that delineates between stable and unstable rock masses. In the Southern Alps, tectonic forces have pervasively fractured intrinsically weak rock to the verge of instability, such that the entire rock mass is susceptible to failure and landslides can potentially extend to great depths. Conversely, in Fiordland, tectonic fracturing of the strong intact rock has produced fracture densities less than the regional stability threshold. Therefore, bedrock failure in Fiordland generally occurs only after geomorphic fracturing has further reduced the rock mass strength. This dependence on geomorphic fracturing limits the depths of bedrock landslides to within this geomorphically weakened zone.
Continent-ocean-transition across a trans-tensional margin segment: off Bear Island, Barents Sea
(2011)
P>A 410 km long Ocean Bottom Seismometer profile spanning from the Bear Island, Barents Sea to oceanic crust formed along the Mohns Ridge has been modelled by use of ray-tracing with regard to observed P-waves. The northeastern part of the model represents typical continental crust, thinned from ca. 30 km thickness beneath the Bear Island to ca. 13 km within the Continent-Ocean-Transition. Between the Hornsund FZ and the Kn circle divide legga Fault, a 3-4 km thick sedimentary basin, dominantly of Permian/Carboniferous age, is modelled beneath the ca. 1.5 km thick layer of volcanics (Vestbakken Volcanic Province). The P-wave velocity in the 3-4 km thick lowermost continental crust is significantly higher than normal (ca. 7.5 km s-1). We interpret this layer as a mixture of mafic intrusions and continental crystalline blocks, dominantly related to the Paleocene-Early Eocene rifting event. The crystalline portion of the crust within the south-western part of the COT consists of a ca. 30 km wide and ca. 6 km thick high-velocity (7.3 km s-1) body. We interpret the body as a ridge of serpentinized peridotites. The magmatic portion of the ocean crust accreted along the Knipovich Ridge from continental break-up at ca. 35 Ma until ca. 20 Ma is 3-5 km thicker than normal. We interpret the increased magmatism as a passive response to the bending of this southernmost part of the Knipovich Ridge. The thickness of the magmatic portion of the crust formed along the Mohns Ridge at ca. 20 Ma decreases to ca. 3 km, which is normal for ultra slow spreading ridges.
We study changes in effective stress (normal stress minus pore pressure) that occurred in the French Alps during the 2003-2004 Ubaye earthquake swarm. Two complementary data sets are used. First, a set of 974 relocated events allows us to finely characterize the shape of the seismogenic area and the spatial migration of seismicity during the crisis. Relocations are performed by a double-difference algorithm. We compute differences in travel times at stations both from absolute picking times and from cross-correlation delays of multiplets. The resulting catalog reveals a swarm alignment along a single planar structure striking N130 degrees E and dipping 80 degrees W. This relocated activity displays migration properties consistent with a triggering by a diffusive fluid overpressure front. This observation argues in favor of a deep-seated fluid circulation responsible for a significant part of the seismic activity in Ubaye. Second, we analyze time series of earthquake detections at a single seismological station located just above the swarm. This time series forms a dense chronicle of +16,000 events. We use it to estimate the history of effective stress changes during this sequence. For this purpose we model the rate of events by a stochastic epidemic-type aftershock sequence model with a nonstationary background seismic rate lambda(0)(t). This background rate is estimated in discrete time windows. Window lengths are determined optimally according to a new change-point method on the basis of the interevent times distribution. We propose that background events are triggered directly by a transient fluid circulation at depth. Then, using rate-and-state constitutive friction laws, we estimate changes in effective stress for the observed rate of background events. We assume that changes in effective stress occurred under constant shear stressing rate conditions. We finally obtain a maximum change in effective stress close to -8 MPa, which corresponds to a maximum fluid overpressure of about 8 MPa under constant normal stress conditions. This estimate is in good agreement with values obtained from numerical modeling of fluid flow at depth, or with direct measurements reported from fluid injection experiments.
Vertical flow filters and vertical flow constructed wetlands are established wastewater treatment systems and have also been proposed for the treatment of contaminated groundwater. This study investigates the removal processes of volatile organic compounds in a pilot-scale vertical flow filter. The filter is intermittently irrigated with contaminated groundwater containing benzene, MTBE and ammonium as the main contaminants. The system is characterized by unsaturated conditions and high contaminant removal efficiency. The aim of the present study is to evaluate the contribution of biodegradation and volatilization to the overall removal of benzene and MTBE. Tracer tests and flow rate measurements showed a highly transient flow and heterogeneous transport regime. Radon-222, naturally occurring in the treated groundwater, was used as a gas tracer and indicated a high volatilization potential. Radon-222 behavior was reproduced by numerical simulations and extrapolated for benzene and MTBE, and indicated these compounds also have a high volatilization potential. In contrast, passive sampler measurements on top of the filter detected only low benzene and MTBE concentrations. Biodegradation potential was evaluated by the analysis of catabolic genes involved in organic compound degradation and a quantitative estimation of biodegradation was derived from stable isotope fractionation analysis. Results suggest that despite the high volatilization potential, biodegradation is the predominant mass removal process in the filter system, which indicates that the volatilized fraction of the contaminants is still subject to subsequent biodegradation. In particular, the upper filter layer located between the injection tubes and the surface of the system might also contribute to biodegradation, and might play a crucial role in avoiding the emission of volatilized contaminants into the atmosphere.
In the Himalaya of Chamba, NW India, a major orographic barrier in front of the Greater Himalayan Range extracts a high proportion of the monsoonal rainfall along its southern slopes and effectively shields the orogen interior from moisture-bearing winds. Along a similar to 100-km-long orogen perpendicular transect, 28 new apatite fission track (AFT) and 30 new zircon (U-Th)/He (ZHe) cooling ages reveal marked variations in age distributions and long-term exhumation rates between the humid frontal range and the semi-arid orogen interior. On the southern topographic front, very young, elevation-invariant AFT ages of <4 Ma have been obtained that are concentrated in a similar to 30-km-wide zone; 1-D-thermal modeling suggests a Plio-Pleistocene mean erosion rate of 0.8-1.9 mm yr(-1). In contrast, AFT and ZHe ages within the orogen interior are older (4-9 and 7-18 Ma, respectively), are positively correlated with sample elevation, and yield lower mean erosion rates (0.3-0.9 mm yr(-1)). Protracted low exhumation rates within the orogen interior over the last similar to 15 Myr prevailed contemporaneously with overall humid conditions and an effective erosional regime within the southern Himalaya. This suggests that the frontal Dhauladar Range was sufficiently high during this time to form an orographic barrier, focusing climatically enhanced erosional processes and tectonic deformation there. Thrusting along the two frontal range-bounding thrust, the Main Central Thrust and the Main Boundary Thrusts, was initiated at least similar to 15 Ma ago and has remained localized since then. The lack of evidence for localized uplift farther north indicates either a rather flat decollement with no ramp or the absence of active duplex systems beneath the interior of Chamba. Exhumational variability within Chamba is best explained as the result of continuous thrusting along a major basal decollement, with a flat beneath the slowly exhuming internal compartments and a steep frontal ramp at the rapidly exhuming frontal range. The pattern in Chamba contrasts with what is observed elsewhere along the Himalaya, where exhumation is focused in a zone similar to 150 km north of the orogenic front. In the NW Himalaya, preserved High Himalayan Crystalline nappes and Lesser Himalayan windows alternate on a relatively small scale of <100 km; these alternations are closely correlated with the pattern of exhumation. Although the spatial distribution of high-exhumation zones varies considerably between individual Himalayan sectors, all of these zones are closely correlated with locally higher rock-uplift rates, sharp topographic discontinuities, and focused orographic precipitation, suggesting strong feedbacks between tectonically driven rock uplift, orographically enhanced precipitation, and erosional processes.
P>Computing the magnitude of an earthquake requires correcting for the propagation effects from the source to the receivers. This is often accomplished by performing numerical simulations using a suitable Earth model. In this work, the energy magnitude M(e) is considered and its determination is performed using theoretical spectral amplitude decay functions over teleseismic distances based on the global Earth model AK135Q. Since the high frequency part (above the corner frequency) of the source spectrum has to be considered in computing M(e), the influence of propagation and site effects may not be negligible and they could bias the single station M(e) estimations. Therefore, in this study we assess the inter- and intrastation distributions of errors by considering the M(e) residuals computed for a large data set of earthquakes recorded at teleseismic distances by seismic stations deployed worldwide. To separate the inter- and intrastation contribution of errors, we apply a maximum likelihood approach to the M(e) residuals. We show that the interstation errors (describing a sort of site effect for a station) are within +/- 0.2 magnitude units for most stations and their spatial distribution reflects the expected lateral variation affecting the velocity and attenuation of the Earth's structure in the uppermost layers, not accounted for by the 1-D AK135Q model. The variance of the intrastation error distribution (describing the record-to-record component of variability) is larger than the interstation one (0.240 against 0.159), and the spatial distribution of the errors is not random but shows specific patterns depending on the source-to-station paths. The set of coefficients empirically determined may be used in the future to account for the heterogeneities of the real Earth not considered in the theoretical calculations of the spectral amplitude decay functions used to correct the recorded data for propagation effects.
The analysis of palaeoclimate time series is usually affected by severe methodological problems, resulting primarily from non-equidistant sampling and uncertain age models. As an alternative to existing methods of time series analysis, in this paper we argue that the statistical properties of recurrence networks - a recently developed approach - are promising candidates for characterising the system's nonlinear dynamics and quantifying structural changes in its reconstructed phase space as time evolves. In a first order approximation, the results of recurrence network analysis are invariant to changes in the age model and are not directly affected by non-equidistant sampling of the data. Specifically, we investigate the behaviour of recurrence network measures for both paradigmatic model systems with non-stationary parameters and four marine records of long-term palaeoclimate variations. We show that the obtained results are qualitatively robust under changes of the relevant parameters of our method, including detrending, size of the running window used for analysis, and embedding delay. We demonstrate that recurrence network analysis is able to detect relevant regime shifts in synthetic data as well as in problematic geoscientific time series. This suggests its application as a general exploratory tool of time series analysis complementing existing methods.