Refine
Year of publication
Document Type
- Article (2817) (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 (2817) (remove)
The rise and fall of the Classic Maya provides a textbook example of human social evolution. It is therefore significant to discover that the history of theMaya was so closely tied to environmental constraints. If Maya civilization could collapse under the weight of natural climate events, it is of more than academic interest to ponder how modem society will fare in the face of an uncertain climate in the years ahead. An understanding of how ancient cultures responded to climatic changes in the past may thus provide important lessons for humanity in the future
The petrology of two distinct granulite types in the Hengshan Mts, China, and tectonic implications
(2005)
The Archean to Proterozoic Hengshan Complex (North China Craton), comprises tonalitic and granodioritic gneisses with subordinate mafic lenses, pegmatites and granites. Amphibolite facies assemblages predominate, although granulite-facies relics are widespread, and greenschist-facies retrogression occurs in km-wide shear zones. Mafic lenses, locally abundant, occur as strongly deformed amphibolite (hornblende + plagioclase) boudins or sheets. In contrast to previously published models we find two series of mafic rocks with distinctly different granulite-facies evolutions. In the north of the complex, relict high-pressure mafic granulites are garnet + clinopyroxene-bearing rocks with a secondary development of orthopyroxene around both garnet (kelyphites) and clinopyroxene (coronas). South of the newly defined central, E-W-trending, Zhujiafang shear zone, numerous mafic boudins and less-deformed dykes exhibit a macroscopically visible magmatic texture with coronitic growth of metamorphic garnet (full of quartz inclusions) between the magmatic plagioclase and pyroxene domains. Additional orthopyroxene (after magmatic augite) and sodic rims to magmatic plagioclase clearly indicate medium-pressure granulite-facies metamorphism. These findings suggest tectonic juxtaposition in this area of three different structural levels of the same Proterozoic-imprinted crust: high-pressure granulite grade in the northern Hengshan, medium-pressure granulite grade in the southern Hengshan and amphibolite- to greenschist-facies grade in the Wutaishan to the SE. (c) 2004 Elsevier Ltd. All rights reserved
Possible causes to explain platform drowning have been hotly debated by carbonate sedimentologists for more than a decade now. In this paper, we present multiple evidence to explain the drowning of a carbonate megabank that covered most of the modem Northern Nicaragua Rise (NNR) during an interval spanning from late Oligocene to early Miocene by the interaction of several environmental factors. The recovery during ODP Leg 165 of late Oligocene to middle Miocene sedimentary sequences in the sub-seafloor of the modern channels and basin, Pedro Channel and Walton Basin, respectively, that dissect the NNR (Site 1000) and south of the rise in the Colombian Basin (Site 999), combined with information from dredged rock samples, allows us to explore in more detail the timing and possible mechanisms responsible for the drowning of the megabank and its relationship to Miocene climate change. The modern system of isolated banks and shelves dissected by a series of intervening seaways and basins on the NNR has evolved from a continuous, shallow-water carbonate "megabank' that extended from the Honduras/Nicaraguan mainland to the modern island of Jamaica. Available information suggests that this megabank broke apart and partially drowned in the late part of the late Oligocene at around 27 Ma and finally foundered during the late early Miocene around 20 Ma, resulting in limited neritic coral growth in the areas where the modern isolated carbonate banks and shelves are occurring today. Available information also suggests that the southern and central parts of Pedro Channel were already a deep-water area before the major episode of platform drowning, and its formation predates the initiation of the Caribbean Current. However, after the partial drowning of the megabank, the channel has become a major pathway for the Caribbean Current. Stratigraphic units identified in deep-water carbonates sampled at ODP Sites 999 and 1000 help to constrain the environmental setting leading to the drowning of the banks. Changes in lithology and mass accumulation rates of both the carbonate and non-carbonate fraction parallel stable isotope shifts and likely indicate regional changes in climate and circulation during the late Oligocene-middle Miocene interval. Carbonate mass accumulation rates (MARS) at Site 999 suggest increased regional productivity during the early Miocene. Terrigenous MARS at both Sites 999 and 1000 show a general increase from the Burdigalian through the Serravallian. The temporal association among episodes of neritic platform deposition, followed by increased productivity as identified by higher carbonate MARs and positive excursion in carbon isotopes, suggests that oceanographic changes such as local upwelling and nutrification have led to the partial drowning of the NNR "megabank". (c) 2005 Elsevier B.V. All rights reserved
Aftershocks rates seem to follow a power law decay, but the question of the aftershock frequency immediately after an earthquake remains open. We estimate an average aftershock decay rate within one day in southern California by stacking in time different sequences triggered by main shocks ranging in magnitude from 2.5 to 4.5. Then we estimate the time delay before the onset of the power law aftershock decay rate. For the last 20 years, we observe that this time delay suddenly increase after large earthquakes, and slowly decreases at a constant rate during periods of low seismicity. In a band-limited power law model such variations can be explained by different patterns of stress distribution at different stages of the seismic cycle. We conclude that, on regional length scales, the brittle upper crust exhibits a collective behavior reflecting to some extent the proximity of a threshold of fracturing
We investigate the controls on the architecture of coarse-grained delta progradational units (PUs) in the Pliocene Loreto basin (Baja California Sur, Mexico), a half-graben located on the western margin of the Gulf of California. Dorsey et al. (1997b) argued that delta progradation and transgression cycles in the basin were driven by episodic fault-controlled subsidence along the basin-bounding Loreto fault. Here we test this hypothesis by a detailed analysis of the sedimentary architecture of 11 exceptionally well-exposed, vertically arranged fluvio-deltaic PUs, each of which shows lateral facies transition from proximal alluvial facies palaeo-seaward into distal pro-delta facies. Of these 11 PUs, seven exhibit a lateral transition from a shoal water to Gilbert-delta facies associations as they are traced palaeo-seaward. This transition is characterised by down-transport development of foresets, which grow in height up to 35 m. Foreset units thicken in a basinward direction, with initially an oblique topset-foreset geometry that becomes increasingly sigmoidal. Each delta is capped by a shell bed that records drowning of the delta top. This systematic transition in delta architecture records increasing water depth through time during individual episodes of progradation. A mechanism that explains this transition is an accelerating rate of fault-controlled subsidence during each PU. During episodes of low slip rate, shoal-water deltas prograde across the submerged topography of the underlying delta unit. As displacement rate accelerates, increasing bathymetry at the delta front leads to steepening of foresets and initiation of Gilbert deltas. Subsequent delta drowning results from sediment starvation at the shoreline at high slip rates because of sediment trapping upstream. The observed delta architecture suggests that the long-term (> 100 kyr) history of slip on the Loreto fault was characterised by repetitive episodes of accelerating displacement accumulation. Such episodic fault behaviour is most likely to be because of variations in temporal and spatial strain partitioning between the Loreto fault and other faults in the Gulf of California. A physical explanation for the acceleration phenomenon involves evolving frictional properties on the episodically active Loreto fault
We report on a receiver function study of the crust and upper mantle within DESERT, a multidisciplinary geophysical project to study the lithosphere across the Dead Sea Transform (DST). A temporary seismic network was operated on both sides of the DST between 2000 April and 2001 June. The depth of the Moho increases smoothly from about 30 to 34-38 km towards the east across the DST, with significant north-south variations east of the DST. These Moho depth estimates from receiver functions are consistent with results from steep-and wide-angle controlled-source techniques. Steep-angle reflections and receiver functions reveal an additional discontinuity in the lower crust, but only east of the DST. This leads to the conclusion that the internal crustal structure east and west of the DST is different. The P to S converted phases from both discontinuities at 410 and 660 km are delayed by 2 s with respect to the IASP91 global reference model. This would indicate that the transition zone is consistent with the global average, but the upper mantle above 410 km is 3-4 per cent slower than the standard earth model
Seismic tomography, imaging of seismic scatterers, and magnetotelluric soundings reveal a sharp lithologic contrast along a similar to 10 km long segment of the Arava Fault (AF), a prominent fault of the southern Dead Sea Transform (DST) in the Middle East. Low seismic velocities and resistivities occur on its western side and higher values east of it, and the boundary between the two units coincides partly with a seismic scattering image. At 1 - 4 km depth the boundary is offset to the east of the AF surface trace, suggesting that at least two fault strands exist, and that slip occurred on multiple strands throughout the margin's history. A westward fault jump, possibly associated with straightening of a fault bend, explains both our observations and the narrow fault zone observed by others
On 26 December 2004, a moment magnitude M-w = 9.3 earthquake occurred along Northern Sumatra, the Nicobar and Andaman islands, resulting in a devastating tsunami in the Indian Ocean region(1). The rapid and accurate estimation of the rupture length and direction of such tsunami-generating earthquakes is crucial for constraining both tsunami wave- height models as well as the seismic moment of the events. Compressional seismic waves generated at the hypocentre of the Sumatra earthquake arrived after about 12 min at the broadband seismic stations of the German Regional Seismic Network (GRSN)(2,3), located approximately 9,000 km from the event. Here we present a modification of a standard array- seismological approach and show that it is possible to track the propagating rupture front of the Sumatra earthquake over a total rupture length of 1,150 km. We estimate the average rupture speed to be 2.3-2.7 km s(-1) and the total duration of rupture to be at least 430 s, and probably between 480 and 500 s.
The Hengshan complex forms part of the central zone of the North China Craton and consists predominantly of ductilely-deformed late Archaean to Palaeoproterozoic high-grade, partly migmatitic, granitoid orthogneisses, intruded by mafic dykes of gabbroic composition. Many highly strained rocks were previously misinterpreted as supracrustal sequences and represent mylonitized granitoids and sheared dykes. Our single zircon dating documents magmatic granitoid emplacement ages between 2.52 Ga and 2.48 Ga, with rare occurrences of 2.7 Ga gneisses, possibly reflecting an older basement. A few granitic gneisses have emplacement ages between 2.35 and 2.1 Ga and show the same structural features as the older rocks, indicating that the main deformation occurred after similar to 2.1 Ga. Intrusion of gabbroic dykes occurred at similar to 1920 Ma, and all Hengshan rocks underwent granulite-facies metamorphism at 1.88-1.85 Ga, followed by retrogression, shearing and uplift. We interpret the Hengshan and adjacent Fuping granitoid gneisses as the lower, plutonic, part of a late Archaean to early Palaeoproterozoic Japan-type magmatic arc, with the upper, volcanic part represented by the nearby Wutai complex. Components of this arc may have evolved at a continental margin as indicated by the 2.7 Ga zircons. Major deformation and HP metamorphism occurred in the late Palaeoproterozoic during the Luliang orogeny when the Eastern and Western blocks of the North China Craton collided to form the Trans-North China orogen. Shear zones in the Hengshan are interpreted as major lower crustal discontinuities post-dating the peak of HP metamorphism, and we suggest that they formed during orogenic collapse and uplift of the Hengshan complex in the late Palaeoproterozoic (< 1.85 Ga)
Numerical thermodynamic modelling of mineral composition and modes for specified pressure-temperature paths reveals the strong influence of fractional garnet crystallisation, as well as water fractionation, on garnet growth histories in high pressure rocks. Disequilibrium element incorporation in garnet due to the development of chemical inhomogeneities around porphyroblasts leads to pronounced episodic growth and may even cause growth interruptions. Discontinuous growth, together with pressure- and temperature-dependent changes in garnet chemistry, cause zonation patterns that are indicative of different degrees of disequilibrium element incorporation. Chemical inhomogeneities in the matrix surrounding garnet porphyroblasts strongly affect garnet growth and lead to compositional discontinuities and steep compositional gradients in the garnet zonation pattern. Further, intergranular diffusion-controlled calcium incorporation can lead to a characteristic rise in grossular and spessartine contents at lower metamorphic conditions. The observation that garnet zonation patterns diagnostic of large and small fractionation effects coexist within the same sample suggests that garnet growth is often controlled by small-scale variations in the bulk rock chemistry. Therefore, the spatial distribution of garnet grains and their zonation patterns, together with numerical growth models of garnet zonation patterns, yield information about the processes limiting garnet growth. These processes include intercrystalline element transport and dissolution of pre-existing grains. Discontinuities in garnet growth induced by limited element supply can mask traces of the thermobarometric history of the rock. Therefore, thermodynamic modelling that considers fractional disequilibrium crystallisation is required to interpret compositional garnet zonation in terms of a quantitative pressure and temperature path of the host rock
Seaward progradation of several kilometers has been documented mostly for leeward margin low-angle carbonate slope systems with a dominant platform top sediment source. However, steep and high-relief margins fronting deep basins can also prograde and as such are somewhat perplexing. Characteristics of two prograding Carboniferous examples provide a model which may apply elsewhere: (1) outcrops in Asturias, northern Spain serve as important analogs for (2) hydrocarbon reservoirs in steep-sided isolated platforms of the North Caspian Basin, Kazakhstan, such as Tengiz. Seismic and well data from Tengiz corroborate outcrop patterns for slope development, showing progradation of up to 5 and more than 10 km, respectively, despite the high-relief (up to 600 m) and steep (similar to 20-35 degrees) nature of these margins. The two examples share a highly productive microbial boundstone factory extending from the platform break down the slope to nearly 300 m (or more) depth and a lower slope dominated by (mega)breccias and grain flow deposits derived from the margin and slope itself. The broad depth range of microbial boundstone increases the potential for production during both lowstands and high stands of sea level and thereby facilitates progradation independent from platform-top- derived sediment. Rapid in situ lithification of the boundstone provides stability to the steep slopes, but also leads to readjustment through shearing and avalanching. What controls the microbial cement boundstone formation remains a debate but its presence is a key factor controlling the progradational geometry of these and possibly other margins. This new model of "slope" shedding has implications for slope readjustment processes and resulting architecture, sequence stratigraphic interpretation, reservoir characterization, and reservoir modeling. Especially the isotropic character of microbial boundstone will reduce the potential for coherent seismic reflections to develop and possibly invoke, under certain stress regimes, shattering and fracturing thereby generating significant non-matrix permeability. Key considerations are the contrasts with the Bahamian high stand shedding depositional model, slope progradation rates that range from 450 to > 1500 m/My, and net growth rates of in situ boundstone of similar to 1000 m/My, comparable to or higher than accretion rates for metazoan skeletal reef growth. (c) 2005 Published by Elsevier B.V
In this paper we explore the relative control of paleoceanography, eustasy, and water temperature over the evolution of a carbonate slope system deposited on the Marion Plateau (Northeastern Australia). Growth of several carbonate platforms started in the early Miocene on this plateau, and although they occurred in low-latitude subtropical waters they are composed mainly of heterozoan organisms. We investigated an upper to distal slope transect drilled during ODP Leg 194 and located close to the Northern Marion Platform. We reconstructed mass accumulation rates of carbonate as well as the evolution in the ratios of carbon and oxygen stable isotopes. Power spectrum analysis of the carbon isotope record revealed the existence of cycles with main frequencies centered around 409 Kyr and 1800 Kyr. We interpret the 409 Kyr cycle as being paced by changes in the eccentricity of the Earth orbit, and we suggest that the 1800 Kyr cycle could be linked to long-term eustatic changes. Finally, on the basis of the timing of changes in mass accumulation rates of carbonate we infer that the strength and direction of oceanic currents affected sedimentation on the Marion Plateau by shifting depocenters of slope sedimentation, a process probably further modulated by sea-level changes. We argue that the evolution and demise of the heterozoan carbonate systems present on the Marion Plateau were controlled mainly by the evolution of strong benthic currents, and that eustasy and water temperature alone did not account for the drowning of the platforms in the mid Miocene
In a series of timed experiments, monazite inclusions are induced to form in the Durango fluorapatite using 1 and 2 N HCl and H2SO4 solutions at temperatures of 300, 600, and 900 degrees C and pressures of 500 and 1,000 MPa. The monazite inclusions form only in reacted areas, i.e. depleted in (Y+REE)+Si+Na+S+Cl. In the HCl experiments, the reaction front between the reacted and unreacted regions is sharp, whereas in the H2SO4 experiments it ranges from sharp to diffuse. In the 1 N HCl experiments, Ostwald ripening of the monazite inclusions took place both as a function of increased reaction time as well as increased temperature and pressure. Monazite growth was more sluggish in the H2SO4 experiments. Transmission electron microscopic (TEM) investigation of foils cut across the reaction boundary in a fluorapatite from the 1 N HCl experiment (600 degrees C and 500 MPa) indicate that the reacted region along the reaction front is characterized by numerous, sub-parallel, 10-20 nm diameter nano-channels. TEM investigation of foils cut from a reacted region in a fluorapatite from the 1 N H2SO4 experiment at 900 degrees C and 1,000 MPa indicates a pervasive nano- porosity, with the monazite inclusions being in direct contact with the surrounding fluorapatite. For either set of experiments, reacted areas in the fluorapatite are interpreted as replacement reactions, which proceed via a moving interface or reaction front associated with what is essentially a simultaneous dissolution-reprecipitation process. The formation of a micro- and nano-porosity in the metasomatised regions of the fluorapatite allows fluids to permeate the reacted areas. This permits rapid mass transfer in the form of fluid-aided diffusion of cations to and from the growing monazite inclusions. Nano-channels and nano-pores also serve as sites for nucleation and the subsequent growth of the monazite inclusions
Intramontane basins may act as important sediment storage areas, serve as recorders of the history of deformation, record syntectonic deposition, and document the evolution of climatic conditions during deposition. We document the timing, cyclicity, and processes that led to the filling and reexcavation of the intramontane Quebrada del Toro basin in NW Argentina. Geomorphic and geologic observations indicate that the basin was filled with sediment that has been subsequently excavated at least two times in the last similar to 8 m.y. The last filling and excavation cycle occurred within the last 0.98 m.y. and has led to the deposition and removal of similar to 61.4 km(3) of material from the basin, leading to a basin-wide averaged minimum denudation rate of 0.16 mm/yr. Aggradation within the basin took place due to channel steepening of the downstream fluvial system that connects the intramontane basin to the foreland. This portion of the fluvial system is actively incising through an uplifting bedrock zone. We use observations within the Toro to test a quasiphysically based model of channel aggradation behind a rising base level that rises due to downstream channel steepening. Our work shows that the bedrock incision rate constant required to reproduce conditions observed within the Toro basin is consistent with values measured independently in similar rock types. Therefore, in intramontane basins that experience similar processes of filling and evacuation, this model may be used to assess the relative importance of tectonic rock uplift, bedrock resistance to fluvial incision, and climate in determining the geomorphic and sedimentologic history of these basins
Rhodoliths (free-living coralline red algae) can thrive under a wide range of temperatures, reduced light, and increased nutrient levels, and often form a distinct so-called rhodalgal lithofacies that is an important component of Cenozoic shallow-water carbonates. Global distributions illustrate that from the late-early to early-late Miocene (Burdigalian-early Tortonian), rhodalgal facies reached peak abundances and commonly replaced coral-reef environments, accompanied by a decline in other carbonate-producing phototrophs. We argue that the dominance of red algae over coral reefs was triggered in the Burdigalian by enhanced trophic resources associated with a global increase in productivity, as evidenced by a long-term shift toward higher carbon isotope values. Rhodalgal lithofacies expanded further in the middle Miocene when strengthened thermal gradients associated with the establishment of the East Antarctic lee Sheet led to enhanced upwelling while climate change generated increased weathering rates, introducing land-derived nutrients into the oceans. Globally cooler temperatures following a climatic optimum in the early-middle Miocene contributed to sustain the dominance of red algae and prevented the recovery of coral reefs. The global shift in nearshore shallow-water carbonate producers to groups tolerant of higher levels of trophic resources provides further evidence for increased nutrient levels during that time interval and shows the sensitivity of shallow-water carbonate facies as indicators of past oceanographic conditions
"Green water"-the water stored in the soil and productively used for plant transpiration-is an important quantity particularly in rainfed agriculture (in contrast to "blue water" available in streams and lakes, on which irrigation relies). This study provides preliminary estimates of contemporary (1961-1990) global green water flows (i.e. plant transpiration), using a well-established, process-based dynamic global vegetation and water balance model, LPJ. Transpiration is analysed with respect to differences between a simulation that accounts for human land cover changes and a simulation under conditions of potential natural vegetation. We found that historic land cover change usually reduced the green water flow to the atmosphere, resulting in a global decrease of similar to 7% in total. To further explore how the biophysical setting influences the green water flow, we analyse the ratio between soil moisture-limited canopy conductance of carbon and water and energy-controlled potential conductance. This plant physiology-based ratio measures the degree to which actual green water flow falls below the potential flow that would occur when the soil is saturated, thus it represents a measure of the water limitation of terrestrial vegetation. We found that plant water limitation is lowest in the wet tropics and at high latitudes, where soil moisture is high enough to meet the atmospheric demand for transpiration. The present results are preliminary, since irrigation is not yet accounted for, and because the model simulations are compromised primarily by the quality of the input data. A more comprehensive and consistent assessment of global green and blue water flows and limitations using an enhanced LPJ model is identified as a prime task for future studies. (c) 2005 Elsevier Ltd. All rights reserved
The uranium deposit at Niederschlema-Alberoda, Germany, contains a rich variety of Bi minerals deposited between the Permian and the Cretaceous; these have been studied for paragenetic relations, composition, and conditions of formation. Particular attention is given to the rare Bi selenides watkinsonite, nevskite, and cuproan bohdanowiczite. Whereas watkinsonite and nevskite only occur intergrown with clausthalite, bohdanowiczite is more widespread and also is associated with Cu selenides. Watkinsonite from this second confirmed locality worldwide has an average composition (Cu1.47Ag0.49)(Sigma 1.96)(Pb1.01Hg0.01 Fe-0.01)(Sigma 1.03)Bi-3.98(Se7.98S0.05)(Sigma 8.03), ideally (Cu,Ag)(2)PbBi4Se8. These findings suggest that the empirical formula of watkinsonite originally proposed for the type specimen from the Otish Mountains uranium deposit in Quebec [CU2+xPb1+xBi4-xSe,S,Te)(8), x approximate to 0.3] requires revision. The composition of nevskite is (Pb0.06Bi0.95)(Sigma 1.01)Se-0.99, on average. Bohdanowiczite from the Cu- selenide assemblage shows extensive substitution of Cu+ for Ag+, expressed by the crystallochemical formula (Ag1.80- 0.94CU0.16-1.05Pb0.00-0.05)(Sigma 1.97-2.07)BiSigma 1.97-2.03SeSigma 3.96-4.04. This observation seems to argue for the natural existence of CU2Bi2Se4, the Se-dominant analogue of emplectite. The Bi selenides were deposited at temperatures of about 100 degrees C, in the Jurassic. The lack of thermodynamic data for all the Bi selenides limits reliable inferences on the fugacities of selenium and sulfur that prevailed during their formation. Other Bi minerals from this locality comprise members of the bismuthinite-aikinite solid-solution series of Permian age and, more importantly, native Bi and Bi sulfides (matildite, bismuthinite, wittichinite), deposited in the Cretaceous
Potential fields are classically represented on the sphere using spherical harmonics. However, this decomposition leads to numerical difficulties when data to be modelled are irregularly distributed or cover a regional zone. To overcome this drawback, we develop a new representation of the magnetic and the gravity fields based on wavelet frames. In this paper, we first describe how to build wavelet frames on the sphere. The chosen frames are based on the Poisson multipole wavelets, which are of special interest for geophysical modelling, since their scaling parameter is linked to the multipole depth (Holschneider et al.). The implementation of wavelet frames results from a discretization of the continuous wavelet transform in space and scale. We also build different frames using two kinds of spherical meshes and various scale sequences. We then validate the mathematical method through simple fits of scalar functions on the sphere, named 'scalar models'. Moreover, we propose magnetic and gravity models, referred to as 'vectorial models', taking into account geophysical constraints. We then discuss the representation of the Earth's magnetic and gravity fields from data regularly or irregularly distributed. Comparisons of the obtained wavelet models with the initial spherical harmonic models point out the advantages of wavelet modelling when the used magnetic or gravity data are sparsely distributed or cover just a very local zone
Basin formation dynamics of the Tertiary Piedmont Basin (TPB) are here investigated by means of cross-section numerical modelling. Previous works hypothesised that basin subsidence occurred due first to extension (Oligocene) and then to subsequent loading due to back-thrusting (Miocene). However, structural evidence shows that the TPB was mainly under contraction from Oligocene until post Pliocene time while extension played a minor role. Furthermore, thermal indicators strongly call for a cold (flexure-induced) mechanism but are strictly inconsistent with a hot (thermally induced) mechanism. Our new modelling shows that the TPB stratigraphic features can be reproduced by flexure of a visco- elastic plate loaded by backthrusts active in the Western Alps in Oligo-Miocene times. Far-field compression contributed to the TPB subsidence and controlled the basin infill geometry by enhancing basin tilting, forebulge uplift and erosion of the southern margin of the basin. These results suggest that the TPB subsidence is the result of a combination of mechanisms including thrust loading and farfield compressional stresses. (c) 2005 Elsevier B.V. All rights reserved
Carpholite+chloritoid+pyrophyllite association occurs widely in the Triassic metaclastic rocks of the Afyon Zone in west-central Turkey. Fe-Mg-carpholite is associated with rare aragonite pseudomorphs and glaucophane in marbles and metabasites, respectively. The Afyon Zone consists stratigraphically of a Pan-African basement and an overlying Mesozoic cover sequence. The Pan-African basement, which shows Barrovian-type amphibolite-facies metamorphism, comprises garnet-mica schists, intruded by sodic amphibole-bearing metagabbros and leucocratic metagranites. It is unconformably overlain by a continuous metasedimentary sequence extending from Triassic to early Palaeocene. This cover sequence begins with metaconglomerates, which pass upwards into phyllites. Fe-Mg-carpholite occurs within this metaclastic sequence as rosette-like crystals in metapelites and fibres in quartz segregations. The metaclastic rocks are succeeded by metamorphosed platform carbonates, grading into Latest Mesozoic metamorphosed pelagic limestones, which in turn progress up to a Late Mesozoic-Early Tertiary olistostrome. This sequence is tectonically overlain by the blueschists of the Tavsanh Zone. Fe-Mg-carpholite-bearing assemblages imply temperatures of about 350 degrees C and minimum pressures of 6-9 kbar, corresponding to burial depths of about 30 km for the Mesozoic passive continental margin sediments and the underlying Pan-African supracrustal metasediments and metaintrusives. The metamorphic rocks of the Afyon Zone are unconformably overlain by Upper Palaeocene-Lower Eocene sedimentary rocks, indicating a Paleocene age for the regional HP/LT metamorphism. This implies continuous younging of HP/LT metamorphism in the Anatolides related to northward subduction of the Anatolide-Tauride platform beneath the Sakarya Zone. From north to south this involved the Tavsanh Zone (Campanian, 80 +/- 5 Ma), the Afyon Zone (Palaeocene?), the Menderes Massif (Middle Eocene) and the Lycian Nappes (Late Cretaceous-Eocene?), all of which were probably derived from the frontal part of the Anatolide-Tauride platform. (c) 2005 Elsevier B.V. All rights reserved
[1] The Puna-Altiplano plateau in South America is a high-elevation, low internal relief landform that is characterized by internal drainage and hyperaridity. Thermochronologic and sedimentologic observations from the Sierra de Calalaste region in the southwestern Puna plateau, Argentina, place new constraints on early plateau evolution by resolving the timing of uplift of mountain ranges that bound present-day basins and the filling pattern of these basins during late Eocene-Miocene time. Paleocurrent indicators, sedimentary provenance analyses, and apatite fission track thermochronology indicate that the original paleodrainage setting was disrupted by exhumation and uplift of the Sierra de Calalaste range between 24 and 29 Ma. This event was responsible for basin reorganization and the disruption of the regional fluvial system that has ultimately led to the formation of internal drainage conditions, which, in the Salar de Antofalla, were established not later than late Miocene. Upper Eocene-Oligocene sedimentary rocks flanking the range contain features that suggest an arid environment existed prior to and during its uplift. Provenance data indicate a common similar source located to the west for both the southern Puna and the Altiplano of Bolivia during the late Eocene- Oligocene with sporadic local sources. This suggests the existence of an extensive, longitudinally oriented foreland basin along the central Andes during this time
An analysis of the shear (S) waves recorded during the wide-angle reflection/refraction (WRR) experiment as part of the DESERT project crossing the Dead Sea Transform (DST) reveals average crustal S-wave velocities of 3.3-3.5 km s(-1) beneath the WRR profile. Together with average crustal P-wave velocities of 5.8-6.1 km s(-1) from an already published study this provides average crustal Poisson's ratios of 0.26-0.27 (V-p/V-s = 1.76-1.78) below the profile. The top two layers consisting predominantly of sedimentary rocks have S- wave velocities of 1.8-2.7kms(-1) and Poisson's ratios of 0.25-0.31 (V-p/V-s = 1.73-1.91). Beneath these two layers the seismic basement has average S- wave velocities of around 3.6 km s(-1) east of the DST and about 3.7 km s(-1) west of the DST and Poisson's ratios of 0.24-0.25 (V-p/V-s = 1.71-1.73). The lower crust has an average S-wave velocity of about 3.75 km s(-1) and an average Poisson's ratio of around 0.27 (V-p/V-s = 1.78). No Sn phase refracted through the uppermost mantle was observed. The results provide for the first time information from controlled source data on the crustal S-wave velocity structure for the region west of the DST in Israel and Palestine and agree with earlier results for the region east of the DST in the Jordanian highlands. A shear wave splitting study using SKS waves has found evidence for crustal anisotropy beneath the WRR profile while a receiver function study has found evidence for a lower crustal, high S-wave velocity layer east of the DST below the profile. Although no evidence was found in the S-wave data for either feature, the S-wave data are not incompatible with crustal anisotropy being present as the WRR profile only lies 30 degrees off the proposed symmetry axis of the anisotropy where the difference in the two S-wave velocities is still very small. In the case of the lower crustal, high S-wave velocity layer, if the velocity change at the top of this layer comprises a small first-order discontinuity underlain by a 2 km thick transition zone, instead of just a large first-order discontinuity, then both the receiver function data and the WRR data presented here can be satisfied. Finally, the S-wave velocities and Poisson's ratios which have been derived in this study are typical of continental crust and do not require extensional processes to explain them
The Sierra del Cuera (Cantabrian Mountains, Asturias, Spain) exhibits an exceptionally well-preserved upper Bashkirian-lower Moscovian succession of a high-relief carbonate platform with steep adjacent slope. During the late Bashkirian, the platform margin prograded rapidly basinward, and during the early Moscovian, platform aggradation increased. Fusulinida distributions have been estimated in three platform study windows representative of the Bashkirian and Moscovian lithofacies. Microbial boundstone accumulations at the platform margin and slope contain scarce fusulinoideans (particularly during the Moscovian aggradation, when water depth was relatively greater). Boundstones are characterized by tetrataxids, lasiodiscids, pseudoammodiscids, biseriamminids, endothyrids and palaeotextulariids. Grainstone units with sediments deposited above wave base are characterized by the predominance of Pseudoendothyra, Eostaffella (Bashkirian) or Profusulinella (Moscovian), Bradyina and palaeotextulariids. The open-marine facies, below effective wave-base but intermittently affected by storm currents, is enriched in endothyrids and palaeotextulariids in the Bashkirian part, and in Profusulinella, Eofusulina and Schubertella in the Moscovian part. The Moscovian facies, of lower energy and deeper water depth, contain some microbially precipitated micrite and have scarce Fusulinida, mostly represented by endothyrids, biseriamminids, and fusulinoideans (Eostaffella, Ozawainella and Profusulinella). Lagoonal open- to restricted-marine facies with beresellid algae contain Profusulinella, Schubertella and biseriamminids. Crinoidal packstones were deposited in moderate-energy settings close to the platform margin, below wave-base, and are characterized by Ozawainella, Eostaffella, palaeotextulariids and tetrataxids. These analyses show that fusulinid foraminifera inhabited mostly shallow and open-marine environments. Their distribution was controlled by depth-related variables, water energy and open-marine vs. restricted-marine conditions. Fusulinoideans increase in the Moscovian strata, and show different trends in the upper Bashkirian and lower Moscovian platforms, possibly because of: a) changes in the paleoecology of Fusulinida at the species level from the Bashkirian to the Moscovian; and/or b) variations in the depositional environments and physicochemical characteristics of seawater from the Bashkirian progradational phase to the Moscovian aggradational phase. Pseudoendothyra, Profusulinella, Schubertella, bradyinids and palacotextulariids showed an affinity for high-energy settings. Profusulinella and Schubertella, however, preferred the shallowest facies and also tolerated restricted marine conditions. Eostaffella's distribution largely differs between the Bashkirian and Moscovian parts because it characterizes the high-energy grainstones of the upper Bashkirian and the open-marine, moderate- to low-energy environments of the lower Moscovian. Endothyrids were excluded from lagoonal environments with restricted circulation and abnormal salinity. Biscriamminids preferred open-marine, moderate- to low-energy settings with muddy substrates, but were not li
The effect of water activity on the oxidation and structural state of Fe in a ferro-basaltic melt
(2005)
Experimental investigations have been performed at T = 1200 degrees C, P = 200 MPa and fH(2) corresponding to H2O-MnO-Mn3O4 and H2O-QFM redox buffers to study the effect of H2O activity on the oxidation and structural state of Fe in an iron-rich basaltic melt. The analysis of Mossbauer and Fe K-edge X-ray absorption nearedge structure (XANES) spectra of the quenched hydrous ferrobasaltic glasses shows that the Fe3+/Sigma Fe ratio of the glass is directly related to aH(2)O in a H-2-buffered system and, consequently, to the prevailing oxygen fugacity (through the reaction of water dissociation H2O <-> H-2 + 1/2 O-2). However, water as a chemical component of the silicate melt has an indistinguishable effect on the redox state of iron at studied conditions. The experimentally obtained relationship between fO(2) and Fe3+/Fe2+ in the hydrous ferrobasaltic melt can be adequately predicted in the investigated range by the existing empiric and thermodynamic models. The ratio of ferric and ferrous Fe is proportional to the oxygen fugacity to the power of similar to 0.25 which agrees with the theoretical value from the stoichiometry of the Fe redox reaction (FeO + 1/4 O-2 = FeO1.5). The mean centre shifts for Fe2+ and Fe3+ absorption doublets in Mossbauer spectra show little change with increasing Fe3+/Sigma Fe, suggesting no significant change in the type of iron coordination. Similarly, XANES preedge spectra indicate a mixed (C3h, Td, and Oh, i.e., 5-, 4-, and sixfold) coordination of Fe in hydrous basaltic glasses. Copyright (c) 2005 Elsevier Ltd
Late Quaternary intensified monsoon phases control landscape evolution in the northwest Himalaya
(2005)
The intensity of the Asian summer-monsoon circulation varies over decadal to millennial time scales and is reflected in changes in surface processes, terrestrial environments, and marine sediment records. However, the mechanisms of long-lived (2-5 k.y.) intensified monsoon phases, the related changes in precipitation distribution, and their effect on landscape evolution and sedimentation rates are not yet well understood. The and high-elevation sectors of the orogen correspond to a climatically sensitive zone that currently receives rain only during abnormal (i.e., strengthened) monsoon seasons. Analogous to present-day rainfall anomalies, enhanced precipitation during an intensified monsoon phase is expected to have penetrated far into these geomorphic threshold regions where hillslopes are close to the angle of failure. We associate landslide triggering during intensified monsoon phases with enhanced precipitation, discharge, and sediment flux leading to an increase in pore-water pressure, lateral scouring of rivers, and over- steepening of hillslopes, eventually resulting in failure of slopes and exceptionally large mass movements. Here we use lacustrine deposits related to spatially and temporally clustered large landslides (>0.5 km(3)) in the Sutlej Valley region of the northwest Himalaya to calculate sedimentation rates and to infer rainfall patterns during late Pleistocene (29-24 ka) and Holocene (10-4 ka) intensified monsoon phases. Compared to present-day sediment-flux measurements, a fivefold increase in sediment-transport rates recorded by sediments in landslide-dammed lakes characterized these episodes of high climatic variability. These changes thus emphasize the pronounced imprint of millennial-scale climate change on surface processes and landscape evolution
The interplay between topography and Indian summer monsoon circulation profoundly controls precipitation distribution, sediment transport, and river discharge along the Southern Himalayan Mountain Front (SHF). The Higher Himalayas form a major orographic barrier that separates humid sectors to the south and and regions to the north. During the Indian summer monsoon, vortices transport moisture from the Bay of Bengal, swirl along the SHF to the northwest, and cause heavy rainfall when colliding with the mountain front. In the eastern and central parts of the Himalaya, precipitation measurements derived from passive microwave analysis (SSM/I) show a strong gradient, with high values at medium elevations and extensive penetration of moisture along major river valleys into the orogen. The end of the monsoonal conveyer belt is near the Sutlej Valley in the NW Himalaya, where precipitation is lower and rainfall maxima move to lower elevations. This region thus comprises a climatic transition zone that is very sensitive to changes in Indian summer monsoon strength. To constrain magnitude, temporal, and spatial distribution of precipitation, we analyzed high-resolution passive microwave data from the last decade and identified an abnormal monsoon year (AMY) in 2002. During the 2002 AMY, violent rainstorms conquered orographic barriers and penetrated far into otherwise and regions in the northwest Himalaya at elevations in excess of 3 km asl. While precipitation in these regions was significantly increased and triggered extensive erosional processes (i.e., debris flows) on sparsely vegetated, steep hillslopes, mean rainfall along the low to medium elevations was not significantly greater in magnitude. This shift may thus play an important role in the overall sediment flux toward the Himalayan foreland. Using extended precipitation and sediment flux records for the last century, we show that these events have a decadal recurrence interval during the present-day monsoon circulation. Hence, episodically occurring AMYs control geomorphic processes primarily in the high-elevation and sectors of the orogen, while annual recurring monsoonal rainfall distribution dominates erosion in the low- to medium- elevation parts along the SHF. (C) 2004 Elsevier B.V. All rights reserved
The up to similar to4 km high southern Patagonian Andes form a pronounced topographic barrier to atmospheric circulation in the southern hemisphere westerlies, and cause one of the most drastic orographic rain shadows on earth. Geologic data imply that this climatic pattern has been established or significantly enhanced during Miocene surface uplift of this Andean segment. We report evidence for important climatic and ecologic changes in the eastern foreland of the Patagonian Andes that appear to be the result of this uplift. To provide constraints on Miocene plant ecosystems and precipitation in the eastern (leeward) foreland of the Patagonian Andes, we determined carbon and oxygen isotope values of pedogenic carbonate nodules from a similar to500 m thick section of the continental Santa Cruz Formation. The age of these deposits was constrained by Ar/Ar dating of intercalated tuffs, which range from similar to22 to 14 Ma. At similar to16.5 Ma, the delta(13)C values increase by similar to3parts per thousand, the delta(18)O values decrease by >2parts per thousand, and the scatter in the oxygen isotope data increases significantly. We interpret these changes as the consequence of >1 km surface uplift in this Andean segment (from the delta(18)O values), and increased aridity to its east (from the delta(13)C values and the increased scatter in the delta(18)O values). Sediments overlying the Santa Cruz Formation are very limited in extent and volume, and dominated by coarse conglomerates related to Pleistocene and older glaciations. It thus seems that, by similar to14 Ma, deposition in the eastern foreland of the Southern Patagonian Andes had essentially ceased as the result of rain shadow formation. (C) 2004 Elsevier B.V. All rights reserved
Composite ground-motion models and logic trees: Methodology, sensitivities, and uncertainties
(2005)
Logic trees have become a popular tool in seismic hazard studies. Commonly, the models corresponding to the end branches of the complete logic tree in a probabalistic seismic hazard analysis (PSHA) are treated separately until the final calculation of the set of hazard curves. This comes at the price that information regarding sensitivities and uncertainties in the ground-motion sections of the logic tree are only obtainable after disaggregation. Furthermore, from this end-branch model perspective even the designers of the logic tree cannot directly tell what ground-motion scenarios most likely would result from their logic trees for a given earthquake at a particular distance, nor how uncertain these scenarios might be or how they would be affected by the choices of the hazard analyst. On the other hand, all this information is already implicitly present in the logic tree. Therefore, with the ground-motion perspective that we propose in the present article, we treat the ground-motion sections of a complete logic tree for seismic hazard as a single composite model representing the complete state-of-knowledge-and-belief of a particular analyst on ground motion in a particular target region. We implement this view by resampling the ground-motion models represented in the ground-motion sections of the logic tree by Monte Carlo simulation (separately for the median values and the sigma values) and then recombining the sets of simulated values in proportion to their logic-tree branch weights. The quantiles of this resampled composite model provide the hazard analyst and the decision maker with a simple, clear, and quantitative representation of the overall physical meaning of the ground-motion section of a logic tree and the accompanying epistemic uncertainty. Quantiles of the composite model also provide an easy way to analyze the sensitivities and uncertainties related to a given logic-tree model. We illustrate this for a composite ground- motion model for central Europe. Further potential fields of applications are seen wherever individual best estimates of ground motion have to be derived from a set of candidate models, for example, for hazard rnaps, sensitivity studies, or for modeling scenario earthquakes
This study concerns the Quantitative Phase Analysis (QPA) of historical bricks coming from the complex of the Great Palace of the Byzantine Emperors in Istanbul. The studied samples are characterised by different chemical compositions (low and high calcium content), variable firing temperatures and different amounts of soluble salts as damage products. In the low-Ca samples, the decrease of the phyllosilicate content (from 23.4 to 6.9 wt%) is associated to the increase of the amorphous fraction (from 24 to 48%). This clear negative correlation between the phyllosilicate content and the amorphous fraction indicates that in low-Ca systems vitrification processes are overwhelming with respect to nucleation and recrystallisation processes. By contrast, high-Ca samples present newly formed Ca(Mg) silicates (diopside from 5.7 to 27.2%; anorthite from 1.4 to 8.7%) and aluminium silicates (gehlenite only in two samples, 6.2 and 7.7%) associated to the decrease of quartz (from 27.7 to 11.5%), phyllosilicate (from 6.5% until complete break down) and amorphous (from 30 to 14%) phase fractions. These findings support the role played by the CaO(MgO) content deriving from carbonates decomposition which reacts with Al2O3 and SiO2 oxides from dehydroxylated clay minerals and quartz grains. The above results have been obtained by X-ray powder diffraction data using the combined Rietveld refinement - internal standard method in order to estimate both the crystalline and the amorphous phase fractions. In addition, the coexistence of two distinct plagioclases in high-Ca samples was modelled as follows: a primary albite, which tends to incorporate Ca during the firing process as demonstrated by the increasing of gamma crystallographic angle, and a newly formed anorthite. Finally, by difference between the X-ray fluorescence data and the chemical compositions inferred by QPA, it proved possible to roughly estimate the residual chemical composition attributable to the amorphous fraction. On the basis of our data, we believe that Rietveld refinement combined with the internal standard method represent a powerful tool to better characterise complex polycrystalline and amorphous mixture as in the case of historical bricks
The PEGASOS project was a major international seismic hazard study, one of the largest ever conducted anywhere in the world, to assess seismic hazard at four nuclear power plant sites in Switzerland. Before the report of this project has become publicly available, a paper attacking both methodology and results has appeared. Since the general scientific readership may have difficulty in assessing this attack in the absence of the report being attacked, we supply a response in the present paper. The bulk of the attack, besides some misconceived arguments about the role of uncertainties in seismic hazard analysis, is carried by some exercises that purport to be validation exercises. In practice, they are no such thing; they are merely independent sets of hazard calculations based on varying assumptions and procedures, often rather questionable, which come up with various different answers which have no particular significance. (C) 2005 Elsevier B.V. All rights reserved
Microseismic data from observatories in Europe, which have been continuously recorded since about 100 years, contain information on the wave-climate in the North Atlantic. They can potentially be used as additional constraints in high-resolution temporal and spatial reconstructions of the storminess and oceanic waveheights in the past. To resolve spatial patterns data from observatories in different regions are needed. While previous recent studies analyzed only few observatory archives and relatively short time ranges, this is a first attempt to process the whole available data archive from different observatories. We correct and compare smoothed microseismic data from different stations and discuss their correlation and possible use for studies of storminess variability. Microseismic amplitudes at four seismic stations in northern Europe show amplitude peaks in 1920 and 1925, a slow decline in amplitudes till the middle of the 1930's followed by a steady increase of amplitudes till about 1990. From 1990 on microseismic amplitudes decrease. We find a good correlation between the average surface wind velocity in the North Atlantic and microseismic amplitudes at inland stations far away from the coast. Coastal stations are more influenced by local swell and are thus potentially useful to recover regional changes in wind and ocean wavefields with time. The study demonstrates that the analysis of microseismic has the potential to assess climate changes during the last 100 years
Annite and Fe-rich siderophyllite constitute the rock-forming micas in the late-Variscan composite granite pluton of Konigshain, Lausitz, Germany. This multiphase pluton is composed of three fractionated, but not chemically specialized monzogranite types, which contain lithophile elements such as Li, Rb, Cs, Sn, and F in average quantities. Abundant miarolitic pegmatites of the NYF family with a broad diversity of rare minerals occur in the apical part of the pluton. These pegmatitic cavities locally contain di- and trioctabedral micas as well as cation-deficient micas. Trioctahedral micas comprise F-rich manganoan lithian siderophyllite to manganoan zinnwaldite, zinnwaldite, and minor lepidolite. The formula [calculated on the basis of 22 anion valencies and 2 (F + OH + Cl)] of the most Mn-rich siderophyllite is (K0.85Rb0.08Na0.04)(0.97)(Al0.99Li0.91Fe0.51Mn0.42Ti0.01Zn0.01)(2.85) (Si3.21Al0.79)(4)O- 10(F1.80OH0.19Cl0.01)(2). This mica constitutes one of the most Mn-rich siderophyllite compositions reported to date. The lithium micas poorer in Mn are distinguished by elevated concentrations of Rb (up to 2.5 wt % Rb2O), CS (UP to 1.2 wt % Cs2O), and F (up to 9.6 wt %). This fluorine content is probably consistent with the maximum possible F occupation of 2 of the (F,OH,Cl)-site. The structural formula of the most Li-rich lepidolite is (K0.83Rb0.07Cs0.03)(0.93) (Li1.62Al1.00Fe0.38)(3.00)(Si3.62Al0.38)(4) O-10(F1.91OH0.09)(2). During hydrothermal alteration, lepidolite and zinnwaldite became partially depleted in K, Li, Rb, Cs, and F and gradually transformed into cation-deficient micas (lithian phengite to illite of phengitic affinity)
The "Bohemians" of New Zealand - an ethnic group? In 1982 a small group of New Zealanders established contacts with the region of origin of their about 200 German-speaking ancestors who had emigrated from Bohemia for economic reasons in the 1860s and 1870s. They had all come from about twenty villages situated west of plzen and founded a rural settlement and participated in the foundation of a second one in New Zealand. Since World War I there had been no further contacts between the emigrants and their descendants on the one hand and their relatives in Bohemia on the other hand. For two reasons new contacts were established after such a long time: (1) the back-to the-roots-movement had spread to New Zealand from the USA, Canada and Australia, and (2) the status of cultural diversity keeps being enhanced in New Zealand since about 1970. These processes also influenced those people in New Zealand who call themselves "Bohemians" because of their ancestors' region of origin. Their total number is estimated at 10, 000 to 15, 000 at present. Up to now hardly any attention was attached to them in New Zealand by academic research and/or the general public. This paper discusses the history and today's situation of the former immigrants' community as well as the New Zealand "Bohemians " in general, raising the questions to what extent they can be defined as an ethnic group now and whether they will retain their status as a specific group in future
Cretaceous magmatism in the Eastern Cordillera of Colombia is related to lithospheric stretching during the late Early-early Late Cretaceous. The small amount of preserved igneous material is represented by small mafic intrusions. This study focuses on three localities, from east to west: Pajarito, Pacho, and Caceres. The investigated igneous bodies are classified as gabbros, pyroxene-hornblende-gabbros, and pyroxene-hornblendites mainly composed of plagioclase, clinopyroxene, and/or amphibole. Although their timing of emplacement and geodynamic position seem similar, significant differences in their geochemical and petrological characteristics rule out simple models of melt genesis. Clinopyroxene and bulk chemistry indicate increasing alkalinity from west to east. Trace element concentrations point to melt sources that range from a slightly enriched mantle in the west to a highly enriched one in the east. In addition, the data reflect a decreasing degree of partial melting from west to east and the decreasing importance of residual garnet in the mantle source. Probable mantle metasomatism in the source region by slab-derived fluids, as displayed by high Ba/Nb and moderate Sr-n/P-n, is clear in the west and very slight to the east. Mantle metasomatism and melt generation probably are processes of different epochs. The lack of large volumes of igneous rocks and the absence of tectonically controlled unconformities in the investigated areas indicate that a mantle plume did not affect the regional tectonics and magmatism. We favor a model of rift-related magmatism in which melt composition is modified from east to west from a highly enriched to a less enriched mantle region, the latter metasomatized by fluids derived from an older subduction phase. (c) 2005 Elsevier Ltd. All rights reserved
Early Ilerdian (Early Eocene, Shallow Benthic Zones 5 and 6) carbonate systems of the Pyrenees shelf were deposited after a time of severe climatic ('Paleocene-Eocene Thermal Maximum, PETM') and phylogenetic ('Larger Foraminifer Turnover') changes. They reflect the radiation of nummulitid, alveolinid, and orbitolitid larger foraminifera after remarkable biotic changes at the end of the Paleocene, and announce their subsequent flourishing in the Middle Eocene. A paleoenvironmental model for tropical carbonate environments of this particular time interval is provided herein. During the Early Ilerdian, the inner and middle ramp deposits from Minerve, Campo and Serraduy revealed the end-member of a tropical carbonate factory with carbonate production dominated by the end-members of biotically (photo-autotrophic skeletal) controlled and biotically induced carbonate precipitation. Inner platform environments are dominated by alveolinids and in part by orbitolitids, middle platform environments are dominated by nummulitids. Corals are present, but they do not form reefs, which is a typical feature for the Eocene. Nummulite shoal complexes, which are well-known from the Middle Eocene are also absent during the studied Early Ilerdian interval, which may reflect the early evolutionary stage of this group
In silicate glasses and melts, water acts according to two main processes. First, it can be dissolved in high temperature/high pressure melts. Second, it constitutes a weathering agent on the glass surface. A number of in-situ x- ray absorption fine structure (XAFS) studies for Fe, Ni, Zr, Th and U show that the more charged cations (Zr, Nb, Mo, Ta, Sn, Th and U) are little affected by the presence of dissolved water in the melt. In contrast, divalent iron and nickel are highly sensitive to the presence of water, which enhance nucleation processes, for example, of phyllosilicates at the angstrom-scale. Such information provides additional constraints on the role of water deep in the Earth, particularly in magmatology. By contrast, the weathering of glass surfaces by water can be studied from a durability perspective. Experimental weathering experiments Of nuclear waste glasses performed in the laboratory show a variety of surface enrichments (carbon, chlorine, alkalis, iron) after exposure to atmospheric fluids and moisture. Mn-, and Fe-surface enrichments of analogous glasses of the XIVth century are related to the formation of Mn and Fe oxy/ hydroxides on the surface. The impact on the glass darkening is considered in terms of urban pollution and mass tourism
Fe in magma : an overview
(2005)
The strong influence of physical conditions during magma formation on Fe equilibria offers a large variety of possibilities to deduce these conditions from Fe-bearing phases and phase assemblages found in magmatic rocks. Conditions of magma genesis and their evolution are of major interest for the understanding of volcanic eruptions. A brief overview on the most common methods used is given together with potential problems and limitations. Fe equilibria are not only sensitive to changes in intensive parameters (especially T and fO(2)) and extensive parameters like composition also have major effects, so that direct application of experimentally calibrated equilibria to natural systems is not always possible. Best estimates for pre-eruptive conditions are certainly achieved by studies that relate field observations directly to experimental observations for the composition of interest using as many constraints as possible (phase stability relations, Fe-Ti oxides, Fe partitioning between phases, Fe oxidation state in glass etc.). Local structural environment of Fe in silicate melts is an important parameter that is needed to understand the relationship between melt transport properties and melt structure. Assignment of Fe co-ordination and its relationship to the oxidation state seems not to be straightforward. In addition, there is considerable evidence that the co- ordination of Fe in glass differs from that in the melt, which has to be taken into account when linking melt structure to physical properties of silicate melts at T and P.
Logic trees are widely used in probabilistic seismic hazard analysis as a tool to capture the epistemic uncertainty associated with the seismogenic sources and the ground-motion prediction models used in estimating the hazard. Combining two or more ground-motion relations within a logic tree will generally require several conversions to be made, because there are several definitions available for both the predicted ground-motion parameters and the explanatory parameters within the predictive ground-motion relations. Procedures for making conversions for each of these factors are presented, using a suite of predictive equations in current use for illustration. The sensitivity of the resulting ground-motion models to these conversions is shown to be pronounced for some of the parameters, especially the measure of source-to-site distance, highlighting the need to take into account any incompatibilities among the selected equations. Procedures are also presented for assigning weights to the branches in the ground-motion section of the logic tree in a transparent fashion, considering both intrinsic merits of the individual equations and their degree of applicability to the particular application
We test the capability of broadband arrays at teleseismic distances to image the spatio-temporal characteristics of the seismic energy release during the Dec 26, 2004 Sumatra earthquake at early observation times. Using a non-plane-wave array location technique previously reported values for rupture length (about 1150 km), duration (about 480 s), and average rupture velocity (2.4-2.7 km/s) are confirmed. Three dominant energy releases are identified: one near the hypocenter, a second at 6 degrees N94 degrees E about 130 s later and a third one after 300 s at 9 degrees N92-93 degrees E. The spatio-temporal distribution of the radiated seismic energy in the source region is calculated from the stacked broadband recordings of two arrays in Germany and Japan and results in rough estimates of the total seismic energy of 0.55.10(18) Nm (GRSN) and 1.53.10(18) Nm (FNET) respectively. Changes in the relative ratio of energy as function of spatio-temporal location indicate a rotation of the focal mechanism during the rupture process
The structure and alterations of subducted oceanic lithosphere ( e. g., thickness and seismic velocity of oceanic crust) can be obtained by analyzing guided seismic waves generated by earthquakes within the slab (Wadati- Benioff zone). In northern Chile prominent secondary phases from intermediate-depth seismicity, observed in the forearc region can be interpreted as guided waves. For the observation of guided waves it is usually required to have stations close to the wave guide, a fact which is not directly given for forearc stations in subduction zone environments. With the help of finite difference simulations we model the decoupling mechanism of guided waves at the contact between the descending oceanic plate and the upper plate crust where the wave guide is opened due to the equalization of seismic velocities. Provided that suited stations are available, this mechanism allows for the use of intermediate depth seismicity to study the shallow subduction zone structure ( <= 100 km depth)
Basement-cored uplift provinces are often characterized by high-angle reverse faulting along preexisting crustal heterogeneities, which may greatly affect the mechanics of deformation and the coupling between erosion and orogenic structure. Herein we construct a coupled deformation-erosion model to understand the mechanics and erosion of mountain belts in which the spatial distribution of deformation is largely influenced by the presence of preexisting high-angle faults. In this case, deformation is accommodated along, and topography is built above, these structures. This topographic loading leads to increasing lithostatic stresses beneath these regions. As a result, active deformation may migrate to frictionally stronger structures in adjacent regions where lithostatic loading is lower. The migration of deformation to such nearby structures depends on the Hubbert-Rubey pore fluid pressure ratio of the crust (lambda), the orientations of the frictionally weaker and stronger preexisting faults (beta(1) and beta(2), respectively), the friction coefficients (mu(b1) and mu(b2)) and Hubbert-Rubey fluid-pressure ratios (lambda(b1) and lambda(b2)) of these faults, and the surface slope of the topography above the frictionally weaker structure (alpha), assuming zero surface slope above the frictionally stronger structure. In general, we found that for a given alpha and beta(1), as mu(b1) increases lambda=lambda(b1)=lambda(b2) increases, and beta(2) decreases, the value of mu(b2) required to force deformation to migrate increases. However, as erosional processes lead to decreasing surface slopes, deformation will be inhibited from migrating to frictionally stronger structures in adjacent regions. Our model results may help to explain some aspects of the deformation observed and the possible correlation between precipitation and the migration of deformation within these tectonic provinces
Stable isotope paleoaltimetry makes use of systematic trends in the distribution and isotopic composition of modern precipitation with climate and topography, and of the potential to estimate the isotopic composition of paleoprecipitation from authigenic (in-situ formed) minerals. To illustrate the usefulness as well as potential limitations of this method, we review (1) processes controlling the isotopic composition of modem precipitation, (2) stable isotope data from modern precipitation across regions of high topography, and (3) stable isotope data from authigenic minerals that have been used to infer paleotopography. From this we conclude that stable isotope studies of authigenic minerals can permit useful inferences on paleotopography, with uncertainties that critically depend on a detailed understanding of local- to global-scale paleoclimate during the time interval of interest
Subduction factory : 1. Theoretical mineralogy, densities, seismic wave speeds, and H2O contents
(2005)
Analysis of contemporary and past gully erosion and infilling processes allowed to reconstruct the long-term evolution of a permanent gully system under cropland. An active and a buried gully under cropland were investigated. The recent sediment deposits within the active gully, adjacent to the buried gully, showed that the recent gully was filling in at a mean rate of 6.4 cm a(-1). In the buried gully, several erosion and deposition phases could be identified and the type of deposited sediments revealed a complex infilling history. Charcoal, pottery and brick fragments of different sizes were found at all depths of the gully infilling. Their age indicates that the first gully incised after the midst of the 17th century, most probably in the second half of the 18th century or the early 19th century. Gully morphology and analogy with the processes in the recent gully indicate that the buried gully filled in rapidly. Overall, five cycles of cut and fill occurred in 350 years or less and four cycles even within little less than a few decades, indicating that gully development and infilling under cropland can be very rapid processes. (c) 2005 Elsevier B.V All rights reserved
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.
In a recent BAMS article, it is argued that community-based Open Source Software (OSS) could foster scientific progress in weather radar research, and make weather radar software more affordable, flexible, transparent, sustainable, and interoperable.
Nevertheless, it can be challenging for potential developers and users to realize these benefits: tools are often cumbersome to install; different operating systems may have particular issues, or may not be supported at all; and many tools have steep learning curves.
To overcome some of these barriers, we present an open, community-based virtual machine (VM). This VM can be run on any operating system, and guarantees reproducibility of results across platforms. It contains a suite of independent OSS weather radar tools (BALTRAD, Py-ART, wradlib, RSL, and Radx), and a scientific Python stack. Furthermore, it features a suite of recipes that work out of the box and provide guidance on how to use the different OSS tools alone and together. The code to build the VM from source is hosted on GitHub, which allows the VM to grow with its community.
We argue that the VM presents another step toward Open (Weather Radar) Science. It can be used as a quick way to get started, for teaching, or for benchmarking and combining different tools. It can foster the idea of reproducible research in scientific publishing. Being scalable and extendable, it might even allow for real-time data processing.
We expect the VM to catalyze progress toward interoperability, and to lower the barrier for new users and developers, thus extending the weather radar community and user base.
The results of streamflow trend studies are often characterized by mostly insignificant trends and inexplicable spatial patterns. In our study region, Western Austria, this applies especially for trends of annually averaged runoff. However, analysing the altitudinal aspect, we found that there is a trend gradient from higher-altitude to lower-altitude stations, i.e. a pattern of mostly positive annual trends at higher stations and negative ones at lower stations. At midaltitudes, the trends are mostly insignificant. Here we hypothesize that the streamflow trends are caused by the following two main processes: on the one hand, melting glaciers produce excess runoff at higher-altitude watersheds. On the other hand, rising temperatures potentially alter hydrological conditions in terms of less snowfall, higher infiltration, enhanced evapotranspiration, etc., which in turn results in decreasing streamflow trends at lower-altitude watersheds. However, these patterns are masked at mid-altitudes because the resulting positive and negative trends balance each other. To support these hypotheses, we attempted to attribute the detected trends to specific causes. For this purpose, we analysed trends of filtered daily streamflow data, as the causes for these changes might be restricted to a smaller temporal scale than the annual one. This allowed for the explicit determination of the exact days of year (DOYs) when certain streamflow trends emerge, which were then linked with the corresponding DOYs of the trends and characteristic dates of other observed variables, e.g. the average DOY when temperature crosses the freezing point in spring. Based on these analyses, an empirical statistical model was derived that was able to simulate daily streamflow trends sufficiently well. Analyses of subdaily streamflow changes provided additional insights. Finally, the present study supports many modelling approaches in the literature which found out that the main drivers of alpine streamflow changes are increased glacial melt, earlier snowmelt and lower snow accumulation in wintertime.
Climate change is likely to impact the seasonality and generation processes of floods in the Nordic countries, which has direct implications for flood risk assessment, design flood estimation, and hydropower production management. Using a multi-model/multi-parameter approach to simulate daily discharge for a reference (1961–1990) and a future (2071–2099) period, we analysed the projected changes in flood seasonality and generation processes in six catchments with mixed snowmelt/rainfall regimes under the current climate in Norway. The multi-model/multi-parameter ensemble consists of (i) eight combinations of global and regional climate models, (ii) two methods for adjusting the climate model output to the catchment scale, and (iii) one conceptual hydrological model with 25 calibrated parameter sets. Results indicate that autumn/winter events become more frequent in all catchments considered, which leads to an intensification of the current autumn/winter flood regime for the coastal catchments, a reduction of the dominance of spring/summer flood regimes in a high-mountain catchment, and a possible systematic shift in the current flood regimes from spring/summer to autumn/winter in the two catchments located in northern and south-eastern Norway. The changes in flood regimes result from increasing event magnitudes or frequencies, or a combination of both during autumn and winter. Changes towards more dominant autumn/winter events correspond to an increasing relevance of rainfall as a flood generating process (FGP) which is most pronounced in those catchments with the largest shifts in flood seasonality. Here, rainfall replaces snowmelt as the dominant FGP primarily due to increasing temperature.We further analysed the ensemble components in contributing to overall uncertainty in the projected changes and found that the climate
projections and the methods for downscaling or bias correction tend to be the largest contributors. The relative role of hydrological parameter uncertainty, however, is highest for those catchments showing the largest changes in flood seasonality, which confirms the lack of robustness in hydrological model parameterization for simulations under transient hydrometeorological conditions.