Refine
Has Fulltext
- no (56)
Year of publication
- 2004 (56) (remove)
Document Type
- Article (56) (remove)
Language
- English (56) (remove)
Is part of the Bibliography
- yes (56) (remove)
Institute
- Institut für Geowissenschaften (56) (remove)
In this paper we present densely sampled fumarole temperature data, recorded continuously at a high-temperature fumarole of Mt. Merapi volcano (Indonesia). These temperature time series are correlated with continuous records of rainfall and seismic waveform data collected at the Indonesian - German multi-parameter monitoring network. The correlation analysis of fumarole temperature and precipitation data shows a clear influence of tropical rain events on fumarole temperature. In addition, there is some evidence that rainfall may influence seismicity rates, indicating interaction of meteoric water with the volcanic system. Knowledge about such interactions is important, as lava dome instabilities caused by heavy-precipitation events may result in pyroclastic flows. Apart from the strong external influences on fumarole temperature and seismicity rate, which may conceal smaller signals caused by volcanic degassing processes, the analysis of fumarole temperature and seismic data indicates a statistically significant correlation between a certain type of seismic activity and an increase in fumarole temperature. This certain type of seismic activity consists of a seismic cluster of several high-frequency transients and an ultra-long-period signal (< 0.002 Hz), which are best observed using a broadband seismometer deployed at a distance of 600 m from the active lava dome. The corresponding change in fumarole temperature starts a few minutes after the ultra-long-period signal and simultaneously with the high-frequency seismic cluster. The change in fumarole temperature, an increase of 5 degreesC on average, resembles a smoothed step. Fifty-four occurrences of simultaneous high-frequency seismic cluster, ultra-long period signal and increase of fumarole temperature have been identified in the data set from August 2000 to January 2001. The observed signals appear to correspond to degassing processes in the summit region of Mt. Merapi. (C) 2004 Elsevier B.V. All rights reserved
We remelted and analyzed crystallized silicate melt inclusions in quartz from a porphyritic albitezinnwaldite microgranite dike to determine the composition of highly evolved, shallowly intruded, Li- and F-rich granitic magma and to investigate the role of crystal fractionation and aqueous fluid exsolution in causing the extreme extent of magma differentiation. This dike is intimately associated with tin- and tungsten-mineralized granites of Zinnwald, Erzgebirge, Germany. Prior research on Zinnwald granite geochemistry was limited by the effects of strong and pervasive greisenization and alkali-feldspar metasomatism of the rocks. These melt inclusions, however, provide important new constraints on magmatic and mineralizing processes in Zinnwald magmas. The mildly peraluminous granitic melt inclusions are strongly depleted in CAFEMIC constituents (e.g., CaO, FeO, MgO, TiO2), highly enriched in lithophile trace elements, and highly but variably enriched in F and Cl. The melt inclusions contain up to several thousand ppm Cl and nearly 3 wt% F, on average; several inclusions contain more than 5 wt% F. The melt inclusions are geochemically similar to the corresponding whole-rock sample, except that the former contain much more F and less CaO, FeO, Zr, Nb, Sr, and Ba. The Sr and Ba abundances are very low implying the melt inclusions represent magma that was more evolved than that represented by the bulk rock. Relationships involving melt constituents reflect increasing lithophile-element and halogen abundances in residual melt with progressive magma differentiation. Modeling demonstrates that differentiation was dominated by crystal fractionation involving quartz and feldspar and significant quantities of topaz and F-rich zinnwaldite. The computed abundances of the latter phases greatly exceed their abundances in the rocks, suggesting that the residual melt was separated physically from phenocrysts during magma movement and evolution. Interactions of aqueous fluids with silicate melt were also critical to magma evolution. To better understand the role of halogen-charged, aqueous fluids in magmatic differentiation and in subsequent mineralization and metasomatism of the Zinnwald granites, Cl-partitioning experiments were conducted with a F-enriched silicate melt and aqueous fluids at 2,000 bar (200 MPa). The results of the experimentally determined partition coefficients for Cl and F, the compositions of fluid inclusions in quartz and other phenocrysts, and associated geochemical modeling point to an important role of magmatic-hydrothermal fluids in influencing magma geochemistry and evolution. The exsolution of halogen-charged fluids from the Li- and F- enriched Zinnwald granitic magma modified the Cl, alkali, and F contents of the residual melt, and may have also sequestered Li, Sri, and W from the melt. Many of these fluids contained strongly elevated F concentrations that were equivalent to or greater than their Cl abundances. The exsolution of F-, Cl-, Li-, +/- W- and Sn-bearing hydrothermal fluids from Zinnwald granite magmas was important in effecting the greisenizing and alkali-feldspathizing metasomatism of the granites and the concomitant mineralization
Fe K-edge X-ray absorption near edge structure (XANES) and Mossbauer spectra were collected on synthetic glasses of basaltic composition and of glasses on the sodium oxide-silica binary to establish a relation between the pre- edge of the XANES at the K-edge and the Fe oxidation state of depolymerised glasses. Charges of sample material were equilibrated at ambient pressure, superliquidus temperatures and oxygen fugacities that were varied over a range of about 15 orders of magnitude. Most experiments were carried out in gas-flow furnaces, either with pure oxygen, air, or different CO/CO2 mixtures. For the most reduced conditions, the samples charges were enclosed together with a pellet of the IQF oxygen buffer in an evacuated silica glass ampoule. Fe3+/SigmaFe x 100 of the samples determined by Mossbauer spectroscopy range between 0% and 100%. Position and intensity of the pre-edge centroid position vary strongly depending on the Fe oxidation state. The pre-edge centroid position and the Fe oxidation state determined by Mossbauer spectroscopy are nonlinearly related and have been fitted by a quadratic polynomial. Alternatively, the ratio of intensities measured at positions sensitive to Fe2+ and Fe3+, respectively, provides an even more sensitive method. Pre- edge intensities of the sample suite indicate average Fe co-ordination between 4 and 6 for all samples regardless of oxidation state. A potential application of the calibration given here opens the possibility of determining Fe oxidation state in glasses of similar compositions with high spatial resolution by use of a Micro-XANES setup (e.g., glass inclusions in natural minerals). (C) 2004 Elsevier B.V. All rights reserved
The use of ground-motion-prediction equations to estimate ground shaking has become a very popular approach for seismic-hazard assessment, especially in the framework of a logic-tree approach. Owing to the large number of existing published ground-motion models, however, the selection and ranking of appropriate models for a particular target area often pose serious practical problems. Here we show how observed around-motion records can help to guide this process in a systematic and comprehensible way. A key element in this context is a new, likelihood based, goodness-of-fit measure that has the property not only to quantify the model fit but also to measure in some degree how well the underlying statistical model assumptions are met. By design, this measure naturally scales between 0 and 1, with a value of 0.5 for a situation in which the model perfectly matches the sample distribution both in terms of mean and standard deviation. We have used it in combination with other goodness-of-fit measures to derive a simple classification scheme to quantify how well a candidate ground-rnotion-prediction equation models a particular set of observed-response spectra. This scheme is demonstrated to perform well in recognizing a number of popular ground-motion models from their rock-site- recording, subsets. This indicates its potential for aiding the assignment of logic-tree weights in a consistent and reproducible way. We have applied our scheme to the border region of France, Germany, and Switzerland where the M-w 4.8 St. Die earthquake of 22 February 2003 in eastern France recently provided a small set of observed-response spectra. These records are best modeled by the ground-motion-prediction equation of Berge-Thierry et al. (2003), which is based on the analysis of predominantly European data. The fact that the Swiss model of Bay et al. (2003) is not able to model the observed records in an acceptable way may indicate general problems arising from the use of weak-motion data for strong-motion prediction
Along the Southern Himalayan Front (SHF), areas with concentrated precipitation coincide with rapid exhumation, as indicated by young mineral cooling ages. Twenty new, young ( < 1-5 Ma) apatite fission track (AFT) ages have been obtained from the Himalayan Crystalline Core along the Sutlej Valley, NW India. The AFT ages correlate with elevation, but show no spatial relationship to tectonic structures, such as the Main Central Thrust or the Southern Tibetan Fault System. Monsoonal precipitation in this region exerts a strong influence on erosional surface processes. Fluvial erosional unloading along the SHF is focused on high mountainous areas, where the orographic barrier forces out > 80% of the annual precipitation. AFT cooling ages reveal a coincidence between rapid erosion and exhumation that is focused in a similar to 50-70-km-wide sector of the Himalaya, rather than encompassing the entire orogen. Assuming simplified constant exhumation rates, the rocks of two age vs. elevation transects were exhumed at similar to 1.4 +/- 0.2 and similar to 1.1 +/- 0.4 mm/a with an average cooling rate of similar to 40-50degreesC/Ma during Pliocene-Quarternary time. Following other recently published hypotheses regarding the relation between tectonics and climate in the Himalaya, we suggest that this concentrated loss of material was accommodated by motion along a back-stepping thrust to the south and a normal fault zone to the north as part of an extruding wedge. Climatically controlled erosional processes focus on this wedge and suggest that climatically controlled surface processes determine tectonic deformation in the internal part of the Himalaya. (C) 2004 Elsevier B.V. All rights reserved
In recent years, H/V measurements have been increasingly used to map the thickness of sediment fill in sedimentary basins in the context of seismic hazard assessment. This parameter is believed to be an important proxy for the site effects in sedimentary basins (e.g. in the Los Angeles basin). Here we present the results of a test using this approach across an active normal fault in a structurally well known situation. Measurements on a 50 km long profile with 1 km station spacing clearly show a change in the frequency of the fundamental peak of H/V ratios with increasing thickness of the sediment layer in the eastern part of the Lower Rhine Embayment. Subsequently, a section of 10 km length across the Erft-Sprung system, a normal fault with ca. 750 m vertical offset, was measured with a station distance of 100 m. Frequencies of the first and second peaks and the first trough in the H/V spectra are used in a simple resonance model to estimate depths of the bedrock. While the frequency of the first peak shows a large scatter for sediment depths larger than ca. 500 m, the frequency of the first trough follows the changing thickness of the sediments across the fault. The lateral resolution is in the range of the station distance of 100 m. A power law for the depth dependence of the S-wave velocity derived from down hole measurements in an earlier study [Budny, 1984] and power laws inverted from dispersion analysis of micro array measurements [Scherbaum et al., 2002] agree with the results from the H/V ratios of this study
We derive a slip rate for a thrust at the central Qilian Shan mountain front by combining structural investigations, satellite imagery, topographic profiling, luminescence dating, and Be-10 exposure dating. The seismically active Zhangye thrust transects late Pleistocene alluvial fan deposits and forms a prominent north facing scarp. The fault consists of two segments that differ in orientation, scarp height, and age. A series of loess-covered terraces records the uplift history of the western thrust segment. Loess accumulation on all terraces started at 8.5 +/- 1.5 kyr and postdates terrace formation. Gravels from the highest terrace yielded a Be-10 exposure age of 90 +/- 11 kyr, which dates the onset of faulting. With a displacement of 55-60 m derived from fault scarp profiles, this yields a vertical slip rate of 0.64 +/- 0.08 mm yr(-1). Along the eastern thrust segment, three Be-10 ages from the uplifted alluvial fan constrain that faulting started at similar to31 +/- 5 kyr. Together with a displacement of 25-30 m this leads to a vertical faulting rate of 0.88 +/- 0.16 mm yr(-1). A dip estimate of 40degrees to 60degrees for the fault plane combined with lower and upper limits of similar to0.6 and similar to0.9 mm yr(-1) for the vertical slip rate gives minimum and maximum horizontal shortening rates of 0.4 and 1.1 mm yr(-1) across the Zhangye thrust. Our results are consistent with geologic and GPS constraints, which suggest that NNE directed shortening across the northeastern Tibetan Plateau is distributed on several active faults with a total shortening rate of 4 to 10 mm yr(-1)
One of the major challenges in engineering seismology is the reliable prediction of site-specific ground motion for particular earthquakes, observed at specific distances. For larger events, a special problem arises, at short distances, with the source-to-site distance measure, because distance metrics based on a point-source model are no longer appropriate. As a consequence, different attenuation relations differ in the distance metric that they use. In addition to being a source of confusion, this causes problems to quantitatively compare or combine different ground- motion models; for example, in the context of Probabilistic Seismic Hazard Assessment, in cases where ground-motion models with different distance metrics occupy neighboring branches of a logic tree. In such a situation, very crude assumptions about source sizes and orientations often have to be used to be able to derive an estimate of the particular metric required. Even if this solves the problem of providing a number to put into the attenuation relation, a serious problem remains. When converting distance measures, the corresponding uncertainties map onto the estimated ground motions according to the laws of error propagation. To make matters worse, conversion of distance metrics can cause the uncertainties of the adapted ground-motion model to become magnitude and distance dependent, even if they are not in the original relation. To be able to treat this problem quantitatively, the variability increase caused by the distance metric conversion has to be quantified. For this purpose, we have used well established scaling laws to determine explicit distance conversion relations using regression analysis on simulated data. We demonstrate that, for all practical purposes, most popular distance metrics can be related to the Joyner-Boore distance using models based on gamma distributions to express the shape of some "residual function." The functional forms are magnitude and distance dependent and are expressed as polynomials. We compare the performance of these relations with manually derived individual distance estimates for the Landers, the Imperial Valley, and the Chi-Chi earthquakes
Anisotropic material properties are usually neglected during inversions for source parameters of earthquakes. In general anisotropic media, however, moment tensors for pure-shear sources can exhibit significant non-double-couple components. Such effects may be erroneously interpreted as an indication for volumetric changes at the source. Here we investigate effects of anisotropy on seismic moment tensors and radiation patterns for pure-shear and tensile-type sources. Anisotropy can significantly influence the interpretation of the source mechanisms. For example, the orientation of the slip within the fault plane may affect the total seismic moment. Also, moment tensors due to pure- shear and tensile faulting can have similar characteristics depending on the orientation of the elastic tensor. Furthermore, the tensile nature of an earthquake can be obscured by near-source anisotropic properties. As an application, we consider effects of inhomogeneous anisotropic properties on the seismic moment tensor and the radiation patterns of a selected type of micro-earthquakes observed in W-Bohemia. The combined effects of near-source and along- path anisotropy cause characteristic amplitude distortions of the P, S1 and S2 waves. However, the modeling suggests that neither homogeneous nor inhomogeneous anisotropic properties alone can explain the observed large non-double-couple components. The results also indicate that a correct analysis of the source mechanism, in principle, is achievable by application of anisotropic moment tensor inversion
Field observations, digital elevation model (DEM) data, and longitudinal profile analysis reveal a perched low- relief upland landscape in the Red River region, Yunnan Province, China, which correlates to an uplifted, regional low- relief landscape preserved over the eastern margin of the Tibetan Plateau. As with other major rivers of the plateau margin, the Red River has deeply incised the low-relief upland landscape, which we interpret to be the remnants of a pre- uplift or relict landscape. We examine longitudinal river profiles for 97 tributaries of the Red River. Most profiles consist of three segments separated by sharp knickpoints: an upper, low-gradient channel segment, a steeper middle channel segment, and a very steep lower channel segment. Upper channel segments correspond to the relict landscape and have not yet experienced river incision. Steeper middle and lower segments indicate onset of rapid, two-phase river incision, on the basis of which changes in external forcings, such as climate or uplift, can be inferred. In terms of two end-member scenarios, two-phase incision could be the result of pulsed plateau growth, in which relatively slow uplift during the first phase is followed by rapid uplift during the second phase, or it could reflect adjustments of the main channel to changing climate conditions against the backdrop of steady plateau growth. Reconstruction of the paleo-Red River indicates 1400 m river incision, 1400-1500 m surface uplift, and a maximum of 750 m vertical displacement across the northern Red River fault, elevating the northern Ailao Shan range above the surrounding relict landscape. On the basis of stratigraphic constraints, incision along the Red River likely began in Pliocene time
A recent K-Ar study elucidated that eruptive style in the eastern Izu peninsula changed from polygenetic to monogenetic volcano at 0.3-0.2 Ma. To narrow down the time of change, we determined 10 K-Ar ages on Togasayama Andesite of Amagi volcano, the youngest polygenetic volcano in the area, and Togasayama Monogenetic Volcano, one of the oldest monogenetic volcanoes in the area, which overlies a part of the Togasayama Andesite. Dating results showed that the Togasayama Andesite effused at least from 0.34 to 0.20 Ma, whereas the Togasayama Monogenetic Volcano erupted at 0.26- 0.29 Ma, suggesting that the northern part of the Togasayama Andesite effused after the eruption of the Togasayama Monogenetic Volcano. Considering previous data, it is therefore inferred that change of eruptive style in the eastern Izu area occurred during the period 0.29-0.20 Ma, with considerable overlap of both polygenetic and monogenetic volcanism
K-Ar ages of 37 samples collected from the Bicol peninsula, the Luzon island, Philippines, were determined by the unspiked sensitivity method in order to constrain the timing of initiation of subduction along the Philippine Trench. The measured K-Ar ages range from 0 to 7 Ma with two old outliers of 27 and 43 Ma. Together with K-Ar ages previously reported on volcanics in Leyte and eastern Mindanao, subduction volcanism has likely propagated from north to south: similar to 6.6 Ma in Bicol and similar to 3.5 Ma in Leyte and its vicinity. The temporal and spatial distribution suggests that the subduction volcanism started earlier in the north than in the south. This is consistent with the southern propagation of subduction along the Philippine Trench from similar to 8 Ma. (C) 2003 Elsevier Ltd. All rights reserved
The ellipticity of Rayleigh surface waves, which is an important parameter characterizing the propagation medium, is studied for several models with increasing complexity. While the main focus lies on theory, practical implications of the use of the horizontal to vertical component ratio (H/V-ratio) to Study the subsurface structure are considered as well. Love's approximation of the ellipticity for an incompressible layer over an incompressible half-space is critically discussed especially concerning its applicability for different impedance contrasts. The main result is an analytically exact formula of H/V for a 2-layer model of compressible media, which is a generalization of Love's formula. It turns out that for a limited range of models Love's approximation can be used also in the general case. (C) 2003 Elsevier B.V. All rights reserved
As part of the international refraction measurements in Central Europe in the year 2000, three profiles traversed the region of earthquake swarms in West-Bohemia/Vogtland. The shots were also recorded at the permanent stations of the local seismic networks. The travel times of P-waves, observed in the West-Bohemian region, are discussed and interpreted in the present paper. In general, significantly lower P-wave velocities were found in the Saxothuringian (northern) part of the studied area than in the adjacent southern parts. The observed travel times are interpreted separately for the individual geological units, in particular for the plutons, crystallinicum, and the Marianske Lazne(Marienbad) Complex. After smoothing the selected data using rational approximations, the Wiechert-Herglotz method was used to compute vertically inhomogeneous velocity models. The characteristic features of the derived models are relatively low P-wave velocities at the surface and prominent velocity increases within the uppermost crust down to a depth of about one kilometer
With controlled seismic sources and specifically designed receiver arrays, we image a subvertical boundary between two lithological blocks at the Arava Fault (AF) in the Middle East. The AF is the main strike-slip fault of the Dead Sea Transform (DST) in the segment between the Dead Sea and the Red Sea. Our imaging (migration) method is based on array beamforming and coherence analysis of P to P scattered seismic phases. We use a 1-D background velocity model and the direct P arrival as a reference phase. Careful resolution testing is necessary, because the target volume is irregularly sampled by rays. A spread function describing energy dispersion at localized point scatterers and synthetic calculations for large planar structures provides estimates of the resolution of the images. We resolve a 7 km long steeply dipping reflector offset roughly 1 km from the surface trace of the AF. The reflector can be imaged from about 1 km down to 4 km depth. Previous and ongoing studies in this region have shown a strong contrast across the fault: low seismic velocities and electrical resistivities to the west and high velocities and resistivities to the east of it. We therefore suggest that the imaged reflector marks the contrast between young sedimentary fill in the west and Precambrian rocks in the east. If correct, the boundary between the two blocks is offset about 1 km east of the current surface trace of the AF
VLT on-axis optical spectroscopy of the z = 0.144 radio-loud quasar HE 1434-1600 is presented. The spatially resolved spectra of the host galaxy are deconvolved and separated from those of the central quasar in order to study the dynamics of the stars and gas as well as the physical conditions of the ISM. We find that the host of HE 1434-1600 is an elliptical galaxy that resides in a group of at least 5 member galaxies, and that most likely experienced a recent collision with its nearest companion. Compared with other quasar host galaxies, HE 1434-1600 has a highly ionized ISM. The ionization state corresponds to that of typical Seyferts, but the ionized regions are not distributed in a homogeneous way around the QSO, and are located preferentially several kiloparsecs away from it. While the stellar absorption lines do not show any significant velocity field, the gas emission lines do. The observed gas velocity field is hard to reconcile with dynamical models involving rotating disk. modified Hubble laws or power laws, that all require extreme central masses (M > 10(9) M-circle dot) to provide only poor fit to the data. Power law models, which best fit the data, provide a total mass of M(<10 kpc) = 9.2 x 10(10) M-&ODOT;. We conclude that the recent interaction between HE 1434-1600 and its closest companion has strongly affected the gas velocity and ionization state, from the center of the galaxy to its most external parts