Refine
Year of publication
Document Type
- Article (2810) (remove)
Keywords
- Holocene (35)
- climate change (30)
- Earthquake source observations (21)
- Climate change (18)
- Pollen (18)
- erosion (18)
- permafrost (18)
- Himalaya (17)
- Tibetan Plateau (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2810) (remove)
Lake sediments in 10 Ethiopian, Kenyan, and Tanzanian rift basins suggest that there were three humid periods at 2.7 to 2.5 million years ago (Ma), 1.9 to 1.7 Ma, and 1.1 to 0.9 Ma, superimposed on the longer-term aridification of East Africa. These humid periods correlate with increased aridity in northwest and northeast Africa and with substantial global climate transitions. These episodes could have had important impacts on the speciation and dispersal of mammals and hominins, because a number of key events, such as the origin of the genus Homo and the evolution of the species Homo erectus, took place in this region during that time.
Permian basins
(2005)
A localized dehydration zone, Sondrum stone quarry, Halmstad, SW Sweden, consists of a central, 1 m wide granitic pegmatoid dyke, on either side of which extends a 2.5-3 m wide dehydration zone (650-700 degrees C; 800 MPa; orthopyroxene-clinopyroxene-biotite-amphibole-garnet) overprinting a local migmatized granitic gneiss (amphibole-biotite- garnet). Whole-rock chemistry indicates that dehydration of the granitic gneiss was predominantly isochemical. Exceptions include [Y + heavy rare earth elements (HREE)], Ba, Sr, and F, which are markedly depleted throughout the dehydration zone. Systematic trends in the silicate and fluorapatite mineral chemistry across the dehydration zone include depletion in Fe, (Y + HREE), Na, K, F, and Cl, and enrichment in Mg, Mn, Ca, and Ti. Fluid inclusion chemistry is similar in all three zones and indicates the presence of a fluid containing CO2, NaCl, and H2O components. Water activities in the dehydration zone average 0.36, or XH2O = 0.25. All lines of evidence suggest that the formation of the dehydration zone was due to advective transport of a CO2-rich fluid with a minor NaCl brine component originating from a tectonic fracture. Fluid infiltration resulted in the localized partial breakdown of biotite and amphiboles to pyroxenes releasing Ti and Ca, which were partitioned into the remaining biotite and amphibole, as well as uniform depletion in (Y + HREE), Ba, Sr, Cl, and F. At some later stage, H2O-rich fluids (H2O activity > 0.8) gave rise to localized partial melting and the probable injection of a granitic melt into the tectonic fracture, which resulted in the biotite and amphibole recording a diffusion profile for F across the dehydration zone into the granitic gneiss as well as a diffusion profile in Fe, Mn, and Mg for all Fe-Mg silicate minerals within 100 cm of the pegmatoid dyke
A confocal set-up is presented that improves micro-XRF and XAFS experiment with high-pressure e diamond-anvil cells (DACs) In this experiment a probing volume is defined by the focus of the incoming synchrotron radiation beam and that of a polycapillary X-ray half-lens with a very long working distance, which is placed in front of the fluorescence detector This set-up enhances the quality of the fluorescence and XAFS spectra, and thus the sensitivity for detecting elements at low concentrations. It efficiently suppresses signal from outside the sample chamber, which stems from elastic and inelastic scattering of the incoming beam by the diamond anvils as well as from excitation of fluorescence from the body of the DAC
To balance the steady decrease of conventional hydrocarbon resources, increased utilization of unconventional and new energy resources, such as shale gas and geothermal energy, is required. Also, the geological sequestration of carbon dioxide is being considered as a technology that may temporarily mitigate the effects of CO2 emission. Sites suitable for shale gas production, geothermal exploration, or CO2 sequestration are commonly characterized by electrical resistivities distinctly different from those of the surrounding rocks. Therefore, electromagnetic methods can be viable tools to help identify target sites suitable for exploration, and to monitor reservoirs during energy production or CO2 injection. Among the wide variety of electromagnetic methods available, controlled-source magnetotelluric (CSMT) may be particularly suitable because of (i) its ability to resolve both electrically resistive and conductive structures, (ii) controlled sources offering noise control and thus facilitating surveys in populated regions, and (iii) the potential of penetration throughout the depth range accessible by drilling. Nevertheless, CSMT has not yet been widely employed because of logistical challenges of field operations and the requirement of complex and highly computer-intensive data processing. With these difficulties gradually being mitigated by recent technological developments, CSMT may now be reconsidered as an exploration tool. Here, we investigate by 1D and 3D numerical simulations the feasibility of detecting gas shales and identifying sites eligible for geothermal exploration or CO2 sequestration from CSMT data. We consider surface-to-surface, borehole-to-surface, and cross-hole configurations of the sources and receivers. Results and conclusions on the detectability of the targets of interest are presented for various exploration and monitoring scenarios, which are roughly representative of the geological setting of the North German Basin.
Characterization of polarization attributes of seismic waves using continuous wavelet transforms
(2006)
Complex-trace analysis is the method of choice for analyzing polarized data. Because particle motion can be represented by instantaneous attributes that show distinct features for waves of different polarization characteristics, it can be used to separate and characterize these waves. Traditional methods of complex-trace analysis only give the instantaneous attributes as a function of time or frequency. However. for transient wave types or seismic events that overlap in time, an estimate of the polarization parameters requires analysis of the time-frequency dependence of these attributes. We propose a method to map instantaneous polarization attributes of seismic signals in the wavelet domain and explicitly relate these attributes with the wavelet-transform coefficients of the analyzed signal. We compare our method with traditional complex-trace analysis using numerical examples. An advantage of our method is its possibility of performing the complete wave-mode separation/ filtering process in the wavelet domain and its ability to provide the frequency dependence of ellipticity, which contains important information on the subsurface structure. Furthermore, using 2-C synthetic and real seismic shot gathers, we show how to use the method to separate different wave types and identify zones of interfering wave modes
What is the most appropriate sampling scheme to estimate event-based average throughfall? A satisfactory answer to this seemingly simple question has yet to be found, a failure which we attribute to previous efforts' dependence on empirical studies. Here we try to answer this question by simulating stochastic throughfall fields based on parameters for statistical models of large monitoring data sets. We subsequently sampled these fields with different sampling designs and variable sample supports. We evaluated the performance of a particular sampling scheme with respect to the uncertainty of possible estimated means of throughfall volumes. Even for a relative error limit of 20%, an impractically large number of small, funnel-type collectors would be required to estimate mean throughfall, particularly for small events. While stratification of the target area is not superior to simple random sampling, cluster random sampling involves the risk of being less efficient. A larger sample support, e.g., the use of trough-type collectors, considerably reduces the necessary sample sizes and eliminates the sensitivity of the mean to outliers. Since the gain in time associated with the manual handling of troughs versus funnels depends on the local precipitation regime, the employment of automatically recording clusters of long troughs emerges as the most promising sampling scheme. Even so, a relative error of less than 5% appears out of reach for throughfall under heterogeneous canopies. We therefore suspect a considerable uncertainty of input parameters for interception models derived from measured throughfall, in particular, for those requiring data of small throughfall events.
In an attempt to map the shallow geometry of the Maleme Fault Zone (North Island, New Zealand) and estimate vertical displacements of selected fault strands, we have collected 2D and 3D georadar data using 100 MHz antennae. The 2D data consisted of three parallel georadar lines recorded perpendicular to the axis of the well-defined graben of the Maleme Fault Zone. These similar to 160 in long lines, which were 7.5 m apart, crossed several fault strands on either side of the graben axis. The processed georadar sections revealed two prominent parallel reflections that originated from the boundaries of Late Pleistocene lacustrine and tephra deposits. Distinct vertical offsets of these reflections allowed us to estimate displacernents at individual fault strands across the entire inner graben. The total displacements represented by these offsets was similar to 10-20% greater than that inferred from geomorphological studies, thus demonstrating the limitations of surface observations for determining cumulative fault movements. The 3D georadar data set, recorded across an area of similar to 70x similar to 20 in to one side of the graben axis, provided key details on individual fault strands. For the 3D visualization of fault-related structures, various spatial attribute analyses based on the cosine of the instantaneous phase proved to be useful
Metabasites were sampled from rock series of the subducted margin of the Indian Plate, the so-called Higher Himalayan Crystalline, in the Upper Kaghan Valley, Pakistan. These vary from corona dolerites, cropping out around Saif- ul-Muluk in the south, to coesite-eclogite close to the suture zone against rocks of the Kohistan arc in the north. Bulk rock major- and trace-element chemistry reveals essentially a single protolith as the source for five different eclogite types, which differ in fabric, modal mineralogy as well as in mineral chemistry. The study of newly-collected samples reveals coesite (confirmed by in situ Raman spectroscopy) in both garnet and omphacite. All eclogites show growth of amphiboles during exhumation. Within some coesite-bearing eclogites the presence of glaucophane cores to barroisite is noted whereas in most samples porphyroblastic sodic-calcic amphiboles are rimmed by more aluminous calcic amphibole (pargasite, tschermakite, and edenite). Eclogite facies rutile is replaced by ilmenite which itself is commonly surrounded by titanite. In addition, some eclogite bodies show leucocratic segregations containing phengite, quartz, zoisite and/or kyanite. The important implication is that the complex exhumation path shows stages of initial cooling during decompression (formation of glaucophane) followed by reheating: a very similar situation to that reported for the coesite-bearing eclogite series of the Tso Morari massif, India, 450 km to the south-east.
Simulation of spatial sensor characteristics in the context of the EnMAP Hyperspectral mission
(2010)
The simulation of remote sensing images is a valuable tool for defining future Earth observation systems, optimizing instrument parameters, and developing and validating data-processing algorithms. A scene simulator for optical Earth observation data has been developed within the Environmental Mapping and Analysis Program (EnMAP) hyperspectral mission. It produces EnMAP-like data following a sequential processing approach consisting of five independent modules referred to as reflectance, atmospheric, spatial, spectral, and radiometric modules. From a modeling viewpoint, the spatial module is the most complex. The spatial simulation process considers the satellite-target geometry, which is adapted to the EnMAP orbit and operating characteristics, the instrument spatial response, and the sources of spatial nonuniformity (keystone, telescope distortion and smile, and detector coregistration). The spatial module of the EnMAP scene simulator is presented in this paper. The EnMAP spatial and geometric characteristics will be described, the simulation methodology will be presented in detail, and the capability of the EnMAP simulator will be shown by illustrative examples.
Map form information on forest biomass is required for estimating bioenergy potentials and monitoring carbon stocks. In Finland, the growing stock of forests is monitored using multi-source forest inventory, where variables are estimated in the form of thematic maps and area statistics by combining information of field measurements, satellite images and other digital map data. In this study, we used the multi-source forest inventory methodology for estimating forest biomass characteristics. The biomass variables were estimated for national forest inventory field plots on the basis of measured tree variables. The plot-level biomass estimates were used as reference data for satellite image interpretation. The estimates produced by satellite image interpretation were tested by cross-validation. The results indicate that the method for producing biomass maps on the basis of biomass models and satellite image interpretation is operationally feasible. Furthermore, the accuracy of the estimates of biomass variables is similar or even higher than that of traditional growing stock volume estimates. The technique presented here can be applied, for example, in estimating biomass resources or in the inventory of greenhouse gases.
The deterministic calculation of earthquake scenarios using complete waveform modelling plays an increasingly important role in estimating shaking hazard in seismically active regions. Here we apply 3-D numerical modelling of seismic wave propagation to M 6+ earthquake scenarios in the area of the Lower Rhine Embayment, one of the seismically most active regions in central Europe. Using a 3-D basin model derived from geology, borehole information and seismic experiments, we aim at demonstrating the strong dependence of ground shaking on hypocentre location and basin structure. The simulations are carried out up to frequencies of ca. 1 Hz. As expected, the basin structure leads to strong lateral variations in peak ground motion, amplification and shaking duration. Depending on source-basin-receiver geometry, the effects correlate with basin depth and the slope of the basin flanks; yet, the basin also affects peak ground motion and estimated shaking hazard thereof outside the basin. Comparison with measured seismograms for one of the earthquakes shows that some of the main characteristics of the wave motion are reproduced. Cumulating the derived seismic intensities from the three modelled earthquake scenarios leads to a predominantly basin correlated intensity distribution for our study area
Two types of electrical conductivity sensors were evaluated to prospect circular ditches surrounding former Bronze Age burial mounds, complementing aerial photography. The first sensor was based on the electrical resistivity (ER) method, while the second sensor was based on frequency-domain electromagnetic induction (FDEM). Both sensors were designed with multiple receivers, which measure several depth sensitivities simultaneously. First, the sensors were tested on an experimental site where a rectangular structure with limited dimensions was dug in a sandy soil. The structure appeared as a higher conductivity anomaly in the low-conductivity sand. Then, both methods were applied on two Bronze Age sites with different soil properties, which were discovered by aerial photography. The first site, in a sandy soil, gave only very weak anomalies. Soil augering revealed that the ditch filling consisted of the same sandy material as the surrounding, therefore this filling was not able to cause a high-conductivity contrast. Due to its lower sensitivity to noise in the low-conductive range, the ER-sensor produced a more pronounced anomaly than the FDEM-sensor. The second site was located on top of a ridge with a shallow substrate of Tertiary, coastal sediments. The ditch was very clearly visible on the sensor maps as a conductive low. At this location, the soil augering revealed that the ditch was dug through an alternating clay-sand layer and subsequently filled up with silty material from the topsoil. Overall, the shallow receiver separation produced anomalies that were both stronger and that corresponded better to the geometry of the ditches. The other receiver separations provided more information on the natural soil layering, and in the case of the ER-array they could be used to obtain a cross-section of the actual electrical conductivity with 2-D inversion modelling. The results of this study proofed that conductivity sensors can detect Bronze Age ditches, with varying contrast depending on the soil geomorphology. Moreover, the sensor maps combined with soil observations by coring provided insight in the environmental conditions that influence the contrast of the anomalies seen on the aerial photographs and the sensor maps.
The strongly peraluminous and P-rich, protolithionite and zinnwaldite leucogranites from Podlesi, western Krusne Hory Mts., Czech Republic, contain accessory zircon with extraordinary enrichment of several elements, which constitute trace elements in common zircon. Elements showing a not yet reported anomalous enrichment include P (up to 20.2 wt.% P2O5; equivalent to 0.60 apfu, formula calculated on the basis of 4 oxygen atoms), Bi (up to 9.0 wt.% Bi2O3; 0.086 apfu), Nb (up to 6.7 wt.% Nb2O5, 0.12 apfu), Sc (up to 3.45 wt.% Sc2O3; 0.10 apfu), U (up to 14.8 wt.% UO2; 0.12 apfu) and F (up to 3.81 wt.% F; 0.42 apfu). Strong enrichment of P preferentially involved the berlinite-type substitution (2 Si4+ double left right arrow P5+ + Al3+) implying that significant Al may enter the Si position in zircon. Incorporation of other exotic elements is primarily governed by the xenotime (Si4++Zr4+ double left right arrow P5++Y3+), pretulite (Sc3++P5+ double left right arrow Zr4++Si4+), brabantite-type (Ca2++(U, Th)(4+)+2P(5+) double left right arrow 2Zr(4+)+2Si(4+)), and ximengite-type (Bi3++P5+double left right arrow Zr4++Si4+) substitution reactions. One part of the anomalous zircons formed late-magmatically, from a strongly peraluminous, P-F-U-rich hydrous residual melt that gave rise to the zinnwaldite granite. Interaction with aggressive residual fluids and metamictization have further aided in element enrichment or depletion, particularly in altered parts of zircon contained in the protolithionite granite. Most of the zircon from F-rich greisens have a composition close to endmember ZrSiO4 and are chemically distinct from zircon in its granite parent. This discrepancy implies that at Podlesi, granitic zircon became unstable and completely dissolved during greisenization. Part of the mobilized elements was reprecipitated in newly grown, hydrothermal zircon.
The correlation of deformation fabrics and metamorphic reactions with geochronologic data of UHP metamorphic rocks demonstrate that the multistage subduction and exhumation history of the Central Dabie Complex requires rapid subduction and rapid initial exhumation. Moreover, these data show that volume diffusion is not the major resetting mechanism of radiogenic isotope systems. Thus, our age data do not simply reflect a thermal/cooling history. In the investigated section, the maximum age for UHP is given by the 244 +/- 3 Ma (2 sigma) U-Pb age of a pre-UHP titanite phenocryst that survived UHP metamorphism and subsequent tectonometamorphic events. A minimum age for UHP is set by the 238 +/- 1 Ma (2 sigma) U-238-Pb-206 mineral isochron age of titanite and cogenetic epidote. These minerals formed from local partial melts during ascent and their age suggests fast exhumation and emplacement in the middle crust. In the period of ca. 238-218 Ma, the UHP terrain records HT metamorphism, local partial melting, and extensive pervasive strain below the eclogite (jd+grt) stability field. Exhumation was polyphase with a first phase of fast exhumation, succeeded by episodes of HT metamorphism and concomitant deformation at deep/mid crustal level between 238 and 218 Ma. Slow exhumation related to the final emplacement of tectonic units along greenschist facies shear zones did not cease before ca. 209-204 Ma. The resetting and homogenization of radiogenic isotope systems were aided by dissolution precipitation creep, which was the dominant deformation mechanism in quartz-feldspar rocks, in combination with fluid influx. (c) 2005 Elsevier B.V. All rights reserved
Intramontane sedimentary basins along the margin of continental plateaus often preserve strata that contain fundamental information regarding the pattern of orogenic growth. The sedimentary record of the elastic Miocene-Pliocene sequence deposited in the Fiambala Basin, at the southern margin of the Puna Plateau (NW Argentina), documents the late Miocene paleodrainage evolution from headwaters to the west, towards headwaters in the ranges that constitute the border of the Puna Plateau to the north. Apatite Fission track (AFT) thermochronology of sedimentary and basement rocks show that the southern Puna Plateau was the source for the youngest, middle Miocene, detrital population detected in late Miocene rocks; and that the margin of the Puna Plateau expressed a high relief, possibly similar to or higher than at present, by late Miocene time. Cooling ages obtained from basement rocks at the southern Puna margin suggest that exhumation started in the Oligocene and continued until the middle Miocene. We interpret the basin reorganization and the creation of a high relief plateau margin to be the direct response of the source-basin system to a wholesale surface uplift event that may have occurred during the late Cenozoic in the Puna-Altiplano region. At this time coeval paleodrainage reorganization is observed not only in the Fiambala Basin, but also in different basins along the southern and eastern Puna margin, suggesting a genetic link between the last stage of plateau formation and basin response. However, this event did not cause sufficient exhumation of basin bounding ranges to be recorded by AFT thermochronology. Our new data thus document a decoupling between late Cenozoic surface uplift and exhumation in the southern Puna Plateau. High relief achieved at the Puna margin by late Miocene time is linked to Oligocene-Miocene exhumation; no significant erosion (< 3 km) has occurred since in this and highland.
A comprehensive survey of the accessory-mineral assemblages in Variscan granites of the German Erzgebirge and Pan-African granites from Jordan revealed the occurrence of intermediate solid solutions of the tetragonal thorite- xenotime-zircon-coffinite mineral group with partially novel compositions. These solid solutions preferentially formed in evolved and metasomatically altered, P-poor leucogranites of either I- or A-type affinity. Thorite from the Erzgebirge contained up to 18-8 Wt-% Y2O3, 16.1 wt.% ZrO2, and 23.3 Wt-% UO2. Xenotime and zircon have incorporated Th in abundances up to 36.3 wt.% and 41.8 wt.% ThO2, respectively. Extended compositional gradation with only minor gaps is confined to hydrated members of this mineral group, and is observed to exist between thorite and xenotime, thorite and coffinite, and Y-HREE-bearing thorite and zircon. Complex, hydrous solid solutions containing elevated abundances of three or more of the endmembers are subordinate. Previously reported intermediate solid solutions between anhydrous zircon and xenotime, and anhydrous zircon and thorite, are not observed and are in conflict with experimental work demonstrating very limited miscibility between anhydrous species of endmember composition. The majority of hydrous intermediate solid solutions in the Th-Y-Zr-U system are likely thermodynamically unstable. Instead, they are probably metastable responses to unusual physico-chemical conditions involving various parameters and conditions, the relative importance of which is incompletely known. Leaching and dissolution of preexisting accessory phases during interaction with F-bearing hydrous fluids enriched in Th, Y(HREE), Zr, and/or U, and common deposition of the various elements at disequilibrium (supersaturation) seems to play a key role, but other processes may be of similar importance. Experimental work involving hydrous conditions and complex systems composed of more than two endmembers are needed to shed light into the stability relations of the chemically uncommon compositions treated in this study.
Metamorphic dome complexes occur within the internal structures of the northern Himalaya and southern Tibet. Their origin, deformation, and fault displacement patterns are poorly constrained. We report new field mapping, structural data, and cooling ages from the western flank of the Leo Pargil dome in the northwestern Himalaya in an attempt to characterize its post-middle Miocene structural development. The western flank of the dome is characterized by shallow, west-dipping pervasive foliation and WNW-ESE mineral lineation. Shear-sense indicators demonstrate that it is affected by east-west normal faulting that facilitated exhumation of high-grade metamorphic rocks in a contractional setting. Sustained top-to-northwest normal faulting during exhumation is observed in a progressive transition from ductile to brittle deformation. Garnet and kyanite indicate that the Leo Pargil dome was exhumed from the mid-crust. Ar- 40/Ar-39 mica and apatite fission track (AFT) ages constrain cooling and exhumation pathways front 350 to 60 degrees C and suggest that the dome cooled in three stages since the middle Miocene. Ar-40/Ar-39 white mica ages of 16-14 Ma suggest a first phase of rapid cooling and provide minimum estimates for the onset of dome exhumation. AFT ages between 10 and 8 Ma suggest that ductile fault displacement had ceased by then, and AFT track-length data from high-elevation samples indicate that the rate of cooling had decreased significantly. We interpret this to indicate decreased fault displacement along the Leo Pargil shear zone and possibly a transition to the Kaurik-Chango normal fault system between 10 and 6 Ma. AFT ages from lower elevations indicate accelerated cooling since the Pliocene that cannot be related to pure fault displacement, and therefore may reflect more pronounced regionally distributed and erosion-driven exhumation
Inversions of an individual geophysical data set can be highly nonunique, and it is generally difficult to determine petrophysical parameters from geophysical data. We show that both issues can be addressed by adopting a statistical multiparameter approach that requires the acquisition, processing, and separate inversion of two or more types of geophysical data. To combine information contained in the physical-property models that result from inverting the individual data sets and to estimate the spatial distribution of petrophysical parameters in regions where they are known at only a few locations. we demonstrate the potential of the fuzzy c-means (FCM) clustering technique. After testing this new approach on synthetic data, we apply it to limited crosshole georadar, crosshole seismic, gamma-log, and slug-test data acquired within a shallow alluvial aquifer. The derived multiparameter model effectively outlines the major sedimentary units observed in numerous boreholes and provides plausible estimates for the spatial distributions of gamma-ray emitters and hydraulic conductivity
The statistics of time delays between successive earthquakes has recently been claimed to be universal and to show the existence of clustering beyond the duration of aftershock bursts. We demonstrate that these claims are unjustified. Stochastic simulations with Poissonian background activity and triggered Omori-type aftershock sequences are shown to reproduce the interevent-time distributions observed on different spatial and magnitude scales in California. Thus the empirical distribution can be explained without any additional long-term clustering. Furthermore, we find that the shape of the interevent-time distribution, which can be approximated by the gamma distribution, is determined by the percentage of main-shocks in the catalog. This percentage can be calculated by the mean and variance of the interevent times and varies between 5% and 90% for different regions in California. Our investigation of stochastic simulations indicates that the interevent-time distribution provides a nonparametric reconstruction of the mainshock magnitude-frequency distribution that is superior to standard declustering algorithm
[1] The Kyrgyz Range, the northernmost portion of the Kyrgyzstan Tien Shan, displays topographic evidence for lateral propagation of surface uplift and exhumation. The highest, most deeply dissected segment lies in the center of the range. To the east, topography and relief decrease, and preserved remnants of a Cretaceous regional erosion surface imply minimal amounts of bedrock exhumation. The timing of exhumation of range segments defines the lateral propagation rate of the range-bounding reverse fault and quantifies the time and erosion depth needed to transform a mountain range from a juvenile to a mature morphology. New multicompositional apatite fission track ( AFT) data from three transects from the eastern Kyrgyz Range, combined with published AFT data, demonstrate that the range has propagated over 110 km eastward over the last similar to 7 - 11 Myr. On the basis of the thermal and topographic evolutionary history, we present a model for a time-varying exhumation rate driven by rock uplift and changes in erodability and the timescale of geomorphic adjustment to surface uplift. Easily eroded, Cenozoic sedimentary rocks overlying resistant basement control early, rapid exhumation and exhibit slow surface uplift rates. As increasing amounts of resistant basement are exposed, exhumation rates decrease while surface uplift rates are sustained or increase, thereby growing topography. As the range becomes high enough to cause ice accumulation and to develop steep river valleys, fluvial and glacial erosion becomes more powerful, and exhumation rates once again increase. Independently determined range-normal shortening rates also varied over time, suggesting a feedback between erosional efficiency and shortening rate
The iron speciation in hydrous haplotonalitic and haplogranitic silicate glasses was studied using XAFS spectroscopy and transmission electron microscopy (TEM). Spectral features occurring at the main crest of the XANES at the iron K-edge of hydrous glasses indicate contributions to the spectra by iron-moieties present in a more ordered structural environment than found in the dry glass. These differences are also suggested by analysis of the EXAFS. These effects are not completely suppressed even for those samples that were quenched with a higher cooling rate. Strongest differences to the dry glass are observed for a sample that was quenched slowly through the temperature of glass transformation. Crystals (60 to 1500 nm in size) of magnetite, maghemite and another unidentified phase were observed in this sample by TEM, whereas no crystals were found in samples quenched with regular or high cooling rates. In-situ XANES measurements up to 700 degrees C and 500 MPa were performed to reveal the origin (i.e., during synthesis or quench) of the structural differences for those hydrous glasses that do not display any detectable crystallization. The comparison of XANES spectra collected on Fe2+ in water-saturated haplogranitic melt at 700 degrees C and 500 MPa and on Fe2+ in dry melt at 1150 degrees C shows that the local structural environment of Fe2+ in both systems is similar. This indicates that there is no detectable and direct influence of water on the local structure around iron in this type of melt. Hence, the differences observed between hydrous and dry glasses can only be related to artefacts formed during the quench process. (c) 2006 Elsevier B.V. All rights reserved
In this paper, two sets of earthquake ground-motion relations to estimate peak ground and response spectral acceleration are developed for sites in southern Spain and in southern Norway using a recently published composite approach. For this purpose seven empirical ground-motion relations developed from recorded strong-motion data from different parts of the world were employed. The different relations were first adjusted based on a number of transformations to convert the differing choices of independent parameters to a single one. After these transformations, which include the scatter introduced, were performed, the equations were modified to account for differences between the host and the target regions using the stochastic method to compute the host-to-target conversion factors. Finally functions were fitted to the derived ground-motion estimates to obtain sets of seven individual equations for use in probabilistic seismic hazard assessment for southern Spain and southern Norway. The relations are compared with local ones published for the two regions. The composite methodology calls for the setting up of independent logic trees for the median values and for the sigma values, in order to properly separate epistemic and aleatory uncertainties after the corrections and the conversions
Receiver functions (RF) are used to investigate the upper mantle structure beneath the Eifel, the youngest volcanic area of Central Europe. Data from 96 teleseismic events recorded by 242 seismological stations from permanent and a temporary network has been analysed. The temporary network operated from 1997 November to 1998 June and covered an area of approximately 400 x 250 km(2) centred on the Eifel volcanic fields. The average Moho depth in the Eifel is approximately 30 km, thinning to ca. 28 km under the Eifel volcanic fields. RF images suggest the existence of a low velocity zone at about 60-90 km depth under the West Eifel. This observation is supported by P- and S-wave tomographic results and absorption (but the array aperture limits the resolution of the tomographic methods to the upper 400 km). There are also indications for a zone of elevated velocities at around 200 km depth, again in agreement with S-wave and absorption tomographic results. This anomaly is not visible in P-wave tomography and could be due to S-wave anisotropy. The RF anomalies at the Moho, at 60-90 km, and near 200 km depth have a lateral extent of about 100 km. The 410 km discontinuity under the Eifel is depressed by 15-25 km, which could be explained by a maximum temperature increase of +200 degrees C to +300 degrees C. In the 3-D RF image of the Eifel Plume we also notice two additional currently unexplained conversions between 410 and 550 km depth. They could represent remnants of previous subduction or anomalies due to delayed phase changes. The lateral extent of these conversions and the depression of the 410 km discontinuity is about 200 km. The 660 km discontinuity does not show any depth deviation from its expected value. Our observations are consistent with interpretation in terms of an upper mantle plume but they do not rule out connections to processes at larger depth
Reaction rims of titanite on ilmenite are described in samples from four terranes of amphibolite-facies metapelites and amphibolites namely the Tamil Nadu area, southern India; the Val Strona, area of the Ivrea-Verbano Zone, northern Italy, the Bamble Sector, southern Norway, and the northwestern Austroalpine Otztal Complex. The titanite rims, and hence the stability of titanite (CaTiSiO4O) and Al-OH titanite, i.e. vuaganatite (hypothetical end-member CaAlSiO4OH), are discussed in the light of fH(2)O- and fO(2)-buffered equilibria involving clinopyroxene, amphibole, biotite, ilmenite, magnetite, and quartz in the systems CaO-FeO/Fe2O3-TiO2-SiO2-H2O-O-2 (CFTSH) and CaO-FeO/Fe2O3-Al2O3- SiO2-H2O-O-2 (CFASH) present in each of the examples. Textural evidence suggests that titanite reaction rims on ilmenite in rocks from Tamil Nadu, Val Strona, and the Bamble Sector originated most likely due to hydration reactions such as clinopyroxene + ilmenite +quartz+ H2O = amphibole +titanite and oxidation reactions such as amphibole + ilmenite + O-2 = titanite + magnetite + quartz + H2O during amphibolite-facies metamorphism, or, as in the case of the Otztal Complex, during a subsequent greenschist-facies overprint. Overstepping of these reactions requires fH(2)O and fO(2) to be high for titanite formation, which is also in accordance with equilibria involving Al-OH titanite. This study shows that, in addition to P, T, bulk-rock composition and composition of the coexisting fluid, fO(2) and fH(2)O also play an important role in the formation of Al-bearing titanite during amphibolite- and greenschist-facies metamorphism.
Several previous studies have addressed the diagenetic evolution of heterozoan carbonate assemblages. Generally it is assumed that early diagenetic processes in heterozoan settings are mainly destructive, including abrasion and dissolution on the sea floor. Constructive diagenesis (cementation) is delayed to later stages in the burial environment, with pressure solution of calcitic grains acting as a cement source. This paper presents a study of Oligo- Miocene inner- to outer-ramp heterozoan carbonates from the Central Mediterranean (Maltese Islands and Sicily) indicating that early diagenetic processes are more important than previously assumed. Four to five different cement types, including fibrous, two types of epitaxial, bladed and blocky cement, are distinguished based on transmitted light microscopy. Cathodoluminescence microscopy allowed a differentiation between primary high-Mg calcitic (fibrous and epitaxial cement I) and primary low-Mg calcitic (epitaxial cement II, bladed and blocky) cements. Stable-isotope data indicate cement precipitation from marine, marine-derived, and meteoric waters. Trace-element analyses point to cementation in an open system (Maltese Islands) and a closed system (Sicily). Our investigations show that the majority of constructive diagenetic processes in these rocks occurs rather early in the shallow, marine burial environment, which is transitional between the marine seafloor and the deep-burial diagenctic environment. The main cement source in this environment is assumed to be aragonite. We suggest careful consideration of the importance of aragonitic components in fossil heterozoan settings, which seem to be more abundant than previously assumed and can act as a major early cement source. Due to the low preservation potential of these components, detailed geochemical studies are necessary to detect aragonite as the cement source. Our findings also have implications when considering the reservoir qualities of these rocks, because primary porosity can be occluded early and secondary porosity is not preserved
[ 1] The Eastern Cordillera of Colombia is key to understanding the role of inherited basement anisotropies in the evolution of active noncollisional mountain belts. In particular, the Rio Blanco-Guatiquia region of the Eastern Cordillera is exemplary in displaying a variety of phenomena that document the importance of the orientation, geometry, and segmentation of preorogenic anisotropies. We document the first unambiguous evidence that extensional basement structures played an important role in determining the locus of deformation during contractional reactivation in the Eastern Cordillera. Detailed structural field mapping and analysis of industry seismic reflection profiles have helped to identify the inherited San Juanito, Naranjal, and Servita normal faults and associated transfer faults as important structures that were inverted during the Cenozoic Andean orogeny. Apparently, the more internal faults in the former rift basin were not properly oriented for an efficient reactivation in contraction. However, these faults have a fundamental role as strain risers, as folding is concentrated west of them. In contrast, reactivated normal faults such as the more external Servita fault are responsible for uplifting the eastern flank of the Eastern Cordillera. In addition, these structures are adjacent and intimately linked to the development of thin-skinned faults farther east. In part, the superimposed compression in this prestrained extensional region is compensated by lateral escape. The dominant presence of basement involved buckling and thrusting, and the restricted development of thin-skinned thrusting in this inversion orogen makes the Eastern Cordillera a close analog to the intraplate Atlas Mountains of Morocco and other inverted sectors of the Andean orogen farther south
In low-seismicity regions, such as France or Germany, the estimation of probabilistic seismic hazard must cope with the difficult identification of active faults and with the low amount of seismic data available. Since the probabilistic hazard method was initiated, most studies assume a Poissonian occurrence of earthquakes. Here we propose a method that enables the inclusion of time and space dependences between earthquakes into the probabilistic estimation of hazard. Combining the seismicity model Epidemic Type Aftershocks-Sequence (ETAS) with a Monte Carlo technique, aftershocks are naturally accounted for in the hazard determination. The method is applied to the Pyrenees region in Southern France. The impact on hazard of declustering and of the usual assumption that earthquakes occur according to a Poisson process is quantified, showing that aftershocks contribute on average less than 5 per cent to the probabilistic hazard, with an upper bound around 18 per cent
Leech et al. [Mary L. Leech, S. Singh, A.K. Jain, Simon L. Klemperer and R.M. Manickavasagam, Earth and Planetary Science Letters 234 (2005) 83-97], present 3 clusters of ages for growth stages in zircon from quartzo- feldspathic gneisses hosting coesite-bearing eclogite from the Tso Morari Complex, NW India. These age clusters, from oldest to youngest, are interpreted to represent the age of ultrahigh-pressure metamorphism, a subsequent eclogite facies overprint and a later amphibolite facies retrogression and require subduction of Indian crust to have started earlier than previously accepted. However, no petrographic evidence, such as inclusions in the zircons relating to particular metamorphic events, is presented to substantiate the proposed sequence of metamorphic stages. Previously published data from eclogites of the same area indicate that coesite-eclogite is not the first but at least the second eclogite facies stage. In addition, the newly proposed time interval between coesite-eclogite and the amphibolite facies overprint is longer than previously indicated by diffusion modelling of natural garnet-garnet couples in eclogite. Neither the age of ultrahigh-pressure metamorphism nor the timing of initiation of subduction is reliably constrained by the presented data
The estimation of minimum-misfit stochastic models from empirical ground-motion prediction equations
(2006)
In areas of moderate to low seismic activity there is commonly a lack of recorded strong ground motion. As a consequence, the prediction of ground motion expected for hypothetical future earthquakes is often performed by employing empirical models from other regions. In this context, Campbell's hybrid empirical approach (Campbell, 2003, 2004) provides a methodological framework to adapt ground-motion prediction equations to arbitrary target regions by using response spectral host-to-target-region-conversion filters. For this purpose, the empirical ground-motion prediction equation has to be quantified in terms of a stochastic model. The problem we address here is how to do this in a systematic way and how to assess the corresponding uncertainties. For the determination of the model parameters we use a genetic algorithm search. The stochastic model spectra were calculated by using a speed-optimized version of SMSIM (Boore, 2000). For most of the empirical ground-motion models, we obtain sets of stochastic models that match the empirical models within the full magnitude and distance ranges of their generating data sets fairly well. The overall quality of fit and the resulting model parameter sets strongly depend on the particular choice of the distance metric used for the stochastic model. We suggest the use of the hypocentral distance metric for the stochastic Simulation of strong ground motion because it provides the lowest-misfit stochastic models for most empirical equations. This is in agreement with the results of two recent studies of hypocenter locations in finite-source models which indicate that hypocenters are often located close to regions of large slip (Mai et al., 2005; Manighetti et al., 2005). Because essentially all empirical ground-motion prediction equations contain data from different geographical regions, the model parameters corresponding to the lowest-misfit stochastic models cannot necessarily be expected to represent single, physically realizable host regions but to model the generating data sets in an average way. In addition, the differences between the lowest-misfit stochastic models and the empirical ground-motion prediction equation are strongly distance, magnitude, and frequency dependent, which, according to the laws of uncertainty propagation, will increase the variance of the corresponding hybrid empirical model predictions (Scherbaum et al., 2005). As a consequence, the selection of empirical ground-motion models for host-to-target-region conversions requires considerable judgment of the ground-motion analyst
In the estimate of dispersion with the help of wavelet analysis considerable emphasis has been put on the extraction of the group velocity using the modulus of the wavelet transform. In this paper we give an asymptotic expression of the full propagator in wavelet space that comprises the phase velocity as well. This operator establishes a relationship between the observed signals at two different stations during wave propagation in a dispersive and attenuating medium. Numerical and experimental examples are presented to show that the method accurately models seismic wave dispersion and attenuation
We report measurements of radial and azimuthal anisotropy in the tipper mantle beneath southern and central Finland, which we obtained by array analysis of fundamental-mode Rayleigh and Love waves. Azimuthally averaged phase velocities were analysed in the period range 15 to 190 s for Rayleigh waves and 15 to 100 s for Love waves. The azimuthal variation of the Rayleigh wave phase velocities was obtained in the period range 20 to 100s. The limited depth resolution of fundamental-mode surface waves necessitated strong damping constraints in the inversion for anisotropic parameters. We investigated the effects of non-unicity on the final model by experimenting with varying model geometries. The radial anisotropy beneath Finland can be explained by a lithosphere at least 200km thick, predominantly (> 50% by volume) composed of olivine crystals having their a-axes randomly distributed in the horizontal plane. On the contrary, the measured lithospheric azimuthal anisotropy is small. This call be reconciled with body-wave observations made in the area that indicate a complex pattern of rapidly varying anisotropy. Below 200-250km depth, that is below the petrologic lithosphere as revealed by xenolith analyses conducted in the area, the magnitude of the azimuthal anisotropy increases and would be compatible with a mantle containing 15-20% by volume of olivine crystals whose a-axes are coherently aligned in the N-NE direction. The alignment of the a-axes is off the direction of present-day absolute plate motion in either the no-net-rotation or hot-spot reference frame, currently N55-N60. We interpret this mismatch as evidence for a complex convective flow pattern of the mantle beneath the shield, which, by inference, is decoupled from the overlying lithosphere.
In the humid tropics, continuing high deforestation rates are seen alongside an increasing expansion of secondary forests. In order to understand and model the consequences of these dynamic land-use changes for regional water cycles, the response of soil hydraulic properties to forest disturbance and recovery has to be quantified.At a site in the Brazilian Amazonia, we annually monitored soil infiltrability and saturated hydraulic conductivity (K-s) at 12.5, 20 cm, and 50 cm soil depth after manual forest conversion to pasture (year zero to four after pasture establishment), and during secondary succession after pasture abandonment (year zero to seven after pasture abandonment). We evaluated the hydrological consequences of the detected changes by comparing the soil hydraulic properties with site-specific rainfall intensities and hydrometric observations. Within one year after grazing started, infiltrability and K-s at 12.5 and 20 cm depth decreased by up to one order of magnitude to levels which are typical for 20-year-old pasture. In the three subsequent monitoring years, infiltrability and K-s remained stable. Land use did not impact on subsoil permeability. Whereas infiltrability values are large enough to allow all rainwater to infiltrate even after the conversion, the sudden decline of near-surface K-s is of hydrological relevance as perched water tables and overland flow occur more often on pastures than in forests at our study site. After pasture abandonment and during secondary succession, seven years of recovery did not suffice to significantly increase infiltrability and K-s at 12.5 depth although a slight recovery is obvious. At 20 cm soil depth, we detected a positive linear increase within the seven-year time frame but annual means did not differ significantly. Although more than a doubling of infiltrability and K-s is still required to achieve pre-disturbance levels, which will presumably take more than a decade, the observed slight increases of K-s might already decrease the probability of perched water table generation and overland flow development well before complete recovery.
Major earthquakes ( M > 8) have repeatedly ruptured the Nazca-South America plate interface of south-central Chile involving meter scale land-level changes. Earthquake recurrence intervals, however, extending beyond limited historical records are virtually unknown, but would provide crucial data on the tectonic behavior of forearcs. We analyzed the spatiotemporal pattern of Holocene earthquakes on Santa Maria Island (SMI; 37 degrees S), located 20 km off the Chilean coast and approximately 70 km east of the trench. SMI hosts a minimum of 21 uplifted beach berms, of which a subset were dated to calculate a mean uplift rate of 2.3 +/- 0.2 m/ky and a tilting rate of 0.022 +/- 0.002 degrees/ky. The inferred recurrence interval of strandline-forming earthquakes is similar to 180 years. Combining coseismic uplift and aseismic subsidence during an earthquake cycle, the net gain in strandline elevation in this environment is similar to 0.4 m per event
Three ODP sites located on the Marion Plateau, Northeast Australian margin, were investigated for clay mineral and bulk mineralogy changes through the early to middle Miocene. Kaolinite to smectite (K/S) ratios, as well as mass accumulation rates of clays, point to a marked decrease in accumulation of smectite associated with an increase in accumulation of kaolinite starting at similar to 15.6 Ma, followed by a, second increase in accumulation of kaolinite at similar to 13.2 Ma. Both of these increases are correlative to an increase in the calcite to detritus ratio. Comparison of our record with published precipitation proxies from continental Queensland indicates that increases in kaolinite did not correspond to more intense tropical-humid conditions, but instead to periods of greater aridity. Three mechanisms are explored to explain the temporal trends in clad, on the Marion Plateau: sea-level changes, changes in oceanic currents, and denudation of the Australian continent followed by reworking and eolian transport of clays. Though low mass accumulation rates of kaolinite are compatible with a possible contribution of eolian material after 14 Ma, when Australia became more arid, the lateral distribution of kaolinite along slope indicates mainly fluvial input for all clays and thus rules out this mechanism as well as oceanic current transport as the main controls behind clay accumulation on the plateau. We propose a model explaining the good correlation between long-term sea-level fall, decrease in smectite accumulation, increase in kaolinite accumulation and increase in carbonate input to the distal slope locations. We hypothesize that during low sea level and thus periods of drier continental climate in Queensland, early Miocene kaolinite-rich lacustrine deposits were being reworked, and that the progradation of the heterozoan carbonate platforms towards the basin center favored input of carbonate to the distal slope sites. The major find of our study is that increase kaolinite fluxes on the Queensland margin during the early and middle Miocene did not reflect the establishment of a tropical climate, and this stresses that care must be taken when reconstructing Australian climate based on deep-sea clay records alone.
We employ P to S converted waveforms to investigate effects of the hot mantle plume on seismic discontinuities of the crust and upper mantle. We observe the Moho at depths between 13 and 17 km, regionally covered by a strong shallow intracrustal converted phase. Coherent phases on the transverse component indicate either dipping interfaces, 3- D heterogeneities or lower crustal anisotropy. We find anomalies related to discontinuities in the upper mantle down to the transition zone evidently related to the hot mantle plume. Lithospheric thinning is confirmed in greater detail than previously reported by Li et al., and we determine the dimensions of the low-velocity zone within the asthenosphere with greater accuracy. Our study mainly focuses on the temperature-pressure dependent discontinuities of the upper mantle transition zone. Effects of the hot diapir on the depths of mineral phase transitions are verified at both major interfaces at 410 and 660 km. We determine a plume radius of about 200 km at the 660 km discontinuity with a core zone of about 120 km radius. The plume conduit is located southwest of Big Island. A conduit tilted in the northeast direction is required in the upper mantle to explain the observations. The determined positions of deflections of the discontinuities support the hypothesis of decoupled upper and lower mantle convection
In June 2013, widespread flooding and consequent damage and losses occurred in Central Europe, especially in Germany. This paper explores what data are available to investigate the adverse impacts of the event, what kind of information can be retrieved from these data and how well data and information fulfil requirements that were recently proposed for disaster reporting on the European and international levels. In accordance with the European Floods Directive (2007/60/EC), impacts on human health, economic activities (and assets), cultural heritage and the environment are described on the national and sub-national scale. Information from governmental reports is complemented by communications on traffic disruptions and surveys of flood-affected residents and companies.
Overall, the impacts of the flood event in 2013 were manifold. The study reveals that flood-affected residents suffered from a large range of impacts, among which mental health and supply problems were perceived more seriously than financial losses. The most frequent damage type among affected companies was business interruption. This demonstrates that the current scientific focus on direct (financial) damage is insufficient to describe the overall impacts and severity of flood events.
The case further demonstrates that procedures and standards for impact data collection in Germany are widely missing. Present impact data in Germany are fragmentary, heterogeneous, incomplete and difficult to access. In order to fulfil, for example, the monitoring and reporting requirements of the Sendai Framework for Disaster Risk Reduction 2015–2030 that was adopted in March 2015 in Sendai, Japan, more efforts on impact data collection are needed.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Widespread flooding in June 2013 caused damage costs of €6 to 8 billion in Germany, and awoke many memories of the floods in August 2002, which resulted in total damage of €11.6 billion and hence was the most expensive natural hazard event in Germany up to now. The event of 2002 does, however, also mark a reorientation toward an integrated flood risk management system in Germany. Therefore, the flood of 2013 offered the opportunity to review how the measures that politics, administration, and civil society have implemented since 2002 helped to cope with the flood and what still needs to be done to achieve effective and more integrated flood risk management. The review highlights considerable improvements on many levels, in particular (1) an increased consideration of flood hazards in spatial planning and urban development, (2) comprehensive property-level mitigation and preparedness measures, (3) more effective flood warnings and improved coordination of disaster response, and (4) a more targeted maintenance of flood defense systems. In 2013, this led to more effective flood management and to a reduction of damage. Nevertheless, important aspects remain unclear and need to be clarified. This particularly holds for balanced and coordinated strategies for reducing and overcoming the impacts of flooding in large catchments, cross-border and interdisciplinary cooperation, the role of the general public in the different phases of flood risk management, as well as a transparent risk transfer system. Recurring flood events reveal that flood risk management is a continuous task. Hence, risk drivers, such as climate change, land-use changes, economic developments, or demographic change and the resultant risks must be investigated at regular intervals, and risk reduction strategies and processes must be reassessed as well as adapted and implemented in a dialogue with all stakeholders.